Backward Euler method: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
local error is O(h^2), see talk
 
Line 1: Line 1:
{| border="1" cellpadding="5" cellspacing="0" align="right" style="width:220px; text-align:justify"
Hi there. Allow me begin by introducing the author, her title is Sophia. Office supervising is where her main earnings comes from but she's already utilized for an additional one. Playing badminton is a factor that he is completely addicted to. Ohio is where her home is.<br><br>My web-site: best psychics ([http://Cartoonkorea.com/ce002/1093612 Cartoonkorea.com])
! style="background:#ffdead;" | [[Feynman diagrams]]
|-
| align="center" | '''Annihilation'''<br>[[Image:Electron-positron-annihilation.svg|220px]]
|-
| align="center" | '''Scattering'''<br>[[Image:Electron-positron-scattering.svg|220px]]
|-
|Note: Time moves forward from the left side of the diagram to the right. '''The arrows are simply markers of particle motion, and are ''not'' the same as the arrows conventionally written into Feynman diagrams'''.
|}
 
In [[quantum electrodynamics]], '''Bhabha scattering''' is the [[electron]]-[[positron]] [[scattering]] process:
::<math>e^+ e^- \rightarrow e^+ e^-</math>
 
There are two leading-order [[Feynman diagram]]s contributing to this interaction: an annihilation process and a scattering process. Bhabha scattering is named after the Indian physicist [[Homi J. Bhabha]].
 
The Bhabha scattering rate is used as a [[luminosity]] monitor in electron-positron colliders.
 
==Differential cross section==
To [[leading-order|leading order]], the spin-averaged [[differential cross section]] for this process is
 
::<math>\frac{\mathrm{d} \sigma}{\mathrm{d} (\cos\theta)} = \frac{\pi \alpha^2}{s} \left( u^2 \left( \frac{1}{s} + \frac{1}{t} \right)^2 + \left( \frac{t}{s} \right)^2 + \left( \frac{s}{t} \right)^2 \right) \,</math>
 
where ''s'',''t'', and ''u'' are the [[Mandelstam variables]], <math>\alpha</math> is the [[fine-structure constant]], and <math>\theta</math> is the scattering angle.
 
This cross section is calculated neglecting the electron mass relative to the collision energy and including only the contribution from photon exchange.  This is a valid approximation at collision energies small compared to the mass scale of the [[Z boson]], about 91 GeV; at higher energies the contribution from Z boson exchange also becomes important.
 
