International Journal of Statistical Mechanics

Volume 2013 (2013), Article ID 317491, 5 pages

http://dx.doi.org/10.1155/2013/317491

Research Article

Spectral Functions and Properties of Nuclear Matter

V. A. Danilenko,¹ K. A. Gridnev,¹ and A. S. Kondratyev²

¹Department of Nuclear Physics, Saint Petersburg State University, Neva River Embankment 7/9, Saint Petersburg 199034, Russia

²Department of Physics, Herzen State Pedagogical University of Russia, Moika River Embankment 48, Saint Petersburg 191168, Russia

Received 14 April 2013; Accepted 19 June 2013

Academic Editor: Adam Lipowski

Copyright © 2013 V. A. Danilenko et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The Green’s function method in the Kadanoff-Baym version provides a basic theory for nuclear dynamics which is applicable also to nonzero temperature and to nonequilibrium systems. At the same time, it maintains the basic many-body techniques of the Brueckner theory that makes reasonable a comparison of the numerical results of the two methods for equilibrium systems. The correct approximation to the spectral function which takes into account the widths of energy levels is offered and discussed, and the comparison of the values of binding energy in the two methods is produced.

1. Introduction

The only microscopic theory which is capable of describing dynamical and statistical properties of quantum many-body systems in a comprehensive way is the Green’s function approach initiated by Martin and Schwinger and later developed by Kadanoff and Baym [1, 2] with specific application to nuclear matter in [3–6]. Traditional many-body theories of nuclei such as Brueckner’s are based on a quasiparticle picture ab initio neglecting the widths of the energy levels associated with the strong interactions when calculating the binding energy. In its original form, the Brueckner theory only considers particle propagation in intermediate states by defining a two-body reaction matrix of “effective” interaction or similarly defined objects [3]. Different approximations to the spectral functions in the Kadanoff-Baym method established a link between the Brueckner and the Green’s function theories [3–6]. Similar approximations to the spectral functions and kinetic equations were considered also in [7–9]. Unfortunately, all the approximations offered in [3–9] turned to be unsatisfactory as it was shown in [10].

The purpose of this paper is to generalize the results of [10] concerning the spectral function and to clarify some of the quantitative relations between the two methods in nuclear physics on the basis of these results. We show that the self-consistent approximations in the two theories lead to very close results, although numerical calculations of the binding energy are still limited by our limited knowledge of nuclear forces and many-body effects.

2. Spectral Functions of the Particle and Quasiparticle States in the Kadanoff-Baym Approach

The KB formalism results in the following general expression for a one-particle spectral function for a system in equilibrium state [1, 2, 10]:

where

and is a one-particle energy in the Hartree-Fock approximation. Real and imaginary parts of the correlation self-energy function are related through the Hilbert transform

Here, refers to a principal value integration.

The spectral function satisfies the exact sum rule:

for all values of momentum . This result follows directly from the commutator relations for field operators and can serve as a keystone for checking all approximations to the spectral function (1).

It was shown in [1, 2] that, in the case of slowly varying in space and time disturbances, after the transition to the Wigner coordinates and Fourier transform with respect to the difference of the initial time and space variables, the spectral function obeys (9–27) in [1, 2] ((10) in [10]) and can be represented by the same equation (1), and only all the quantities entering (1) become also the functions of space and time variables and . This property is one more keystone for checking the validity of different approximations to the spectral function.

We consider a system in equilibrium state and define the quasiparticle energy as a solution of [1, 2]

In the case of a small width of the energy level, when , we get

where is the Dirac delta function. It is shown in [1, 2] that the argument of the delta function is a monotonically increasing function of for all near the Fermi level where . It then follows that, in this case for every , there exists just one root of (5).

Now we expand as a function of in Taylor series near the value and save only linear terms:

Substituting (7) into (1), we get the formula for called the spectral function of the quasiparticle state:

where is a renormalizing factor defined by the expression

It is easy to prove that for all values of when correlation energy is taken into account [1, 2]. The quantity can be considered the width of the quasiparticle state, and (8) can be rewritten as

A severe drawback of is that it normalizes to rather than to 1 as in (4):

This deficiency can be remedied by better approximations to the spectral function (1).

If we neglect the width of the quasiparticle state , (10) reads

