Calculating asymptotic behavior of a correlation function

In summary: This gives\begin{align}\int_{BZ} d^2q \exp{[i \mathbf{q}\cdot\mathbf{r}]} \Big( \ldots \Big) &\rightarrow \frac{1}{N^2} \sum_{q,q'} \cos{[(\mathbf{q}-\mathbf{q}')\cdot(\mathbf{r}_1-\mathbf{r}_2)]} \frac{\Delta_q \Delta_{q'} - \epsilon_q \epsilon_{q'}}{E_q E
  • #1
conana
23
0
I was hoping to get some assistance in reproducing a calculation from https://arxiv.org/abs/0803.1292 (https://journals.aps.org/pra/abstract/10.1103/PhysRevA.78.012304 for the published version).

A certain four-spin correlation function for Kitaev's spin-1/2 model on the honeycomb lattice is defined in Eq. 28, and from Eq.'s 29 and 30 one arrives at the final expression for the correlation function as

[tex] C(\mathbf{r}_1, \mathbf{r}_2) = \frac{1}{N^2} \sum_{q,q'} \cos{[(\mathbf{q}-\mathbf{q}')\cdot(\mathbf{r}_1-\mathbf{r}_2)]} \frac{\Delta_q \Delta_{q'} - \epsilon_q \epsilon_{q'}}{E_q E_{q'}},[/tex]

where N is the number of sites in the lattice, the summations are over momenta in the first Brillouin zone, [itex]\mathbf{r}_1, \mathbf{r}_2[/itex] denote unit cell positions, [itex] \Delta_q = \sin{q_1} + \sin{q_2} [/itex], [itex] \epsilon_q = \cos{q_1} + \cos{q_2} + 1 [/itex], and [itex] E_{q} = \sqrt{\epsilon_q^2 + \Delta_q^2} [/itex] with [itex] q_1, q_2 [/itex] being the components of [itex] \mathbf{q} [/itex] in the basis of reciprocal lattice vectors.

Arriving at this expression is not a problem. However, the authors go on to argue that for large [itex] |\mathbf{r}_1 - \mathbf{r}_2| [/itex], the correlation function decays algebraically as
[tex] C(\mathbf{r}_1, \mathbf{r}_2) \sim \frac{1}{|\mathbf{r}_1 - \mathbf{r}_2|^4}. [/tex]
In arriving at this expression, the authors state "[...] the denominator in Eq. (30) has two zero points, which are of order 1/N in the large N limit. Their contribution causes the summation to be finite in the thermodynamic limit. Then using the stationary phase method, we can evaluate the exponents of the correlation function at long distance to be 4." Eq. 30 refers to (essentially) the above expression for the correlation function involving the summations over the Brillouin zone, and the two zero points of the denominator correspond to the two Dirac points in the fermionic quasiparticle spectrum.

It is the calculation of this asymptotic behavior which I am trying to reproduce, but I seem to be unable to do this with the information given. It seems there are no momenta at which the phase is stationary due to the phase being a linear function of momentum. So then I thought, perhaps the authors intend to imply that the phase is never stationary and, thus, the summation would vanish at all points in the Brillouin zone EXCEPT for at the Dirac points where the denominator vanishes, yielding a singularity in the summand. However, the numerator also vanishes at these points and, in fact, the limit of the summand at these points is finite.

Having run into only dead-ends with the stationary phase method, I though about approximating the sums as integrals in the [itex] N \rightarrow \infty [/itex] limit and performing successive integration by parts to extract the asymptotic behavior. (For sake of simplicity, in the following I take [itex]\mathbf{r} = \mathbf{r}_1 - \mathbf{r}_2[/itex] to be parallel to one of the lattice vectors). However, I seem to be getting terms proportional to
[tex] \frac{\sin^2{[\pi r]}}{r^{2n}}, \qquad n\in \mathbb{Z}^+,[/tex]
where the lattice spacing is taken to be unity. These terms all vanish as [itex] r [/itex] is an integer (it denotes a distance along one direction of the lattice).

It seems clear to me that I am not understanding the spirit of the calculation/am making math errors as all of my attempts suggest that the correlation function vanishes. The numerics (of the authors and of my own) show, however, that the power law behavior described above is correct.

Any help would be greatly appreciated. Thank you in advance.
 
Physics news on Phys.org
  • #2
I think they're just using an integral trick in order to approximate the Correlation function over a BZ. The function has a sum over q, which can be written as an integral in some limit. You could then rewrite the trigonometric terms as exponential functions and use this method: https://en.wikipedia.org/wiki/Stationary_phase_approximation. I think the integral will just produce a term proportional to x^-4. In this method you only evaluate the integral in an approximation, not exactly as it looks like you have done.
 
  • #3
Thank you for your response, DeathbyGreen.

Despite the authors claiming to have used the stationary phase approximation, I don't believe it applies as the cosine (or complex exponential) does not possesses a stationary phase. This is why I believe one should develop an asymptotic series with the boundary terms from successive integration by parts as outlined in, e.g., "Advanced Mathematical Methods for Scientists and Engineers" by Carl M. Bender and Steven A. Orszag.

Regarding the two zero points of the denominator mentioned in my original post, although the integrand remains finite at these points, I failed to realize initially that the integrand is discontinuous at these points and, thus, I was not correctly taking them into account when evaluating the boundary terms of the integration by parts.

The idea I then pursued was to express the integrand in cylindrical coordinates relative to the ill-behaved points and integrate by parts as (keeping in mind I also have to integrate everything with respect to q' as well)

[tex]
\begin{align}
\int_{BZ} d^2q \exp{[i \mathbf{q}\cdot\mathbf{r}]} \Big( \ldots \Big) &\rightarrow \lim_{\epsilon\rightarrow 0^+} \int_0^{2\pi} d\phi \int_\epsilon^\infty dq~q \exp{[i q r \cos{\phi}]} \Big( \ldots \Big) \nonumber\\
&= \lim_{\epsilon\rightarrow 0^+} \frac{1}{i r} \int_0^{2\pi} \frac{d\phi}{\cos{\phi}} ~q \exp{[i q r \cos{\phi}]} \Big(\ldots\Big) \Bigg|_{q=\epsilon}^\infty + \int d^2q \ldots \nonumber
\end{align}
[/tex]

I extended my q-integration to infinity and figured I would just ignore that part of the boundary term due to the rapidly oscillating phase. Then my hope was to take the [itex]\epsilon\rightarrow 0^+[/itex] limit before performing the angular integral, yielding a simpler integrand. Here the discontinuity of the integrand at [itex]\epsilon = 0[/itex] should now show up as a [itex]\phi[/itex]-dependent limit. Following this reasoning (and assuming my math beyond this point holds), the boundary term due to this first integration by parts vanishes in the limit due to an overall prefactor of [itex]\epsilon[/itex]. I would then integrate by parts once more yielding two boundary terms, one of which vanishes when I take the limit, the other of which does not. I then perform this integration by parts procedure twice more for the integral over q' with it's phase factor [itex]\exp{[-i\mathbf{q}'\cdot\mathbf{r}]}[/itex].

Having done all that I wind up with
[tex]
C(\mathbf{r}) \sim \frac{1}{r^4} \int_0^{2\pi} \frac{d\phi}{\cos^2{\phi}} \int_0^{2\pi} \frac{d\phi'}{\cos^2{\phi'}} \lim_{\epsilon\rightarrow 0^+} \lim_{\epsilon'\rightarrow 0^+} \Big( ... \Big)
[/tex]

At this stage I'd be happy if I could show that the remaining integration neither vanishes nor diverges. Unfortunately, it does vanish by anti-symmetry under [itex]\phi\rightarrow\phi+\pi[/itex].

Ok, if I'm making some simple mistake during all of my integrating by parts and limiting etc. then I need to be careful and find it. I don't expect anyone to go through all of that math for me hunting down my mistake. My question at this stage is more, is my method for approximating this integral conceptually sound? That is, integrating by parts, extending limits of integration and ignoring [itex]q\rightarrow\infty[/itex] contribution, etc.

Thanks again for reading and any help you all may have.
 
Last edited:

1. What is a correlation function?

A correlation function is a mathematical representation of the relationship between two variables. It measures how closely the values of one variable correspond to the values of another variable.

2. How is the asymptotic behavior of a correlation function calculated?

The asymptotic behavior of a correlation function is calculated by taking the limit as the number of data points approaches infinity. This is done by determining the trend of the correlation function as the number of data points increases, and then extrapolating this trend to an infinite number of data points.

3. Why is it important to calculate the asymptotic behavior of a correlation function?

Calculating the asymptotic behavior of a correlation function allows us to understand the long-term relationship between two variables. This can help us make predictions and identify patterns that may not be apparent in smaller data sets.

4. What factors can affect the asymptotic behavior of a correlation function?

The asymptotic behavior of a correlation function can be affected by the sample size, the distribution of the data, and the strength of the relationship between the two variables. It may also be influenced by outliers or other anomalies in the data.

5. Can the asymptotic behavior of a correlation function change over time?

Yes, the asymptotic behavior of a correlation function can change over time. This may be due to changes in the data, such as shifts in the relationship between the variables or the introduction of new data points. It is important to regularly recalculate the asymptotic behavior to ensure the accuracy of the correlation function.

Similar threads

  • Classical Physics
Replies
1
Views
784
  • Atomic and Condensed Matter
Replies
3
Views
733
  • Advanced Physics Homework Help
Replies
1
Views
2K
  • Electromagnetism
Replies
1
Views
1K
Replies
10
Views
1K
  • Atomic and Condensed Matter
Replies
1
Views
3K
  • Thermodynamics
Replies
7
Views
1K
  • Math Proof Training and Practice
Replies
28
Views
5K
  • Math Proof Training and Practice
Replies
16
Views
5K
  • High Energy, Nuclear, Particle Physics
Replies
5
Views
3K
Back
Top