Relativistic Breit–Wigner distribution

<h2 id="usage">Usage</h2>
The probability of producing the resonance at a given energy E is proportional to f (E), so that a plot of the production rate of the unstable particle as a function of energy traces out the shape of the relativistic Breit–Wigner distribution. Note that for values of E off the maximum at M such that   | E 2 − M 2| = M Γ   , (hence   | E − M | = Γ / 2   for   M ≫ Γ  ), the distribution f has attenuated to half its maximum value, which justifies the name for Γ , width at half-maximum.
In the limit of vanishing width,   Γ → 0   , the particle becomes stable as the Lorentzian distribution f sharpens infinitely to   2 M δ(E 2 − M 2)   , where   δ   is the <a href="/facts/Dirac_delta_function/V8PjRP8T">Dirac delta function</a> (point impulse).
In general, Γ can also be a function of E ; this dependence is typically only important when Γ is not small compared to M and the <a href="/facts/Phase_space/qPe4Fq57">phase space</a>-dependence of the width needs to be taken into account. (For example, in the decay of the <a href="/facts/Rho_meson/55vNvssd">rho meson</a> into a pair of <a href="/facts/Pion/x5HamzG3">pions</a>.) The factor of   M 2   that multiplies  Γ2   should also be replaced with E 2 
(or   E 4 / M 2   , etc.) when the resonance is wide.<a class="footnote-ref" id="fnref:3" href="#fn:3">3</a>
The form of the relativistic Breit–Wigner distribution arises from the <a href="/facts/Propagator/FER8tNKP">propagator</a> of an unstable particle,<a class="footnote-ref" id="fnref:4" href="#fn:4">4</a> which has a denominator of the form   p 2 − M 2 + i M Γ   . (Here, p2 is the square of the <a href="/facts/Four-momentum/t1F5Lmbr">four-momentum</a> carried by that particle in the tree Feynman diagram involved.) The propagator in its rest frame then is proportional to the <a href="/facts/Probability_amplitude/ryPgJvs8">quantum-mechanical amplitude</a> for the decay utilized to reconstruct that resonance,

k
               
            
            
               
              
                (
                
                  
                    E
                    
                      2
                    
                  
                  −
                  
                    M
                    
                      2
                    
                  
                
                )
              
              +
              i
               
              M
               
              Γ
               
            
          
        
         
        .
      
    
    {\displaystyle \ {\frac {\sqrt {k\ }}{\ \left(E^{2}-M^{2}\right)+i\ M\ \Gamma \ }}~.}

The resulting probability distribution is proportional to the absolute square of the amplitude, so then the above relativistic Breit–Wigner distribution for the probability density function.
The form of this distribution is similar to the amplitude of the solution to the classical equation of motion for a <a href="/facts/Harmonic_oscillator/7IQNToGt">driven harmonic oscillator</a> damped and driven by a <a href="/facts/Sinusoidal/mD1CB09M">sinusoidal</a> external force. It has the standard <a href="/facts/Resonance/q7VLALUz">resonance</a> form of the Lorentz, or <a href="/facts/Cauchy_distribution/3HS6D3Xk">Cauchy distribution</a>, but involves relativistic variables   s = p 2   , here   = E 2   . The distribution is the solution of the differential equation for the amplitude squared w.r.t. the energy energy (frequency), in such a classical forced oscillator,

