QuaCCAToo is a Python library for simulating and analyzing spin dynamics of color centers for quantum technology applications.
The systems' time evolution under pulsed experiments are calculated through quantum master equations based on the provided Hamiltonian, with realistic pulses in the laboratory frame.
The software is built on top of QuTip, inheriting its object-oriented framework and the Qobj
class.
For learning more about the package, we recommend first checking the tutorials section.
If you used QuaCCAToo in your work, please cite this repository and 10.1103/PhysRevA.111.022606.
- Repository: https://github.com/QISS-HZB/QuaCCAToo
- Documentation: https://qiss-hzb.github.io/QuaCCAToo/
- PyPI: https://pypi.org/project/QuaCCAToo/
We strongly recommend using a virtual environment (use whichever tool like venv
/conda
/uv
that you prefer) so that the system Python remains untouched.
pip install quaccatoo
Check here for detailed installation instructions.
- L. Tsunaki, A. Singh, K. Volkova, S. Trofimov, T. Pregnolato, T. Schröder, & B. Naydenov. (2025). Ambiguous resonances in multipulse quantum sensing with nitrogen-vacancy centers. Physical Review A, 111(2), 022606. doi: 10.1103/PhysRevA.111.022606
- S. Trofimov, C. Thessalonikios, V. Deinhart, A. Spyrantis, L. Tsunaki, K. Volkova, K. Höflich, & B. Naydenov. (2025). Local nanoscale probing of electron spins using NV centers in diamond. arXiv:2507.13295 quant-ph. 2507.13295.
QuaCCAToo is an object-oriented package organized with the following classes:
-
QSys
defines the quantum system of the problem. It has an obligatory intrinsic internal Hamiltonian$H_0$ , optional initial state, observable and a set of collapse operators. OnQSys
, calculates the eigenstates and eigenvalues of the system and has methods for truncating the systems and adding other spins. QuaCCAToo providesNV
as a predefined system for nitrogen vacancy centers in diamonds, more systems will be provided soon. -
PulsedSim
contains the logic for performing the simulation of pulsed experiments upon aQSys
object. It has attributes of a pulse sequence containing a set of pulses and free evolutions, control Hamiltonian$H_1$ , experiment variable and simulation results. Many predefined common pulse sequences are given inpredef_seqs
andpredef_dd_seqs
modules. Different pulse shapes are predefined in thepulse_shapes
module. -
ExpData
is a class to load experimental data and perform basic data processing, such as rescaling, subtracting columns or performing polynomial baseline corrections. -
Analysis
can be used either on simulation or experimental results, with a series of methods like for fitting (based onlmfit
), Fourier transforms and data comparison. The class can also used for plotting the results in multiple forms, including density matrix histograms and Bloch spheres. Several fit models and functions relevant for analysis of color centers are provided in thefit_functions
module.
Any contribution or bug report are welcome.
- To contribute, fork the main branch and make a pull request.
- We use
hatch/hatchling
as the build backend. The other development dependencies includepytest
andruff
. They can be installed by runningpip install -e '.[dev]'
from within the cloned repository. See here for details. - Properly document everything in details following the
numpy
docstring format. - Test your branch by running
pytest
and the tutorial notebooks. Feel free to add more tests. - Please pay attention to linter warnings (
ruff check
) and format your code withruff format
. - Module level refactors require corresponding changes in the
sphinx
setup, too. - Use US-English, not British-English. Eg: analyze instead of analyse, color instead of colour, center instead of centre.