Entanglement and firewalls in quantum circuit model of black hole evaporation

Da-Shin Lee; Tien Hsieh; Yi-Chen Lu

REVIEW 3 major objections 3 minor 23 references

A one-parameter scrambling of the infalling-matter qubits changes where and when firewall entanglement appears in a tripartite quantum-circuit model of black-hole evaporation, and for a derived range of the scrambling strength the Hawking r

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 10:40 UTC pith:IMEVUTFX

load-bearing objection A genuinely new θ-parameterized extension of the 2018 circuit model, internally coherent, but the two-level Fock truncation is unchecked at exactly the small-Mω regime where the firewall and all-Mω claims live. the 3 major comments →

arxiv 2606.21181 v2 pith:IMEVUTFX submitted 2026-06-19 gr-qc

Entanglement and firewalls in quantum circuit model of black hole evaporation

Yi-Chen Lu , Tien Hsieh , Da-Shin Lee This is my paper

classification gr-qc

keywords black hole information paradoxfirewallquantum circuit modelquantum scramblingentanglement negativityHawking radiationquantum monogamyPage curve

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish that the initial state of the qubits falling into a black hole—not just the evaporation dynamics—controls whether a firewall forms. Applying a one-parameter scrambling unitary to that initial state enlarges the range of Mω (black-hole mass times radiation frequency) for which a firewall appears, and maximal scrambling makes it appear earlier. Within an analytically determined interval of the scrambling parameter, the model also allows all information initially in the black hole to be carried off by the late radiation, in contrast to the unscrambled case. If correct, this means that in these toy models quantum monogamy, not the gate dynamics alone, decides whether an infalling observer hits a firewall. The result matters because it isolates a clean, qubit-level mechanism connecting scrambling, entropy, and information recovery.

Core claim

The central claim is that one-parameter scrambling of the initial black-hole qubit state (θ = 0 no scrambling, θ = π/2 maximal scrambling) changes the entanglement structure enough to shift the firewall transition. At θ = π/2 a firewall between BH and JR emerges earlier in the circuit evolution than at θ = 0. For θ lying in the range derived in Eq. (53), information is carried away by radiation for all Mω, and the initial black-hole qubit state can be reconstructed from its imprint on the final radiation, even though it was initially behind the horizon. The paper presents analytic negativity and mutual-information results from step 4 through step 7 supporting this picture.

What carries the argument

The machinery is a one-parameter scrambling unitary acting on the ground state of infalling-matter qubits, interpolating between no scrambling and maximal scrambling, followed by the circuit's CNOT-U gates that create Hawking pairs. The U gate is parameterized by tan γ = exp(−4πMω), so γ encodes how strongly a Hawking pair is entangled; small Mω gives strong pair entanglement. Firewalls are diagnosed by entanglement negativity and mutual information between BH and JR, and the analytic reduced density matrices (Appendices A–B) carry the argument. Quantum monogamy is the mechanism: the BH qubit cannot remain maximally entangled with both its interior partner and the early radiation, so when th

Load-bearing premise

The load-bearing assumption is that each Hawking mode can be represented by only its vacuum and one-particle states, with tan γ = exp(−4πMω); because e^(−4πMω) is order one exactly in the small-Mω regime where firewall claims are strongest, higher modes could change the entanglement and negativity results.

What would settle it

Compute the same BH–JR negativity and mutual information keeping the next Fock levels (|0⟩, |1⟩, |2⟩, ...) of the Hawking state with tan γ = exp(−4πMω), at θ = π/2 and small Mω, and compare the firewall onset and the Eq. (53) information-carrying range; if the firewall transitions shift or disappear, the truncated two-level computation is the decisive approximation. A numerical scan over θ and Mω with truncation level N → ∞ would settle it.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