with the same normalization (11).

Another deficiency associated with the quasiparticle spectral function (10) is the change of in the width for , which contradicts the dispersion relation (3). As a result, the following approximations to the spectral function (1) were introduced in [3–9].

Quasiclassical approximation [3, 4]:

Quasiparticle approximation in the form (12) [3–5]:

Extended quasiparticle approximation [3–6]:

Improved quasiparticle approximation was offered in [4] in the form

but this approximation was never used in numerical calculations.

The approximations (13) and (14) correspond to the cases when the widths of the energy levels are neglected, and the approximation (13) is twice contradictory: the correlation part of interaction is taken into account, but the normalization caused by this interaction is neglected. The approximations (15) and (16) correspond to the cases when the widths of the levels are taken into account, but they were obtained by a nonadequate expansion of the expression (1) in powers of starting with the Dirac delta function. It was shown in [10] that these approximations do not obey (9–27) in [1, 2] in the case of slowly varying disturbances and consequently cannot be correct; as it was stated previously, equilibrium and nonequilibrium approximations are genetically linked with one another [1, 2].

The way of obtaining a correct approximation to the spectral function which takes into account the widths of the levels and obeys (10) in [10] can be associated with the integral Fourier transform. Such an approximation looks like

The proof of (17) starts with the formula

We expand the first exponent in the left side of (18) in powers of and integrate term by term using the relations [11–13]

where is the 2th derivative of the Dirac delta function. We should mention that the uniform convergence of a series is only a sufficient, but not a necessary, condition for the term-by-term integration validity [14, 15]. Although the uniform convergence of the Taylor series resulting from the expansion of the exponent in (18) is absent, the term-by-term integration, as we will see, leads to the reasonable result.

The terms of the expansion containing the odd powers of in correspondence with (19) generate a geometric sequence with the common ratio :

The sum of the sequence equals to

The expression (22) coincides with the general formula (1) for the spectral function.

The terms of the expansion containing the even powers of can be considered on the basis of the formula [15, page 92]

It is easy to prove by mathematical induction that

In correspondence with (20) and (24) the terms even in powers of generate a geometric sequence with the same common ratio (21). The sum of the sequence is equal to

The contribution of to the sum rule (4) is equal to zero:

The obtained results prove the validity of the produced expansion and term-by-term integration and testify to the approximation of the spectral function with the precision to the linear in terms in the form (17).

3. Binding Energy of Nuclear Matter

The Brueckner theory of nuclear matter [16] is based on the matrix of effective interaction with the diagonal elements defined by formula (1) in [4]. One-particle energies entering the expression for the diagonal elements are given by the equation

The final expression for the quasiparticle energy in the Brueckner theory looks like [4, 17]

The expressions for the potentials, , and are given correspondingly by the formulas (6), (7), and (10) in [4]. The results of numerical calculations presented in [3–6] were based on these expressions with different approximate transformations. The discussion of these transformations can be found in [4, 6].

Numerical calculations on the basis of spectral functions were presented with a comparison of the results obtained through the two approaches, based on the Brueckner theory and on the spectral functions [3–6]. The results obtained with the different methods are as follows [4].

The binding energy (per nucleon) in the quasiclassical approximation (13) is

The binding energy in the quasiparticle approximation (14) is

The binding energy in the extended quasiparticle approximation (15) is

In the Brueckner theory, when the formula (27) for the one-particle energy is used in (28) for the quasiparticle energy , the binding energy is

When in the formula (28) for the quasiparticle energy is changed for defined by (5) ( in the notation of [4] with exactly the same definition),

the binding energy turns to be equal

The idea of the change (33) was considered in [3, 4] by the comparison of the formulas in the Brueckner theory and in the Green’s function method with the diagram expansions for the interaction energy which differ in the two approaches. Equation (33) corresponds to several energy diagrams in the perturbation theory. But the solution of (5) for corresponds to the completely renormalized energy when formally all diagrams of infinite series are taken into account. Thus, here, one has to be very careful with overcounting, the fact also mentioned in [4]. As a result, it is reasonable to compare the value of the binding energy in the spectral function method with its value in the Brueckner theory  MeV/nucleon but not  MeV/nucleon. In this comparison, formula (17) for the approximation is to be used.

