NOTE: This project is still in a work-in-progress state.
This software aims to simulate quantum lattice gas algorithms using graphics processing units.
We have included copies of PhD dissertations from Steven Smith and Jasper Taylor in the References folder. These can be used as references for both the modeling methodology and inplementation.
Hardware:
Main software:
Python Packages:
After downloading all the packages, look in the base directory you will find two files "Master.py" and "SimulationMaster.py". These are the main files that will need to be used to run simulations. "SimulationMaster.py" offers more flexibility, such as setting use of specific GPUs or setting up loops to run multiple simulations in succession. "Master.py" is an early version of a simple GUI, to facilitate use by inexperienced coders. Be aware that you may have to use "SimulationMaster.py" first, to generate .pyc files, before using "Master.py" GUI.
Here is a list of some of the built-in models included in the current version of the software. The code is modular, so that a user can also add their own models.
Make sure that the particle number indicator (1P,2P) at the end of the model name matches the number of particles you are simulating. Options described using the format:
Function(Keyword options)
Model and intial condition keywords can both be input into the "Experiment Keywords" field using python dictionary syntax, i.e. {'key':value}. Function based methods, such as Function_1P take lists created by the user (see Additional info). Some methods may works across different models, but no guarentee.
Models
Free_1P(None)
Free_2P(None)
Initial Conditions
Function_1P(func, px, py, pz)
- func is a user specified function written in CUDA/C++ syntax using variables vectorSize, X, Y, Z, Lx, Ly, and Lz.
- px, py, & pz are integers that give momentum kicks in the x, y, & z directions respectively.
Gaussians_1D_1P(momentums, shifts, sigmas)
- momentums are a list of N momentum kicks for N gaussians.
- shifts are a list of N discribution centers for N gaussians.
- sigmas are a list of N widths for N gaussians.
Gaussians_1D_2P(momentums, shifts, sigmas)
- momentums are a list of N momentum kicks for N gaussians.
- shifts are a list of N discribution centers for N gaussians.
- sigmas are a list of N widths for N gaussians.
Gaussians_2D_2P(momentums_x, momentums_y, shifts_x, shifts_y, sigmas_x, sigmas_y)
- momentums_x/y are a list of N momentum kicks along x/y for N gaussians.
- shifts_x/y are a list of N discribution centers along _x/y for N gaussians.
- sigmas_x/y are a list of N widths along x/y for N gaussians.
Visualization
Density_1D_1P(None)
Density_1D_2P(save_density, save_entanglement, save_expectation_value, exp_range)
- save_density saves the density if set to True.
- save_entanglement saves the entanglement if set to True.
- save_expectation value saves the expectation value over the range exp_range=[lower bound, upper bound]if set to True.
Density_Phase_2D_1PNone)
Projection_1D_2P(save_density, save_entanglement, save_expectation_value, exp_range)
- save_density saves the density if set to True.
- save_entanglement saves the entanglement if set to True.
- save_expectation value saves the expectation value over the range exp_range=[lower bound, upper bound]if set to True.
Projection_2D_2P(save_density, save_entanglement, save_expectation_value, exp_range)
- save_density saves the density if set to True.
- save_entanglement saves the entanglement if set to True.
- save_expectation value saves the expectation value over the range exp_range=[lower bound, upper bound]if set to True.
Models
SpinHalf_FC_1P(G0,G1,scaling)
- G0 & G1 are the effective coupling constants in the GP-equation
- scaling rescales the dimensions of the lattice
Initial Conditions
Function_1P(func, px, py, pz)
- func is a user specified function written in CUDA/C++ syntax using variables vectorSize, X, Y, Z, Lx, Ly, and Lz.
- px, py, & pz are integers that give momentum kicks in the x, y, & z directions respectively.
Skyrmion_Quandrupole(G0, G1, K, scaling, p1x, p1y)
- G0 & G1 are the effective coupling constants in the GP-equation.
- K/scaling scale the size of the Skyrmion/lattice respectively.
- p1x & p1y give momentum kicks in the x & y directions respectively.
Skyrmion_Octopole(G0, G1, K, scaling, p1x, p1y)
- G0 & G1 are the effective coupling constants in the GP-equation.
- K/scaling scale the size of the Skyrmion/lattice respectively.
- p1x & p1y give momentum kicks in the x & y directions respectively.
Visualization
3D_FC(None)
Models
Spin2_BEC_1P(G0, G1, G2, scaling)
- G0, G1, & G2 are the effective coupling constants in the GP-equation.
- scaling rescales the dimensions of the lattice.
Initial Conditions
Function_1P(func, px, py, pz)
- func is a user specified function written in CUDA/C++ syntax using variables vectorSize, X, Y, Z, Lx, Ly, and Lz.
- px, py, & pz are integers that give momentum kicks in the x, y, & z directions respectively.
Bright_Solitons(center1, center2, G0, G1, G2, scaling, state_1, state_2, p1x, p2x)
- center1 & center2 give locations of centers for soliton 1/2 respectively.
- G0, G1, & G2 are the effective coupling constants in the GP-equation
- p1x & p2x give momentum kicks for soliton 1/2 respectively.
- state_1, state_2 can be integers between 0-16, where 1-16 are different bright soliton solutions and 0 is no soliton.
- scaling rescales the dimensions of the lattice.
Dark_Soliton(center1, center2, G0, G1, G2, scaling, state_1, p1x)
- center1 & center2 give locations of centers for soliton kinks respectively.
- G0, G1, & G2 are the effective coupling constants in the GP-equation
- p1x give momentum kick for the dark soliton.
- state_1 can be integers between 0-16, where 1-16 are different dark soliton solutions and 0 is no soliton.
- scaling rescales the dimensions of the lattice.
Quadrupoles(G0, G1, G2, scaling, solution1, solution2, p1x, p1y, p1z, p2x, p2y, p2z, lr_shift1, lr_shift2, ud_shift1, ud_shift2, separation, orientation1, orientation2)
- G0, G1, & G2 are the effective coupling constants in the GP-equation.
- p1x & p2x (or y,z) give momentum kicks for quadrupole 1/2 respectively.
- lr_shift and ud_shift shifts the quadrupoles left/right and up/down.
- solution1/2 can be integers between 0-8, where 1-16 are different Pade approximates and 0 is no quadrupole.
- orientation can be either 'x','y', or 'z' and speficies the orientation.
- scaling rescales the dimensions of the lattice.
Visualization
1D_BEC(None)
2D_BEC(None)
3D_BEC(None)
Models
Function_1P(cond_list,func_list)
Function_2P(cond_list,func_list)
- cond_list is a list of conditions that determine which function to use in func_list (like a piecewise function).
- func_list is a list of CUDA/C++ syntax functions using variables vectorSize, X, Y, Z, T, Lx, Ly, and Lz.
External potentials and measurements are also utilized when running "SimulationMaster.py". To do this you change the "POTENTIAL" & "MEASUREMENT" parameters to match the name of the file of the desired method located in "/Code/Initializations/CUDA_physics_models/". Then change the kwargs to desired values, found in the method file.
- Steven Smith
- Jasper Taylor
- Anthony Gasbarro
This project is licensed under the MIT License - see the LICENSE.md file for details
- This work could not have been done without the guidance of Jeffrey Yepez.