pith. sign in

arxiv: 2502.14207 · v2 · submitted 2025-02-20 · 🪐 quant-ph · cond-mat.mes-hall

Quantum stick-slip motion in nanoscaled friction

Pith reviewed 2026-05-23 03:06 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.mes-hall
keywords quantum frictionPrandtl-Tomlinson modelLandau-Zener tunnelingstick-slip motionnanoscaled frictiondensity matrixfrictional dissipationMarkov approximation
0
0 comments X

The pith

Landau-Zener tunneling reduces frictional dissipation in the quantum Prandtl-Tomlinson model compared to classical motion.

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

This paper develops the quantum version of the Prandtl-Tomlinson model for friction, in which a nanoparticle moves over a periodic potential while dissipating energy to a heat bath. The system's evolution is described by the Liouville-von Neumann equation for the density matrix under the Markov approximation. By contrasting quantum and classical behaviors, the authors show that Landau-Zener tunneling between energy levels lowers the frictional dissipation. The study maps how velocity, interaction strength, and temperature influence the frictional force, providing a framework for understanding quantum effects in nanoscale friction experiments.

Core claim

In the quantum mechanical Prandtl-Tomlinson model, the Landau-Zener tunneling is a key factor in the overall reduction of the frictional dissipation when compared to the classical motion in which such tunneling is absent. Other regimes of motion controlled by the corrugation parameter are also identified.

What carries the argument

Landau-Zener tunneling in the quantum Prandtl-Tomlinson model with Markovian dissipation to a bath, which enables transitions that reduce energy loss.

If this is right

  • The frictional force decreases due to quantum tunneling effects.
  • Regimes of motion depend on the corrugation parameter and other system properties.
  • Interplay of velocity, interaction strength, and temperature controls the frictional behavior.
  • Guidelines are provided for interpreting experimental data on nanoscaled friction.

Where Pith is reading between the lines

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

  • This could imply that quantum corrections are necessary for accurate modeling of friction in devices operating at low temperatures or small scales.
  • Connections might exist to other quantum transport phenomena where tunneling affects dissipation.
  • Testing in systems with tunable corrugation could reveal the predicted motion regimes.

Load-bearing premise

The dissipation process is assumed to occur via release of heat to an external bath under the Markov approximation for the Liouville-von Neumann equation, which treats the coupling as memoryless.

What would settle it

An experiment measuring the frictional force in a nanosystem where Landau-Zener tunneling is expected but finding dissipation rates identical to the classical prediction without reduction would challenge the claim.

Figures

