Thouless Pump Simulation with Local Uniform Bounded Disorder

This Python package simulates the Thouless pump of the Aubry-André model in the presence of local uniform bounded disorder. The simulation supports both single-particle and many-body initial states of non-interacting particles.

Features

• Simulation of Thouless pumping with configurable parameters
• Support for local uniform bounded disorder
• Parallel computation capabilities for multiple disorder strengths

Requirements

• Python 3.x
• NumPy
• SciPy
• Matplotlib
• joblib (for parallel processing)

Usage

Run the simulation using command line arguments:

python thouless_pump.py  

Arguments:

• log_pumping_rate: Logarithm of the pumping rate
• potential_strength: Strength of the potential (V)

Example:

python thouless_pump.py -2 1.0

Key Functions

• generate_hamiltonian(): Generates the Hamiltonian matrix for the Aubry-André model
• compute_pumping(): Computes the pumping dynamics for given parameters
• plot_time_evolution(): Visualizes the time evolution of the system
• run_simulations(): Handles parallel execution of multiple simulations
• plot_final_profile(): Creates visualizations of the final density profiles

Output

Results are saved in a directory structure organized by simulation parameters.

Parameters

Key configurable parameters include:

• beta: Incommensurate parameter
• V: Potential strength
• J: Hopping amplitude
• phi_0: Initial phase of the quasiperiodic modulation
• phi_end: Final phase of the quasiperiodic modulation
• n_sites: Number of lattice sites
• noise_amp_list: Range of disorder strengths to simulate

Author

Emmanuel Gottlob
Contact: emm.gottlob@gmail.com