249 Downloads (Pure)


A numerical method is set out which efficiently computes stationary (z-independent) two- and three-dimensional spatiotemporal solitons in second-harmonic-generating media. The method relies on a Chebyshev decomposition with an infinite mapping, bunching the collocation points near the soliton core. Known results for the type-I interaction are extended and a stability boundary is found by two-parameter continuation as defined by the Vakhitov-Kolokolov criteria. The validity of this criteria is demonstrated in (2+1)-dimensions by simulation and direct calculation of the linear spectrum. The method has wider applicability for general soliton-bearing equations in (2+1)- and (3+1)-dimensions.
Original languageEnglish
Publication statusUnpublished - 2004

Bibliographical note

Sponsorship: The authors would like to thank John P. Boyd, Boris A. Malomed and Dimitru Mihalache for their useful remarks. The authors would also like to thank Richard R. Kerswell for initially suggesting the numerical discretisation.

Structured keywords

  • Engineering Mathematics Research Group


  • optical soliton
  • second-harmonic-generation
  • type-I interaction
  • continuation


Dive into the research topics of 'Efficient numerical continuation and stability analysis of spatiotemporal quadratic optical solitons'. Together they form a unique fingerprint.

Cite this