Skip to content

Commit 366237c

Browse files
authored
Merge pull request #22 from numericalEFT/interaction
Interaction
2 parents bb09510 + 83632b1 commit 366237c

15 files changed

+580
-268
lines changed

Project.toml

+1-1
Original file line numberDiff line numberDiff line change
@@ -12,7 +12,7 @@ Lehmann = "95bf888a-8996-4655-9f35-1c0506bdfefe"
1212
Parameters = "d96e819e-fc66-5662-9728-84c9c7592b0a"
1313

1414
[compat]
15-
CompositeGrids = "0.0.4"
15+
CompositeGrids = ">= 0.0.4"
1616
Cuba = "2"
1717
GreenFunc = "0.1"
1818
LegendrePolynomials = "0.3"

docs/make.jl

+15-13
Original file line numberDiff line numberDiff line change
@@ -1,28 +1,30 @@
11
using ElectronGas
22
using Documenter
33

4-
DocMeta.setdocmeta!(ElectronGas, :DocTestSetup, :(using ElectronGas); recursive=true)
4+
DocMeta.setdocmeta!(ElectronGas, :DocTestSetup, :(using ElectronGas); recursive = true)
55

66
makedocs(;
7-
modules=[ElectronGas],
8-
authors="Kun Chen",
9-
repo="https://github.com/numericalEFT/ElectronGas.jl/blob/{commit}{path}#{line}",
10-
sitename="ElectronGas.jl",
11-
format=Documenter.HTML(;
12-
prettyurls=get(ENV, "CI", "false") == "true",
13-
canonical="https://numericalEFT.github.io/ElectronGas.jl",
14-
assets=String[],
7+
modules = [ElectronGas],
8+
authors = "Kun Chen",
9+
repo = "https://github.com/numericalEFT/ElectronGas.jl/blob/{commit}{path}#{line}",
10+
sitename = "ElectronGas.jl",
11+
format = Documenter.HTML(;
12+
prettyurls = get(ENV, "CI", "false") == "true",
13+
canonical = "https://numericalEFT.github.io/ElectronGas.jl",
14+
assets = String[]
1515
),
16-
pages=[
16+
pages = [
1717
"Home" => "index.md",
1818
"Manual" => [
1919
"manual/polarization.md",
20+
"manual/polarization_2D.md",
2021
"manual/legendreinteraction.md",
22+
"manual/legendreinteraction_2D.md",
2123
],
22-
],
24+
]
2325
)
2426

2527
deploydocs(;
26-
repo="github.com/numericalEFT/ElectronGas.jl",
27-
devbranch="master",
28+
repo = "github.com/numericalEFT/ElectronGas.jl",
29+
devbranch = "master"
2830
)

docs/src/manual/legendreinteraction.md

+2-2
Original file line numberDiff line numberDiff line change
@@ -53,13 +53,13 @@ could be decomposed with ``l``.
5353
For gap-function equation we have
5454
```math
5555
\begin{aligned}
56-
\Delta_l(k) = \int \frac{p^2dp}{{2\pi}^2} w_l(k, p) f_l(p),
56+
\Delta_l(k) = \int \frac{p^2dp}{{4\pi}^2} w_l(k, p) f_l(p),
5757
\end{aligned}
5858
```
5959
and for self-energy we have
6060
```math
6161
\begin{aligned}
62-
\Sigma(k) = \int \frac{p^2dp}{{2\pi}^2} w_0(k, p) G(p).
62+
\Sigma(k) = \int \frac{p^2dp}{{4\pi}^2} w_0(k, p) G(p).
6363
\end{aligned}
6464
```
6565

+79
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,79 @@
1+
# Decomposition of Interaction in two dimensions
2+
3+
In GW-approximation, we calculation self-energy as
4+
```math
5+
\Sigma(\mathbf{k},\omega_n)=-T\int \frac{{\rm d}^d \mathbf{q}}{(2\pi)^d} \sum_m G(\mathbf{p},\omega_m)W(\mathbf{k-p},\omega_n-\omega_m) \,,
6+
\tag{1}
7+
```
8+
where ``G`` is the Green's function and W is the effective interaction. Here, we suppress spin index.
9+
10+
## Spherical harmonic representation
11+
We first express the ``W(q,\tau)`` function as an expansion in Legendre polynomials ``P_\ell(\chi)``
12+
```math
13+
\begin{gathered}
14+
W(|\mathbf{k}-\mathbf{p}|, \tau)=\sum_{\ell=0}^{\infty} \bar{w}_{\ell}(k, p, \tau) P_{\ell}(\hat{k p}) \,, \\
15+
\bar{w}_{\ell}(k, p, \tau)=\frac{N(d,\ell)}{2} \int_{-1}^{1} d \chi P_{\ell}(\chi) W\left(\sqrt{k^{2}+p^{2}-2 k p \chi} ,\tau\right)\,.
16+
\end{gathered}
17+
```
18+
Since the Legendre polynomials of a scalar product of unit vectors can be expanded with spherical harmonics using
19+
```math
20+
P_{\ell}(\hat{k p})=\frac{\Omega_{d}}{N(d,\ell)} \sum_{m=1}^{N(d,\ell)} Y_{\ell m}(\hat{k}) Y_{\ell m}^{*}(\hat{p})\,,
21+
```
22+
where ``\Omega_{d}`` is the solid angle in ``d`` dimensions, and
23+
```math
24+
N(d, \ell)=\frac{2 \ell+d-2}{\ell}\left(
25+
\begin{array}{c}
26+
\ell+d-3 \\
27+
\ell-1
28+
\end{array}\right)
29+
```
30+
denotes the number of linearly independent homogeneous harmonic polynomials of degree ``\ell`` in ``d`` dimensions. The spherical harmonics are orthonormal as
31+
```math
32+
\int {\rm d}\Omega_{\hat k} Y_{\ell m}(\hat k) Y_{\ell^\prime m^\prime}(\hat k) = \delta_{\ell \ell^\prime} \delta_{mm^\prime} \,.
33+
```
34+
Hence, the ``W(q,\tau)`` function can expressed further on as
35+
```math
36+
W(|\mathbf{k}-\mathbf{p}|, \tau)=\sum_{\ell} \frac{\Omega_{d}}{N(d,\ell)} \bar{w}_{\ell}(k, p, \tau) \sum_{m} Y_{\ell m}(\hat{k}) Y_{\ell m}^{*}(\hat{p})
37+
```
38+
39+
or
40+
```math
41+
W(|\mathbf{k}-\mathbf{p}|, \tau)= \frac{\Omega_{d}}{2}\sum_{\ell} w_{\ell}(k, p, \tau) \sum_{m} Y_{\ell m}(\hat{k}) Y_{\ell m}^{*}(\hat{p} )\,.
42+
```
43+
with
44+
```math
45+
w_{\ell}(k, p, \tau)=\int_{-1}^{1} d \chi P_{\ell}(\chi) W\left(\sqrt{k^{2}+p^{2}-2 k p \chi} ,\tau\right)
46+
```
47+
In addition, the Green's function ``G(\mathbf p, \tau)`` has
48+
```math
49+
\begin{aligned}
50+
G(\mathbf p, \tau)= \sum_{\ell=0}^{\infty}\sum_{m=1}^{N(d,\ell)} G_{\ell m}(p,\tau) Y_{\ell m}(\hat p)\,, \\
51+
G_{\ell m}(p,\tau) = \int {\rm d}\Omega_{\hat k} G(\mathbf p,\tau) Y^*_{\ell m}(\hat p)
52+
\end{aligned}
53+
```
54+
### Decouple with channels
55+
By the sphereical harmonic expansion of Eq.(1), the self-energy
56+
57+
```math
58+
\begin{aligned}
59+
\sum_{\ell m} \Sigma_{\ell m }(k,\tau)Y_{\ell m}(\hat k) &= \frac{\Omega_d}{2} \int \frac{{\rm d}\mathbf p}{(2\pi)^d} \sum_{\ell m}G_{\ell m}(p, \tau) Y_{\ell m}(\hat p) \sum_{\ell^\prime m^\prime} w_{\ell^\prime}(k,p,\tau) Y_{\ell^\prime m^\prime}(\hat{k}) Y_{\ell^\prime m^\prime}^{*}(\hat{p} ) \\
60+
&=\frac{\Omega_d}{2} \sum_{\ell m} \int \frac{p^{d-1}dp}{(2\pi)^d} G_{\ell m}(p, \tau) w_{\ell}(k,p,\tau) Y_{\ell m}(\hat k)
61+
\end{aligned}
62+
```
63+
Since self-energy is symmertric with ``\hat k``, we just need to project on the s-wave channel, namely
64+
```math
65+
\Sigma(k,\tau) =\frac{\Omega_d}{2} \int \frac{p^{d-1}dp}{(2\pi)^d} G(p,\tau) w_0(k,p,\tau)
66+
```
67+
68+
## Two dimensions
69+
```math
70+
\begin{aligned}
71+
N(2,\ell) &= 2 \\
72+
Y_{\ell 1}(\hat k) &= \cos(\ell \theta),\; Y_{\ell 2}(\hat k) = \sin(\ell \theta) \\
73+
P_{\ell}(\hat{kp}) &= \pi \cos[\ell(\theta_{\hat k} - \theta_{\hat p})]
74+
\end{aligned}
75+
```
76+
Hence, the self-energy is
77+
```math
78+
\Sigma(k,\tau) = \int \frac{pdp}{4\pi} G(p,\tau)w_0(k,p,\tau)
79+
```

docs/src/manual/polarization.md

+31-34
Original file line numberDiff line numberDiff line change
@@ -1,67 +1,64 @@
11

22
# Polarization of free electron
33

4-
## Definition
4+
## Generic formalism
55

66
The bare polarization is defined as
77

88
```math
99
\begin{aligned}
10-
\Pi_0(\omega_n, \vec{q})=
11-
N_ST\sum_{m}\int\frac{d^Dk}{{(2\pi)}^D}G_0(\omega_m, \vec{k})G_0(\omega_{m+n}, \vec{k}+\vec{q})
10+
\Pi_0(\Omega_n, \vec{q})=
11+
-S T\sum_{m}\int\frac{d^d k}{{(2\pi)}^d}G_0(\omega_m, \vec{k})G_0(\omega_m+\Omega_n, \vec{k}+\vec{q})
1212
\end{aligned}
1313
```
14-
15-
where ``N_S`` is spin number, ``G_0`` is bare Green's function.
14+
in the matsubara representation, where ``S`` is spin number, ``G_0`` is bare Green's function.
1615
We have
1716

1817
```math
1918
\begin{aligned}
20-
G_0(\omega_m, \vec{k})=\frac{1}{i\omega_m-\varepsilon_{\vec{k}}+\mu}
19+
G_0(\omega_m, \vec{k})=\frac{1}{i\omega_m-\epsilon_{\vec{k}}}
2120
\end{aligned}
2221
```
2322

24-
with dispersion given by ``\varepsilon_{\vec{k}}`` and chemical potential ``\mu``.
25-
26-
27-
## Free electron in 3D
28-
29-
From now on we consider free electron in 3D, where ``D=3``, ``N_S=2`` and ``\varepsilon_{\vec{k}}=\frac{k^2}{2m}``.
30-
Then by summing over frequency we have
23+
with bare electron dispersion given by ``\epsilon_{\vec{k}}``. By summing over frequency we have
3124

3225
```math
3326
\begin{aligned}
34-
\Pi_0(\omega_n, \vec{q})&=2T\sum_{m}\int\frac{d^3k}{{(2\pi)}^3}
35-
\frac{1}{(i\omega_m-\varepsilon_{\vec{k}}+\mu)(i\omega_{m+n}-\varepsilon_{\vec{k}+\vec{q}}+\mu)}\\
36-
&=-2\int\frac{d^3k}{{(2\pi)}^3}
37-
\frac{n(\vec{k}+\vec{q})-n(\vec{k})}{i\omega_n-\varepsilon_{\vec{k}+\vec{q}}+\varepsilon_{\vec{k}}}\\
27+
\Pi_0(\Omega_n, \vec{q}) &= -S \int\frac{d^d \vec k}{{(2\pi)}^d}
28+
\frac{n(\epsilon_{\vec{k}+\vec{q}})-n(\epsilon_{\vec{k}})}{i\omega_n-\epsilon_{\vec{k}+\vec{q}}+\epsilon_{\vec{k}}} \\
29+
&=-S\int \frac{d^d \vec k}{{(2\pi)}^d} n(\epsilon_{\vec k}) \left[ \frac{1}{i\Omega+\epsilon_{\vec k}-\epsilon_{\vec k+\vec q}}-\frac{1}{i\Omega+\epsilon_{\vec k-\vec q}-\epsilon_{\vec k}}\right]
3830
\end{aligned}
3931
```
40-
41-
with
32+
with the fermi distribution function
4233

4334
```math
44-
\begin{aligned}
45-
n(\vec{k})=\frac{1}{e^{\beta(\varepsilon_{\vec{k}}-\mu)}+1}.
46-
\end{aligned}
35+
n(\epsilon_{\vec k}) =\frac{1}{e^{\beta\varepsilon_{\vec{k}}}+1}
4736
```
4837

49-
Now plug in the dispersion and we have
38+
## Free electron in 3D
39+
40+
From now on we consider free electron in 3D, where ``d=3`` and dispersion ``\epsilon_{\vec{k}}=k^2/2m-\mu``. We have
5041

5142
```math
5243
\begin{aligned}
53-
\Pi_0(\omega_n, \vec{q})
54-
&=-2\int\frac{d^3k}{{(2\pi)}^3}
55-
[\frac{n(\vec{k})}{i\omega_n-\varepsilon_{\vec{k}}+\varepsilon_{\vec{k}+\vec{q}}}
56-
-\frac{n(\vec{k})}{i\omega_n-\varepsilon_{\vec{k}+\vec{q}}+\varepsilon_{\vec{k}}}
57-
]\\
58-
&=-2\int\frac{d^3k}{{(2\pi)}^3}
59-
\frac{4mn(\vec{k})(q^2+2kq\cos(\theta))}
60-
{4m^2\omega_n^2+{(q^2+2kq\cos(\theta))}^2}\\
61-
&=-\int\frac{kmdk}{2\pi^2q}
62-
n(\vec{k})\ln(\frac{4m^2\omega_n^2+{(q^2+2kq)}^2}{4m^2\omega_n^2+{(q^2-2kq)}^2})\\
44+
\Pi_0(\Omega, \vec{q})&=-S\int_0^{\infty} \frac{k^2dk}{4\pi^2} n(\epsilon_k) \int_{-1}^{1} d(\cos \theta) \left[ \frac{1}{i\Omega+\epsilon_{\vec k}-\epsilon_{\vec k+\vec q}}-\frac{1}{i\Omega+\epsilon_{\vec k-\vec q}-\epsilon_{\vec k}}\right]\\
45+
&=-S\int_0^{\infty} \frac{k^2dk}{4\pi^2} n(\epsilon_k) \frac{m}{kq}\ln\frac{4m^2\Omega^2+(q^2-2kq)^2}{4m^2\Omega^2+(q^2+2kq)^2} \,,
6346
\end{aligned}
6447
```
65-
6648
which could be handled with one dimensional integral of ``k``.
6749

50+
- In the limit ``q^2+2k_F q \ll 2m\Omega_n ``, the intergrand of ``\Pi_0`` is expanded as
51+
```math
52+
\frac{m}{kq}\ln\frac{4m^2\Omega^2+(q^2-2kq)^2}{4m^2\Omega^2+(q^2+2kq)^2}=-\frac{2q^2}{m\Omega^2}+\frac{2k^2q^4}{m^3\Omega^4}+\frac{(-4k^2+m^2\Omega^2)q^6}{2m^5\Omega^6}+...
53+
```
54+
- Zero temperature polarization can be calculated explicitly
55+
```math
56+
\Pi_0(\Omega,q) = -\frac{N_F}{2}\left[1-\frac{1}{8 k_F q}\left\{ \left[\frac{(i2m\Omega-q^2)^2}{q^2}-4 k_F^2\right]\log\left(\frac{i2m\Omega-q^2-2 k_F q}{i2m\Omega-q^2+2 k_F q}\right)+\left[\frac{(i2m\Omega+q^2)^2}{q^2}-4 k_F^2\right]\log\left(\frac{i2m\Omega+q^2+2 k_F q}{i2m\Omega+q^2-2 k_F q}\right)\right\}\right]
57+
```
58+
- In the static limit ``\Omega=0``,
59+
```math
60+
\Pi_0(0, q) = -N_F F(q/2k_F) \,,
61+
```
62+
where ``N_F=Smk_F/(2\pi^2)`` is the density of states, and ``F(x)=\frac{1}{2}-\frac{x^2-1}{4x}\ln \left|\frac{1+x}{1-x}\right|`` is the Lindhard function.
63+
64+
The weak logarithmic singularity near ``2k_F`` is the cause of the Friedel oscillation and Kohn-Luttinger superconductivity.

docs/src/manual/polarization_2D.md

+35
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,35 @@
1+
# Polarization of free electron in two dimensions
2+
3+
We consider the polarziation of free electron gas in 2D,
4+
```math
5+
\Pi_0(q, \Omega_n)
6+
=-S\int\frac{d^2 \vec k}{{(2\pi)}^2}
7+
\frac{n(\epsilon_{\vec{k}})-n(\epsilon_{\vec{k}+\vec{q}})}{i\Omega_n+\epsilon_{\vec{k}}-\epsilon_{\vec{k}+\vec{q}}}\\
8+
```
9+
where ``n(\vec k) = 1/(e^{\beta(\epsilon_{\vec k}-\mu)}+1)``, ``\epsilon_{\vec k}=k^2/(2m)-\mu``, and ``S`` is the spin factor. It is expressed as
10+
```math
11+
\Pi_0(q, \Omega_n)=-S\int_0^{\infty} \frac{mkdk}{2\pi^2} n(\epsilon_k) \int_{0}^{2\pi} d\theta \left[ \frac{1}{i2m\Omega_n-2kq \cos\theta-q^2}-\frac{1}{i2m\Omega_n-2kq\cos \theta+q^2}\right]
12+
```
13+
14+
## Static limit ``\Omega_n=0``
15+
16+
For real ``a,b``, the integral ``\int_0^{2\pi}\frac{1}{a-b\cos \theta}=\frac{2\pi}{\sqrt{a^2-b^2}}`` and ``\int_0^{2\pi}\frac{1}{a+b\cos \theta}=C\frac{2\pi}{\sqrt{a^2-b^2}}`` where ``C=1`` for ``a>b`` and ``C=-1`` for ``a<b``. Hence, we have
17+
```math
18+
\Pi_0(q, 0)= -S\int_0^{q/2}\frac{mkdk}{2\pi^2} n(\epsilon_k) \frac{4\pi}{\sqrt{q^4-4k^2q^2}}
19+
```
20+
21+
At zero temperature,
22+
- for ``q<2k_F``,
23+
```math
24+
\Pi_0(q, 0)=-S\int_0^{q/2} \frac{m}{\pi q} \frac{dk^2}{\sqrt{q^2-4k^2}}=-\int_0^1\frac{mS}{4\pi}\frac{dx}{\sqrt{1-x}}=-\frac{mS}{2\pi} \,;
25+
```
26+
- for ``q>2k_F``,
27+
```math
28+
\Pi_0(q, 0)=-S\int_0^{k_F} \frac{m}{\pi q} \frac{dk^2}{\sqrt{q^2-4k^2}}=-\int_0^{4k_F^2/q^2}\frac{mS}{4\pi}\frac{dx}{\sqrt{1-x}}=-\frac{mS}{2\pi}\left( 1-\sqrt{1-\frac{4k_F^2}{q^2}}\right) \,.
29+
```
30+
31+
## Zero temperature
32+
```math
33+
\Pi_0(q, \Omega_n)=-\frac{mS}{2\pi} \left[1-\frac{2k_F}{q}{\rm Re}\sqrt{\left(\frac{q}{2k_F}+ i \frac{m\Omega_n}{qk_F}\right)^2-1} \right] \,.
34+
```
35+
For ``q\to 0`` and ``\Omega_n \neq 0``, ``\Pi_0(q, \Omega_n) \to 0``.

example/spinfermion.jl

+41
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,41 @@
1+
# Self-energy for spin-fermion model
2+
3+
using ElectronGas
4+
using Test, Printf, DelimitedFiles
5+
6+
dim = 2
7+
beta, rs = 1e2, 1.0
8+
espin = 1.0
9+
param = Interaction.Parameter.defaultUnit(1 / beta, rs, dim)
10+
11+
Λa = Polarization.Polarization0_ZeroTemp(0.0, 0, param) * (-espin^2) / param.ϵ0
12+
println(Λa)
13+
14+
param = Parameter.Para(param, e0 = 0.0, espin = espin, Λa = Λa)
15+
16+
Euv, rtol = 100 * param.EF, 1e-10
17+
# Nk, order, minK = 8, 4, 1e-7
18+
Nk, order, minK = 11, 4, 1e-8
19+
20+
Σ = SelfEnergy.G0W0(param, Euv, rtol, Nk, 10 * param.kF, minK * param.kF, order, :rpa)
21+
Σ = SelfEnergy.GreenFunc.toMatFreq(Σ)
22+
23+
kgrid = Σ.spaceGrid
24+
kF = kgrid.panel[3]
25+
kF_label = searchsortedfirst(kgrid.grid, kF)
26+
println(kF_label)
27+
ωgrid = Σ.dlrGrid
28+
29+
ΣR = real.dynamic)
30+
ΣI = imag.dynamic)
31+
println.instant[1, 1, :])
32+
println(ΣR[1, 1, kF_label, :])
33+
34+
# f = open("./data/Nk$Nk/SigmaIm_b1e2_m0.txt", "w")
35+
for (n, sigma) in enumerate(ΣI[1, 1, kF_label, :])
36+
println(n, ' ', ωgrid.n[n], ' ', ωgrid.ωn[n], ' ', sigma)
37+
writedlm(f, [ωgrid.n[n] ωgrid.ωn[n] sigma])
38+
end
39+
# close(f)
40+
41+
println(SelfEnergy.zfactor(Σ))

