Parametric interactions in periodically poled lithium niobate are studied by using 2D non-paraxial finite-difference beam propagation method based on Pade (1,1) approximant of propagation operator and implicit quasi-Crank-Nicholson scheme. The nonlinear perfectly matched layer boundary conditions for non-paraxial beam propagation method are developed and applied to the simulation of parametric spatial solitons. Its effectiveness is verified by way of comparison of different realization of transparent boundary conditions.
"Modeling of nonparaxial propagation of parametric spatial solitons", Proc. SPIE 4415, Optical Organic and Inorganic Materials, (10 April 2001); doi: 10.1117/12.425498; https://doi.org/10.1117/12.425498