===Mandelstam variables===
In this article, the [[Mandelstam variables]] are defined by
::{|
|align="right"|<math>s= \,</math>
|align="right"|<math>(k+p)^2= \,</math>
|align="right"|<math>(k'+p')^2 \approx \,</math>
|align="right"|<math>2 k \cdot p \approx\,</math>
|align="right"|<math> 2 k' \cdot p' \,</math>
|rowspan="3"|&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[[Image:Mandelstam01.png]]
|-
|align="right"|<math>t= \,</math>
|align="right"|<math>(k-k')^2= \,</math>
|align="right"|<math>(p-p')^2\approx  \,</math>
|<math> -2 k \cdot k' \approx \,</math>
|<math> -2 p \cdot p' \,</math>
|-
|align="right"|<math>u= \,</math>
|align="right"|<math>(k-p')^2= \,</math>
|align="right"|<math>(p-k')^2\approx \,</math>
|<math> -2 k \cdot p' \approx \,</math>
|<math> -2 k' \cdot p \,</math>
|}
 
where the approximations are for the high-energy (relativistic) limit.
 
==Deriving unpolarized cross section==
 
===Matrix elements===
Both the scattering and annihilation diagrams contribute to the transition matrix element. By letting ''k'' and ''k' '' represent the four-momentum of the positron, while letting ''p'' and ''p' '' represent the four-momentum of the electron, and by using [[Feynman rules]] one can show the following diagrams give these matrix elements:
:{| border="0" cellpadding="5" cellspacing="0"
|
|align="center" | [[Image:Feynman-scattering.svg|160px]]
|align="center" | [[Image:Feynman-annihilation.svg|160px]]
|Where we use:<br><math>\gamma^\mu \,</math> are the [[Gamma matrices]],<br><math>u, \ \mathrm{and} \ \bar{u}\,</math> are the four-component spinors for fermions, while<br><math>v, \ \mathrm{and} \ \bar{v}\,</math> are the four-component spinors for anti-fermions (see [[Dirac equation#Four spinor|Four spinors]]).
|-
|
|align="center" | (scattering)
|align="center" | (annihilation)
|
|-
|<math>\mathcal{M} = \,</math>
|<math>-e^2 \left( \bar{v}_{k} \gamma^\mu v_{k'} \right) \frac{1}{(k-k')^2} \left( \bar{u}_{p'} \gamma_\mu u_p \right) </math>
|<math>+e^2 \left( \bar{v}_{k} \gamma^\nu u_p \right) \frac{1}{(k+p)^2} \left( \bar{u}_{p'} \gamma_\nu v_{k'} \right) </math>
|
|}
 
Notice that there is a relative sign difference between the two diagrams.
 
===Square of matrix element===
To calculate the unpolarized [[Cross section (physics)|cross section]], one must ''average'' over the spins of the incoming particles (''s''<sub>e-</sub> and ''s''<sub>e+</sub> possible values) and ''sum'' over the spins of the outgoing particles.  That is,
::{|
|<math>\overline{|\mathcal{M}|^2} \,</math>
|<math> = \frac{1}{(2s_{e-} + 1)(2 s_{e+} + 1)} \sum_{\mathrm{spins}} |\mathcal{M}|^2 \,</math>
|-
|
|<math>= \frac{1}{4} \sum_{s=1}^2 \sum_{s'=1}^2 \sum_{r=1}^2 \sum_{r'=1}^2 |\mathcal{M}|^2 \,</math>
|}
 
First, calculate <math>|\mathcal{M}|^2 \,</math>:
::{| cellpadding=4
|<math>|\mathcal{M}|^2 \,</math>=
|<math> e^4 \left| \frac{(\bar{v}_{k} \gamma^\mu v_{k'} )( \bar{u}_{p'} \gamma_\mu u_p)}{(k-k')^2} \right|^2 \,</math>
| align="center" | (scattering)
|-
|
|<math>{}-  e^4 \left( \frac{ (\bar{v}_{k} \gamma^\mu v_{k'} )( \bar{u}_{p'} \gamma_\mu u_p)}{(k-k')^2} \right)^* \left( \frac{ (\bar{v}_{k} \gamma^\nu u_p )( \bar{u}_{p'} \gamma_\nu v_{k'}) }{(k+p)^2} \right)  \,</math>
| align="center" | (interference)
|-
|
|<math>{}-  e^4 \left( \frac{ (\bar{v}_{k} \gamma^\mu v_{k'} )( \bar{u}_{p'} \gamma_\mu u_p)}{(k-k')^2} \right) \left( \frac{ (\bar{v}_{k} \gamma^\nu u_p )( \bar{u}_{p'} \gamma_\nu v_{k'}) }{(k+p)^2} \right)^*  \,</math>
| align="center" | (interference)
|-
|
|<math>{}+ e^4 \left| \frac{(\bar{v}_{k} \gamma^\nu u_p )( \bar{u}_{p'} \gamma_\nu v_{k'} )}{(k+p)^2} \right|^2 \,</math>
| align="center" | (annihilation)
|}
 
===Scattering term (t-channel)===
 
====Magnitude squared of M====
::{|
|<math>|\mathcal{M}|^2 \,</math>
|<math>= \frac{e^4}{(k-k')^4} \Big( (\bar{v}_{k} \gamma^\mu v_{k'} )( \bar{u}_{p'} \gamma_\mu u_p) \Big)^* \Big( (\bar{v}_{k} \gamma^\nu v_{k'})( \bar{u}_{p'} \gamma_\nu u_p) \Big) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(1) \,</math>
|-
|
|<math>= \frac{e^4}{(k-k')^4} \Big( (\bar{v}_{k} \gamma^\mu v_{k'} )^* ( \bar{u}_{p'} \gamma_\mu u_p)^* \Big) \Big( (\bar{v}_{k} \gamma^\nu v_{k'})( \bar{u}_{p'} \gamma_\nu u_p) \Big) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(2) \,</math>
|-
|
|align="center" | (complex conjugate will flip order)
|
|-
|
|<math>= \frac{e^4}{(k-k')^4} \Big( \left(\bar{v}_{k'} \gamma^\mu v_{k} \right) \left( \bar{u}_{p} \gamma_\mu u_{p'} \right) \Big) \Big( \left( \bar{v}_{k} \gamma^\nu v_{k'} \right) \left( \bar{u}_{p'} \gamma_\nu u_p \right) \Big) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(3) \,</math>
|-
|
|align="center" | (move terms that depend on same momentum to be next to each other)
|
|-
|
|<math>= \frac{e^4}{(k-k')^4} \left( \bar{v}_{k'} \gamma^\mu v_{k} \right) \left( \bar{v}_{k} \gamma^\nu v_{k'} \right) \left( \bar{u}_{p} \gamma_\mu u_{p'} \right) \left(  \bar{u}_{p'} \gamma_\nu u_p \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(4) \,</math>
|}
 
====Sum over spins====
Next, we'd like to sum over spins of all four particles. Let ''s'' and ''s' '' be the spin of the electron and ''r'' and ''r' '' be the spin of the positron.
::{|
|<math> \sum_{\mathrm{spins}} |\mathcal{M}|^2 \,</math>
|<math>= \frac{e^4}{(k-k')^4}\left(\sum_{r'} \bar{v}_{k'} \gamma^\mu (\sum_{r}v_{k} \bar{v}_{k}) \gamma^\nu v_{k'} \right) \left(\sum_{s} \bar{u}_{p} \gamma_\mu (\sum_{s'}{u_{p'} \bar{u}_{p'}}) \gamma_\nu u_p \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(5) \,</math>
|-
|
|<math>= \frac{e^4}{(k-k')^4}\operatorname{Tr}\left( \Big(\sum_{r'} v_{k'} \bar{v}_{k'} \Big) \gamma^\mu \Big(\sum_{r}v_{k} \bar{v}_{k} \Big) \gamma^\nu \right) \operatorname{Tr} \left( \Big(\sum_{s} u_p \bar{u}_{p} \Big) \gamma_\mu \Big( \sum_{s'}{u_{p'} \bar{u}_{p'}} \Big) \gamma_\nu \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(6) \,</math>
|-
|
|align="center" | (now use [[Bhabha scattering#Completeness relations|Completeness relations]])
|
|-
|
|<math>=\frac{e^4}{(k-k')^4}\operatorname{Tr}\left( (k\!\!\!/' - m) \gamma^\mu (k\!\!\!/ - m) \gamma^\nu \right) \cdot \operatorname{Tr}\left( (p\!\!\!/' + m) \gamma_\mu (p\!\!\!/ + m) \gamma_\nu \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(7) \,</math>
|-
|
|align="center" | (now use [[Bhabha scattering#Trace identities|Trace identities]])
|
|-
|
|<math>=\frac{e^4}{(k-k')^4}\left(4 \left( {k'}^\mu k^\nu - (k' \cdot k)\eta^{\mu\nu} + k'^\nu k^\mu \right) + 4 m^2 \eta^{\mu\nu}  \right) \left( 4 \left( {p'}_\mu p_\nu - (p' \cdot p)\eta_{\mu\nu} + p'_\nu p_\mu \right) + 4 m^2 \eta_{\mu\nu} \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(8) \,</math>
|-
|
|<math>=\frac{32{e^4}}{(k-k')^4}\left( (k' \cdot p') (k \cdot p) + (k' \cdot p) (k \cdot p') -m^2 p' \cdot p - m^2 k' \cdot k + 2m^4 \right) \,</math>
|&nbsp;&nbsp;&nbsp;&nbsp; <math>(9) \,</math>
|}
 
Now that is the exact form, in the case of electrons one is usually interested in energy scales that far exceed the electron mass.  Neglecting the electron mass yields the simplified form:
::{|
|<math>\frac{1}{4} \sum_{\mathrm{spins}} |\mathcal{M}|^2 \,</math>
|<math> = \frac{32e^4}{4(k-k')^4} \left( (k' \cdot p') (k \cdot p) + (k' \cdot p) (k \cdot p') \right) \,</math>
|-
|
|(use the [[Bhabha scattering#Mandelstam variables|Mandelstam variables]] in this relativistic limit)
|-
|
|<math>=\frac{8e^4}{t^2} \left(\tfrac{1}{2} s \tfrac{1}{2}s + \tfrac{1}{2}u \tfrac{1}{2} u \right) \,</math>
|-
|
|<math>= 2 e^4 \frac{s^2 +u^2}{t^2} \,</math>
|}
 
===Annihilation term (s-channel)===
The process for finding the annihilation term is similar to the above. Since the two diagrams are related by [[crossing symmetry]], and the initial and final state particles are the same, it is sufficient to permute the momenta, yielding
::{|
|<math>\frac{1}{4} \sum_{\mathrm{spins}} |\mathcal{M}|^2 \,</math>
|<math> = \frac{32e^4}{4(k+p)^4} \left( (k \cdot k') (p \cdot p') + (k' \cdot p) (k \cdot p') \right) \,</math>
|-
|
|<math>=\frac{8e^4}{s^2} \left(\tfrac{1}{2} t \tfrac{1}{2}t + \tfrac{1}{2}u \tfrac{1}{2} u \right) \,</math>
|-
|
|<math>= 2 e^4 \frac{t^2 +u^2}{s^2} \,</math>
|}
(This is proportional to
<math>(1 + \cos^2\theta)</math>
where <math>\theta</math> is the scattering angle in the center-of-mass frame.)
 
===Solution===
Evaluating the interference term along the same lines and adding the three terms yields the final result
::<math>\frac{\overline{|\mathcal{M}|^2}}{2e^4} = \frac{u^2 + s^2}{t^2} + \frac{2 u^2}{st} + \frac{u^2 + t^2}{s^2} \,</math>
 
==Simplifying steps==
 
===Completeness relations===
The completeness relations for the [[Dirac_spinor#Four-spinor_for_particles|four-spinors]] ''u'' and ''v'' are
::<math>\sum_{s=1,2}{u^{(s)}_p \bar{u}^{(s)}_p} = p\!\!\!/ + m \,</math>
::<math>\sum_{s=1,2}{v^{(s)}_p \bar{v}^{(s)}_p} = p\!\!\!/ - m \,</math>
:where
::<math>p\!\!\!/ = \gamma^\mu p_\mu  \,</math> &nbsp;&nbsp;&nbsp;&nbsp; (see [[Feynman slash notation#With four-momentum|Feynman slash notation]])
::<math>\bar{u} = u^{\dagger} \gamma^0 \,</math>
 
===Trace identities===
:''Main article:'' [[Gamma matrix#Identities|Trace identities]]
To simplify the trace of the [[Gamma matrices|Dirac gamma matrices]], one must use trace identities. Three used in this article are:
#The Trace of any product of an ''odd number'' of <math>\gamma_\mu \,</math>'s is zero
#<math>\operatorname{Tr} (\gamma^\mu\gamma^\nu) = 4\eta^{\mu\nu}</math>
#<math>\operatorname{Tr}\left( \gamma_\rho \gamma_\mu \gamma_\sigma \gamma_\nu \right) = 4 \left( \eta_{\rho\mu}\eta_{\sigma\nu}-\eta_{\rho\sigma}\eta_{\mu\nu}+\eta_{\rho\nu}\eta_{\mu\sigma} \right) \,</math>
 
Using these two one finds that, for example,
::{|
|<math>\operatorname{Tr}\left( (p\!\!\!/' + m) \gamma_\mu (p\!\!\!/ + m) \gamma_\nu \right) \,</math>
|<math> = \operatorname{Tr}\left( p\!\!\!/' \gamma_\mu p\!\!\!/ \gamma_\nu \right) + \operatorname{Tr}\left(m \gamma_\mu p\!\!\!/ \gamma_\nu \right)  \,</math>
|-
|
|&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <math>+ \operatorname{Tr}\left( p\!\!\!/' \gamma_\mu m \gamma_\nu \right) + \operatorname{Tr}\left(m^2 \gamma_\mu \gamma_\nu \right) \,</math>
|-
|
|align="center"|(the two middle terms are zero because of (1))
|-
|
|<math> = \operatorname{Tr}\left( p\!\!\!/' \gamma_\mu p\!\!\!/ \gamma_\nu \right) + m^2 \operatorname{Tr}\left(\gamma_\mu \gamma_\nu \right)  \,</math>
|-
|
|align="center"|(use identity (2) for the term on the right)
|-
|
|<math>= {p'}^{\rho} p^\sigma \operatorname{Tr}\left( \gamma_\rho \gamma_\mu \gamma_\sigma \gamma_\nu \right) + m^2 \cdot 4\eta_{\mu\nu} \,</math>
|-
|
|align="center"|(now use identity (3) for the term on the left)
|-
|
|<math>= {p'}^{\rho} p^\sigma 4 \left( \eta_{\rho\mu}\eta_{\sigma\nu}-\eta_{\rho\sigma}\eta_{\mu\nu}+\eta_{\rho\nu}\eta_{\mu\sigma} \right) + 4 m^2 \eta_{\mu\nu} \,</math>
|-
|
|<math>=4 \left( {p'}_\mu p_\nu - \mathbf{p' \cdot p}\eta_{\mu\nu} + p'_\nu p_\mu \right) + 4 m^2 \eta_{\mu\nu} \,</math>
|}
 
==Uses==
Bhabha scattering has been used as a [[luminosity]] monitor in a number of e+e- collider physics experiments.  The accurate measurement of luminosity is necessary for accurate measurements of cross sections. 
* Small-angle Bhabha scattering was used to measure the luminosity of the 1993 run of the [[Stanford Large Detector]] (SLD), with a relative uncertainty of less than 0.5%.<ref>[http://adsabs.harvard.edu/abs/1995PhDT.......160W A Study of Small Angle Radiative Bhabha Scattering and Measurement of the Luminosity at SLD]</ref>
* Electron-positron colliders operating in the region of the low-lying hadronic resonances (about 1 GeV to 10 GeV), such as the [[Beijing Electron Synchrotron]] (BES) and the [[Belle experiment|Belle]] and [[BaBar experiment|BaBar]] "B-factory" experiments, use large-angle Bhabha scattering as a luminosity monitor.  To achieve the desired precision at the 0.1% level, the experimental measurements must be compared to a theoretical calculation including next-to-leading-order [[renormalization|radiative corrections]].[http://arxiv.org/abs/hep-ph/0003268]  The high-precision measurement of the total hadronic cross section at these low energies is a crucial input into the theoretical calculation of the [[anomalous magnetic dipole moment]] of the [[muon]], which is used to constrain [[supersymmetry]] and other models of physics [[beyond the Standard Model]].
 
==References==
<references/>
 
*{{cite book | author=Halzen, Francis; Martin, Alan | title=Quarks & Leptons: An Introductory Course in Modern Particle Physics | publisher=John Wiley & Sons | year=1984 | isbn=0-471-88741-2}}
*{{cite book | author=Peskin, Michael E.; Schroeder, Daniel V. | title=An Introduction to Quantum Field Theory | publisher=Perseus Publishing | year=1994 | isbn=0-201-50397-2}}
*[http://xstructure.inr.ac.ru/x-bin/theme3.py?level=1&index1=270563 Bhabha scattering on arxiv.org]
{{QED}}
 
[[Category:Quantum field theory]]
[[Category:Quantum electrodynamics]]

Latest revision as of 22:04, 20 April 2014

Hi there. Allow me begin by introducing the author, her title is Sophia. Office supervising is where her main earnings comes from but she's already utilized for an additional one. Playing badminton is a factor that he is completely addicted to. Ohio is where her home is.

My web-site: best psychics (Cartoonkorea.com)