Compare the expressions (13) and (14) with the values of and to get the normalizing factor :

which consequently gives .

Compare the expressions (14) and (15) with the values of and to find the contribution of the second term in the right side of (15) to the value of the binding energy:

This corresponds to the contribution of the second term in the right side of (17) to the value of the binding energy:

Finally, the binding energy corresponding to the approximation (17) to the spectral function turns to be equal

which is very close to the value  MeV/nucleon in the Brueckner theory.

4. Conclusions

Self-consistent approximations used in the Brueckner theory and in the Green’s function method lead to almost coinciding values of the binding energy provided the correct formula (17) for the spectral function is used. The Brueckner formalism as well as its applications is found in many publications using various approaches of presentation [4]. This theory is strictly applicable to zero-temperature case. On the contrary, the Green’s function method in Kadanoff-Baym version has certain advantage going beyond these restrictions. As mentioned in [6], “it can be now conclusively stated that self-consistent spectral functions should (and perhaps equally important now can) be included in any serious consideration of nuclear properties being it in ground state, excited state, or in a state of nonequilibrium.” For the study of nonequilibrium phenomena, the method can be used with constraints to define an initial condition. It should be stated that, in the case of nonequilibrium phenomena, different forms of the Kadanoff-Baym equations can be preferable for different physical systems [18, 19].

References

L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics, Benjamin, New York, NY, USA, 1962.

L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics, Perseus Books, Cambridge, Mass, USA, 1989.

H. S. Köhler, “Spectral functions in nuclear matter,” Physical Review C, vol. 46, no. 5, pp. 1687–1696, 1992.

H. S. Köhler and R. Malfliet, “Extended quasiparticle approximation and Brueckner theory,” Physical Review C, vol. 48, no. 3, pp. 1034–1045, 1993.

H. S. Köhler and K. Morawetz, “Correlations in many-body systems with two-time Green’s functions,” Physical Review C, vol. 64, no. 2, Article ID 024613, 2001.

H. S. Köhler, “Beyond the quasi-particle picture in nuclear matter calculations using Green's function techniques,” Journal of Physics, vol. 35, article 384, 2006.

W. D. Kraeft, D. Kremp, W. Ebeling, and G. Ropke, Quantum Statistics of Charged Particle Systems, Springer, Berlin, Germany, 1986.

V. Spicka and P. Lipavsky, “Quasiparticle Boltzmann equation in semiconductors,” Physical Review Letters, vol. 73, no. 25, pp. 3439–3442, 1994.

V. Spicka and P. Lipavsky, “Quasiparticle Boltzmann equation in semiconductors,” Physical Review B, vol. 52, no. 20, pp. 14615–14635, 1995.

M. Arshad, A. S. Kondratyev, and I. Siddique, “Spectral function and kinetic equation for a normal Fermi liquid,” Physical Review B, vol. 76, no. 5, Article ID 54306, 2007.

E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals, Clarendon, Oxford, UK, 1975.

M. J. Lighthill, Introduction to Fourier Analysis and Generalized Functions, Cambridge University Press, Cambridge, UK, 1960.

Y. A. Brychkov and A. P. Prudnikov, Integral Transform of Generalized Functions, Gordon and Breach, New York, NY, USA, 1989.

W. Appel, Mathematics for Physics and Physisists, Princeton University Press, Princeton, NJ, USA; Oxford, UK, 2007.

G. B. Arfken and H. J. Weber, Mathematical Methods For Physisists, Academic Press; Elsevier, 2005.

K. A. Brueckner and D. T. Goldman, “Single particle energies in the theory of nuclear matter,” Physical Review, vol. 117, no. 1, pp. 207–213, 1960.

H. S. Köhler, “Theory of finite nuclei,” Physical Review, vol. 137, no. 5, pp. B1145–B1157, 1965.

A. S. Kondratyev and I. Siddique, “Kinetic equations in the theory of normal Fermi liquid,” Low Temperature Physics, vol. 34, no. 2, pp. 137–141, 2008.

A. S. Kondratyev and N. Shahid, “Different forms of the Kadanoff-Baym equations in quantum statistical mechanics,” Low Temperature Physics, vol. 37, article 777, no. 10, 2011.