Figures reproduced from arXiv: 2502.14207 by Dai-Nam Le, Lilia M. Woods, Pablo Rodriguez-Lopez.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematics of a particle with mass [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The average energy [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Kinematic properties of the moving particle as function of [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Dynamic properties as a function of [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
read the original abstract

Friction in atomistic systems is usually described by the classical Prandtl-Tomlinson model suitable for capturing the dragging force of a nanoparticle in a periodic potential. Here we consider the quantum mechanical version of this model in which the dissipation is facilitated by releasing heat to an external bath reservoir. The time evolution of the system is captured with the Liouville-von Neumann equation through the density matrix of the system in the Markov approximation. We examine several kinetic and dissipative properties of the nanoparticle by delineating classical vs quantum mechanical effects. We find that the Landau-Zener tunneling is a key factor in the overall reduction of the frictional dissipation when compared to the classical motion in which such tunneling is absent. Other regimes of motion, controlled by the corrugation parameter and other properties, are also found. This in-depth study analyzes the interplay between velocity, strength of interaction, and temperature to control the frictional {force} and provide useful guidelines for experimental data interpretation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript extends the classical Prandtl-Tomlinson model to a quantum version of a nanoparticle in a periodic potential, with dissipation modeled via the Markovian Liouville-von Neumann equation coupled to a bath. It compares classical and quantum regimes, claiming that Landau-Zener tunneling is the dominant mechanism reducing frictional dissipation relative to the classical stick-slip case, while also exploring regimes controlled by corrugation parameter, velocity, interaction strength, and temperature.

Significance. If the central attribution to Landau-Zener tunneling is rigorously isolated and the comparison to the classical limit is placed on equal footing, the work could offer concrete guidelines for interpreting nanoscale friction experiments where quantum effects become relevant.

major comments (2)
  1. [Abstract] Abstract and results sections: the claim that 'Landau-Zener tunneling is a key factor in the overall reduction of the frictional dissipation' is not supported by any reported control (e.g., artificially suppressed LZ probability at fixed ħ via slower velocity or larger gap, or explicit correlation between computed LZ transition probability and observed force reduction). Without such isolation, the reduction cannot be attributed specifically to LZ transitions rather than other quantum or bath effects.
  2. Methods and results: the classical comparison is performed separately rather than recovered as the ħ→0 limit of the same open-system Lindblad/LvN dynamics; this leaves open whether the reported dissipation reduction survives a consistent semiclassical limit or arises from differences in how the bath is coupled.
minor comments (1)
  1. [Abstract] The abstract states that 'several kinetic and dissipative properties' are examined but does not list them; a brief enumeration would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and indicate where revisions will be made to strengthen the presentation and support for our claims.

read point-by-point responses
  1. Referee: [Abstract] Abstract and results sections: the claim that 'Landau-Zener tunneling is a key factor in the overall reduction of the frictional dissipation' is not supported by any reported control (e.g., artificially suppressed LZ probability at fixed ħ via slower velocity or larger gap, or explicit correlation between computed LZ transition probability and observed force reduction). Without such isolation, the reduction cannot be attributed specifically to LZ transitions rather than other quantum or bath effects.

    Authors: We agree that an explicit isolation of the Landau-Zener contribution would make the attribution more rigorous. Our existing parameter scans (velocity, corrugation) exhibit trends consistent with LZ expectations, but we did not compute direct correlations or perform controls such as gap variation at fixed ħ. In the revised manuscript we will add a dedicated analysis that extracts LZ transition probabilities from the dynamics and correlates them with the observed frictional force reduction; the abstract claim will be revised accordingly to reflect the new evidence. revision: yes

  2. Referee: Methods and results: the classical comparison is performed separately rather than recovered as the ħ→0 limit of the same open-system Lindblad/LvN dynamics; this leaves open whether the reported dissipation reduction survives a consistent semiclassical limit or arises from differences in how the bath is coupled.

    Authors: The referee correctly identifies that our classical results use the deterministic Prandtl-Tomlinson equations while the quantum results employ the open-system LvN equation. We will add, where numerically tractable, a direct ħ→0 limit study within the same Lindblad framework and compare the resulting dissipation to the separate classical runs. Any residual differences in bath coupling will be discussed explicitly; if the semiclassical limit is not fully recoverable we will state the limitation and clarify the modeling choices. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation rests on explicit solution of open quantum dynamics

full rationale

The paper solves the Markovian Liouville-von Neumann equation for the quantum Prandtl-Tomlinson model and compares the resulting frictional force and dissipation to the classical limit. The central claim that Landau-Zener tunneling contributes to reduced dissipation is presented as an outcome of those dynamics rather than a re-expression of fitted parameters or a self-citation chain. No equations are shown to reduce to their own inputs by construction, no uniqueness theorem is imported from the authors' prior work, and no ansatz is smuggled via citation. The Markov approximation and bath coupling are stated assumptions whose validity is external to the derivation itself. The comparison to classical motion is performed separately and does not rely on re-labeling a fitted quantity as a prediction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review supplies minimal detail on parameters or assumptions; the Markov approximation is the only explicitly invoked modeling choice.

free parameters (1)
  • corrugation parameter
    Controls regimes of motion as stated in the abstract; value not specified.
axioms (1)
  • domain assumption Markov approximation for bath coupling in the Liouville-von Neumann equation
    Invoked to describe dissipation to the external reservoir.

pith-pipeline@v0.9.0 · 5694 in / 1200 out tokens · 31082 ms · 2026-05-23T03:06:23.380115+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Tuning of quantum nanoscaled friction within the Prandtl-Tomlinson model

    quant-ph 2026-04 unverdicted novelty 5.0

    Within the Prandtl-Tomlinson framework, frictional dynamics in nanoscaled systems are controlled by corrugation and characteristic length ratio parameters, revealing multiple motion types and the importance of quantum...

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages · cited by 1 Pith paper

  1. [6]

    Wolloch, G

    M. Wolloch, G. Levita, P. Restuccia, and M. C. Righi, Interfacial Charge Density and Its Connection to Adhesion and Frictional Forces, Phys. Rev. Lett. 121, 026804 (2018)

  2. [22]

    Meyer, T

    E. Meyer, T. Gyalog, R. M. Overney, and K. Dransfeld, Nanoscience: Friction and Rheology on the Nanometer Scale (World Scientific, 1998). 12

  3. [42]

    S. M. Rytov, Y. A. Kravtsov, and V. I. Tatarskii, Principles of Statistical Radiophysics 3 (Springer, Berlin, Germany, 1989)

  4. [49]

    Wilkie and M

    J. Wilkie and M. C ¸ etinba¸ s, Variable-stepsize runge–kutta methods for stochastic schr¨ odinger equations, Physics Letters A 337, 166 (2005). S-1 Supplemental materials for: Quantum stick-slip motion in nanoscaled friction Dai-Nam Le 1,∗ , Pablo Rodriguez-Lopez 2,† , Lilia M. Woods 1,‡ 1Department of Physics, University of South Florida, Tampa, Florida...

  5. [50]

    Quantum theory of open system in weak coupling and Born-Markov approximations S-7

  6. [51]

    Caldeira–Leggett model of thermal harmonic bath S-10

  7. [52]

    A r 2πddc a2 e− d dc − 3πC6 8d5a # + 2

    Choice of F(ˆx) and Matrix representations in the basis set of moving harmonic oscillator n(0)(t) S-12 S-XI. Classical dynamical evolution with enviroment: the case of Caldeira-Leggett model S-13 S-XII. Geometry phase S-15 S-XIII. Tunneling and slipping of quantum quasi-static eigenstates and classical trajectory S-16 S-XIV. Animation of quantum and class...

  8. [53]

    (S-32) In (S-31), dimensionless factor λ is introduced to manipulate the magnitude of the interaction from non-interacting case λ = 0 to the full interacting case λ = 1

    Quantum theory of open system in weak coupling and Born-Markov approximations Let us consider the general Hamiltonian describing a system ( S) interacting with a surrounding environment i.e bath (B) [25] ˆH(t) = ˆHS(t) ⊗ ˆ1B + ˆ1S ⊗ ˆHB + λ ˆHSB(t), (S-31) where the interacting Hamiltonian ˆHSB(t) generally takes the following form ˆHSB(t) = X α Aα(t) ⊗ B...

  9. [54]

    Caldeira–Leggett model of thermal harmonic bath The Caldeira–Leggett or harmonic bath model is commonly used to demonstrate the quantum dissipation [33, 34] ˆHint+B = X i " ˆp2 i 2mi + 1 2 miω2 i ˆxi − ci miω2 i F (ˆx) 2# , (S-56) where the choice of F (ˆx) strongly depends on the global sysmetry of the system ’s Hamiltonian ˆHS. The bath is characterized...

  10. [55]

    xhomo i (t) − X i ci miω2 i sin 2πℓ a x(0) cos ωit #) cos 2πℓ a x(t) −4π2ℓ2 a2 cos 2πℓ a x(t) tZ 0

    Choice of F(ˆx) and Matrix representations in the basis set of moving harmonic oscillator n(0)(t) E The Hamiltonian of the system ˆHS can be re-written in another form, ˆHS(t) = −1 2 ∂2 ∂x2 + 1 2 x − vt 2 + 1 2 u0 − 1 2 u0 cos 2πℓ a x , (S-76) in order to illustrate its periodicity of 2πℓ a x . Then it is obvious to choose F(ˆx) as a smooth periodic funct...

  11. [56]

    A. I. Volokitin and B. N. J. Persson, Near-field radiative heat transfer and noncontact friction, Rev. Mod. Phys. 79, 1291 (2007)

  12. [57]

    Vakis, V

    A. Vakis, V. Yastrebov, J. Scheibert, L. Nicola, D. Dini, C. Minfray, A. Almqvist, M. Paggi, S. Lee, G. Limbert, J. Molinari, G. Anciaux, R. Aghababaei, S. Echeverri Restrepo, A. Papangelo, A. Cammarata, P. Nicolini, C. Putignano, G. Carbone, S. Stupkiewicz, J. Lengiewicz, G. Costagliola, F. Bosia, R. Guarino, N. Pugno, M. M¨ user, and M. Ciavarella, Mode...

  13. [58]

    Wolloch, G

    M. Wolloch, G. Losi, M. Ferrario, and M. C. Righi, High-throughput screening of the static friction and ideal cleavage strength of solid interfaces, Sci. Rep. 9, 17062 (2019)

  14. [59]

    Cahangirov, C

    S. Cahangirov, C. Ataca, M. Topsakal, H. Sahin, and S. Ciraci, Frictional Figures of Merit for Single Layered Nanostruc- tures, Phys. Rev. Lett. 108, 126103 (2012)

  15. [60]

    R. K. Barik and L. M. Woods, Frictional Properties of Two-Dimensional Materials: Data-Driven Machine Learning Pre- dictive Modeling, ACS Applied Materials & Interfaces 16, 40149 (2024), https://doi.org/10.1021/acsami.4c05532

  16. [61]

    Wolloch, G

    M. Wolloch, G. Levita, P. Restuccia, and M. C. Righi, Interfacial Charge Density and Its Connection to Adhesion and Frictional Forces, Phys. Rev. Lett. 121, 026804 (2018). S-19

  17. [62]

    P. C. Torche, A. Silva, D. Kramer, T. Polcar, and O. Hovorka, Multi-scale model predicting friction of crystalline materials, Adv. Mater. Interfaces 9, 2100914 (2022)

  18. [63]

    Szlufarska, M

    I. Szlufarska, M. Chandross, and R. W. Carpick, Recent advances in single-asperity nanotribology, J. Phys. D: Appl. Phys. 41, 123001 (2008)

  19. [64]

    Tomlinson, CVI

    G. Tomlinson, CVI. A molecular theory of friction , London Edinb. Dublin Philos. Mag. J. Sci. 7, 905 (1929)

  20. [65]

    Prandtl, Ein Gedankenmodell zur kinetischen Theorie der festen K¨ orper, J

    L. Prandtl, Ein Gedankenmodell zur kinetischen Theorie der festen K¨ orper, J. Appl. Math. Mech. 8, 8, 85 (1928)

  21. [66]

    F. R. Krajewski and M. H. M¨ user, Quantum Creep and Quantum-Creep Transitions in 1D Sine-Gordon Chains, Phys. Rev. Lett. 92, 030601 (2004)

  22. [67]

    F. R. Krajewski and M. H. M¨ user, Quantum dynamics in the highly discrete, commensurate Frenkel Kontorova model: A path-integral molecular dynamics study, J. Chem. Phys. 122, 124711 (2005)

  23. [68]

    H. Xu, W. Chen, and Y. Zhu, Influence of the bond defect in driven Frenkel-Kontorova chains, Phys. Rev. B 75, 224303 (2007)

  24. [69]

    Landau, Zur Theorie der Energieubertragung

    L. Landau, Zur Theorie der Energieubertragung. II, Physikalische Zeitschrift der Sowjetunion 2, 46–51 (1932)

  25. [70]

    Zener, Non-Adiabatic Crossing of Energy Levels, Proc

    C. Zener, Non-Adiabatic Crossing of Energy Levels, Proc. R. Soc. London Ser. A 137, 696 (1932)

  26. [71]

    V. L. Pokrovsky and N. A. Sinitsyn, Landau-Zener transitions in a linear chain, Phys. Rev. B 65, 153105 (2002)

  27. [72]

    Huang and Y

    Z. Huang and Y. Zhao, Dynamics of dissipative Landau-Zener transitions, Phys. Rev. A 97, 013803 (2018)

  28. [73]

    Zanca, F

    T. Zanca, F. Pellegrini, G. E. Santoro, and E. Tosatti, Frictional lubricity enhanced by quantum mechanics, Proc. Natl. Acad. Sci. U.S.A. 115, 3547 (2018)

  29. [74]

    Arceci, S

    L. Arceci, S. Barbarino, R. Fazio, and G. E. Santoro, Dissipative Landau-Zener problem and thermally assisted Quantum Annealing, Phys. Rev. B 96, 054301 (2017)

  30. [75]

    Stone, The Theory of Intermolecular Forces (Oxford University Press, 2013)

    A. Stone, The Theory of Intermolecular Forces (Oxford University Press, 2013)

  31. [76]

    [48, 49] therein

    See Supplementary Information for detail derivations of coupling potential between the nanoparticle and the atomic chain, matrix elements of quantum Prandtl-Tomlison Hamiltonian in moving harmonic oscillator basis set, numerical calculation for quasi-static eigenenergies and eigenstates, quantum evolution with and without external bath reservoir, classica...

  32. [77]

    Meyer, T

    E. Meyer, T. Gyalog, R. M. Overney, and K. Dransfeld, Nanoscience: Friction and Rheology on the Nanometer Scale (World Scientific, 1998)

  33. [78]

    M. H. M¨ user, Velocity dependence of kinetic friction in the Prandtl-Tomlinson model, Phys. Rev. B 84, 125419 (2011)

  34. [79]

    J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed. (Cambridge University Press, 2020)

  35. [80]

    Yamaguchi, T

    M. Yamaguchi, T. Yuge, and T. Ogawa, Markovian quantum master equation beyond adiabatic regime, Phys. Rev. E 95, 012136 (2017)

  36. [81]

    The considered time domain for the moving particle (from tmin = 0 to tmax = 3 T in this study) is divided evenly by infinitesimal small intervals ∆ t who must be smaller than the smallest period of Bohr oscillation between eigenstates ∆t < max(En(t) − Em(t)) −1

  37. [82]

    N. A. Peters, T.-C. Wei, and P. G. Kwiat, Mixed-state sensitivity of several quantum-information benchmarks, Phys. Rev. A 70, 052309 (2004)

  38. [83]

    D. M. Tong, E. Sj¨ oqvist, L. C. Kwek, and C. H. Oh, Kinematic Approach to the Mixed State Geometric Phase in Nonunitary Evolution, Phys. Rev. Lett. 93, 080405 (2004)

  39. [84]

    Socoliuc, R

    A. Socoliuc, R. Bennewitz, E. Gnecco, and E. Meyer, Transition from Stick-Slip to Continuous Sliding in Atomic Friction: Entering a New Regime of Ultralow Friction, Phys. Rev. Lett. 92, 134301 (2004)

  40. [85]

    Gangloff, A

    D. Gangloff, A. Bylinskii, I. Counts, W. Jhe, and V. Vuleti´ c, Velocity tuning of friction with two trapped atoms, Nat. Phys. 11, 915–919 (2015)

  41. [86]

    Counts, D

    I. Counts, D. Gangloff, A. Bylinskii, J. Hur, R. Islam, and V. Vuleti´ c, Multislip Friction with a Single Ion, Phys. Rev. Lett. 119, 043601 (2017)

  42. [87]

    Le, N.-T

    D.-N. Le, N.-T. D. Hoang, and V.-H. Le, Exact analytical solutions of the schr¨ odinger equation for a two dimensional purely sextic double-well potential, Journal of Mathematical Physics 59, 032101 (2018), https://doi.org/10.1063/1.4997532

  43. [88]

    A. O. Caldeira and A. J. Leggett, Influence of Dissipation on Quantum Tunneling in Macroscopic Systems, Phys. Rev. Lett. 46, 211 (1981)

  44. [89]

    Weiss, Quantum Dissipative Systems , 4th ed

    U. Weiss, Quantum Dissipative Systems , 4th ed. (World Scientific, Singapore, 2012)

  45. [90]

    Dittrich, P

    T. Dittrich, P. H¨ aanggi, G.-L. Ingold, B. Krammer, G. Sch¨ on, and W. Zwerger,Quantum transport and dissipation (Wiley- VCH, Weinheim, Germany, 1998)

  46. [91]

    Shevchenko, S

    S. Shevchenko, S. Ashhab, and F. Nori, Landau–zener–st¨ uckelberg interferometry, Physics Reports492, 1 (2010)

  47. [92]

    Bouwmeester, G

    D. Bouwmeester, G. P. Karman, C. A. Schrama, and J. P. Woerdman, Observation of interference in transitions due to local geometric phases, Phys. Rev. A 53, 985 (1996)

  48. [93]

    Gasparinetti, P

    S. Gasparinetti, P. Solinas, and J. P. Pekola, Geometric landau-zener interferometry, Phys. Rev. Lett. 107, 207002 (2011)

  49. [94]

    Viotti, F

    L. Viotti, F. C. Lombardo, and P. I. Villar, Geometric phase in a dissipative jaynes-cummings model: Theoretical expla- nation for resonance robustness, Phys. Rev. A 105, 022218 (2022)

  50. [95]

    F. C. Lombardo and P. I. Villar, Geometric phases in open systems: A model to study how they are corrected by decoherence, Phys. Rev. A 74, 042311 (2006)

  51. [96]

    Viotti, A

    L. Viotti, A. L. Gramajo, P. I. Villar, F. C. Lombardo, and R. Fazio, Geometric phases along quantum trajectories, Quantum 7, 1029 (2023)

  52. [97]

    S. M. Rytov, Y. A. Kravtsov, and V. I. Tatarskii, Principles of Statistical Radiophysics 3 (Springer, Berlin, Germany, S-20 1989)

  53. [98]

    Yang, K.-t

    C.-W. Yang, K.-t. Leung, R.-F. Ding, H.-C. Ko, Y.-H. Lu, C.-K. Fang, and I.-S. Hwang, Lateral force microscopy of interfacial nanobubbles: Friction reduction and novel frictional behavior, Scientific Reports 8, 3125 (2018)

  54. [99]

    Sch¨ affer, S

    E. Sch¨ affer, S. F. Nørrelykke, and J. Howard, Surface forces and drag coefficients of microspheres near a plane surface measured with optical tweezers, Langmuir 23, 3654 (2007)

  55. [100]

    Gnecco, R

    E. Gnecco, R. Roth, and A. Baratoff, Analytical expressions for the kinetic friction in the Prandtl-Tomlinson model, Phys. Rev. B 86, 035443 (2012)

  56. [101]

    Bonvin, L

    E. Bonvin, L. Devaud, M. Rossi, A. Militaru, L. Dania, D. S. Bykov, M. Teller, T. E. Northup, L. Novotny, and M. Frimmer, Hybrid paul-optical trap with large optical access for levitated optomechanics, Phys. Rev. Res. 6, 043129 (2024)

  57. [102]

    Zenesini, H

    A. Zenesini, H. Lignier, G. Tayebirad, J. Radogostowicz, D. Ciampini, R. Mannella, S. Wimberger, O. Morsch, and E. Arimondo, Time-resolved measurement of landau-zener tunneling in periodic potentials, Phys. Rev. Lett. 103, 090403 (2009)

  58. [103]

    I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products (Academic press, 2014)

  59. [104]

    Wilkie and M

    J. Wilkie and M. C ¸ etinba¸ s, Variable-stepsize runge–kutta methods for stochastic schr¨ odinger equations, Physics Letters A 337, 166 (2005)