Two-dimensional solitons in saturable media with a quasi-one-dimensional lattice potential

Boris A. Malomed; Thawatchai Mayteevarunyoo

REVIEW

Two-dimensional solitons in saturable media with a quasi-one-dimensional lattice potential

Thawatchai Mayteevarunyoo, Boris A. Malomed

Pattern Formation and Solitons · Submitted 25 Feb 2006 · 16 references

Pith's verdict: Sound numerical study of Q1D lattice solitons in saturable media, but the 'stable in the entire existence region' claim is overreach given the finite-time evidence.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDnlin/0602052 (version 1)

Pith IDTE3EO3EM, Pith's permanent record of this paper and its review

CitationBoris A. Malomed, Thawatchai Mayteevarunyoo. "Two-dimensional solitons in saturable media with a quasi-one-dimensional lattice potential." arXiv:nlin/0602052 (2006). https://pith.science/paper/TE3EO3EM

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Colleague,

The paper combines a quasi-1D lattice with saturable nonlinearity in 2D and reports soliton families, thresholds, mobility, and collisions. That combination is new, and the work is a reasonable extension of the authors' own program. The numerical method is standard, and the main qualitative findings—finite minimum power, side lobes for broad solitons, a critical boost for cross-lattice motion, and inelastic collisions—are plausible and worth knowing.

The real soft spot is the stability claim. The paper says solitons are stable in the entire existence region, but the evidence is direct simulation over z ≤ 200 with random noise, with only two examples plotted (I0=5, Ip=56.76 and Ip=1.1). The authors themselves note that the VK criterion is not a full guarantee and even cite a counter-example in a related model. In the near-threshold case the peak intensity slowly grows over the whole run; no saturation is shown, so a weak instability with growth rate <1/200 would be invisible. A linear-stability eigenvalue calculation or a convergence study in z would have settled this. Without it, the absolute claim outruns the evidence.

There are also minor gaps: no error bars on the threshold values or critical boosts, and no code or data. That is typical for the era and not disqualifying.

The threshold and existence results, plus the collision phenomenology, look solid enough to stand on their own. The paper is honest about the VK limitation in the text, which makes the overstatement in the abstract and conclusion more of a framing problem than a fabrication.

For you: if you work on lattice solitons or photorefractive spatial solitons, this is a citable numerical study. I would send it to peer review; a referee should ask for a linear-stability analysis or at least a more systematic finite-time convergence check before the full-stability claim is accepted.

Referee report

Significance

If the claims hold, this is a useful numerical contribution to the physics of photorefractive lattice solitons, a topic with active experimental relevance (the model and lattice geometry match experiments in Ref. [1,2]). The paper provides concrete numerical evidence for threshold behavior, the shape transition from single-lobe to multi-lobe structures, and a mobility limit for solitons across a lattice, along with a first study of collisions in this model. The threshold curves in Figs. 2 and 3 and the mobility data in Fig. 5 are plausible and will be of interest to the soliton community. However, the most strongly worded claim—that solitons are stable in the entire existence region—is backed only by two finite-time simulations, and the authors themselves note that the VK criterion is not a complete stability guarantee. The paper also contains a specific technical error in the Galilean transformation used to construct moving solitons along y, which undermines part of the collision analysis. These issues need to be addressed before the paper can be recommended for publication.

Major objections (2)

