S-matrix formula from Lippmann-schwinger eqn

  • Context: Graduate 
  • Thread starter Thread starter muppet
  • Start date Start date
  • Tags Tags
    Formula S-matrix
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
muppet
Messages
602
Reaction score
0
Hi all,
In chapter 3.2 of Weinberg's QFT text he asserts that one can derive the expression
[tex]S_{\beta\alpha}=\delta(\beta-\alpha)-2\pi i\delta(E_\alpha -E_\beta)T_{\beta\alpha}[/tex]
for the S-matrix in terms of the matrix elements [tex]T^{\pm}_{\beta \alpha}=(\Phi_\beta , V\Psi^{\pm}_\alpha)[/tex], where I'm using Weinberg's "linear algebraic" notation for the inner product of state vectors, V is a scattering potential, [tex]\Phi_{\alpha}[/tex] is an eigenstate of a free hamiltonian [tex]H_0[/tex] with eigenvalue [tex]E_\alpha[/tex], and the [tex]\Psi^{\pm}[/tex] are in/out states, from the Lippmann -Schwinger equations
[tex]\Psi^{\pm}_\alpha=\Phi_{\alpha}+(E_\alpha-H_0 \pm i\epsilon )^{-1} V\Psi^{\pm}_\alpha[/tex].

However, he never actually gets around to it. Here's my attempt, which seems to get me close but contains either a mistake or a pathologically strange representation of the delta function.

Consider the matrix element [tex](\Psi_{\beta}^-,V\Psi_\alpha^+)[/tex]. We apply the Lippman-Schwinger equations to both the 'in' and 'out' states and equate the results.
Note that I'll use [tex]\Psi^{\pm}=\int d\gamma(\Phi_\gamma,\Psi^{\pm})\Phi_\gamma[/tex] a few times.
Starting with the in state:
[tex](\Psi_{\beta}^-,V\Psi_\alpha^{\plus})=(\Psi_{\beta}^-,V\Phi_{\alpha}+V(E_\alpha-H_0 + i\epsilon )^{-1} V\Psi^{+}_\alpha)[/tex]
[tex]=(\Psi_{\beta}^-,V\Phi_{\alpha})+(\Psi_{\beta}^-,V(E_\alpha-H_0 + i\epsilon )^{-1} V\Psi^{+}_\alpha)[/tex]
[tex]=(T^{-}_{\alpha \beta})^*+\int d \gamma (\Psi_{\beta}^-,V(E_\alpha-H_0 + i\epsilon )^{-1} T^+_{\gamma\alpha}\Phi_\gamma)[/tex]
[tex]=T^{-}_{ \beta \alpha}+\int d \gamma (E_\alpha-E_\gamma + i\epsilon )^{-1} T^+_{\gamma\alpha} (T^-_{\gamma \beta})^*=T^{-}_{ \beta \alpha}+\int d \gamma (E_\alpha-E_\gamma + i\epsilon )^{-1} T^-_{\beta \gamma}T^+_{\gamma\alpha}[/tex]

Now the out state:
[tex](\Psi_{\beta}^-,V\Psi_\alpha^{+})=(\Phi_{\beta}+(E_\blpha-H_0 - i\epsilon )^{-1} V\Psi^-_\beta,V\Psi_\alpha^{+})[/tex]
[tex]=(\Phi_{\beta},V\Psi_\alpha^{+})+(\Psi^-_\beta,V(E_\beta-H_0 + i\epsilon )^{-1} V\Psi_\alpha^{+})[/tex]
[tex]=T^{+}_{\beta \alpha}+ \int d \gamma (\Psi_{\beta}^-,V(E_\alpha-H_0 + i\epsilon )^{-1} T^+_{\gamma\alpha}\Phi_\gamma)[/tex]
[tex]=T^{+}_{ \beta \alpha}+\int d \gamma (E_\beta-E_\gamma + i\epsilon )^{-1} T^-_{\beta \gamma}T^+_{\gamma\alpha}[/tex]

Now we can get a relationship between [tex]T^{+}_{ \beta \alpha}[/tex] and [tex]T^{-}_{ \beta \alpha}[/tex] by expanding the defining matrix element for in states using the completeness relation for the out states:
[tex]T^{+}_{ \beta \alpha}=(\Phi_\beta,V\Psi_\alpha^+)=\int d\gamma (\Phi_\beta,V(\Psi^-_\gamma,\Psi_\alpha^+)\Psi^-_\gamma)=\int d\gamma S_{\gamma \alpha}T^-_{\beta\gamma}[/tex]

So equating the two expansions of [tex](\Psi_{\beta}^-,V\Psi_\alpha^{\plus})[/tex] and rearranging for [tex]T^+_{\beta\alpha}[/tex] leads to
[tex]T^+_{\beta\alpha}=\int d\gamma S_{\gamma \alpha}T^-_{\beta\gamma}[/tex]
[tex]=T^{-}_{ \beta \alpha}+\int d \gamma (E_\alpha-E_\gamma + i\epsilon )^{-1} T^-_{\beta \gamma}T^+_{\gamma\alpha} - (E_\beta-E_\gamma + i\epsilon )^{-1} T^-_{\beta \gamma}T^+_{\gamma\alpha}[/tex]
[tex]=\int d \gamma T^-_{\beta \gamma}(\delta(\alpha-\gamma)+\{(E_\alpha-E_\gamma + i\epsilon )^{-1} - (E_\beta-E_\gamma + i\epsilon )^{-1} \}T^+_{\gamma\alpha} )[/tex]

Now this argument is only correct if
[tex]\{(E_\alpha-E_\gamma + i\epsilon )^{-1} - (E_\beta-E_\gamma + i\epsilon )^{-1} \} = -2\pi i\delta(E_\alpha -E_\beta)[/tex]
inside an integral, which really doesn't appear to be true.

Any help would be greatly appreciated, I've spent most of today deriving this and then trying to correct my error. Thanks in advance.
 
Last edited:
Physics news on Phys.org
A:I think the problem is that you are not using the right expression for the delta function.$\delta(E_\alpha-E_\beta)$ is a Dirac delta function, which means that it is equal to $0$ everywhere expect for $E_\alpha=E_\beta$, where it is equal to $+\infty$. This means that it does not appear as part of an integral.The expression $\frac{1}{2\pi i}\bigg(\frac{1}{E_\alpha-E_\beta+i\epsilon}-\frac{1}{E_\alpha-E_\beta-i\epsilon}\bigg)$ is a distribution, and it is equal to $\delta(E_\alpha-E_\beta)$ when it appears as part of an integral.That is why Weinberg's expression for the S-matrix contains this term.