Recent content by dieon

I'm currently working on reproducing a graph from the Gutzwiller-RVB theory of high-temperature superconductivity using Python. The graph plots χ_BBxy versus Δ, where Δ ranges from 16 to 26. I've implemented the theory's iterative parameter update scheme but encountered discrepancies in the...

I was supposed to use thee formulas: g_t = 2 * x / (1 + x) g_s = 4 / (1 + x)**2 After that the issue got resolved

You're right, the steps do get smaller. However, at the same time, the value of \Delta at x=0 also gets smaller, which should not happen since its value should be near 0.35 at x = 0. When N_1 and N_2 are 4x, \Delta is nearing 0.05. Is there something else I can try?

I tried solving for a single value of x but got the same result. And, no, I don't think the iteration is converging, but I can't find a reason why. I think it has something to do with the \mu values, but after so many trials I still can't correct it. I think this curve is supposed to be smooth...

I’ve been working on solving gap equations related to superconductivity using Python, in the context of RMFT using the t-J model of the Hubbard Hamiltonian. My goal is to calculate and plot the variables (\Delta_x), (\Delta_y), (\xi_x), and (\xi_y) against the doping concentration (x). However...