|
|
| Line 1: |
Line 1: |
| In [[many-body theory]], the term '''Green's function''' (or '''Green function''') is sometimes used interchangeably with [[Correlation function (quantum field theory)|correlation function]], but refers specifically to correlators of [[field operator]]s or [[creation and annihilation operators]].
| | It's easy to make money online. There is truth to the fact that you can start making money on the Internet as soon as you're done with this article. After all, so many others are making money online, why not you? |
|
| |
|
| The name comes from the [[Green's functions]] used to solve inhomogeneous [[differential equations]], to which they are loosely related. (Specifically, only two-point 'Green's functions' in the case of a non-interacting system are Green's functions in the mathematical sense; the linear operator that they invert is the [[Hamiltonian (quantum mechanics)|Hamiltonian operator]], which in the non-interacting case is quadratic in the fields.)
| | Keep your mind open and you can make a lot of money. As you [http://www.comoganhardinheiro101.com/opcoes-binarias/ ganhando dinheiro na internet] can see, there are many ways to approach the world of online income. With various streams of income available, you are sure to find one, or two, that can help you with your income needs. Take this como ganhar dinheiro na internet [http://comoconseguirdinheiro.comoganhardinheiro101.com como conseguir dinheiro] information to heart, put it to use and build your own online success story. Here's more information [http://comoganhardinheironainternet.comoganhardinheiro101.com ganhando dinheiro na internet] in regards to [http://comoganhardinheiropelainternet.comoganhardinheiro101.com/ como ganhar dinheiro] look into comoganhardinheiropelainternet.comoganhardinheiro101.com/ <br><br><br>fee to watch your webinar at their convenience. Once it is in place, [http://www.comoganhardinheiro101.com/?p=66 como ganhar dinheiro] promotion and possibly answering questions will be your only tasks.<br><br>Getting [http://ganhardinheironainternet.comoganhardinheiro101.com/ paid money] to work online isn't the easiest thing to do in the world, but it is possible. If this is something you wish to work with, then the tips presented above should have helped you. |
|
| |
|
| ==Spatially uniform case== | | Take some time, do things the right way and then you can succeed. Start your online [http://www.comoganhardinheiro101.com/ferramentas/ ganhe dinheiro na internet] earning today by following the great advice discussed in this article. Earning money is not as hard as it may seem, you just need to know how to get started. By choosing to put your right foot forward, you are heading off to a great start earning money to make ends meet.<br><br>Make money online by selling your talents. Good music is always in demand and with today's technological advances, anyone with musical talent can make music and offer it for sale [http://www.comoganhardinheiro101.com/?p=73 como ganhar dinheiro na internet] to a broad audience. |
|
| |
|
| ===Basic definitions===
| | By setting up your own website and using social media for promotion, you can share your music with others and sell downloads with a free PayPal account. |
| | |
| We consider a many-body theory with field operator (annihilation operator written in the position basis) <math>\psi(\mathbf{x})</math>.
| |
| | |
| The [[Heisenberg picture|Heisenberg operators]] can be written in terms of [[Schrödinger picture|Schrödinger operators]] as
| |
| :<math>
| |
| \psi(\mathbf{x},t) = \mathrm{e}^{\mathrm{i} K t} \psi(\mathbf{x}) \mathrm{e}^{-\mathrm{i} K t},
| |
| </math>
| |
| and <math>\bar\psi(\mathbf{x},t) = [\psi(\mathbf{x},t)]^\dagger</math>, where <math>K = H - \mu N</math> is the [[Grand canonical ensemble|grand-canonical]] Hamiltonian.
| |
| | |
| Similarly, for the [[Imaginary time|imaginary-time]] operators,
| |
| :<math>
| |
| \psi(\mathbf{x},\tau) = \mathrm{e}^{K \tau} \psi(\mathbf{x}) \mathrm{e}^{-K\tau}
| |
| </math>
| |
| :<math>
| |
| \bar\psi(\mathbf{x},\tau) = \mathrm{e}^{K \tau} \psi^\dagger(\mathbf{x}) \mathrm{e}^{-K\tau}.
| |
| </math>
| |
| (Note that the imaginary-time creation operator <math>\bar\psi(\mathbf{x},\tau)</math> is not the [[Hermitian conjugate]] of the annihilation operator <math>\psi(\mathbf{x},\tau)</math>.)
| |
| | |
| In real time, the <math>2n</math>-point Green function is defined by
| |
| :<math>
| |
| G^{(n)}(1 \ldots n | 1' \ldots n')
| |
| = \mathrm{i}^n \langle T\psi(1)\ldots\psi(n)\bar\psi(n')\ldots\bar\psi(1')\rangle,
| |
| </math>
| |
| where we have used a condensed notation in which <math>j</math> signifies <math>\mathbf{x}_j, t_j</math> and <math>j'</math> signifies <math>\mathbf{x}_j', t_j'</math>. The operator <math>T</math> denotes [[time ordering]], and indicates that the field operators that follow it are to be ordered so that the their time arguments increase from right to left.
| |
| | |
| In imaginary time, the corresponding definition is
| |
| :<math>
| |
| \mathcal{G}^{(n)}(1 \ldots n | 1' \ldots n')
| |
| = \langle T\psi(1)\ldots\psi(n)\bar\psi(n')\ldots\bar\psi(1')\rangle,
| |
| </math>
| |
| where <math>j</math> signifies <math>\mathbf{x}_j, \tau_j</math>. (The imaginary-time variables <math>\tau_j</math> are restricted to the range from <math>0</math> to the inverse temperature <math>\beta=\frac{1}{k_{B}T}</math>.)
| |
| | |
| '''Note''' regarding signs and normalization used in these definitions: The signs of the Green functions have been chosen so that [[Fourier transform]] of the two-point (<math>n=1</math>) thermal Green function for a free particle is
| |
| :<math>
| |
| \mathcal{G}(\mathbf{k},\omega_n) = \frac{1}{-\mathrm{i}\omega_n + \xi_\mathbf{k}},
| |
| </math>
| |
| and the retarded Green function is
| |
| :<math>
| |
| G^{\mathrm{R}}(\mathbf{k},\omega) = \frac{1}{-(\omega+\mathrm{i}\eta) + \xi_\mathbf{k}},
| |
| </math>
| |
| where
| |
| :<math>
| |
| \omega_n = {[2n+\theta(-\zeta)]\pi}/{\beta}
| |
| </math>
| |
| is the [[Matsubara frequency]].
| |
| | |
| Throughout, <math>\zeta</math> is <math>+1</math> for [[boson]]s and <math>-1</math> for [[fermion]]s and <math>[\ldots,\ldots]=[\ldots,\ldots]_{-\zeta}</math> denotes either a [[commutator]] or anticommutator as appropriate.
| |
| | |
| (See [[Green's function (many-body theory)#Non-interacting case|below]] for details.)
| |
| | |
| ===Two-point functions===
| |
| | |
| The Green function with a single pair of arguments (<math>n=1</math>) is referred to as the two-point function, or propagator. In the presence of both spatial and temporal translational symmetry, it depends only on the difference of its arguments. Taking the Fourier transform with respect to both space and time gives
| |
| :<math>
| |
| \mathcal{G}(\mathbf{x}\tau|\mathbf{x}'\tau') = \int_\mathbf{k} d\mathbf{k} \frac{1}{\beta}\sum_{\omega_n} \mathcal{G}(\mathbf{k},\omega_n) \mathrm{e}^{\mathrm{i} \mathbf{k}\cdot(\mathbf{x}-\mathbf{x}')-\mathrm{i}\omega_n (\tau-\tau')},
| |
| </math>
| |
| where the sum is over the appropriate [[Matsubara frequency|Matsubara frequencies]] (and the integral involves an implicit factor of <math>(L/2\pi)^{d}</math>, as usual).
| |
| | |
| In real time, we will explicitly indicate the time-ordered function with a superscript T:
| |
| :<math>
| |
| G^{\mathrm{T}}(\mathbf{x} t|\mathbf{x}' t') = \int_\mathbf{k} d \mathbf{k} \int \frac{\mathrm{d}\omega}{2\pi} G^{\mathrm{T}}(\mathbf{k},\omega) \mathrm{e}^{\mathrm{i} \mathbf{k}\cdot(\mathbf{x} -\mathbf{x} ')-\mathrm{i}\omega(t-t')}.
| |
| </math>
| |
| | |
| The real-time two-point Green function can be written in terms of `retarded' and `advanced' Green functions, which will turn out to have simpler analyticity properties. The retarded and advanced Green functions are defined by
| |
| :<math>
| |
| G^{\mathrm{R}}(\mathbf{x} t|\mathbf{x} 't') = \mathrm{i}\langle[\psi(\mathbf{x} ,t),\bar\psi(\mathbf{x} ',t')]\rangle\Theta(t-t')
| |
| </math>
| |
| and
| |
| :<math>
| |
| G^{\mathrm{A}}(\mathbf{x} t|\mathbf{x} 't') = -\mathrm{i}\langle[\psi(\mathbf{x} ,t),\bar\psi(\mathbf{x} ',t')]\rangle\Theta(t'-t),
| |
| </math>
| |
| respectively.
| |
| | |
| They are related to the time-ordered Green function by
| |
| :<math>
| |
| G^{\mathrm{T}}(\mathbf{k},\omega) = [1+\zeta n(\omega)]G^{\mathrm{R}}(\mathbf{k},\omega) - \zeta n(\omega) G^{\mathrm{A}}(\mathbf{k},\omega),
| |
| </math>
| |
| where
| |
| :<math>
| |
| n(\omega) = \frac{1}{\mathrm{e}^{\beta \omega}-\zeta}
| |
| </math>
| |
| is the [[Bose-Einstein statistics|Bose-Einstein]] or [[Fermi–Dirac statistics|Fermi-Dirac]] distribution function.
| |
| | |
| ====Imaginary-time ordering and <math>\beta</math>-periodicity====
| |
| | |
| The thermal Green functions are defined only when both imaginary-time arguments are within the range <math>0</math> to <math>\beta</math>. The two-point Green function has the following properties. (The position or momentum arguments are suppressed in this section.)
| |
| | |
| Firstly, it depends only on the difference of the imaginary times:
| |
| :<math>
| |
| \mathcal{G}(\tau,\tau') = \mathcal{G}(\tau - \tau').
| |
| </math>
| |
| The argument <math>\tau - \tau'</math> is allowed to run from <math>-\beta</math> to <math>\beta</math>.
| |
| | |
| Secondly, <math>\mathcal{G}(\tau)</math> is periodic under shifts of <math>\beta</math>. Because of the small domain within which the function is defined, this means just
| |
| :<math>
| |
| \mathcal{G}(\tau - \beta) = \zeta \mathcal{G}(\tau),
| |
| </math>
| |
| for <math>0 < \tau < \beta</math>. (Note that the function is antiperiodic for fermions.) Time ordering is crucial for this property, which can be proved straightforwardly, using the cyclicity of the trace operation.
| |
| | |
| These two properties allow for the Fourier transform representation and its inverse,
| |
| :<math>
| |
| \mathcal{G}(\omega_n) = \int_0^\beta \mathrm{d}\tau \, \mathcal{G}(\tau)\, \mathrm{e}^{\mathrm{i}\omega_n \tau}.
| |
| </math>
| |
| | |
| Finally, note that <math>\mathcal{G}(\tau)</math> has a discontinuity at <math>\tau = 0</math>; this is consistent with a long-distance behaviour of <math>\mathcal{G}(\omega_n) \sim 1/|\omega_n|</math>.
| |
| | |
| ===Spectral representation===
| |
| | |
| The propagators in real and imaginary time can both be related to the spectral density (or spectral weight), given by
| |
| :<math>
| |
| \rho(\mathbf{k},\omega) = \frac{1}{\mathcal{Z}}\sum_{\alpha,\alpha'} 2\pi \delta(E_\alpha-E_{\alpha'}-\omega)\;
| |
| |\langle\alpha|\psi_\mathbf{k}^\dagger|\alpha'\rangle|^2\left(\mathrm{e}^{-\beta E_{\alpha'}}-\zeta\mathrm{e}^{-\beta E_{\alpha}}\right),
| |
| </math>
| |
| where <math>|\alpha \rangle</math> refers to a (many-body) eigenstate of the grand-canonical Hamiltonian <math>H-\mu N</math>, with eigenvalue <math>E_\alpha</math>.
| |
| | |
| The imaginary-time propagator is then given by
| |
| :<math>
| |
| \mathcal{G}(\mathbf{k},\omega_n) = \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega'}{2\pi}
| |
| \frac{\rho(\mathbf{k},\omega')}{-\mathrm{i}\omega_n+\omega'}.
| |
| </math>
| |
| and the retarded propagator by | |
| :<math>
| |
| G^{\mathrm{R}}(\mathbf{k},\omega) = \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega'}{2\pi}
| |
| \frac{\rho(\mathbf{k},\omega')}{-(\omega+\mathrm{i}\eta)+\omega'},
| |
| </math>
| |
| where the limit as <math>\eta\rightarrow 0^+</math> is implied.
| |
| | |
| The advanced propagator is given by the same expression, but with <math>-\mathrm{i}\eta</math> in the denominator. The time-ordered function can be found in terms of <math>G^{\mathrm{R}}</math> and <math>G^{\mathrm{A}}</math>. As claimed above, <math>G^{\mathrm{R}}(\omega)</math> and <math>G^{\mathrm{A}}(\omega)</math> have simple analyticity properties: the former (latter) has all its poles and discontinuities in the lower (upper) half-plane. The thermal propagator <math>\mathcal{G}(\omega_n)</math> has all its poles and discontinuities on the imaginary <math>\omega_n</math> axis.
| |
| | |
| The spectral density can be found very straightforwardly from <math>G^{\mathrm{R}}</math>, using the [[Sokhatsky–Weierstrass theorem]]
| |
| :<math>
| |
| \lim_{\eta\rightarrow 0^+}\frac{1}{x\pm\mathrm{i}\eta} = {P}\frac{1}{x}\mp i\pi\delta(x),
| |
| </math>
| |
| where <math>P</math> denotes the [[Cauchy principal part]].
| |
| This gives
| |
| :<math>
| |
| \rho(\mathbf{k},\omega) = 2\mathrm{Im}\, G^{\mathrm{R}}(\mathbf{k},\omega).
| |
| </math>
| |
| | |
| This furthermore implies that <math>G^{\mathrm{R}}(\mathbf{k},\omega)</math> obeys the following relationship between its real and imaginary parts:
| |
| :<math>
| |
| \mathrm{Re}\, G^{\mathrm{R}}(\mathbf{k},\omega) = -2 P \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega'}{2\pi}
| |
| \frac{\mathrm{Im}\, G^{\mathrm{R}}(\mathbf{k},\omega')}{\omega-\omega'},
| |
| </math>
| |
| where <math>P</math> denotes the principal value of the integral.
| |
| | |
| The spectral density obeys a sum rule:
| |
| :<math>
| |
| \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega}{2\pi} \rho(\mathbf{k},\omega) = 1,
| |
| </math>
| |
| which gives
| |
| :<math>
| |
| G^{\mathrm{R}}(\omega)\sim\frac{1}{|\omega|}
| |
| </math>
| |
| as <math>|\omega| \rightarrow \infty</math>.
| |
| | |
| ====Hilbert transform====
| |
| | |
| The similarity of the spectral representations of the imaginary- and real-time Green functions allows us to define the function
| |
| :<math>
| |
| G(\mathbf{k},z) = \int_{-\infty}^\infty \frac{\mathrm{d} x}{2\pi} \frac{\rho(\mathbf{k},x)}{-z+x},
| |
| </math>
| |
| which is related to <math>\mathcal{G}</math> and <math>G^{\mathrm{R}}</math> by
| |
| :<math>
| |
| \mathcal{G}(\mathbf{k},\omega_n) = G(\mathbf{k},\mathrm{i}\omega_n)
| |
| </math>
| |
| and
| |
| :<math>
| |
| G^{\mathrm{R}}(\mathbf{k},\omega) = G(\mathbf{k},\omega + \mathrm{i}\eta).
| |
| </math>
| |
| A similar expression obviously holds for <math>G^{\mathrm{A}}</math>.
| |
| | |
| The relation between <math>G(\mathbf{k},z)</math> and <math>\rho(\mathbf{k},x)</math> is referred to as a [[Hilbert transform]].
| |
| | |
| ====Proof of spectral representation====
| |
| | |
| We demonstrate the proof of the spectral representation of the propagator in the case of the thermal Green function, defined as
| |
| :<math>
| |
| \mathcal{G}(\mathbf{x} ,\tau|\mathbf{x} ',\tau') = \langle T\psi(\mathbf{x} ,\tau)\bar\psi(\mathbf{x} ',\tau')\rangle.
| |
| </math>
| |
| | |
| Due to translational symmetry, it is only necessary to consider <math>\mathcal{G}(\mathbf{x} ,\tau|\mathbf{0},0)</math> for <math>\tau > 0</math>, given by
| |
| :<math>
| |
| \mathcal{G}(\mathbf{x},\tau|\mathbf{0},0) = \frac{1}{\mathcal{Z}}\sum_{\alpha'} \mathrm{e}^{-\beta E_{\alpha'}}
| |
| \langle\alpha' | \psi(\mathbf{x},\tau)\bar\psi(\mathbf{0},0) |\alpha' \rangle.
| |
| </math>
| |
| Inserting a complete set of eigenstates gives
| |
| :<math>
| |
| \mathcal{G}(\mathbf{x} ,\tau|\mathbf{0},0) = \frac{1}{\mathcal{Z}}\sum_{\alpha,\alpha'} \mathrm{e}^{-\beta E_{\alpha'}}
| |
| \langle\alpha' | \psi(\mathbf{x} ,\tau)|\alpha \rangle\langle\alpha |\bar\psi(\mathbf{0},0) |\alpha' \rangle.
| |
| </math>
| |
| | |
| Since <math>|\alpha \rangle</math> and <math>|\alpha' \rangle</math> are eigenstates of <math>H-\mu N</math>, the Heisenberg operators can be rewritten in terms of Schrödinger operators, giving
| |
| :<math>
| |
| \mathcal{G}(\mathbf{x} ,\tau|\mathbf{0},0) = \frac{1}{\mathcal{Z}}\sum_{\alpha,\alpha'} \mathrm{e}^{-\beta E_{\alpha'}}
| |
| \mathrm{e}^{\tau(E_{\alpha'} - E_\alpha)}\langle\alpha' | \psi(\mathbf{x} )|\alpha \rangle\langle\alpha |\psi^\dagger(\mathbf{0}) |\alpha' \rangle.
| |
| </math>
| |
| Performing the Fourier transform then gives
| |
| :<math>
| |
| \mathcal{G}(\mathbf{k},\omega_n) = \frac{1}{\mathcal{Z}} \sum_{\alpha,\alpha'} \mathrm{e}^{-\beta E_{\alpha'}}
| |
| \frac{1-\zeta \mathrm{e}^{\beta(E_{\alpha'} - E_\alpha)}}{-\mathrm{i}\omega_n + E_\alpha - E_{\alpha'}} \int_{\mathbf{k}'} d\mathbf{k}' \langle\alpha |\psi(\mathbf{k})|\alpha' \rangle\langle\alpha' |\psi^\dagger(\mathbf{k}')|\alpha \rangle.
| |
| </math>
| |
| | |
| Momentum conservation allows the final term to be written as (up to possible factors of the volume)
| |
| :<math>
| |
| |\langle\alpha' |\psi^\dagger(\mathbf{k})|\alpha \rangle|^2,
| |
| </math>
| |
| which confirms the expressions for the Green functions in the spectral representation.
| |
| | |
| The sum rule can be proved by considering the expectation value of the commutator,
| |
| :<math>
| |
| 1 = \frac{1}{\mathcal{Z}} \sum_\alpha \langle\alpha |\mathrm{e}^{-\beta(H-\mu N)}[\psi_\mathbf{k},\psi_\mathbf{k}^\dagger]_{-\zeta}|\alpha \rangle,
| |
| </math>
| |
| and then inserting a complete set of eigenstates into both terms of the commutator:
| |
| :<math>
| |
| 1 = \frac{1}{\mathcal{Z}} \sum_{\alpha,\alpha'} \mathrm{e}^{-\beta E_\alpha} \left(
| |
| \langle\alpha |\psi_\mathbf{k} |\alpha' \rangle\langle\alpha' | \psi_\mathbf{k}^\dagger|\alpha \rangle - \zeta \langle\alpha |\psi_\mathbf{k}^\dagger |\alpha' \rangle\langle\alpha' | \psi_\mathbf{k}|\alpha \rangle
| |
| \right).
| |
| </math>
| |
| | |
| Swapping the labels in the first term then gives
| |
| :<math>
| |
| 1 = \frac{1}{\mathcal{Z}} \sum_{\alpha,\alpha'}
| |
| \left(\mathrm{e}^{-\beta E_{\alpha'}} - \zeta \mathrm{e}^{-\beta E_\alpha} \right)
| |
| |\langle\alpha | \psi_\mathbf{k}^\dagger|\alpha' \rangle|^2
| |
| ,
| |
| </math>
| |
| which is exactly the result of the integration of <math>\rho</math>.
| |
| | |
| ====Non-interacting case====
| |
| | |
| In the non-interacting case, <math>\psi_\mathbf{k}^\dagger|\alpha' \rangle</math> is an eigenstate with (grand-canonical) energy <math>E_{\alpha'} + \xi_\mathbf{k}</math>, where <math>\xi_\mathbf{k} = \epsilon_\mathbf{k} - \mu</math> is the single-particle dispersion relation measured with respect to the chemical potential. The spectral density therefore becomes
| |
| :<math>
| |
| \rho_0(\mathbf{k},\omega) = \frac{1}{\mathcal{Z}}\,2\pi\delta(\xi_\mathbf{k} - \omega) \sum_{\alpha'}\langle\alpha' |\psi_\mathbf{k}\psi_\mathbf{k}^\dagger|\alpha' \rangle(1-\zeta \mathrm{e}^{-\beta\xi_\mathbf{k}})\mathrm{e}^{-\beta E_{\alpha'}}.
| |
| </math>
| |
| | |
| From the commutation relations,
| |
| :<math>
| |
| \langle\alpha' |\psi_\mathbf{k}\psi_\mathbf{k}^\dagger|\alpha' \rangle =
| |
| \langle\alpha' |(1+\zeta\psi_\mathbf{k}^\dagger\psi_\mathbf{k})|\alpha' \rangle,
| |
| </math>
| |
| with possible factors of the volume again. The sum, which involves the thermal average of the number operator, then gives simply <math>[1 + \zeta n(\xi_\mathbf{k})]\mathcal{Z}</math>, leaving
| |
| :<math>
| |
| \rho_0(\mathbf{k},\omega) = 2\pi\delta(\xi_\mathbf{k} - \omega).
| |
| </math>
| |
| | |
| The imaginary-time propagator is thus
| |
| :<math>
| |
| \mathcal{G}_0(\mathbf{k},\omega) = \frac{1}{-\mathrm{i}\omega_n + \xi_\mathbf{k}}
| |
| </math>
| |
| and the retarded propagator is
| |
| :<math>
| |
| G_0^{\mathrm{R}}(\mathbf{k},\omega) = \frac{1}{-(\omega+\mathrm{i} \eta) + \xi_\mathbf{k}}.
| |
| </math>
| |
| | |
| ====Zero-temperature limit====
| |
| | |
| As <math>\beta\rightarrow\infty</math>, the spectral density becomes
| |
| :<math>
| |
| \rho(\mathbf{k},\omega) = 2\pi\sum_{\alpha} \left[ \delta(E_\alpha - E_0 - \omega)
| |
| |\langle\alpha |\psi_\mathbf{k}^\dagger|0 \rangle|^2
| |
| - \zeta \delta(E_0 - E_{\alpha} - \omega)
| |
| |\langle0 |\psi_\mathbf{k}^\dagger|\alpha \rangle|^2\right]
| |
| </math>
| |
| where <math>\alpha = 0</math> corresponds to the ground state. Note that only the first (second) term contributes when <math>\omega</math> is positive (negative).
| |
| | |
| ==General case==
| |
| | |
| ===Basic definitions===
| |
| | |
| We can use `field operators' as above, or creation and annihilation operators associated with other single-particle states, perhaps eigenstates of the (noninteracting) kinetic energy. We then use
| |
| :<math>
| |
| \psi(\mathbf{x} ,\tau) = \varphi_\alpha(\mathbf{x} ) \psi_\alpha(\tau),
| |
| </math>
| |
| where <math>\psi_\alpha</math> is the annihilation operator for the single-particle state <math>\alpha</math> and <math>\varphi_\alpha(\mathbf{x} )</math> is that state's wavefunction in the position basis. This gives
| |
| :<math>
| |
| \mathcal{G}^{(n)}_{\alpha_1\ldots\alpha_n|\beta_1\ldots\beta_n}(\tau_1 \ldots \tau_n | \tau_1' \ldots \tau_n')
| |
| = \langle T\psi_{\alpha_1}(\tau_1)\ldots\psi_{\alpha_n}(\tau_n)\bar\psi_{\beta_n}(\tau_n')\ldots\bar\psi_{\beta_1}(\tau_1')\rangle
| |
| </math>
| |
| with a similar expression for <math>G^{(n)}</math>.
| |
| | |
| ===Two-point functions===
| |
| | |
| These depend only on the difference of their time arguments, so that
| |
| :<math>
| |
| \mathcal{G}_{\alpha\beta}(\tau|\tau') = \frac{1}{\beta}\sum_{\omega_n}
| |
| \mathcal{G}_{\alpha\beta}(\omega_n)\,\mathrm{e}^{-\mathrm{i}\omega_n (\tau-\tau')}
| |
| </math>
| |
| and
| |
| :<math>
| |
| G_{\alpha\beta}(t|t') = \int_{-\infty}^{\infty}\frac{\mathrm{d}\omega}{2\pi}\,
| |
| G_{\alpha\beta}(\omega)\,\mathrm{e}^{-\mathrm{i}\omega(t-t')}.
| |
| </math>
| |
| | |
| We can again define retarded and advanced functions in the obvious way; these are related to the time-ordered function in the same way as above.
| |
| | |
| The same periodicity properties as described in above apply to <math>\mathcal{G}_{\alpha\beta}</math>. Specifically,
| |
| :<math>
| |
| \mathcal{G}_{\alpha\beta}(\tau|\tau') = \mathcal{G}_{\alpha\beta}(\tau-\tau')
| |
| </math>
| |
| and
| |
| :<math>
| |
| \mathcal{G}_{\alpha\beta}(\tau) = \mathcal{G}_{\alpha\beta}(\tau + \beta),
| |
| </math>
| |
| for <math>\tau < 0</math>.
| |
| | |
| ===Spectral representation===
| |
| | |
| In this case,
| |
| :<math>
| |
| \rho_{\alpha\beta}(\omega) = \frac{1}{\mathcal{Z}}\sum_{m,n} 2\pi \delta(E_n-E_m-\omega)\;
| |
| \langle m |\psi_\alpha|n \rangle\langle n |\psi_\beta^\dagger|m \rangle
| |
| \left(\mathrm{e}^{-\beta E_m} - \zeta \mathrm{e}^{-\beta E_n}\right) ,
| |
| </math>
| |
| where <math>m</math> and <math>n</math> are many-body states.
| |
| | |
| The expressions for the Green functions are modified in the obvious ways:
| |
| :<math>
| |
| \mathcal{G}_{\alpha\beta}(\omega_n) = \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega'}{2\pi}
| |
| \frac{\rho_{\alpha\beta}(\omega')}{-\mathrm{i}\omega_n+\omega'}
| |
| </math>
| |
| and
| |
| :<math>
| |
| G^{\mathrm{R}}_{\alpha\beta}(\omega) = \int_{-\infty}^{\infty} \frac{\mathrm{d}\omega'}{2\pi}
| |
| \frac{\rho_{\alpha\beta}(\omega')}{-(\omega+\mathrm{i}\eta)+\omega'}.
| |
| </math>
| |
| | |
| Their analyticity properties are identical. The proof follows exactly the same steps, except that the two matrix elements are no longer complex conjugates.
| |
| | |
| ====Noninteracting case====
| |
| | |
| If the particular single-particle states that are chosen are `single-particle energy eigenstates', i.e.
| |
| :<math>
| |
| [H-\mu N,\psi_\alpha^\dagger] = \xi_\alpha\psi_\alpha^\dagger,
| |
| </math>
| |
| then for <math>|n \rangle</math> an eigenstate:
| |
| :<math>
| |
| (H-\mu N)|n \rangle = E_n |n \rangle,
| |
| </math>
| |
| so is <math>\psi_\alpha |n \rangle</math>:
| |
| :<math>
| |
| (H-\mu N)\psi_\alpha|n \rangle = (E_n - \xi_\alpha) \psi_\alpha |n \rangle,
| |
| </math>
| |
| and so is <math>\psi_\alpha^\dagger|n \rangle</math>:
| |
| :<math>
| |
| (H-\mu N)\psi_\alpha^\dagger|n \rangle = (E_n + \xi_\alpha) \psi_\alpha^\dagger |n \rangle.
| |
| </math>
| |
| | |
| We therefore have
| |
| :<math>
| |
| \langle m |\psi_\alpha|n \rangle\langle n |\psi_\beta^\dagger|m \rangle =
| |
| \delta_{\xi_\alpha,\xi_\beta}\delta_{E_n,E_m+\xi_\alpha}\langle m |\psi_\alpha|n \rangle\langle n |\psi_\beta^\dagger|m \rangle.
| |
| </math>
| |
| | |
| We then rewrite
| |
| :<math>
| |
| \rho_{\alpha\beta}(\omega) = \frac{1}{\mathcal{Z}}\sum_{m,n} 2\pi \delta(\xi_\alpha-\omega)
| |
| \delta_{\xi_\alpha,\xi_\beta}\langle m |\psi_\alpha|n \rangle\langle n |\psi_\beta^\dagger|m \rangle
| |
| \mathrm{e}^{-\beta E_m}\left(1 - \zeta \mathrm{e}^{-\beta \xi_\alpha}\right),
| |
| </math>
| |
| therefore
| |
| :<math>
| |
| \rho_{\alpha\beta}(\omega) = \frac{1}{\mathcal{Z}}\sum_m 2\pi \delta(\xi_\alpha-\omega)
| |
| \delta_{\xi_\alpha,\xi_\beta}\langle m |\psi_\alpha\psi_\beta^\dagger\mathrm{e}^{-\beta (H-\mu N)}|m \rangle
| |
| \left(1 - \zeta \mathrm{e}^{-\beta \xi_\alpha}\right),
| |
| </math>
| |
| use
| |
| :<math>
| |
| \langle m |\psi_\alpha \psi_\beta^\dagger|m \rangle = \delta_{\alpha,\beta}\langle m |\zeta \psi_\alpha^\dagger \psi_\alpha + 1|m \rangle
| |
| </math>
| |
| and the fact that the thermal average of the number operator gives the Bose-Einstein or Fermi-Dirac distribution function.
| |
| | |
| Finally, the spectral density simplifies to give
| |
| :<math>
| |
| \rho_{\alpha\beta} = 2\pi \delta(\xi_\alpha - \omega)\delta_{\alpha\beta},
| |
| </math>
| |
| so that the thermal Green function is
| |
| :<math>
| |
| \mathcal{G}_{\alpha\beta}(\omega_n) = \frac{\delta_{\alpha\beta}}{-\mathrm{i}\omega_n + \xi_\beta}
| |
| </math>
| |
| and the retarded Green function is
| |
| :<math>
| |
| G_{\alpha\beta}(\omega) = \frac{\delta_{\alpha\beta}}{-(\omega+\mathrm{i}\eta) + \xi_\beta}.
| |
| </math>
| |
| Note that the noninteracting Green function is diagonal, but this will not be true in the interacting case.
| |
| | |
| ==See also==
| |
| *[[Fluctuation theorem]]
| |
| *[[Green-Kubo relations]]
| |
| *[[Linear response function]]
| |
| *[[Lindblad equation]]
| |
| | |
| == References ==
| |
| | |
| ===Books===
| |
| *Bonch-Bruevich V. L., [[Sergei Tyablikov|Tyablikov S. V.]] (1962): ''The Green Function Method in Statistical Mechanics.'' North Holland Publishing Co.
| |
| *Abrikosov, A. A., Gorkov, L. P. and Dzyaloshinski, I. E. (1963): ''Methods of Quantum Field Theory in Statistical Physics'' Englewood Cliffs: Prentice-Hall.
| |
| *Negele, J. W. and Orland, H. (1988): ''Quantum Many-Particle Systems'' AddisonWesley.
| |
| *[[Dmitry Zubarev|Zubarev D. N.]], Morozov V., Ropke G. (1996): ''Statistical Mechanics of Nonequilibrium Processes: Basic Concepts, Kinetic Theory'' (Vol. 1). John Wiley & Sons. ISBN 3-05-501708-0.
| |
| *Mattuck Richard D. (1992), ''A Guide to Feynman Diagrams in the Many-Body Problem'', Dover Publications, ISBN 0-486-67047-3.
| |
| | |
| ===Papers===
| |
| *[[Nikolay Bogoliubov|Bogolyubov N. N.]], [[Sergei Tyablikov|Tyablikov S. V.]] Retarded and advanced Green functions in statistical physics, Soviet Physics Doklady, Vol. '''4''', p.589 (1959).
| |
| *[[Dmitry Zubarev|Zubarev D. N.]], [http://dx.doi.org/10.1070/PU1960v003n03ABEH003275 Double-time Green functions in statistical physics], Soviet Physics Uspekhi '''3'''(3), 320—345 (1960).
| |
| | |
| [[Category:Quantum field theory]]
| |
| [[Category:Statistical mechanics]]
| |