SciPost Phys. 10, 091 (2021) ·
published 26 April 2021
|
· pdf
We propose a method to improve the computational and memory efficiency of numerical solvers for the nonequilibrium Dyson equation in the Keldysh formalism. It is based on the empirical observation that the nonequilibrium Green's functions and self energies arising in many problems of physical interest, discretized as matrices, have low rank off-diagonal blocks, and can therefore be compressed using a hierarchical low rank data structure. We describe an efficient algorithm to build this compressed representation on the fly during the course of time stepping, and use the representation to reduce the cost of computing history integrals, which is the main computational bottleneck. For systems with the hierarchical low rank property, our method reduces the computational complexity of solving the nonequilibrium Dyson equation from cubic to near quadratic, and the memory complexity from quadratic to near linear. We demonstrate the full solver for the Falicov-Kimball model exposed to a rapid ramp and Floquet driving of system parameters, and are able to increase feasible propagation times substantially. We present examples with 262144 time steps, which would require approximately five months of computing time and 2.2 TB of memory using the direct time stepping method, but can be completed in just over a day on a laptop with less than 4 GB of memory using our method. We also confirm the hierarchical low rank property for the driven Hubbard model in the weak coupling regime within the GW approximation, and in the strong coupling regime within dynamical mean-field theory.
SciPost Phys. 10, 087 (2021) ·
published 23 April 2021
|
· pdf
We solve the quantum-mechanical antiferromagnetic Heisenberg model with spins positioned on vertices of the truncated icosahedron using the density-matrix renormalization group (DMRG). This describes magnetic properties of the undoped C$_{60}$ fullerene at half filling in the limit of strong on-site interaction $U$. We calculate the ground state and correlation functions for all possible distances, the lowest singlet and triplet excited states, as well as thermodynamic properties, namely the specific heat and spin susceptibility. We find that unlike smaller C$_{20}$ or C$_{32}$ that are solvable by exact diagonalization, the lowest excited state is a triplet rather than a singlet, indicating a reduced frustration due to the presence of many hexagon faces and the separation of the pentagonal faces, similar to what is found for the truncated tetrahedron. This implies that frustration may be tuneable within the fullerenes by changing their size. The spin-spin correlations are much stronger along the hexagon bonds and exponentially decrease with distance, so that the molecule is large enough not to be correlated across its whole extent. The specific heat shows a high-temperature peak and a low-temperature shoulder reminiscent of the kagomé lattice, while the spin susceptibility shows a single broad peak and is very close to the one of C$_{20}$.
SciPost Phys. 10, 083 (2021) ·
published 15 April 2021
|
· pdf
Fascinating structures have arisen from the study of the fractional quantum Hall effect (FQHE) at the even denominator fraction of $5/2$. We consider the FQHE at another even denominator fraction, namely $\nu=2+3/8$, where a well-developed and quantized Hall plateau has been observed in experiments. We examine the non-Abelian state described by the "$\bar{3}\bar{2}^{2}1^{4}$" parton wave function and numerically demonstrate it to be a feasible candidate for the ground state at $\nu=2+3/8$. We make predictions for experimentally measurable properties of the $\bar{3}\bar{2}^{2}1^{4}$ state that can reveal its underlying topological structure.