[II. Two-dimensional solitons and their stability] The central claim that all solitons are stable in the entire existence region is not supported by the presented evidence. Stability is tested by direct simulation for only two parameter points (Fig. 4: I0=5, Ip=56.76 and Ip=1.1), integrated to z=200. The authors themselves state in Section II that the VK criterion does not provide full stability and that its applicability is not proven in this context, citing a counterexample [14]. In Fig. 4(b) the peak intensity grows monotonically over the whole run, and no saturation or return to a stationary value is demonstrated; a weak exponential instability with a growth rate below roughly 1/200 would be indistinguishable from this behavior. To support the 'entire existence region' statement, a linear-stability analysis (eigenvalues of the perturbation operator) or a systematic scan over I0, d, and µ with convergence in z is required. Alternatively, the claim should be narrowed to the tested parameter points. As it stands, the statement 'they are stable in the entire existence region' is a load-bearing overclaim.
[II.A, Eq. (9)] Equation (9) is not a Galilean transformation of Eq. (1). A genuine boost along the free direction y would read U(x, y - 2Q_y z, z) exp(i Q_y y - i Q_y^2 z). The formula as written shifts the argument in y but puts a phase linear in x. Substituting this into Eq. (1) and passing to the moving frame v=y-2Q_yz leaves a residual term 2 i Q_y (∂_x - ∂_v)U that does not vanish for a generic soliton. Consequently, the initial conditions used for the y-directed and diagonal moving-soliton simulations are not exact moving soliton solutions, and the measured critical boosts and collision outcomes in Section II.A may be affected. Please correct the formula (or state explicitly that a shear-type transformation is intended and justify its validity), and re-examine the affected numerical results.

Minor objections (5)

[II] The statement that the soliton families 'belong to the semi-infinite bandgap' is asserted without showing the linear spectrum of Eq. (1); please display the spectrum or cite a reference that establishes this for the present lattice potential.
[II, Eqs. (4)-(7)] The numerical method is described only schematically; for reproducibility, please provide the computational grid (number of points, domain size, resolution), the iteration stopping criterion, and some evidence of convergence (e.g., residual behavior or a grid-independence test).
[Fig. 5] The mobility results Qmax(Ip) are shown only for I0=5 and 25, while Fig. 2 presents soliton families also for I0=10 and 15; please clarify whether the same qualitative behavior holds for intermediate lattice strengths or restrict the claim.
[II.A] The critical boost for collision, |Q|=0.05, is reported only for I0=5; please state whether this value is generic or an example, and indicate how the collision outcome was determined (e.g., threshold resolution in Q).
[Global] There are several typographical and readability issues, including 'quasi-o ne-dimensional' in the title, the use of 'power' and 'intensity' for Ip and P, and inconsistent notation for Fourier variables (e.g., 'k2x' should read k_x^2). Please edit throughout.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Two-dimensional solitons in a saturable photorefractive medium with a quasi-1D lattice exist only above finite power thresholds, are stable throughout their existence region, move up to a critical boost, and collide inelastically.

The reading

The paper studies a two-dimensional model of light propagation in a photorefractive crystal whose refractive index is modulated by a quasi-one-dimensional periodic lattice, under saturable nonlinearity. It establishes that two-dimensional solitons exist in this setting only above finite minimum values of peak and total power, with narrow solitons single-peaked and broader ones carrying side lobes along the lattice direction. Direct numerical simulations of the evolution equation with random noise indicate that the soliton families are stable across their entire existence region, although the paper notes the Vakhitov-Kolokolov criterion does not by itself prove full stability. Because the lattice leaves one transverse direction free, the solitons can be tilted to move across the lattice up to a critical boost, and their collisions are always inelastic: fast collisions destroy them, slow ones merge them into a single stable pulse. The model is directly realizable with a single pair of counterpropagating beams, which is why the results matter for experiments.

What carries the argument

The central objects are the stationary 2D soliton solutions $u(x,y)$ of Eq. (3), computed with the modified Fourier-space iteration (7) whose ratio factors $\alpha_n/\beta_n$ enforce the fixed point; the quasi-1D lattice breaks translation symmetry along $x$ but leaves $y$ free, which is what permits tilted motion. Stability is interpreted through the Vakhitov-Kolokolov criterion via the slope of $P(\mu)$ and then checked by direct propagation, while moving solitons are generated by the Galilean transformation (9) along $y$ and by the exponential boost (10) along $x$.

Extended reading notes

Core claim