For the maximally scrambled initial state θ = π/2, the firewall transition occurs earlier in the circuit than with θ = 0; the paper locates this shift in the negativity between BH and JR.
Introducing initial-state scrambling enlarges the range of Mω for which a firewall forms, so firewall behavior is not fixed solely by mass and frequency but by initial qubit entanglement structure.
When θ is inside the analytically derived interval (Eq. 53), information is carried away by radiation for every Mω, not just the large-Mω regime.
The initial black hole qubit state is recoverable from the final radiation state; unitary gate dynamics preserves the information behind the horizon as an imprint on radiation.
Quantum monogamy—not a modification of the gate dynamics—carries the firewall structure between BH and JR in this model.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Editorial inference: the two-Fock-mode truncation is least reliable at small Mω, where e^(−4πMω) is order one and discarded higher modes are populated at order one; if those modes are included, the precise Mω boundary of the firewall window could move or dissolve.
Editorial inference: scrambling can be read as distributing the initial information so that no single BH–JR pair can carry it; this suggests a direct tradeoff between scrambling strength and the time at which information exits, which the model can test by varying θ continuously.
Editorial inference: a natural next calculation is finite-N versions of the Hawking state with N > 2 Fock levels; if the θ = π/2 firewall advance survives as N grows, the mechanism is robust, and if not, the firewall effect is an artifact of truncation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

A genuinely new θ-parameterized extension of the 2018 circuit model, internally coherent, but the two-level Fock truncation is unchecked at exactly the small-Mω regime where the firewall and all-Mω claims live.

read the letter

Worth your time, but read with the truncation in view. The paper extends the 2018 circuit model [1] by scrambling the infalling-matter qubits with a one-parameter unitary (θ = 0 no scrambling, θ = π/2 maximal), and claims an earlier firewall at maximal scrambling plus an analytic θ-window, Eq. (53), where radiation carries the information away for all values of Mω. What is actually new is the θ knob and the θ-dependent results: the earlier-firewall claim, the enlarged Mω range, and the analytic interval (53). The θ = 0 limit reduces to [1], which is a clean consistency check. The setup is coherent — CNOT-U gate, tan γ = e^{−4πMω}, all following [1] — and the paper promises analytic appendices (wave functions and reduced density matrices for arbitrary θ); I could not inspect those sections, but the visible setup is clean. No parameters are fitted; θ is a hand-tuned knob and γ is a physical input. Within the truncated model, the firewall diagnostic (vanishing negativity between BH and JR) is applied consistently.

The soft spot is the two-level Fock truncation, and it is load-bearing. Each Hawking mode is reduced to |0⟩ and |1⟩ with λ = e^{−4πMω}; the discarded weight is λ⁴ = e^{−16πMω}. That is not small where the claims live: at Mω = 0.05 it is 8%, at 0.02 it is 37%, and at 0.01 it is 61%. The firewall boundary and the Eq. (53) window are determined exactly in this regime, and all negativity and mutual-information results sit inside the truncated subspace. Retaining |2⟩ can plausibly flip the sign of the BH–JR negativity — the firewall criterion — and shift the θ-window. The paper calls these the 'two dominant contributions,' which is not right at small Mω. The truncation is inherited from [1], so the original model shares it, but the new universal claim makes control of it unavoidable. A second, lesser concern: the information-retrieval statement is near-tautological for a closed unitary circuit unless shown for the radiation subsystem alone via an explicit recovery map; the visible text does not show this.

Who this is for: people working on qubit-circuit toy models of evaporation and the monogamy/firewall debate in that restricted setting. For that audience the θ-dependent extension is genuinely useful. I would send it to peer review, but the referee should insist on either a quantitative check of the truncation (retain n ≥ 2 levels at a few small-Mω values and show the negativity sign and Eq. (53) survive) or a restriction of the claims to the regime where the truncation is controlled, plus a sharper retrieval statement. As it stands, the internal story holds up; the extrapolation to the Hawking state does not.

Referee Report

3 major / 3 minor

Summary. This paper revisits the 2018 quantum circuit model of black hole evaporation of Ref. [1], in which BH, JR, and ER are represented by qubits and evolved through CNOT-U gates with a squeezing angle gamma related to Mω by tanγ = e^{-4πMω}. The authors apply a one-parameter scrambling unitary U(θ) to the infalling-matter ground state (θ = 0 no scrambling, θ = π/2 maximal scrambling). They claim that maximal scrambling makes firewalls appear earlier, that scrambling enlarges the Mω range for firewall formation, and that for θ in an analytically determined interval (53) Hawking radiation carries away the information for all Mω, with the initial BH qubit state recoverable from the final radiation. The paper reports analytic wave functions and reduced density matrices for arbitrary θ and computes entropies, mutual information, and negativity.

