Anderson Model HEOM Implementation - Code and Data

Théodore GOUMAÏ VEDEKOÏ, Jean-Pierre Tchapet Njafa, S. G. Nana Engo

Published Jan 11, 2026• 0 citations

Abstract

Complete computational implementation for the manuscript “Comparative Study of the Anderson Model inWeak and Strong Interaction Regimes: Implementations in Julia (HierarchicalEOM.jl) and Python (QuTiP)”. This repository contains: Python implementation using QuTiP Julia implementation using HierarchicalEOM.jl All Jupyter notebooks reproducing manuscript figures Temperature study scripts Complete documentation and installation guides All simulations are fully reproducible using the provided code.

DOI: 10.5281/zenodo.18213731

Related papers

An efficient Julia framework for hierarchical equations of motion in open quantum systems

Yi-Te Huang, Po-Chen Kuo, Neill Lambert +5 more

Abstract The hierarchical equations of motion (HEOM) approach can describe the reduced dynamics of a system simultaneously coupled to multiple bosonic and fermionic environments. The complexity of exactly describing the system-environment interaction with the HEOM method usually results in time-consuming calculations and a large memory cost. Here, we introduce an open-source software package called HierarchicalEOM.jl: a Julia framework integrating the HEOM approach. HierarchicalEOM.jl features a collection of methods to compute bosonic and fermionic spectra, stationary states, and the full dynamics in the extended space of all auxiliary density operators (ADOs). The required handling of the ADOs multi-indexes is achieved through a user-friendly interface. We exemplify the functionalities of the package by analyzing a single impurity Anderson model, and an ultra-strongly coupled charge-cavity system interacting with bosonic and fermionic reservoirs. HierarchicalEOM.jl achieves a significant speedup with respect to the corresponding method in the Quantum Toolbox in Python (QuTiP), upon which this package is founded.

arXiv:2306.07522Communications PhysicsOct 23, 202326 citations

Nonlinear optical spectra from Rydberg-mediated photon-photon interactions

Xinghan Wang, Yupeng Wang, Aishik Panja +1 more

This repository contains the data and simulation codes supporting the results presented in the manuscript: “Nonlinear optical spectra from Rydberg-mediated photon-photon interactions.” All simulations are implemented in Python and designed to be executed using Jupyter Notebook. Running all cells sequentially in each notebook will reproduce the figures reported in the manuscript.

Feb 11, 20260 citations

Code and datasets for "Shortcuts to adiabaticity with a quantum control field"

Emma C. King, Giovanna Morigi, Raphaël Menu

Shortcuts to adiabaticity with a quantum control field — Code & data========================================================== This repository contains a Jupyter notebook and supporting data files (csv) to reproduce the numerical results and figures in our manuscript (arXiv 2508.06304), which investigates adiabatic shortcuts mediated by a quantum control field.The code implements a Landau–Zener two-level system coupled to a spectator qubit, time evolution of the full system (LZ qubit + spectator) via QuTiP, parameter sweeps (parallelized with joblib), and plotting utilities for the manuscript figures. Related publication----------------------------This code is complementary to the preprint: https://arxiv.org/abs/2508.06304. Requirements---------------------The simulations were developed for and tested on Python 3.9+. Recommended packages and minimum versions: - Python >= 3.9- numpy >= 1.22- scipy >= 1.8- matplotlib >= 3.5- qutip >= 4.6 (install a recent stable qutip; some installations require `pip install qutip` or `conda` channels)- joblib >= 1.2- tqdm >= 4.60- pandas (optional, for CSV handling)- jupyterlab / notebook (for interactive use) Running the notebook (interactive) Open the notebook with JupyterLab or Jupyter Notebook. Execute cells sequentially in the order indicated by the notebook outline.The sections “Definitions” and “Figure imports and setup” must be executed before generating any data or figures, as they define all operators, and helper functions. The notebook is organized into separate sections for each figure of the manuscript, allowing you to re-run only the figure(s) of interest. Note: Figure 4 (Density plot) is configured to load precomputed data files to reduce computational cost. To reproduce Figure 4 without recomputing the full dataset, ensure the required .csv files are downloaded and located in the same directory as the notebook. To explore different parameter regimes: Modify the relevant parameter definitions (e.g., anneal_window, g, epsilon, omega_c, x0, x0_vals, q_vals). To include dissipation of the spectator (ancilla) qubit, set the parameter kappa> 0. For questions about reproducing the figures or about parameter choices used in the paper, contact the corresponding authors listed in the preprint. If you reuse this data/code, please cite the associated preprint and this code archive.

Feb 22, 20260 citations

Engineered non-Gaussian Coherence as a Thermodynamic Resource for Quantum Batteries

Kingshuk Adhikary

This repository contains the numerical data and computational codes supporting the results presented in the manuscript: "Engineered non-Gaussian Coherence as a Thermodynamic Resource for Quantum Batteries". The archive (data_QB.zip) includes all datasets, Jupyter notebooks, and figure files required to reproduce the results reported in both the main text and the Supplemental Material. Each figure is organised in a separate directory containing the corresponding data (.txt), code (.ipynb), and generated plots (.png/.pdf). The provided materials are sufficient to reproduce all figures and results without restriction. Detailed instructions and directory structure are described in the accompanying README.txt file. This dataset is made openly available to ensure transparency and reproducibility of the reported results.

Apr 11, 20260 citations

Controlling emergent dynamical behavior via phase-engineered strong symmetries

Marc Nairn

This record contains code, notebooks, data, and figure scripts associated with the project. The description below lists the contents of the upload and their roles. strongsym.zip: Zip archive used for the Zenodo upload. It contains the simul/, finsize_data/, and dependencies/ directories. simul/ : Simulation and analysis notebooks and scripts. MF_*.ipynb - Mean-field analysis notebooks. finite_size_spins.ipynb - Finite-size dynamics for 2 level system alt_finite_size_spins.ipynb - Alternative finite-size with a stress on Liouvullian spectrum analysis in the two level system. liouv_spec_spin1.ipynb -Finite-size spectrum analysis for three level model. find_symmetry.ipynb - Symbolic symmetry-finding notebook. finsize_data/ : Data files used to generate figures. fig_script.ipynb: Notebook to generate the figures overview_sketch.ipynb: Notebook to generate schmatics and diagrams (no raw data) dependencies/ : Environment specification files, for both python (.yml) and julia (.toml) Usage: Unzip strongsym.zip, set up environments using the files in dependencies/, and run the notebooks in simul/. Data required for figures is in finsize_data/; figure scripts are at the repository root. For a more detailed walkthrough and documentation, see the GitHub repository: https://github.com/MarcNairn/spin_strong_symmetries_phase

Feb 1, 20261 citation

Code for Quantum Simulation of Nuclear Dynamics in First Quantization

Luca Spagnoli, Chiara Lissoni, Alessandro Roggero

We provide a direct implementation in Python of the resource estimates discussed in the paper "Quantum Simulation of Nuclear Dynamics in First Quantization". There are two files: the python script contains routines that evaluate the specific gate and qubit counts described in the text while the Jupyter notebook provides an example on how to use these routines to produce the plots present in the paper.

May 5, 20261 citation