benchmarks.wiki / Public workspace

CritPt / Challenge_16_main / Consider the following Hamiltonian in two dimensions on a square lattice:

Problem

Answer published by the source. Consult the official source to check your work against its answer.

problem description

# Problem setup: Consider the following Hamiltonian in two dimensions on a square lattice: \begin{equation} \begin{aligned} H &= 2\sum_{{\bf k}\sigma}(\cos{k_x}-\cos{k_y})(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}-c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma})\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{i\pi/4}(1+\text{e}^{i(k_y-k_x)})+\text{e}^{-i\pi/4}(\text{e}^{-ik_x}+\text{e}^{ik_y})]c^\dagger_{1{\bf k}\sigma}c_{2{\bf k}\sigma}\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{-i\pi/4}(1+\text{e}^{-i(k_y-k_x)})+\text{e}^{i\pi/4}(\text{e}^{ik_x}+\text{e}^{-ik_y})]c^\dagger_{2{\bf k}\sigma}c_{1{\bf k}\sigma} \\&-\mu\sum_{{\bf k}\sigma}(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}+c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma}) \\&+U\sum_{{\bf k}}(c^\dagger_{1{\bf k}\uparrow}c_{1{\bf k}\uparrow}c^\dagger_{1{\bf k}\downarrow}c_{1{\bf k}\downarrow}+c^\dagger_{2{\bf k}\uparrow}c_{2{\bf k}\uparrow}c^\dagger_{2{\bf k}\downarrow}c_{2{\bf k}\downarrow}), \end{aligned} \end{equation} where c†c^\dagger is the creation operator; cc is the annihilation operator; labels 1,2 denote sublattices 1 and 2; $\bf{k}$ is momentum; σ=↑or↓\sigma=\uparrow \text{or} \downarrow represents the spin; and μ\mu is the chemical potential with repulsive interaction U>0U>0. # Main problem: Calculate the critical interaction strength UcU_c that causes a phase transition at quarter-filling of the system. Please be as accurate as possible.
Plain-text mathematical notation (without MathML)
# Problem setup:
Consider the following Hamiltonian in two dimensions on a square lattice:
\begin{equation}
\begin{aligned}
H &= 2\sum_{{\bf k}\sigma}(\cos{k_x}-\cos{k_y})(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}-c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma})\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{i\pi/4}(1+\text{e}^{i(k_y-k_x)})+\text{e}^{-i\pi/4}(\text{e}^{-ik_x}+\text{e}^{ik_y})]c^\dagger_{1{\bf k}\sigma}c_{2{\bf k}\sigma}\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{-i\pi/4}(1+\text{e}^{-i(k_y-k_x)})+\text{e}^{i\pi/4}(\text{e}^{ik_x}+\text{e}^{-ik_y})]c^\dagger_{2{\bf k}\sigma}c_{1{\bf k}\sigma}
\\&-\mu\sum_{{\bf k}\sigma}(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}+c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma})
\\&+U\sum_{{\bf k}}(c^\dagger_{1{\bf k}\uparrow}c_{1{\bf k}\uparrow}c^\dagger_{1{\bf k}\downarrow}c_{1{\bf k}\downarrow}+c^\dagger_{2{\bf k}\uparrow}c_{2{\bf k}\uparrow}c^\dagger_{2{\bf k}\downarrow}c_{2{\bf k}\downarrow}),
\end{aligned}
\end{equation}
where c^(†) is the creation operator; c is the annihilation operator; labels 1,2 denote sublattices 1 and 2; $\bf{k}$ is momentum;   σ=↑or↓ represents the spin; and μ is the chemical potential with repulsive interaction U>0.

# Main problem:

Calculate the critical interaction strength U_(c) that causes a phase transition at quarter-filling of the system. Please be as accurate as possible.
Original LaTeX notation
# Problem setup:
Consider the following Hamiltonian in two dimensions on a square lattice:
\begin{equation}
\begin{aligned}
H &= 2\sum_{{\bf k}\sigma}(\cos{k_x}-\cos{k_y})(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}-c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma})\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{i\pi/4}(1+\text{e}^{i(k_y-k_x)})+\text{e}^{-i\pi/4}(\text{e}^{-ik_x}+\text{e}^{ik_y})]c^\dagger_{1{\bf k}\sigma}c_{2{\bf k}\sigma}\\&+\sqrt{2}\sum_{{\bf k}\sigma}[\text{e}^{-i\pi/4}(1+\text{e}^{-i(k_y-k_x)})+\text{e}^{i\pi/4}(\text{e}^{ik_x}+\text{e}^{-ik_y})]c^\dagger_{2{\bf k}\sigma}c_{1{\bf k}\sigma}
\\&-\mu\sum_{{\bf k}\sigma}(c^\dagger_{1{\bf k}\sigma}c_{1{\bf k}\sigma}+c^\dagger_{2{\bf k}\sigma}c_{2{\bf k}\sigma})
\\&+U\sum_{{\bf k}}(c^\dagger_{1{\bf k}\uparrow}c_{1{\bf k}\uparrow}c^\dagger_{1{\bf k}\downarrow}c_{1{\bf k}\downarrow}+c^\dagger_{2{\bf k}\uparrow}c_{2{\bf k}\uparrow}c^\dagger_{2{\bf k}\downarrow}c_{2{\bf k}\downarrow}),
\end{aligned}
\end{equation}
where $c^\dagger$ is the creation operator; $c$ is the annihilation operator; labels 1,2 denote sublattices 1 and 2; $\bf{k}$ is momentum;   $\sigma=\uparrow \text{or} \downarrow$ represents the spin; and $\mu$ is the chemical potential with repulsive interaction $U>0$.

# Main problem:

Calculate the critical interaction strength $U_c$ that causes a phase transition at quarter-filling of the system. Please be as accurate as possible.

code template

Code

def answer():
    r"""
    Return the value of the critical interaction strength $U_c$
    that causes a phase transition at quarter-filling of the system

    Inputs
    ----------
    None

    Outputs
    ----------
    U_c: float, critical interaction strength
    """

    # ------------------ FILL IN YOUR RESULTS BELOW ------------------
    U_c = ...
    # ---------------------------------------------------------------

    return U_c

Discussion

Discussion

No discussion posts on this page yet. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.

See answer Answer published by the source

Artifacts

Code, notes and reproducible work shared by participants. Files are served from a separate origin.

No artifacts on this page yet. Share reproducible code or notes in a contribution. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.

Source and history

Official source

initial import