Significance. If correct, the paper gives a compact analytic illustration of how initial-state scrambling and quantum monogamy can shift the firewall in a toy circuit model and allow information transfer across all Mω. Its strengths are the closed-form wave functions and reduced density matrices for arbitrary θ, the physically motivated relation γ(Mω), and the absence of fitted parameters: θ is a tunable model parameter rather than a data fit. The main limitation is that all results are obtained in a two-level Fock truncation of the Hawking state, and the claims are made precisely in the regime where this truncation is least controlled. The paper is a potentially useful toy-model study, but its broader conclusions are not yet established.

major comments (3)

[Sec. 2, Eq. (3); Apps. A–B] The two-level Fock truncation is uncontrolled in the regime where the firewall and all-Mω claims live. The full two-mode squeezed state is sqrt(1−λ²)∑_{n≥0}λ^n|n,n> with λ=e^{−4πMω}; keeping only |0>,|1> discards probability λ^4=e^{−16πMω}. This is not small at the small Mω values emphasized in the Introduction: at Mω=0.05, λ^4≈8%; at 0.02, ≈37%; at 0.01, ≈61%. The negativity and mutual-information results in Apps. A–B are computed entirely in this truncated subspace, and the sign of the BH–JR negativity (the firewall diagnostic) and the endpoints of window (53) can shift when |2> is retained. Please provide a truncation-error estimate, extend to a finite cutoff with a convergence check, or restrict the claims to a range where e^{−16πMω} is negligible.
[Abstract; Intro; Eq. (53)] The statement that 'information can be carried away by radiation for all values of Mω' for θ in (53) is broader than the evidence. The interval (53) is derived in Appendix B within the two-qubit truncated Hilbert space; this does not by itself justify an all-Mω statement that includes Mω→0, where the discarded Fock levels have O(1) weight. Either prove that higher Fock contributions leave the interval invariant, or replace 'all Mω' with a qualified range and state the restriction explicitly.
[Abstract; conclusion] The claim that the initial BH qubit state can be retrieved from its imprint on the final radiation state is not demonstrated. The text describes unitary evolution and gives reduced density matrices, but no reconstruction map, fidelity, or decoding procedure is exhibited. Retrieval from radiation alone is nontrivial, since knowing the full final state makes inversion trivial by unitarity. Please specify the subsystem and the operation that recovers the initial BH state, or rephrase the claim as a possibility statement consistent with unitary evolution.

minor comments (3)

[Sec. 2, p. 2] The model description says 'a 6-qubit system for BH, a 1-qubit for ER, and a 6-qubit for ER.' The second 'ER' should presumably be 'JR'. Please correct.
[Sec. 2, θ] The physical interpretation of θ as a scrambling parameter is clear from Ref. [21], but the text should state whether the same θ is applied to all infalling-matter qubits and how the ER degrees of freedom are initialized, since these choices affect the entanglement structure.
[Throughout] The word 'firewall' is used for a negativity-based criterion. Please state explicitly that this is a toy-model diagnostic, not a literal AMPS firewall, to avoid overinterpretation.

Circularity Check

0 steps flagged

No significant circularity: all results are exact analytical consequences of the stated toy-model unitary circuit; no fitted parameter is renamed as a prediction, and no self-citation chain is load-bearing.

full rationale

The derivation chain is self-contained. The only inputs are (i) the six-qubit circuit from [1], (ii) the standard Boltzmann-weighted two-mode squeezed-state relation tan γ = exp(−4πMω), and (iii) a one-parameter scrambling unitary U(θ). No parameter is fitted to data; θ is a scanned knob and γ is a fixed physical input. The firewall diagnostic is operational: negativity of the BH–JR reduced state is computed at each step, and the paper's Mω-window and the Eq. (53) θ-range follow from the analytical reduced density matrices in Appendices A–B. The information-retrieval statement is a consequence of the unitarity of the closed circuit, not an input disguised as an output. The two-level Fock truncation |0⟩,|1⟩ is a modeling approximation with uncontrolled weight e^{−16πMω} at small Mω; that is a correctness/robustness risk, not a circularity, because the discarded levels are not used to define any fitted quantity or to choose θ. No load-bearing self-uniqueness theorem or ansatz-by-citation appears.