src/convention.jl

+4-4
Original file line numberDiff line numberDiff line change
@@ -13,10 +13,10 @@ const INL, OUTL, INR, OUTR = 1, 2, 3, 4
1313
const EPS = 1e-16
1414

1515
# export all conventions
16-
for n in names(@__MODULE__; all=true)
17-
if Base.isidentifier(n) && n (Symbol(@__MODULE__), :eval, :include)
18-
@eval export $n
19-
end
16+
for n in names(@__MODULE__; all = true)
17+
if Base.isidentifier(n) && n (Symbol(@__MODULE__), :eval, :include)
18+
@eval export $n
19+
end
2020
end
2121

2222
end

src/interaction.jl

+8-6
Original file line numberDiff line numberDiff line change
@@ -52,7 +52,7 @@ Bare interaction in momentum space. Coulomb interaction if Λs=0, Yukawa otherwi
5252
"""
5353
function coulomb(q, param)
5454
@unpack me, kF, rs, e0s, e0a, β, Λs, Λa, ϵ0 = param
55-
if (q^2+Λs)*(q^2+Λa) 0.0
55+
if (q^2 + Λs) * (q^2 + Λa) 0.0
5656
return 0.0, 0.0
5757
else
5858
return e0s^2 / ϵ0 / (q^2 + Λs), e0a^2 / ϵ0 / (q^2 + Λa)
@@ -71,12 +71,12 @@ function bubbledyson(V::Float64, F::Float64, Π::Float64, n::Int)
7171
end
7272
if n == 0
7373
if F == 0
74-
K = -V * Π * (1)^2 / (1.0 / V + Π * (1))
74+
K = V * Π * (1)^2 / (1.0 / V - Π * (1))
7575
else
76-
K = -V * Π * (1 - F / V)^2 / (1.0 / V + Π * (1 - F / V))
76+
K = V * Π * (1 - F / V)^2 / (1.0 / V - Π * (1 - F / V))
7777
end
7878
else
79-
K = -(Π) * (V - F)^2 / (1.0 + (Π) * (V - F))
79+
K = Π * (V - F)^2 / (1.0 - (Π) * (V - F))
8080
end
8181
@assert !isnan(K) "nan at V=$V, F=$F, Π=, n=$n"
8282
return K
@@ -116,7 +116,8 @@ end
116116

117117
function RPAwrapped(Euv, rtol, sgrid::SGT, param;
118118
pifunc = Polarization0_ZeroTemp, landaufunc = landauParameterTakada, V_Bare = coulomb) where {SGT}
119-
@unpack me, kF, rs, e0, β, Λs, ϵ0 = param
119+
120+
@unpack β = param
120121
gs = GreenFunc.Green2DLR{Float64}(:rpa, GreenFunc.IMFREQ, β, false, Euv, sgrid, 1; timeSymmetry = :ph, rtol = rtol)
121122
ga = GreenFunc.Green2DLR{Float64}(:rpa, GreenFunc.IMFREQ, β, false, Euv, sgrid, 1; timeSymmetry = :ph, rtol = rtol)
122123
green_dyn_s = zeros(Float64, (gs.color, gs.color, gs.spaceGrid.size, gs.timeGrid.size))
@@ -186,7 +187,8 @@ end
186187

187188
function KOwrapped(Euv, rtol, sgrid::SGT, param;
188189
pifunc = Polarization0_ZeroTemp, landaufunc = landauParameterTakada, V_Bare = coulomb) where {SGT}
189-
@unpack me, kF, rs, e0, β, Λs, ϵ0 = param
190+
191+
@unpack β = param
190192
gs = GreenFunc.Green2DLR{Float64}(:ko, GreenFunc.IMFREQ, β, false, Euv, sgrid, 1; timeSymmetry = :ph, rtol = rtol)
191193
ga = GreenFunc.Green2DLR{Float64}(:ko, GreenFunc.IMFREQ, β, false, Euv, sgrid, 1; timeSymmetry = :ph, rtol = rtol)
192194
green_dyn_s = zeros(Float64, (gs.color, gs.color, gs.spaceGrid.size, gs.timeGrid.size))

0 commit comments

Comments
 (0)