fkoe FresnelKirchhoff

fkoe

fkoe is stand-allone program which does a numerical solution of the Fresnel- Kirchhoff integral (see SPIE article bahrdt:2014). The program is coded in Fortran. It provides basic functionality. It is controlled by simple text files and command line. The calculation can be time consuming (the example with default settings took 3 hours). The program is good to have a reference and benchmark for other codes. See also the mpi version fkoempi.

Besides integration fkoe can generate a Gaussian source, an optical element and evaluate the optical path. The power in the Gaussian source is normalized to 1W. Output files in wave_file format and phase_hdf5 format are provided.

Examples

Gaussian Source

parameters: lambda: 10nm, waist: 70mum, size: 1.6mm^2, points: 301^2, drift: 5m

statistics: fwhm: 280 mum, rms: 119 mum, intensity: 0.5W, max: 5.6e6 W/m^2

hdf5 file

.

Result (mode 3, source to OE to image)

parameter: source to image, toroid 10m:0.5m, 88 deg, y: -10..20mum, 101 points, z: -10..10mum 51points, calculation time: 3 hours

statistics: zfwhm: 3.7 mum, zrms: 1.6 mum, yfwhm: 5.4 mum, yrms: 2.3 mum, intensity: 0.47W, max: 1.4e10 W/m^2

hdf5 file

intensity distribution on image plane (contour plot) intensity distribution on image plane (surface plot) phase distribution on image plane (contour plot, unwraped) .

Result (mode 1, source to OE)

parameters: same as mode 3, calculation time: 4 hours

statistics: zfwhm: 806 mum, zrms: 343 mum, yfwhm: 23 mm, yrms: 9.8 mm, intensity: 0.5W, max: 2.4e4 W/m^2

hdf5 file

intensity distribution on OE (contour plot) phase distribution on OE (contour plot, phunwraped) .

Result (mode 2, OE to image)

parameters: same as mode 3, outputput files from mode 2 as input, calculation time: 10 minutes.

statistics: almost the same as for mode 1

hdf5 file

we skip the plots- since they look identical to mode 1

Result optical path

to be done

Result phaseqt (for comparison)

parameter: source to image, toroid 10m:0.5m, 88 deg, y: -10..20mum, 101 points, z: -10..10mum 51points, 601 points angle grid, -2.2<dy<2.2, -3.6<dz<3.6 fkoe.phase

statistics: zfwhm: 3.9 mum, zrms: 1.7 mum, yfwhm: 5.6 mum, yrms: 2.4 mum, intensity: 0.53W, max: 1.15e10 W/m^2

hdf5 file

intensity distribution on image plane (contour plot)

!! there is a discrepancy in the sign of the output- in phaseqt we had to say deflection downwards to reproduce the result of fkoe !!