Axiom & Free-Parameter Ledger

1 free parameters · 5 axioms · 0 invented entities

The central claims rest on [1]'s circuit architecture and on a two-level truncation of the squeezed Hawking state. θ is the only hand-set knob in the visible text; γ is a physical input from prior literature (Hawking Boltzmann factor). No new entities are introduced and no data are fitted.

free parameters (1)

θ (scrambling angle) = 0 ≤ θ ≤ π/2 (θ=0 no scrambling; θ=π/2 maximal scrambling)
Hand-chosen model parameter that interpolates the initial black-hole state between the [1] configuration and a maximally scrambled one. Not fitted to data, but introduced as a free knob; all new claims are functions of θ.

axioms (5)

domain assumption Two-level truncation of the Hawking squeezed state: keep only |0⟩ and |1⟩ with tan γ = exp(−4πMω).
Stated in Sec. 2 as 'the two dominant contributions.' Uncontrolled where e^(−4πMω) ~ O(1), i.e., the small-Mω regime where the firewall claim lives.
domain assumption The circuit of [1] with CNOT-U gates and the before/after-Page-time causal ordering is a faithful toy representation of black-hole evaporation.
Adopted from [1] and [12,22]; the firewall conclusions inherit this model dependence.
domain assumption Firewall formation is identified with the vanishing of BH–JR entanglement measures (negativity) at late circuit steps.
The paper equates 'entanglement between BH and JR not sustained' with firewall emergence, following [1]'s criterion; this is a definitional rather than derived identification.
standard math Standard entanglement toolkit: von Neumann entropy, mutual information, negativity in finite-dimensional systems.
Used without derivation per [13–16]; standard and uncontroversial.
standard math Unitarity of the overall circuit dynamics.
Gate-based model; unitarity guarantees reversibility, which underlies the information-retrieval claim.

pith-pipeline@v1.3.0-alltime-deepseek · 131826 in / 16104 out tokens · 154222 ms · 2026-08-02T10:40:26.159541+00:00 · methodology

0 comments

read the original abstract

We reexamine the quantum circuit model of black hole evaporation proposed in a previous work (Class. Quantum Grav. 35, 235013, 2018) [1]. This tripartite model incorporates the following systems: black hole ($\mathbf{BH}$), just radiation ($\mathbf{JR}$), and early radiation ($\mathbf{ER}$). We apply a scrambling unitary matrix with a single parameter $\theta$ to the ground state of the qubits in infalling matter toward a black hole in order to generate initial qubit states of the black hole that are more general than those in [1]. Specifically, the scrambling unitary matrix reduces to no scrambling and maximum scrambling when $\theta=0$ and $\theta=\pi/2$, respectively. Our aim is to explore the role of quantum monogamy in the firewall formation between the black hole and radiation. In this model, entanglement and firewall formation depend on the black hole mass $M$ and the frequency of Hawking radiation $\omega$. For the initial state with $\theta=\pi/2$, a firewall emerges at an earlier stage of the evolution than with $\theta=0$. We also find that a firewall structure emerges between $\mathbf{BH}$ and $\mathbf{JR}$, and that the information is carried away by radiation for all values of $M\omega$, provided that $\theta$ lies within a certain analytically determined range. Following unitary gate dynamics, the initial black hole qubit state can be retrieved from its imprint on the final radiation state, which was originally hidden behind the black hole's horizon. These results may provide insight into the properties of multipartite entanglement due to the different initial states in the evolution of a quantum circuit model for black hole evaporation.

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

23 extracted references · 2 linked inside Pith

[1]

Quantum Grav

Tomoro Tokusumi, Akira Matsumura and Yasusada Nambu, Quantum circuit model of black hole evapora- tion, Class. Quantum Grav. 35, 235013 (2018)

2018
[2]

