Response function

The ResponseFunction object in Quanty defines (linear) response functions. For Hamiltonian $H$, ground-state $| \psi_0 \rangle$ and a list of transition operators $T_i$ with $i \in [1,n]$ we define the response function $G(\omega,\Gamma)$ as giving an $n$ by $n$ matrix for each value of $\omega$ and $\Gamma$. The elements of this matrix are given as $$ G_{i,j}(\omega,\Gamma) = \left\langle \psi_0 \middle| T^{\dagger}_i \frac{1}{\omega - H + \mathrm{i} \Gamma/2 + E_0} T^{\phantom{\dagger}}_j \middle| \psi_0 \right\rangle, $$ with $E_0 = \left\lange \psi_0 \middle| H \middle| \psi_0 \right\rangle$. One can calculate response functions using the Quanty function CreateSpectra. This function returns two objects. At the first position a spectra object that contains the intensity for given values of $\omega$ and one specific value of $\Gamma$ on a grid. At the second position CreateSpectra returns a ResponseFunction object.

ResponseFunctions are objects that can be evaluated at any frequency or imaginary onset. For example:

Example.Quanty

H   = Matrix.ToOperator( Matrix.Diagonal({1,2,3,4,5}) )
psi = NewWavefunction(5,0,{{"00000",1}})
T = {}
for i=0,4 do
  T[i+1] = NewOperator(5,0,{{i,1}})
end
S, G = CreateSpectra(H,T,psi)
omega = 1.1
gamma = 0.01
print(G[1](omega,gamma))

returns

(9.9750623441396 - 0.49875311720698 I)

i.e. the value of the response function for the first transition operator defined at $\omega=1.1$ and $\Gamma=0.01$.