Numerical codes developed in our group

MATLAB implementation of the statistical analysis code that performs proper averaging of the temporal traces acquired using THz-TDS setups. The algorithms aim at denoising the THz (or any pump-probe) temporal traces by extracting "slow" noise variations (parameter drifts) from the data.

Related publication:
M. Skorobogatiy, J. Sadasivan, and H. Guerboukha, “Statistical models for averaging of the pump-probe time traces: Example of denoising in terahertz time domain spectroscopy,” arXiv:1711.03566 [physics.ins-det], (2017).

MATLAB implementation of the boundary integral method for calculating leaky and guided modes of microstructured optical fibers is presented. The method can handle a large number of inclusions (hundreds) of arbitrary geometries. Both, solid and hollow core photonic crystal fibers can be treated efficiently. For large systems featuring closely spaced inclusions the computational intensity of the boundary integral method is significantly smaller than that of the multipole method. This is of particular importance in the case of hollow core band gap guiding fibers.

Related publication:
E. Pone, A. Hassani, S. Lacroix, A. Kabashin and M. Skorobogatiy “Boundary integral method for the challenging problems in bandgap guiding, plasmonics and sensing,” Optics Express, vol. 15, p.10231 (2007).

MATLAB implementation of the statistical image analysis code to characterize imperfections of planar photonic crystals. The code is designed to treat high resolution images of planar photonic crystals (PC). Our original motivation was to attempt an intuitive, while rigorous statistical description of fabrication imperfections to provide a realistic input into theoretical modelling of PC device performance.

Related publication:
M. Skorobogatiy, G. Bégin, and A. Talneau, “Statistical analysis of imperfections in experimental planar 2D photonic crystals, 1Mb” Optics Express, vol. 13, p. 2487 (2005).