S. W. Hawking, Particle creation by black holes, Comm. Math. Phys. 43, 199 (1975)

1975
[3]

S. W. Hawking, Breakdown of predictability in gravita- tional collapse, Phys. Rev. D 14, 2460 (1976)

1976
[4]

D. N. Page, Information in black hole radiation, Phys. Rev. Lett. 71, (1993) 3743

1993
[5]

D. N. Page, Average entropy of a subsystem, Phys. Rev. Lett. 71, (1993) 1291

1993
[6]

Susskind, L

L. Susskind, L. Thorlacius, and J. Uglum, The Stretched horizon and black hole complementarity, Phys. Rev. D 48, 3743 (1993)

1993
[7]

Almheiri, D

A. Almheiri, D. Marolf, J. Polchinski, and J. Sully, Black Holes: Complementarity or Firewalls?, JHEP 02, 062 (2013)

2013
[8]

S. Luo, H. Stoltenberg, and A. Albrecht, Multipartite Entanglement and Firewalls, Phys. Rev. D 95, 064039 (2017)

2017
[9]

Hwang, D

J. Hwang, D. S. Lee, D. Nho, J. Oh, H. Park, D.-h. Yeom, and H. Zoe, Page curves for tripartite systems, Class. Quant. Grav. 34, 145004 (2017)

2017
[10]

S. G. Avery, Qubit models of black hole evaporation, JHEP 01 (2013) 176

2013
[11]

Osuga and D

K. Osuga and D. N. Page, Qubit Transport Model for Unitary Black Hole Evaporation without Firewalls, Phys. Rev. D 97, 066023 (2018)

2018
[12]

Broda, Causal unitary qubit model of black hole evap- oration, Phys

B. Broda, Causal unitary qubit model of black hole evap- oration, Phys. Lett. B 820 (2021) 136564

2021
[13]

Peres, Separability Criterion for Density Matrices, Phys

A. Peres, Separability Criterion for Density Matrices, Phys. Rev. Lett. 77, 1413 (1996)

1996
[14]

Wei-Can Syu, Da-Shin Lee, and Chen-Pin Yeh, Entan- glement of quantum oscillators coupled to different heat baths, J. Phys. B 54, 055501 (2021)

2021
[15]

Vidal and R

G. Vidal and R. F. Werner, Computable measure of en- tanglement, Phys. Rev. A 65, 032314 (2002)

2002
[16]

He and G

H. He and G. Vidal, Disentangling theorem and monogamy for entanglement negativity, Phys. Rev. A 91, 012339 (2015)

2015
[17]

Hayden and J

P. Hayden and J. Preskill, Black holes as mirrors: quan- tum information in random subsystems, JHEP, 09 (2007) 120

2007
[18]

K. A. Landsman, C. Figgatt, T. Schuster, N. M. Linke, B. Yoshida, N. Y. Yao, and C. Monroe, Verified Quantum Information Scrambling, Nature, 567 (2019) 61. 19

2019
[19]

Beni Yoshida and Alexei Kitaev, Efficient decoding for the Hayden-Preskill protocol, arXiv:1710.03363 (2017)

Pith/arXiv arXiv 2017
[20]

Yao, Disentangling scrambling and decoherence via quantum teleportation, arXiv:1803.10772 (2018)

Beni Yoshida and Norman Y. Yao, Disentangling scrambling and decoherence via quantum teleportation, arXiv:1803.10772 (2018)

Pith/arXiv arXiv 2018
[21]

MuSeong Kim, Mi-Ra Hwang, Eylee Jung, DaeKil Park, Scrambling and Quantum Teleportation, Quant. Inf. Proc. 22, 176 (2023)

2023
[22]

Broda, Unitary toy qubit transport model for black hole evaporation, Eur

B. Broda, Unitary toy qubit transport model for black hole evaporation, Eur. Phys. J. C 80 (2020) 5, 418

2020
[23]

Horodecki, P

M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996)

1996

[1] [1]

Quantum Grav

Tomoro Tokusumi, Akira Matsumura and Yasusada Nambu, Quantum circuit model of black hole evapora- tion, Class. Quantum Grav. 35, 235013 (2018)

