pkua007 / mbl_ed Goto Github PK

We need a class which is able to generate general hamiltonian for the problem - where hopping constants and diagonal elements should be programmable by classes inheriting after this base class.

There will be different analyzing modes. It should be organized. Analyzer will be a collection of tasks, some of them to be printed "inline", some to separate files.

We need a simple class which can compute level statistics for eigenenergies.

The code should have the whole necessary documentation.

Because we want to go higher than 10x10 ;)

One may don't want to obtain eigenvectors if they are now used or stored. The switch can save memory and enable to diagonalize larger matrices without overflowing int indices :')

We need additional terms in hamiltonian, which are higher order corrections.

Now the first specialized hamiltonian - the one from https://arxiv.org/pdf/1902.00357.pdf for bosons to test the program.

We need separate mode to analyze simulation results afterwards.

The last term after Z^2 + YZ + ZY to be added for cavity higher order considerations.

Make a working lexicographical Fock base generator.

After writing and testing all pieces it's time for first operable version of the program which will compute level statistics for a number of disorder realizations.

The elements of fock base should be searched very fast using hashing.

We need mode where we can look at inverse participation ratio of individual eigenvectors (so as a function of eigenenergy).

We need analyzer mode, where we compute the mean inverse of sum of 4th powers of eigenvector elements for eigenvalues from a given band (like in mean gap ratio).

We need an AnalyzerTask for performing time evolution for some observables: correlations, occupations, fluctuations, for some initial product states.

Let's add second simulation mode, where subsequent simulations take phi0 values from 0 to 2pi.