3. Gaussian Beams > 2. Photodetectors
Beat coefficients
Author: Daniel Brown
If Gaussian modes are detected by a simple photo detector which is much larger than the beam size, the total power is the sum of the powers in each modes because the modes are orthogonal.
However, if the photo detector has some features (or is to small) the modes couple to each other so that cross terms must be computed. This feature is used in so-called composite photo detectors that consist of seprate elements whose signals is are added or subtracted to generate the detector output signal.
The ideal split photo detector consists of two half-infinite planes. The signal from the two planes are subtracte from each other. For a beam which is centered on the line seperating the two planes, the detector will only show the beat between symmetric and asymmetric modes, which is useful for generating alignment signals.
For a beam: $$E~=~~\sum_{n,m}~a_{nm}~u_{nm}(x,y,z)$$ the power is proportional to $\int E*E^*$
For a normal detector the integral would be over the entire x-y plane. For the split-photodectors (here assuming split along the vertical y-axis) we get $$P\sim \int_0^\infty E*E^* - \int_0^{-\infty} E*E^* =2 \int_0^\infty E*E^*$$
In order to simplify the calculation we can combute so called beat coefficents $c_{n1,n2}$ so that we can write $$P~\sim~\sum_{m}\sum_{n1,n2}~ c_{n1,n2} ~a_{n1,m}~a_{n2,m}^*$$ These can be computed as follows:
import numpy as np
import math
from math import factorial
from scipy.integrate import quad
from scipy.special import hermite
def HG_split_diode_coefficient_numerical(n1,n2):
# return zero if n1,n2 are both even or both odd
if (n1+n2) % 2 == 0:
return 0.0
# extract normlisation coefficients
A = 2.0 * np.sqrt(2.0/np.pi) * np.sqrt(1.0 / (2.0**(n1+n2) * factorial(n1) * factorial(n2)))
# define rest of u_n1 * u_n2^* as integrand
f = lambda x: hermite(n1)(np.sqrt(2.0)*x) * math.exp(-2.0*x*x) * hermite(n2)(np.sqrt(2.0)*x)
# perform numerical integration
val, res= quad(f, 0.0, np.Inf, epsabs=1e-10, epsrel=1e-10, limit=200)
return A * val
print(HG_split_diode_coefficient_numerical(0,0))
print(HG_split_diode_coefficient_numerical(0,1))
If both n1 and n2 are even or odd then the beat coefficient is zero, if one is odd and one even, the value is computed by numerical integration.
This function has been implemented in PyKat and can be used as follows:
import pykat.optics.pdtype as pdtype
pdtype.HG_split_diode_coefficient_numerical(0,1)
The beat coefficients can be used to define a new photodiode type with the 'pdtype' command and the respective entry in the kat.ini file. To creat the coefficients in the format for the kat.ini file please use pdtype.finesse_split_photodiode(maxtem,'x' or 'y'), for example:
pdtype.finesse_split_photodiode(5,'x')