2018

[2] [2]

S. W. Hawking, Particle creation by black holes, Comm. Math. Phys. 43, 199 (1975)

1975

[3] [3]

S. W. Hawking, Breakdown of predictability in gravita- tional collapse, Phys. Rev. D 14, 2460 (1976)

1976

[4] [4]

D. N. Page, Information in black hole radiation, Phys. Rev. Lett. 71, (1993) 3743

1993

[5] [5]

D. N. Page, Average entropy of a subsystem, Phys. Rev. Lett. 71, (1993) 1291

1993

[6] [6]

Susskind, L

L. Susskind, L. Thorlacius, and J. Uglum, The Stretched horizon and black hole complementarity, Phys. Rev. D 48, 3743 (1993)

1993

[7] [7]

Almheiri, D

A. Almheiri, D. Marolf, J. Polchinski, and J. Sully, Black Holes: Complementarity or Firewalls?, JHEP 02, 062 (2013)

2013

[8] [8]

S. Luo, H. Stoltenberg, and A. Albrecht, Multipartite Entanglement and Firewalls, Phys. Rev. D 95, 064039 (2017)

2017

[9] [9]

Hwang, D

J. Hwang, D. S. Lee, D. Nho, J. Oh, H. Park, D.-h. Yeom, and H. Zoe, Page curves for tripartite systems, Class. Quant. Grav. 34, 145004 (2017)

2017

[10] [10]

S. G. Avery, Qubit models of black hole evaporation, JHEP 01 (2013) 176

2013

[11] [11]

Osuga and D

K. Osuga and D. N. Page, Qubit Transport Model for Unitary Black Hole Evaporation without Firewalls, Phys. Rev. D 97, 066023 (2018)

2018

[12] [12]

Broda, Causal unitary qubit model of black hole evap- oration, Phys

B. Broda, Causal unitary qubit model of black hole evap- oration, Phys. Lett. B 820 (2021) 136564

2021

[13] [13]

Peres, Separability Criterion for Density Matrices, Phys

A. Peres, Separability Criterion for Density Matrices, Phys. Rev. Lett. 77, 1413 (1996)

1996

[14] [14]

Wei-Can Syu, Da-Shin Lee, and Chen-Pin Yeh, Entan- glement of quantum oscillators coupled to different heat baths, J. Phys. B 54, 055501 (2021)

2021

[15] [15]

Vidal and R

G. Vidal and R. F. Werner, Computable measure of en- tanglement, Phys. Rev. A 65, 032314 (2002)

2002

[16] [16]

He and G

H. He and G. Vidal, Disentangling theorem and monogamy for entanglement negativity, Phys. Rev. A 91, 012339 (2015)

2015

[17] [17]

Hayden and J

P. Hayden and J. Preskill, Black holes as mirrors: quan- tum information in random subsystems, JHEP, 09 (2007) 120

2007

[18] [18]

K. A. Landsman, C. Figgatt, T. Schuster, N. M. Linke, B. Yoshida, N. Y. Yao, and C. Monroe, Verified Quantum Information Scrambling, Nature, 567 (2019) 61. 19

2019

[19] [19]

Beni Yoshida and Alexei Kitaev, Efficient decoding for the Hayden-Preskill protocol, arXiv:1710.03363 (2017)

Pith/arXiv arXiv 2017

[20] [20]

Yao, Disentangling scrambling and decoherence via quantum teleportation, arXiv:1803.10772 (2018)

Beni Yoshida and Norman Y. Yao, Disentangling scrambling and decoherence via quantum teleportation, arXiv:1803.10772 (2018)

Pith/arXiv arXiv 2018

[21] [21]

MuSeong Kim, Mi-Ra Hwang, Eylee Jung, DaeKil Park, Scrambling and Quantum Teleportation, Quant. Inf. Proc. 22, 176 (2023)

2023

[22] [22]

Broda, Unitary toy qubit transport model for black hole evaporation, Eur

B. Broda, Unitary toy qubit transport model for black hole evaporation, Eur. Phys. J. C 80 (2020) 5, 418

2020

[23] [23]

Horodecki, P

M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996)

1996