REVIEW 3 major objections 4 minor 33 references
Classical Simulation of an All-Optical Toffoli Gate using Soliton Scattering through Asymmetric Potential Wells
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Two asymmetric Pöschl–Teller wells can implement an all-optical Toffoli gate by reflecting solitons only when both control bits are set.
desk verdict A plausible incremental numerical demonstration of a soliton-based Toffoli gate, undermined by an unverified gap between single-soliton transport thresholds and two-soliton ordering in the operational window. 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 operating mechanism is the velocity-selective scattering of bright solitons by Pöschl–Teller wells, $V(x)=-\sum_i V_i\,\mathrm{sech}^2((x-x_i)/w_i)$, chosen near the reflectionless linear limit $V_i w_i^2=1$. A soliton slows when it enters the first well; if its reduced velocity stays below the critical velocity of the second well, it reflects there and travels back through the first well to emerge on the incident side. The gate is read out through reflection/transmission coefficients $R_j,L_j,T_j$ defined by integrals of $|\psi_j|^2$ in the reflected, trapped, and transmitted regions; the control-bit encoding is the presence or absence of the wells, and the target-bit encoding is the relative order of the two soliton peaks. Asymmetry $V_1\neq V_2$ separates the critical velocities so that both single-well cases transmit while the double-well case reflects.
What would settle it
Directly simulate the four target-$|1\rangle$ rows, for example by starting with $\psi_1$ left of $\psi_2$, using $u=1.4$, $v=0.51$, and $\alpha=1.08$, and record which soliton peak is on the right at the final time. The central claim fails if any configuration other than $|11\rangle$ exchanges the order, or if $|11\rangle$ fails to exchange it, even when the reported reflection and transmission thresholds are satisfied.
Extended reading notes
Core claim
The central discovery is a parameter regime in which the conditional swap of two solitons reproduces the Toffoli truth table. With soliton amplitude $u=1.4$ and initial velocity in the window $0.4925\le v\le0.5350$, a single asymmetric well at either location transmits the soliton with $T\ge0.97$, while both wells together reflect it with $R_{11}\ge0.94$; the soliton that first reaches the wells is slowed, reflected by the second well, and passes back through the first, so the pair exits with its spatial order exchanged. This exchange occurs only in the double-well configuration, which is exactly the CCNOT rule that the target flips only when both controls are $|1\rangle$. The authors map how the operational region in the $(v,u)$ plane depends on the asymmetry parameter $\alpha$ and find that $\alpha\approx1.04$ (for $1.0<u<1.2$) and $\alpha\approx1.08$ (for $1.2<u<1.6$) give broader windows than the symmetric $\alpha=1.00$ case. They state that the four input rows with target initially $|1\rangle$ are inferred from the symmetry of the uncoupled equations, and that cross-coupling above roughly $g_{12}=0.3$ eliminates the operational region.
Load-bearing premise
The load-bearing premise is that the thresholds $R_{11}>0.9$, $T_{10}>0.9$, and $T_{01}>0.9$ are enough to guarantee the two solitons end up in the correct final order for all eight input rows, even though the final ordering is directly simulated only for the four target-$|0\rangle$ cases.
Editorial extensions
If this is right
- If the central claim is correct, a classical optical circuit element whose transport coefficients match the CCNOT truth table exists for the tested inputs, without requiring quantum superposition or entanglement.
- The reported operational window at $u=1.4$ is $0.4925\le v\le0.5350$, within which the double-well reflectance stays above 0.94 and the single-well transmittances stay above 0.97.
- Tuning the asymmetry parameter broadens the operational window: roughly $\alpha=1.04$ works best for $1.0<u<1.2$ and $\alpha=1.08$ for $1.2<u<1.6$, compared with the narrower symmetric case.
- Weak inter-component coupling ($g_{12}$ up to about 0.3) preserves the gate, while stronger coupling shrinks the operational region and eventually removes it.
- Because the Toffoli gate is universal for reversible classical computation, a working all-optical version would support scaling the same soliton-scattering principle to more complex logic circuits.
Reading between the lines
- The four target-$|1\rangle$ rows are inferred rather than simulated; a direct numerical check of those rows would settle whether the truth table holds under weak coupling.
- The same critical-velocity ordering mechanism could be recreated with other potential shapes, so the scheme points to a general design rule for soliton-based conditional logic rather than a special property of Pöschl–Teller wells.
- Because the gate is not reversible and relies on dissipation-free but direction-dependent scattering, practical use would most naturally be in unidirectional optical circuits or as a classical simulator for benchmarking quantum circuits, not inside a reversible computation loop.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes an all-optical implementation of a Toffoli (CCNOT) gate using the scattering of two distinguishable spatial solitons by up to two Pöschl-Teller wells. The target qubit is encoded in the relative left-right ordering of the two solitons; the control bits are encoded by the presence or absence of the two wells. Solving the uncoupled/coupled nonlinear Schrödinger equations numerically, the authors scan soliton amplitude and velocity and identify parameter regions in which a single soliton is mostly reflected by the double-well configuration (R11>0.9) and mostly transmitted by each single-well configuration (T10>0.9, T01>0.9). They report that asymmetric well depths broaden this operational window relative to symmetric wells, and that weak inter-component coupling progressively degrades the window. The paper is explicitly a classical-analog simulation of the truth table, not a quantum implementation.
Significance. If the central claim were fully established, the work would extend the known soliton-based optical CNOT analog (Ref. [29]) to a three-qubit universal classical reversible gate, with a useful design observation that asymmetry broadens the usable velocity/amplitude region. The paper is honest about its scope: Sec. VI explicitly lists the absence of quantum features, the difficulty of dynamically engineering wells, and the non-reversibility of the implementation. The systematic parameter scans in Figs. 3-4 and the clear statement of thresholds are commendable. However, the current evidence is an existence scan for single-soliton transport coefficients rather than a verified demonstration of clean two-soliton ordering across the claimed operational windows.
major comments (3)
- [Sec. V, Eq. (4), Figs. 3-4] The operational regions are defined solely by the single-soliton transport thresholds R11>0.9, T10>0.9, T01>0.9, but the Toffoli output is the relative ordering of two solitons. A reflectance of 0.94 at the window edge allows about 6% of the norm to be transmitted or trapped, so the final two-component state can contain both a swapped and an unswapped ordering, making the truth-table output ambiguous. The only direct two-soliton validation is the single representative point in Fig. 1; the shaded regions of Figs. 3-4 are built from single-soliton coefficients, and no ordering fidelity (e.g., the fraction of the total norm for which ψ1 ends up left of ψ2) or trapped-mass bound L_j is reported anywhere across the window. I therefore do not think the central claim that these shaded regions are Toffoli operational windows is yet established. Please compute a direct two-soliton output measure for at least representative points, including the window boundaries, and report the resulting ordering fidelity and trapped mass.
- [Sec. II/V] The manuscript gives no numerical convergence information: no time step, grid spacing, domain size, or convergence test for the split-step Fourier method is reported, and the operational window boundaries are quoted to four decimal places (e.g., 0.4925≤v≤0.5350). Because every conclusion rests on these numerical simulations, the paper should include a convergence study in u and v, report the accuracy of the transport coefficients, and use that to set error bars on the claimed windows.
- [Sec. V A, Fig. 3] The headline result that asymmetry significantly broadens the operational parameter window is supported only by visual inspection of shaded areas in Fig. 3. Please define a quantitative measure of the operational window (e.g., area or velocity width at fixed u), report it for each α, and show its sensitivity to the chosen 0.9 thresholds. This is necessary to compare α=1.00, 1.04, 1.08, and 1.12 and to determine whether the broadening persists when the two-soliton ordering fidelity of Major Comment 1 is used instead of single-soliton thresholds.
minor comments (4)
- [Sec. III] The statement that the remaining four input cases can be inferred directly is correct under exact label-exchange symmetry, but the text should explicitly note that this symmetry holds for g12≠0 as well whenever g11=g22 and g12=g21, and that the inference assumes the numerical scheme preserves the symmetry.
- [Eq. (4)] The integration limits l1 and l2 are never given; please specify the values used and confirm that they are chosen after the reflected and transmitted parts have separated from the trapped region.
- [Sec. IV/Fig. 1] The relation between the encoding (target |0⟩ when ψ1 is to the right of ψ2) and the initialization in Eq. (2) (with δ=-10, ψ1 at x0=-30 and ψ2 at x0-10=-40) should be stated explicitly.
- [Sec. IV vs Sec. V] The qualitative discussion in Sec. IV refers to total reflection and complete transmission, while Sec. V defines the operational region by thresholds of 0.9; please reconcile these descriptions so the reader knows that partial transport is expected and accepted.
Circularity Check
No significant circularity: the numerical scan is self-contained, and the R/T threshold proxy is an assumption about gate validity, not a circular reduction.
full rationale
The paper's derivation chain is self-contained. The Toffoli gate is defined by the relative spatial ordering of two soliton components, and the control bits by well presence/absence. The core demonstration in Fig. 1 is a direct two-soliton simulation for the four target-|0> input configurations, showing the ordering flips only for control |11>. The four target-|1> cases are inferred from label-exchange symmetry: in the uncoupled case (g12=0) the two components evolve identically, so swapping which soliton starts at which position simply swaps the final labels, making the inference exact. The operational windows in Sec. V are defined by single-soliton transport thresholds (R11>0.9, T10>0.9, T01>0.9), and the paper explicitly states these are the conditions used to map the regions in Figs. 3 and 4. This is a proxy for correct Toffoli ordering, not a circular definition of it: the thresholds are not the truth-table condition, and the truth-table condition is not fitted from them. The asymmetry-broadening claim is a comparison of the sizes of these threshold-defined regions across alpha, which is an independent numerical observation rather than a prediction derived from its own input. All cited prior work (Refs. [23,24,29]) is by other authors, so no self-citation chain is load-bearing. The main scientific weakness is that the shaded regions are not verified by direct two-soliton ordering checks, but that is a gap in evidence, not circularity.
Assumptions & free parameters
free parameters (6)
- V0 (base well depth of second well) =
4
- alpha (asymmetry parameter, V1 = alpha V0) =
scanned 1.00 to 1.12; optimal 1.04 or 1.08 depending on u
- well width w =
0.5
- soliton amplitude u and velocity v =
e.g., u=1.4, v in [0.4925, 0.5350]
- initial separation delta =
-10
- cross-coupling g12 =
0 (ideal); scanned 0 to 0.3
assumptions (5)
- domain assumption The system is governed by the coupled nonlinear Schrödinger (Manakov) equations with attractive nonlinearity.
- standard math Initial wavefunctions are exact bright solitons of the homogeneous system, well separated from the potential region.
- domain assumption In the uncoupled case (g12=0), the two components evolve independently, so single-component transport coefficients fully determine the order flip.
- ad hoc to paper Transport coefficients R, L, T at final time, with boundaries chosen to avoid re-entry, faithfully represent reflected and transmitted fractions.
- ad hoc to paper Satisfying R11>0.9, T10>0.9, T01>0.9 suffices for correct Toffoli truth table behavior, including the inferred target-|1> cases.
Cite this review
Pith. "Pith review of Classical Simulation of an All-Optical Toffoli Gate using Soliton Scattering through Asymmetric Potential Wells." pith.science (2026). https://pith.science/paper/MLKYV3WO
@misc{pith2026250606883,
author = {Pith},
title = {Pith review of: Classical Simulation of an All-Optical Toffoli Gate using Soliton Scattering through Asymmetric Potential Wells},
year = {2026},
howpublished = {\url{https://pith.science/paper/MLKYV3WO}},
note = {Machine review of arXiv:2506.06883}
}
read the original abstract
We propose and numerically simulate an all-optical Toffoli (controlled-controlled-NOT) gate based on the scattering of spatial solitons by asymmetric P\"oschl-Teller potential wells. In our scheme, the logical state of the target bit is encoded in the relative spatial ordering of two distinguishable soliton components, while the control bits are represented by the presence or absence of external potential wells. We solve the nonlinear Schr\"odinger equations governing the soliton dynamics and systematically scan soliton amplitude and velocity, analyzing reflection and transmission coefficients to identify the operational conditions for Toffoli gate behavior. Our results demonstrate that introducing asymmetry in the potential wells significantly broadens the operational parameter window compared to symmetric configurations. We also investigate the impacts of varying degrees of asymmetry and soliton amplitude on gate performance. Furthermore, we examine the influence of weak inter-component coupling and confirm that it is not essential for gate operation. These findings generalize earlier soliton-based CNOT simulations and support the broader feasibility of classical analog modeling of multi-qubit logic gates in nonlinear optical systems.
Figures
Reference graph
Works this paper leans on
-
[29]
Amaria Javed, T. Uthayakumar, and U. Al Khawaja, “Simulat- ing an all-optical quantum controlled-NOT gate using soliton scattering by a reflectionless potential well,” Physics Letters A 429, 127949 (2022)
work page 2022
-
[1]
The role of optics in computing,
Rodney S. Tucker, “The role of optics in computing,” Nature Photonics4, 405–405 (2010)
work page 2010
-
[2]
G. P. Agrawal,Nonlinear Fiber Optics, Optics and Photonics (Elsevier Science, 2013)
work page 2013
-
[3]
U. Al Khawaja, S. M. Al-Marzoug, and H. Bahlouli, “All- optical switches, unidirectional flow, and logic gates with dis- crete solitons in waveguide arrays,” Optics Express24, 11062– 11074 (2016)
work page 2016
-
[4]
Amplify- ing optical signals with discrete solitons in waveguide arrays,
Amaria Javed, Alaa Shaheen, and U. Al Khawaja, “Amplify- ing optical signals with discrete solitons in waveguide arrays,” Physics Letters A384, 126654 (2020)
work page 2020
-
[5]
Two- Dimensional Discrete Cavity Solitons: Switching and All- Optical Gates,
K. M. Aghdami, M. Golshani, and R. Kheradmand, “Two- Dimensional Discrete Cavity Solitons: Switching and All- Optical Gates,” IEEE Photonics Journal4, 1147–1154 (2012)
work page 2012
-
[6]
All-optical switch based on the local nonlin- ear mach-zehnder interferometer,
Yaw-Dong Wu, Min-Lun Huang, Mao-Hsiung Chen, and Rong-Zhan Tasy, “All-optical switch based on the local nonlin- ear mach-zehnder interferometer,” Opt. Express15, 9883–9892 (2007)
work page 2007
-
[7]
Blocking and rout- ing discrete solitons in two-dimensional networks of nonlinear waveguide arrays,
D.N. Christodoulides and E.D. Eugenieva, “Blocking and rout- ing discrete solitons in two-dimensional networks of nonlinear waveguide arrays,” Physical Review Letters87, 233901 (2001)
work page 2001
Show all 33 references
-
[8]
Adding bi- nary numbers with discrete solitons in waveguide arrays,
Alaa Shaheen, Amaria Javed, and U Al Khawaja, “Adding bi- nary numbers with discrete solitons in waveguide arrays,” Phys- ica Scripta95, 085107 (2020)
2020
-
[9]
Stable dis- crete soliton molecules in two-dimensional waveguide arrays,
P.S. Vinayagam, A. Javed, and U. Al Khawaja, “Stable dis- crete soliton molecules in two-dimensional waveguide arrays,” Physical Review A98, 063839 (2018)
2018
-
[10]
All-optical sequential logic using coupled-resonators-based optical switches,
Jacob Scheuer and Meir Orenstein, “All-optical sequential logic using coupled-resonators-based optical switches,” Journal of the Optical Society of America B22, 1053 (2005)
2005
-
[11]
All-optical switching in nonlinear X junctions,
J.P. Sabini, N. Finlayson, and G.I. Stegeman, “All-optical switching in nonlinear X junctions,” Applied Physics Letters 55, 1176–1178 (1989)
1989
-
[12]
Simulation of universal optical logic gates under energy sharing collisions of Manakov solitons and fulfillment of practical optical logic criteria,
M. Vijayajayanthi, T. Kanna, and M. Lakshmanan, “Simulation of universal optical logic gates under energy sharing collisions of Manakov solitons and fulfillment of practical optical logic criteria,” Physical Review E108, 054213 (2023)
2023
-
[13]
State transformations of colliding optical solitons and possi- ble application to computation in bulk media,
Mariusz H. Jakubowski, Ken Steiglitz, and Richard Squier, “State transformations of colliding optical solitons and possi- ble application to computation in bulk media,” Physical Review E58, 6752–6758 (1998)
1998
-
[14]
Reversible logic elements as a new field of application of optical solitons,
Andrei I. Maimistov, “Reversible logic elements as a new field of application of optical solitons,” Quantum Electronics25, 1009 (1995)
1995
-
[15]
Energy-Exchange Interactions between Col- liding Vector Solitons,
Charalambos Anastassiou, Mordechai Segev, Kenneth Steiglitz, JA Giordmaine, Matthew Mitchell, Ming-Feng Shih, Song Lan, and John Martin, “Energy-Exchange Interactions between Col- liding Vector Solitons,” Physical Review Letters83, 2332–2335 (1999)
1999
-
[16]
All-optical switching and logic gating with spatial solitons in liquid crystals,
Marco Peccianti, Claudio Conti, Gaetano Assanto, Antonio De Luca, and Cesare Umeton, “All-optical switching and logic gating with spatial solitons in liquid crystals,” Applied Physics Letters81, 3335–3337 (2002)
2002
-
[17]
Information transfer between solitary waves in the saturable Schr¨odinger equation,
Mariusz H. Jakubowski, Ken Steiglitz, and Richard Squier, “Information transfer between solitary waves in the saturable Schr¨odinger equation,” Physical Review E56, 7267–7272 (1997)
1997
-
[18]
Bistable soliton switching dynamics in a PT-symmetric coupler with saturable nonlinearity,
Ambaresh Sahoo, Dipti Kanika Mahato, A. Govindarajan, and Amarendra K. Sarma, “Bistable soliton switching dynamics in a PT-symmetric coupler with saturable nonlinearity,” Physical Review A105, 063503 (2022)
2022
-
[19]
Varia- tional analysis of soliton scattering by external potentials,
B A Umarov, A Messikh, N Regaa, and B B Baizakov, “Varia- tional analysis of soliton scattering by external potentials,” Jour- nal of Physics: Conference Series435, 012024 (2013)
2013
-
[20]
Trapping of discrete solitons by defects in nonlinear waveguide arrays,
L. Morales-Molina and R.A. Vicencio, “Trapping of discrete solitons by defects in nonlinear waveguide arrays,” Optics Let- ters31, 966–968 (2006)
2006
-
[21]
Unidirectional soliton flows in PT-symmetric poten- tials,
U. Al Khawaja, S. M. Al-Marzoug, H. Bahlouli, and Yuri S. Kivshar, “Unidirectional soliton flows in PT-symmetric poten- tials,” Physical Review A88, 023830 (2013)
2013
-
[22]
Critical soliton speed for quantum reflec- tion by a reflectionless potential well,
U. Al Khawaja, “Critical soliton speed for quantum reflec- tion by a reflectionless potential well,” Physical Review E103, 062202 (2021)
2021
-
[23]
Unidirectional flow of composite bright-bright solitons through asymmetric double potential barriers and wells,
Amaria Javed, T. Uthayakumar, M. O. D. Alotaibi, S. M. Al- Marzoug, H. Bahlouli, and U. Al Khawaja, “Unidirectional flow of composite bright-bright solitons through asymmetric double potential barriers and wells,” Communications in Non- linear Science and Numerical Simulation1...
2021
-
[24]
Unidirectional flow of symbiotic solitons and nonlinear modes of the Schr¨odinger equation with an external potential,
Amaria Javed, H. Susanto, Rudy Kusdiantara, and I. Kourakis, “Unidirectional flow of symbiotic solitons and nonlinear modes of the Schr¨odinger equation with an external potential,” The Eu- ropean Physical Journal Plus137, 1170 (2022)
2022
-
[25]
Reflectionless potentials and resonant scattering of flat-top and thin-top solitons,
L. Al Sakkaf and U. Al Khawaja, “Reflectionless potentials and resonant scattering of flat-top and thin-top solitons,” Physical 7 Review E107, 014202 (2023)
2023
-
[26]
Soliton-based optical switching in waveguide arrays,
W. Kr ´olikowski and Y .S. Kivshar, “Soliton-based optical switching in waveguide arrays,” Journal of the Optical Society of America B: Optical Physics13, 876–887 (1996)
1996
-
[27]
Directional flow of solitons with asymmetric potential wells: Soliton diode,
M. Asad-uz-zaman and U. Al Khawaja, “Directional flow of solitons with asymmetric potential wells: Soliton diode,” EPL (Europhysics Letters)101, 50008 (2013)
2013
-
[28]
Unidirectional flow of discrete solitons in waveguide arrays,
U. Al Khawaja and Andrey A. Sukhorukov, “Unidirectional flow of discrete solitons in waveguide arrays,” Optics Letters 40, 2719 (2015)
2015
-
[30]
Kivshar and Govind P
Yuri S. Kivshar and Govind P. Agrawal,Optical Solitons: From Fibers to Photonic Crystals(Academic Press, 2003)
2003
-
[31]
Bright solitons in coupled defocus- ing NLS equation supported by coupling: Application to Bose–Einstein condensation,
Sadhan K. Adhikari, “Bright solitons in coupled defocus- ing NLS equation supported by coupling: Application to Bose–Einstein condensation,” Physics Letters A346, 179–185 (2005)
2005
-
[32]
Bemerkungen zur Quantenmechanik des anharmonischen Oszillators,
G. P ¨oschl and E. Teller, “Bemerkungen zur Quantenmechanik des anharmonischen Oszillators,” Zeitschrift f ¨ur Physik83, 143–151 (1933)
1933
-
[33]
Active meta-optics and nanophotonics with halide perovskites,
Dragomir Neshev and Igor Aharonovich, “Active meta-optics and nanophotonics with halide perovskites,” Light: Science & Applications7, 58 (2018)
2018
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.