f
          ′
        
        
        (
        
          E
        
        )
         
        
          (
          
            
              
                (
                
                   
                  
                    
                      E
                    
                    
                      2
                    
                  
                  −
                  
                    M
                    
                      2
                    
                  
                
                )
              
              
                2
              
            
            +
            
              Γ
              
                2
              
            
             
            
              M
              
                2
              
            
          
          )
        
        −
        4
         
        
          E
        
         
        
          (
          
            
              M
              
                2
              
            
            −
            
              
                E
              
              
                2
              
            
          
          )
        
         
        f
        (
        
          E
        
        )
        =
        0
         
        ,
      
    
    {\displaystyle f'\!(\mathrm {E} )\ \left(\left(\ \mathrm {E} ^{2}-M^{2}\right)^{2}+\Gamma ^{2}\ M^{2}\right)-4\ \mathrm {E} \ \left(M^{2}-\mathrm {E} ^{2}\right)\ f(\mathrm {E} )=0\ ,}

or rather

f
                ′
              
              
              (
              
                E
              
              )
               
            
            
               
              f
              (
              
                E
              
              )
               
            
          
        
         
        =
         
        
          
            
               
              4
               
              
                (
                
                  
                    M
                    
                      2
                    
                  
                  −
                  
                    
                      E
                    
                    
                      2
                    
                  
                
                )
              
               
              
                E
              
               
            
            
               
              
                
                  (
                  
                     
                    
                      
                        E
                      
                      
                        2
                      
                    
                    −
                    
                      M
                      
                        2
                      
                    
                  
                  )
                
                
                  2
                
              
              +
              
                Γ
                
                  2
                
              
               
              
                M
                
                  2
                
              
               
            
          
        
         
        ,
      
    
    {\displaystyle {\frac {\ f'\!(\mathrm {E} )\ }{\ f(\mathrm {E} )\ }}~=~{\frac {\ 4\ \left(M^{2}-\mathrm {E} ^{2}\right)\ \mathrm {E} \ }{\ \left(\ \mathrm {E} ^{2}-M^{2}\right)^{2}+\Gamma ^{2}\ M^{2}\ }}\ ,}

with

f
        (
        M
        )
        =
        
          
            k
            
               
              
                Γ
                
                  2
                
              
               
              
                M
                
                  2
                
              
               
            
          
        
         
        .
      
    
    {\displaystyle f(M)={\frac {k}{~\Gamma ^{2}\ M^{2}\ }}~.}

Further information: <a href="/facts/Cauchy_distribution/3HS6D3Xk">Cauchy distribution</a>
<h2 id="gaussian-broadening">Gaussian broadening</h2>
In experiment, the incident beam that produces resonance always has some spread of energy around a central value. Usually, that is a <a href="/facts/Normal_distribution/UapjjPyQ">Gaussian/normal distribution</a>. The resulting resonance shape in this case is given by the <a href="/facts/Convolution/PVPrdz9J">convolution</a> of the Breit–Wigner and the Gaussian distribution,

V
          
            2
          
        
        (
        E
        ;
        M
        ,
        Γ
        ,
        k
        ,
        σ
        )
        =
        
          ∫
          
            −
            ∞
          
          
            ∞
          
        
        
          
            k
            
               
              (
              
                E
                
                  ′
                  
                    2
                  
                
              
              −
              
                M
                
                  2
                
              
              
                )
                
                  2
                
              
              +
              (
              M
               
              Γ
              
                )
                
                  2
                
              
               
            
          
        
         
        
          
            1
            
               
              σ
              
                
                  2
                  π
                   
                
              
               
            
          
        
         
        
          e
          
            −
            
              
                
                   
                  (
                  
                    E
                    ′
                  
                  −
                  E
                  
                    )
                    
                      2
                    
                  
                
                
                   
                  2
                  
                    σ
                    
                      2
                    
                  
                   
                
              
            
          
        
        d
        ⁡
        
          E
          ′
        
         
        .
      
    
    {\displaystyle V_{2}(E;M,\Gamma ,k,\sigma )=\int _{-\infty }^{\infty }{\frac {k}{\ (E'^{2}-M^{2})^{2}+(M\ \Gamma )^{2}\ }}\ {\frac {1}{\ \sigma {\sqrt {2\pi \ }}\ }}\ e^{-{\frac {\ (E'-E)^{2}}{\ 2\sigma ^{2}\ }}}\operatorname {d} E'~.}

This function can be simplified <a class="footnote-ref" id="fnref:5" href="#fn:5">5</a> by introducing new variables,

t
        =
        
          
            
               
              E
              −
              
                E
                ′
              
               
            
            
               
              
                
                  2
                   
                
              
              σ
               
            
          
        
         
        ,
        
        
          u
          
            1
          
        
        =
        
          
            
               
              E
              −
              M
               
            
            
               
              
                
                  2
                   
                
              
              σ
               
            
          
        
         
        ,
        
        
          u
          
            2
          
        
        =
        
          
            
               
              E
              +
              M
               
            
            
               
              
                
                  2
                   
                
              
              σ
               
            
          
        
         
        ,
        
        a
        =
        
          
            
               
              k
               
              π
               
            
            
               
              2
               
              
                σ
                
                  2
                
              
               
            
          
        
         
        ,
      
    
    {\displaystyle \ t={\frac {\ E-E'\ }{\ {\sqrt {2\ }}\sigma \ }}\ ,\quad u_{1}={\frac {\ E-M\ }{\ {\sqrt {2\ }}\sigma \ }}\ ,\quad u_{2}={\frac {\ E+M\ }{\ {\sqrt {2\ }}\sigma \ }}\ ,\quad a={\frac {\ k\ \pi \ }{~2\ \sigma ^{2}\ }}\ ,}

to obtain

V
          
            2
          
        
        (
        E
        ;
        M
        ,
        Γ
        ,
        k
        ,
        σ
        )
        =
        
          
            
               
              
                H
                
                  2
                
              
              (
              a
              ,
              
                u
                
                  1
                
              
              ,
              
                u
                
                  2
                
              
              )
               
            
            
               
              
                σ
                
                  2
                
              
              2
              
                
                  π
                   
                
              
               
            
          
        
         
        ,
      
    
    {\displaystyle \ V_{2}(E;M,\Gamma ,k,\sigma )={\frac {\ H_{2}(a,u_{1},u_{2})\ }{\ \sigma ^{2}2{\sqrt {\pi \ }}\ }}\ ,}

where the relativistic line broadening function <a class="footnote-ref" id="fnref:6" href="#fn:6">6</a> has the following definition,

H
          
            2
          
        
        (
        a
        ,
        
          u
          
            1
          
        
        ,
        
          u
          
            2
          
        
        )
        =
        
          
            
               
              a
               
            
            
               
              π
               
            
          
        
        
          ∫
          
            −
            ∞
          
          
            ∞
          
        
        
          
            
               
              
                e
                
                  −
                  
                    t
                    
                      2
                    
                  
                
              
            
            
               
              (
              
                u
                
                  1
                
              
              −
              t
              
                )
                
                  2
                
              
              (
              
                u
                
                  2
                
              
              −
              t
              
                )
                
                  2
                
              
              +
              
                a
                
                  2
                
              
               
            
          
        
        d
        ⁡
        t
         
        .
      
    
    {\displaystyle \ H_{2}(a,u_{1},u_{2})={\frac {\ a\ }{\ \pi \ }}\int _{-\infty }^{\infty }{\frac {~e^{-t^{2}}}{~(u_{1}-t)^{2}(u_{2}-t)^{2}+a^{2}\ }}\operatorname {d} t~.}

H
 
 2
 
 
 
 
 {\displaystyle H_{2}}
 
 is the relativistic counterpart of the similar line-broadening function <a class="footnote-ref" id="fnref:7" href="#fn:7">7</a> for the <a href="/facts/Voigt_profile/uwQXjc4k">Voigt profile</a> used in spectroscopy (see also § 7.19 of <a class="footnote-ref" id="fnref:8" href="#fn:8">8</a>).

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Breit, G.; Wigner, E. (1936). "Capture of slow neutrons". Physical Review. 49 (7): 519. Bibcode:1936PhRv...49..519B. doi:10.1103/PhysRev.49.519. <a href="/wiki/Physical_Review" target="_blank">/wiki/Physical_Review</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">See Pythia 6.4 Physics and Manual (page 98 onwards) for a discussion of the widths of particles in the PYTHIA manual. Note that this distribution is usually represented as a function of the squared energy. <a href="https://ftp.fau.de/gentoo/distfiles/lutp0613man2.pdf" target="_blank">https://ftp.fau.de/gentoo/distfiles/lutp0613man2.pdf</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Bohm, A.; Sato, Y. (2005). "Relativistic resonances: Their masses, widths, lifetimes, superposition, and causal evolution". Physical Review D. 71 (8): 085018. arXiv:hep-ph/0412106. Bibcode:2005PhRvD..71h5018B. doi:10.1103/PhysRevD.71.085018. S2CID 119417992. <a href="/wiki/ArXiv_(identifier)" target="_blank">/wiki/ArXiv_(identifier)</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">Brown, L.S. (1994). Quantum Field Theory. Cambridge University Press. § 6.3. ISBN 978-0521469463. <a href="978-0521469463" target="_blank">978-0521469463</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">Kycia, Radosław A.; Jadach, Stanisław (15 July 2018). "Relativistic Voigt profile for unstable particles in high energy physics". Journal of Mathematical Analysis and Applications. 463 (2): 1040–1051. arXiv:1711.09304. doi:10.1016/j.jmaa.2018.03.065. ISSN 0022-247X. S2CID 78086748. <a href="http://www.sciencedirect.com/science/article/pii/S0022247X18302774" target="_blank">http://www.sciencedirect.com/science/article/pii/S0022247X18302774</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">Kycia, Radosław A.; Jadach, Stanisław (15 July 2018). "Relativistic Voigt profile for unstable particles in high energy physics". Journal of Mathematical Analysis and Applications. 463 (2): 1040–1051. arXiv:1711.09304. doi:10.1016/j.jmaa.2018.03.065. ISSN 0022-247X. S2CID 78086748. <a href="http://www.sciencedirect.com/science/article/pii/S0022247X18302774" target="_blank">http://www.sciencedirect.com/science/article/pii/S0022247X18302774</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Finn, G.D.; Mugglestone, D. (1 February 1965). "Tables of the line broadening function H(a,v)". Monthly Notices of the Royal Astronomical Society. 129 (2): 221–235. doi:10.1093/mnras/129.2.221. ISSN 0035-8711. <a href="https://academic.oup.com/mnras/article/129/2/221/2604178" target="_blank">https://academic.oup.com/mnras/article/129/2/221/2604178</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">NIST Handbook of Mathematical Functions. U.S. National Institute of Standards and Technology. Olver, Frank W.J. Cambridge, UK: Cambridge University Press. 2010. ISBN 978-0-521-19225-5. OCLC 502037224.{{cite book}}: CS1 maint: others (link) <a href="978-0-521-19225-5" target="_blank">978-0-521-19225-5</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
</ol>

Relativistic Breit–Wigner distribution open-in-new

Relativistic Breit–Wigner distribution