Flux balance analysis: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
Added new citation for R Software package sybil.
en>BG19bot
m WP:CHECKWIKI error fix for #61. Punctuation goes before References. Do general fixes if a problem exists. - using AWB (9916)
Line 1: Line 1:
The '''Lippmann-Schwinger equation''' (named after [[Bernard A. Lippmann]] and [[Julian Schwinger]]<ref>{{cite journal | last=Lippmann|first=Bernard A.|coauthors=Schwinger, Julian|journal=Physical Review Letters|year=1950|volume=79|pages=469|doi = 10.1103/PhysRev.79.469}}</ref>) is one of the most used equations to describe particle colisions - or, more precisely, [[scattering]] - in [[quantum mechanics]]. It may be used in scattering of molecules, atoms, neutrons, photons or any other particles and is important mainly in [[atomic, molecular, and optical physics]], [[nuclear physics]] and [[particle physics]]. It relates the scattered wave function with the interaction that produces the scattering (the scattering potential) and allows to calculate the relevant experimental parameters ([[scattering amplitude]] and [[cross section (physics)|cross section]]s).
Hello, my title is Andrew and my wife doesn't like it at all. For years he's been living in Alaska and he doesn't plan on changing it. Office supervising is where her primary earnings arrives from. The preferred pastime for him and his kids is style and he'll [http://www.edmposts.com/build-a-beautiful-organic-garden-using-these-ideas/ good psychic] be beginning some thing else alongside with it.<br><br>Also visit my web blog free psychic ([http://www.january-yjm.com/xe/index.php?mid=video&document_srl=158289 january-yjm.com]) :: authentic psychic readings ([http://c045.danah.co.kr/home/index.php?document_srl=1356970&mid=qna c045.danah.co.kr])
 
The most fundamental equation to describe any quantum phenomenon, including scattering, is the [[Schrödinger equation]]. In physical problems, this [[differential equation]] must be solved with the input of an additional set of initial and/or [[boundary condition]]s for the specific physical system studied. The Lippmann–Schwinger equation is equivalent to the Schrödinger equation plus the typical boundary conditions for scattering problems. In order to embed the boundary conditions, the Lippmann–Schwinger equation must be written as an [[integral equation]].<ref>[[#Joachain|Joachain, Charles J., 1983]] page 112</ref> For scattering problems, the Lippmann–Schwinger equation has been proved to be mathematically and intuitively more convenient than the original Schrödinger equation.
 
The Lippmann-Schwinger equation general shape is (in reality, two equations are shown below, one for the <math> + \,</math> sign and other for the <math> - \,</math> sign):
:<math> | \psi^{(\pm)} \rangle = | \phi \rangle + \frac{1}{E - H_0 \pm i \epsilon} V |\psi^{(\pm)} \rangle. \,</math>
 
In the equations above, <math> | \psi^{(+)} \rangle \,</math> is the wave function of the whole system (the two colliding systems considered as a whole) at an infinite time before the interaction; and <math> | \psi^{(-)} \rangle \,</math>, at an infinite time after the interaction (the "scattered wave fuction"). The potential energy <math> V \,</math> describes the interaction between the two colliding systems. The [[hamiltonian]] <math> H_0 \,</math> describes the situation in which the two systems are infinitely far apart and do not interact. Its [[eigenfunction]]s are <math> | \phi \rangle \,</math> and its [[eigenvalue]]s are the energies <math> E \,</math>. Finally, <math> i \epsilon \,</math> is a mathematical technicality necessary for the calculation of the integrals needed to solve the equation and has no physical meaning.
 
==Usage==
The Lippmann-Schwinger equation is useful in a very large number of situations involving two-body scattering. For three or more colliding bodies it does not work well because of mathematical limitations; Faddeev equations may be used instead.<ref>[[#Joachain|Joachain, Charles J., 1983]] page 517</ref> However, there are approximations that can reduce a [[many-body problem]] to a set of [[two-body problem]]s in a variety of cases. For example, in a collision between electrons and molecules, there may be tens or hundreds of particles involved. But the phenomenum may be reduced to a two-body problem by describing all the molecule constituent particle potentials together with a [[pseudopotential]].<ref>[[#Joachain|Joachain, Charles J., 1983]] page 576</ref> In these cases, the Lippmann-Schwinger equations may be used. Of course, the main motivations of these approaches are also the possibility of doing the calculations with much lower computational efforts.
 
==Derivation==
We will assume that the [[Hamiltonian (quantum mechanics)|Hamiltonian]] may be written as
:<math>H = H_0 + V \,</math>
where ''H'' and ''H''<sub>0</sub> have the same [[eigenvalues]] and ''H''<sub>0</sub> is a free Hamiltonian.  For example in nonrelativistic quantum mechanics ''H''<sub>0</sub> may be
:<math>H_0 = \frac{p^2}{2m}. \,</math>
 
Intuitively <math>V \,</math> is the interaction energy of the system.  This analogy is somewhat misleading, as interactions generically change the energy levels ''E'' of steady states of the system, but ''H'' and ''H''<sub>0</sub> have identical spectra ''E''<sub>α</sub>.  This means that, for example, a bound state that is an eigenstate of the interacting Hamiltonian will also be an eigenstate of the free Hamiltonian. This is in contrast with the Hamiltonian obtained by turning off all interactions, in which case there would be no bound states.  Thus one may think of ''H''<sub>0</sub> as the free Hamiltonian for the boundstates with effective parameters that are determined by the interactions.
 
Let there be an [[eigenstate]] of <math>H_0 \,</math>:
:<math>H_0 | \phi \rangle = E | \phi \rangle. \,</math>
 
Now if we add the interaction <math> V \,</math> into the mix, we need to solve
:<math>\left( H_0 + V \right) | \psi \rangle = E | \psi \rangle. \,</math>
 
Because of the continuity of the energy eigenvalues, we wish that <math> | \psi \rangle \to | \phi \rangle \,</math> as <math> V \to 0 \,</math>.
 
A potential solution to this situation is
:<math> |\psi \rangle = | \phi \rangle + \frac{1}{E - H_0} V | \psi \rangle. \,</math>
 
However <math>(E-H_0)</math>  is [[Mathematical singularity|singular]] since <math> E </math> is an eigenvalue of <math> H_0 </math>.
 
As is described below, this singularity is eliminated in two distinct ways by making the denominator slightly complex:
:<math> | \psi^{(\pm)} \rangle = | \phi \rangle + \frac{1}{E - H_0 \pm i \epsilon} V |\psi^{(\pm)} \rangle. \,</math>
 
==Interpretation as in and out states==
 
===The S-matrix paradigm===
 
In the [[S-matrix]] formulation of [[particle physics]], which was pioneered by [[John Archibald Wheeler]] among others, all physical processes are modeled according to the following paradigm.
 
One begins with a non-interacting multiparticle state in the distant past.  Non-interacting does not mean that all of the forces have been turned off, in which case for example [[proton]]s would fall apart, but rather that there exists an interaction-free [[Hamiltonian (quantum mechanics)|Hamiltonian]] ''H''<sub>0</sub>, for which the bound states have the same energy level spectrum as the actual Hamiltonian ''H''.  This initial state is referred to as the '''in''' state.  Intuitively, it consists of bound states that are sufficiently well separated that their interactions with each other are ignored.
 
The idea is that whatever physical process one is trying to study may be modeled as a [[scattering]] process of these well separated bound states.  This process is described by the full Hamiltonian ''H'', but once it's over, all of the new bound states separate again and one finds a new noninteracting state called the '''out''' state. The S-matrix is more symmetric under relativity than the Hamiltonian, because it does not require a choice of time slices to define.
 
This paradigm allows one to calculate the probabilities of all of the processes that we have observed in 70 years of particle collider experiments with remarkable accuracy.  But many interesting physical phenomena do not obviously fit into this paradigm. For example, if one wishes to consider the dynamics inside of a neutron star sometimes one wants to know more than what it will finally decay into.  In other words, one may be interested in measurements that are not in the asymptotic future.  Sometimes an asymptotic past or future is not even available.  For example, it is very possible that there is no past before the [[big bang]].
 
In the 1960s, the S-matrix paradigm was elevated by many physicists to a fundamental law of nature. In [[S-matrix theory]], it was stated that any quantity that one could measure should be found in the S-matrix for some process. This idea was inspired by the physical interpretation that S-matrix techniques could give to [[Feynman diagrams]] restricted to the [[mass-shell]], and led to the construction of [[dual resonance model]]s. But it was very controversial, because it denied the validity of [[quantum field theory]] based on local fields and Hamiltonians.
 
===The connection to Lippmann–Schwinger===
 
Intuitively, the slightly deformed eigenfunctions <math> \psi^{(\pm)}</math> of the full Hamiltonian ''H'' are the in and out states. The <math>\phi</math> are noninteracting states that resemble the '''in''' and '''out''' states in the infinite past and infinite future.
 
===Creating wavepackets===
 
This intuitive picture is not quite right, because  <math> \psi^{(\pm)}</math> is an eigenfunction of the Hamiltonian and so at different times only differs by a phase. Thus, in particular, the physical state does not evolve and so it cannot become noninteracting.  This problem is easily circumvented by assembling <math> \psi^{(\pm)}</math> and <math> \phi</math> into wavepackets with some distribution <math>g(E)</math> of energies <math>E</math> over a characteristic scale <math>\Delta E</math>.  The [[uncertainty principle]] now allows the interactions of the asymptotic states to occur over a timescale <math>\hbar/\Delta E</math> and in particular it is no longer inconceivable that the interactions may turn off outside of this interval.  The following argument suggests that this is indeed the case.
 
Plugging the Lippmann–Schwinger equations into the definitions
 
::<math> \psi^{(\pm)}_g(t)=\int dE\, e^{-iEt} g(E)\psi^{(\pm)}</math>
 
and
 
::<math> \phi_g(t)=\int dE\, e^{-iEt} g(E)\phi</math>
 
of the wavepackets we see that, at a given time, the difference between the <math>\psi_g(t)</math> and <math>\phi_g(t)</math> wavepackets is given by an integral over the energy ''E''.
 
===A contour integral===
 
This integral may be evaluated by defining the wave function over the complex ''E'' plane and closing the ''E'' contour using a semicircle on which the wavefunctions vanish.  The integral over the closed contour may then be evaluated, using the [[Cauchy integral theorem]], as a sum of the residues at the various poles.  We will now argue that the residues of <math> \psi^{(\pm)}</math> approach those of <math> \phi</math> at time <math>t\rightarrow\mp\infty</math> and so the corresponding wavepackets are equal at temporal infinity.
 
In fact, for very positive times ''t'' the <math>e^{-iEt}</math> factor in a [[Schrödinger picture]] state forces one to close the contour on the lower half-plane. The pole in the ''V'' from the Lippmann–Schwinger equation reflects the time-uncertainty of the interaction, while that in the wavepackets weight function reflects the duration of the interaction.  Both of these varieties of poles occur at finite imaginary energies and so are suppressed at very large times.  The pole in the energy difference in the denominator is on the upper half-plane in the case of <math> \psi^{-}</math>, and so does not lie inside the integral contour and does not contribute to the <math> \psi^{-}</math> integral. The remainder is equal to the <math>\phi</math> wavepacket.  Thus, at very late times <math> \psi^{-}=\phi</math>, identifying <math> \psi^{-}</math> as the asymptotic noninteracting '''out''' state.
 
Similarly one may integrate the wavepacket corresponding to <math> \psi^{+}</math> at very negative times.  In this case the contour needs to be closed over the upper half-plane, which therefore misses the energy pole of <math> \psi^{+}</math>, which is in the lower half-plane.  One then finds that the <math> \psi^{+}</math>
and <math> \phi</math> wavepackets are equal in the asymptotic past, identifying <math> \psi^{+}</math> as the asymptotic noninteracting '''in''' state.
 
===The complex denominator of Lippmann–Schwinger===
 
This identification of the <math>\psi</math>'s as asymptotic states is the justification for the <math>\pm\epsilon</math> in the denominator of the Lippmann–Schwinger equations.
 
==A formula for the S-matrix==
 
The [[S-matrix]] ''S'' is defined to be the inner product
 
::<math> S_{ab}=(\psi^-_a,\psi^+_b)</math>
 
of the ''a''th and ''b''th [[Heisenberg picture]] asymptotic states.  One may obtain a formula relating the ''S''-matrix to the potential ''V'' using the above contour integral strategy, but this time switching the roles of <math> \psi^+</math> and <math> \psi^-</math>.  As a result, the contour now does pick up the energy pole.  This can be related to the <math>\phi</math>'s if one uses the S-matrix to swap the two <math>\psi</math>'s.  Identifying the coefficients of the <math>\phi</math>'s on both sides of the equation one finds the desired formula relating ''S'' to the potential
 
:<math> S_{ab}=\delta(a-b)-2i\pi\delta(E_a-E_b)(\phi_a,V\psi^+_b).</math>
 
In the [[Born approximation]], corresponding to first order [[perturbation theory]], one replaces this last <math> \psi^+</math> with the corresponding eigenfunction <math> \phi</math> of the free Hamiltonian ''H''<sub>0</sub>, yielding
 
:<math> S_{ab}=\delta(a-b)-2i\pi\delta(E_a-E_b)(\phi_a,V\phi_b)\,</math>
which expresses the S-matrix entirely in terms of ''V'' and free Hamiltonian eigenfunctions.
 
These formulas may in turn be used to calculate the reaction rate of the process <math>b\rightarrow a</math>, which is equal to <math>|S_{ab}-\delta_{ab}|^2.\,</math>
 
==Homogenization==
With the use of Green's function, the Lippmann–Schwinger equation has counterparts in homogenization theory (e.g. mechanics, conductivity, permittivity).
 
== See also ==
*[[Bethe-Salpeter equation]]
 
==References==
 
{{reflist|30em}}
 
==Bibliography==
 
*{{cite book | author=Joachain, Charles J. | title=Quantum collision theory | publisher=North Holland | year=1983 | url=http://pt.scribd.com/doc/76844441/Charles-J-Joachain-Quantum-Collision-Theory | isbn=0-7204-0294-8 | ref=Joachain}}
 
*{{cite book | author=Sakurai, J. J. | title=Modern Quantum Mechanics | publisher=Addison Wesley | year=1994 | isbn=0-201-53929-2}}
 
*{{cite book|authorlink=Steven Weinberg|author=Weinberg, S.|title=The Quantum Theory of Fields|publisher=Cambridge University Press|year=1995|isbn=0-521-67053-5}}
 
{{DEFAULTSORT:Lippmann-Schwinger equation}}
[[Category:Scattering]]

Revision as of 06:53, 6 February 2014

Hello, my title is Andrew and my wife doesn't like it at all. For years he's been living in Alaska and he doesn't plan on changing it. Office supervising is where her primary earnings arrives from. The preferred pastime for him and his kids is style and he'll good psychic be beginning some thing else alongside with it.

Also visit my web blog free psychic (january-yjm.com) :: authentic psychic readings (c045.danah.co.kr)