In model (1), combining a photorefractive saturable nonlinearity with the quasi-1D lattice potential $I_0\cos^2(\pi x/d)$ yields families of 2D solitons in the semi-infinite bandgap. For each lattice strength $I_0$ the families terminate at finite thresholds $P_{\min}$ and $I_{p,\min}$; the $P(\mu)$ curves have negative slope, consistent with VK stability, and direct runs with random noise show stable propagation for all tested family members. Solitons can be boosted across the lattice up to a critical $Q_{\max}$, which is approximately isotropic in the tested case, and collisions below that threshold merge into a stable vibrating pulse, while above it the pulses are destroyed.

Load-bearing premise

The all-stable conclusion rests on direct simulations over a finite propagation distance rather than on a proof, a limitation the paper itself acknowledges when it notes the VK criterion cannot guarantee full stability.

Editorial extensions

If this is right

Any experimental realization must exceed the finite thresholds $P_{\min}$ and $I_{p,\min}$; below them no 2D soliton can form for a given lattice strength.
Across-lattice motion is bounded: boosts above $Q_{\max}$ disperse the input into small-amplitude waves, so steering is only possible within that window.
Head-on or offset collisions are always inelastic: fast collisions destroy both pulses, while slow collisions fuse them into a single vibrating but stable pulse, regardless of collision direction.
The absence of collapse at merger, unlike the cubic-nonlinearity model, means the saturable lattice system offers a collapse-free route to multi-soliton interaction.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the near-isotropy of the critical boost $|Q|=0.05$ holds beyond the tested parameters, then the effective lattice potential acts almost equally on $x$- and $y$-directed motion despite the strong anisotropy of the lattice itself; this could be tested by mapping $Q_{\max}$ across $I_0$ and $d$.
The persistent vibrations of merger products suggest collision outcomes might be controllable by relative phase or impact parameter, possibly producing bound soliton trains; the paper does not test phase-controlled collisions.
A natural three-dimensional analogue, saturable nonlinearity with a quasi-2D lattice, would likely suppress collapse the same way, but the paper does not address 3D.

Original abstract

We study families of solitons in a two-dimensional (2D) model of the light transmission through a photorefractive medium equipped with a (quasi-)one-dimensional photonic lattice. The soliton families are bounded from below by finite minimum values of the peak and total power. Narrow solitons have a single maximum, while broader ones feature side lobes. Stability of the solitons is checked by direct simulations. The solitons can be set in motion across the lattice (actually, made tilted in the spatial domain), provided that the respective boost parameter does not exceed a critical value. Collisions between moving solitons are studied too. Collisions destroy the solitons, unless their velocities are sufficiently small. In the latter case, the colliding solitons merge into a single stable pulse.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,047 input tokens.

What the abstract itself states, as probabilities: claims first 0.02 · states limitation 0.71 · human participants 0.01 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.24 · compares with prior results 0.05 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

We study families of solitons in a two-dimensional (2D) model of the light transmission through a photorefractive medium equipped with a (quasi-)one-dimensional photonic lattice. The soliton families are bounded from below by finite minimum values of the peak and total power. Narrow solitons have a single maximum, while broader ones feature side lobes. Stability of the solitons is checked by direct simulations. The solitons can be set in motion across the lattice (actually, made tilted in the spatial domain), provided that the respective boost parameter does not exceed a critical value. Collisions between moving solitons are studied too. Collisions destroy the solitons, unless their velocities are sufficiently small. In the latter case, the colliding solitons merge into a single stable pulse.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	We study families of solitons in a two-dimensional (2D) model of the light transmission through a photorefractive medium equipped with a (quasi-)one-dimensional photonic lattice.	contribution 0.98	0.90	0.05	asserted 0.94
s1	The soliton families are bounded from below by finite minimum values of the peak and total power.	result 1.00	0.85	0.15	asserted 1.00
s2	Narrow solitons have a single maximum, while broader ones feature side lobes.	result 1.00	0.84	0.08	asserted 1.00
s3	Stability of the solitons is checked by direct simulations.	method 0.99	0.93	0.04	asserted 0.99
s4	The solitons can be set in motion across the lattice (actually, made tilted in the spatial domain), provided that the respective boost parameter does not exceed a critical value.	result 0.95	0.85	0.11	asserted 0.55
s5	Collisions between moving solitons are studied too.	contribution 0.98	0.87	0.04	asserted 0.93
s6	Collisions destroy the solitons, unless their velocities are sufficiently small.	result 1.00	0.89	0.16	asserted 1.00
s7	In the latter case, the colliding solitons merge into a single stable pulse.	result 1.00	0.89	0.16	asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.93 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** Typical examples of two-dimensional solitons found

