piqs
Permutational Invariant Quantum Solver (PIQS) 
PIQS is not actively maintained anymore as an independent package, but it’s been ported to QuTiP and you can use its functionalities as qutip.piqs. Find out more in the related QuTiP documentation and tutorials on the permutational invariant Lindbald solver.
PIQS is an open-source Python library that allows to study the exact Lindbladian dynamics of open quantum systems consisting of identical qubits. The documentation for the package can be found in piqs.readthedocs.io. Example notebooks on how to use the library can be found here.
Exponential reduction
In the case where local processes are included in the model of a system’s dynamics, numerical simulation requires dealing with density matrices of exponential sizes. This becomes infeasible for a large number of qubits. We can simplify the calculations by exploiting the permutational invariance of indistinguishable quantum particles, which allows the user to study hundreds of qubits.
Integrated with QuTiP
A major feature of PIQS is that it allows to build the Liouvillian of the system in an optimal way. It uses Cython to optimize performance and by defining sparse matrices it can deal with large systems. Since it is compatible with the quantum object class of QuTiP one can take full advantage of existing features of this excellent open-source library.
A wide range of applications
- The time evolution of the total density matrix of quantum optics and cavity QED systems for permutationally symmetric initial states (such as the GHZ state, Dicke states, coherent spin states).
- Phase transitions of driven-dissipative out-of-equilibrium quantum systems.
- Correlation functions of collective systems in quantum optics experiments, such as the spectral density and second-order correlation functions.
- Various quantum optics phenomena such as resonance fluorescence, steady-state superradiance, superradiant light emission.
- Spin squeezing for quantum metrology, long-range interaction in noisy spin models, decoherence in quantum information processing.
- Nonlinearities of dissipative cavity QED systems up to the ultrastrong coupling regime.
Installation
PIQS is integrated inside QuTiP, from QuTiP version 4.3.1, as the qutip.piqs module. If QuTiP is installed, no additional installation is required. The PIQS functions can be called as
from qutip.piqs import *
PIQS can also be installed as a stand-alone library and it is still maintained.
Linux and MacOS
We have made the package available in conda-forge. Download and install Anaconda or Miniconda from https://www.anaconda.com/download/ and then install piqs from the conda-forge repository with,
conda config --add channels conda-forge
conda install -c conda-forge piqs
Then you can fire up a Jupyter notebook and run piqs out of your browser.
Windows
We will add a Windows version soon but if you are on Windows, you can build piqs from source by first installing conda. You can add the conda-forge channel with conda config --add channels conda-forge. Then, please install cython, numpy, scipy and qutip as piqs depends on these packages. The command is simply,
conda install cython numpy scipy qutip
Finally, you can download the source code, unzip it and run the setup file inside with,
python setup.py install
If you have any problems installing the tool, please open an issue or write to us.
Use
from piqs import Dicke
from qutip import steadystate
N = 10
system = Dicke(N, emission = 1, pumping = 3)
L = system.liouvillian()
steady = steadystate(L)
For more details and examples on the use of PIQS see the notebooks folder.
License
PIQS is licensed under the terms of the BSD license.
Documentation
PIQS documentation can be found at http://piqs.readthedocs.io/.
Notebooks
PIQS Notebooks include Jupyter notebooks for the paper https://arxiv.org/abs/1805.05129
- Superradiant light emission
- Steady-state superradiance
- Superradiant phase transition out of equilibrium
- Spin squeezing
- Ultrastrong light-matter coupling
- Multiple ensembles of qubits
- Boundary time crystals
- Performance of PIQS
A collection of Jupyter notebooks can be found at https://github.com/nathanshammah/notebooks and at http://qutip.org/tutorials.html under the “Permutational invariant Lindblad dynamics” section.
Citation
DOI:10.5281/zenodo.1212802 and N. Shammah, S. Ahmed, N. Lambert, S. De Liberato, and F. Nori, Open quantum systems with local and collective incoherent processes: Efficient numerical simulation using permutational invariance,
Phys. Rev. A 98, 063815 (2018)
https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.063815
DOI: 10.1103/PhysRevA.98.063815
https://arxiv.org/abs/1805.05129
Disclaimer
The Dicke.liouvillian() object is not to be used directly to compute its full spectrum or the Liouvillian gap.
Resources
Theoretical aspects and applications are in Ref. [1]. Other open-source codes using permutational invariance to study open quantum systems and related research papers can be found in [2-3].
[1] N. Shammah, S. Ahmed, N. Lambert, S. De Liberato, and F. Nori, Phys. Rev. A 98, 063815 (2018)
https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.063815
[2] https://github.com/peterkirton/permutations P. Kirton and J. Keeling Phys. Rev. Lett. 118, 123602 (2017)
[3] https://github.com/modmido/psiquasp M. Gegg and M. Richter, Sci. Rep. 7, 16304 (2017)