**Figure 2.** The total intensity (power) (a) and peak intensity (b

**Figure 3.** The minimum peak intensity, Ip,min, and total power, Pmin, necessary for the existence of the two-dimensional solitons, vs. the strength, I0, of the underlying quasi-one-dimensional lattice

**Figure 4.** (a) Stable evolution of a high-intensity soliton, wi

**Figure 5.** The maximum value of the boost parameter,

**Figure 6.** A typical example of destructive collisions between

**Figure 7.** FIG. 7: Merger of colliding solitons which were set in motion

**Figure 8.** FIG. 8: Merger as a result of head-on collisions between slow

**Figure 9.** FIG. 9: Merger of broad (multi-lobed) solitons caused by the

Figures reproduced from arXiv:nlin/0602052 by the authors.

References

Works cited (16)

[14]

B. B. Baizakov, B. A. Malomed, and M. Salerno, Europhys. Lett. 63, 642 (2003)

View on Pith 2003
[1]

N. K. Efremidis, S. Sears, D. N. Christodoulides, J. W. Fl eischer, and M. Segev, Phys. Rev. E 66, 046602 (2002); N. K. Efremidis, J. Hudock, D. N. Christodoulides, J. W. Fleische r, O. Cohen, and M. Segev, Phys. Rev. Lett. 91, 213906 (2003)

View on Pith 2002
[2]

J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. Chris todoulides, Nature 422, 147 (2003); J. W. Fleischer, G. Bartal, O. Cohen, O. Manela, M. Segev, J. Hudock, and D. N. Christodou lides, Phys. Rev. Lett. 92, 123904 (2004)

View on Pith 2003
[3]

D N. Neshev, J. Alexander, E. A. Ostrovskaya, Y. S. Kivsha r, H. Martin, I. Makasyuk, and Z. Chen, Phys. Rev. Lett. 92, 123903 (2004); J. W. Fleischer, G. Bartal, O. Cohen, O. Manel a, M. Segev, J. Hudock, and D. N. Christodoulides, Phys. Rev. Lett. 92, 123904 (2004)

View on Pith 2004
[4]

O. Manela, O. Cohen, G. Bartal, J. W. Fleischer, and M. Seg ev, Opt. Lett. 29, 2049 (2004)

View on Pith 2004
[5]

H. Sakaguchi and B. A. Malomed, Higher-order vortex solitons, multipoles, and supervorti ces on a square optical lattice , Europhys. Lett., in press

View on Pith
[6]

J. Yang, I. Makasyuk, A. Bezryadina, and Z. Chen, Stud. Ap pl. Math. 113, 389 (2004)

View on Pith 2004
[7]

Further, the mergers caused by head-on collisions in the diagona l or vertical ( y) direction are presented in Fig. 8. As seen from the ﬁgures, the pulse generated by the merger perf orms intrinsic vibrations, but remains stable. The above examples displayed collisions of single-lobe solitons. Collisions between their broader counterparts, which feature s...

View on Pith
[8]

J. Yang, I. Makasyuk, P. G. Kevrekidis, H. Martin, B. A. Ma lomed, D. J. Frantzeskakis, and Z. Chen, Phys. Rev. Lett. 94, 113902 (2005)

View on Pith 2005
[9]

J. W. Fleischer, G. Bartal, O. Cohen, T. Schwartz, O. Mane la, B. Freedman, M. Segev, H. Buljan, N. K. Efremidis, Opt. Exp. 13, 1780 (2005)

View on Pith 2005
[10]

B. B. Baizakov, B. A. Malomed and M. Salerno, Phys. Rev. A 70, 053613 (2004); in: Nonlinear Waves: Classical and Quantum Aspects, ed. by F. Kh. Abdullaev and V. V. Konotop, p. 61 (Kluwer Acade mic Publishers: Dordrecht, 2004; also available at http://rsphy2.anu.edu.au/˜asd124/Baizakov 2004 61 NonlinearWaves.pdf)

View on Pith 2004
[11]

D. Mihalache, D. Mazilu, F. Lederer, Y. V. Kartashov, L. -C. Crasovan, and L. Torner, Phys. Rev. E 70, 055603(R) (2004 )

View on Pith 2004
[12]

R. Fischer, D. Traeger, D. N. Neshev, A. A. Sukhorukov, W . Krolikowski, C. Denz, and Y. S. Kivshar, e-print physics/0509255 (2005)

View on Pith Pith review arXiv 2005
[13]

B. A. Malomed, T. Mayteevarunyoo, E. A. Ostrovskaya, an d Y. S. Kivshar, Phys. Rev. E 71, 056616 (2005)

View on Pith 2005
[15]

I. M. Merhasin, B. V. Gisin, R. Driben, and B. A. Malomed, Phys. Rev. E 71, 016613 (2005). 7 Fig. 1: Typical examples of two-dimensional solitons found from Eq. (3) with a moderately strong quasi-one-dimensiona l photonic lattice, corresponding to I0 = 5. Panels (a), (b), and (c) show the intensity distribution in the lattice, Iphl(x) = I0 cos2 (πx/d ), f...

View on Pith 2005
[16]

M. G. Vakhitov and A. A. Kolokolov, Izv. Vuz. Radioﬁz. 16, 1020 (1973) [in Russian; English translation: Sov. J. Radio phys. Quantum Electr. 16, 783 (1973)]; see also L. Berg´ e, Phys. Rep. 303, 260 (1998). 8 0 0.05 0.1 0.15 0.2 0.25 0.3 0.3510 2 10 3 10 4 10 5 10 6 µ P I0 = 5 I0 = 10 I0 = 15 (a) 0 0.05 0.1 0.15 0.2 0.25 0.3 0.3510 −1 10 0 10 1 10 2 10 3 ...

View on Pith 1973

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [14]

B. B. Baizakov, B. A. Malomed, and M. Salerno, Europhys. Lett. 63, 642 (2003)

View on Pith 2003

[2] [1]

N. K. Efremidis, S. Sears, D. N. Christodoulides, J. W. Fl eischer, and M. Segev, Phys. Rev. E 66, 046602 (2002); N. K. Efremidis, J. Hudock, D. N. Christodoulides, J. W. Fleische r, O. Cohen, and M. Segev, Phys. Rev. Lett. 91, 213906 (2003)

View on Pith 2002

[3] [2]

J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. Chris todoulides, Nature 422, 147 (2003); J. W. Fleischer, G. Bartal, O. Cohen, O. Manela, M. Segev, J. Hudock, and D. N. Christodou lides, Phys. Rev. Lett. 92, 123904 (2004)

View on Pith 2003

[4] [3]

D N. Neshev, J. Alexander, E. A. Ostrovskaya, Y. S. Kivsha r, H. Martin, I. Makasyuk, and Z. Chen, Phys. Rev. Lett. 92, 123903 (2004); J. W. Fleischer, G. Bartal, O. Cohen, O. Manel a, M. Segev, J. Hudock, and D. N. Christodoulides, Phys. Rev. Lett. 92, 123904 (2004)

View on Pith 2004

[5] [4]

O. Manela, O. Cohen, G. Bartal, J. W. Fleischer, and M. Seg ev, Opt. Lett. 29, 2049 (2004)

View on Pith 2004

[6] [5]

H. Sakaguchi and B. A. Malomed, Higher-order vortex solitons, multipoles, and supervorti ces on a square optical lattice , Europhys. Lett., in press

View on Pith

[7] [6]

J. Yang, I. Makasyuk, A. Bezryadina, and Z. Chen, Stud. Ap pl. Math. 113, 389 (2004)

View on Pith 2004

[8] [7]

Further, the mergers caused by head-on collisions in the diagona l or vertical ( y) direction are presented in Fig. 8. As seen from the ﬁgures, the pulse generated by the merger perf orms intrinsic vibrations, but remains stable. The above examples displayed collisions of single-lobe solitons. Collisions between their broader counterparts, which feature s...

View on Pith

[9] [8]

J. Yang, I. Makasyuk, P. G. Kevrekidis, H. Martin, B. A. Ma lomed, D. J. Frantzeskakis, and Z. Chen, Phys. Rev. Lett. 94, 113902 (2005)

View on Pith 2005

[10] [9]

J. W. Fleischer, G. Bartal, O. Cohen, T. Schwartz, O. Mane la, B. Freedman, M. Segev, H. Buljan, N. K. Efremidis, Opt. Exp. 13, 1780 (2005)

View on Pith 2005

[11] [10]

B. B. Baizakov, B. A. Malomed and M. Salerno, Phys. Rev. A 70, 053613 (2004); in: Nonlinear Waves: Classical and Quantum Aspects, ed. by F. Kh. Abdullaev and V. V. Konotop, p. 61 (Kluwer Acade mic Publishers: Dordrecht, 2004; also available at http://rsphy2.anu.edu.au/˜asd124/Baizakov 2004 61 NonlinearWaves.pdf)

View on Pith 2004

[12] [11]

D. Mihalache, D. Mazilu, F. Lederer, Y. V. Kartashov, L. -C. Crasovan, and L. Torner, Phys. Rev. E 70, 055603(R) (2004 )

View on Pith 2004

[13] [12]

R. Fischer, D. Traeger, D. N. Neshev, A. A. Sukhorukov, W . Krolikowski, C. Denz, and Y. S. Kivshar, e-print physics/0509255 (2005)

View on Pith Pith review arXiv 2005

[14] [13]

B. A. Malomed, T. Mayteevarunyoo, E. A. Ostrovskaya, an d Y. S. Kivshar, Phys. Rev. E 71, 056616 (2005)

View on Pith 2005

[15] [15]

I. M. Merhasin, B. V. Gisin, R. Driben, and B. A. Malomed, Phys. Rev. E 71, 016613 (2005). 7 Fig. 1: Typical examples of two-dimensional solitons found from Eq. (3) with a moderately strong quasi-one-dimensiona l photonic lattice, corresponding to I0 = 5. Panels (a), (b), and (c) show the intensity distribution in the lattice, Iphl(x) = I0 cos2 (πx/d ), f...

View on Pith 2005

[16] [16]

M. G. Vakhitov and A. A. Kolokolov, Izv. Vuz. Radioﬁz. 16, 1020 (1973) [in Russian; English translation: Sov. J. Radio phys. Quantum Electr. 16, 783 (1973)]; see also L. Berg´ e, Phys. Rep. 303, 260 (1998). 8 0 0.05 0.1 0.15 0.2 0.25 0.3 0.3510 2 10 3 10 4 10 5 10 6 µ P I0 = 5 I0 = 10 I0 = 15 (a) 0 0.05 0.1 0.15 0.2 0.25 0.3 0.3510 −1 10 0 10 1 10 2 10 3 ...

View on Pith 1973