Skip to main content

Simulations of point-like particles in any dimension with any pair potential

Project description

The dompap simulation package

The dompap package focuses on simulations of point-like particles in any dimension with any pair potential. The package uses NumPy and Numba for efficient calculations and SymPy to implement any pair potentials. The user is not expected to be familiar with these packages but only basic Python syntax.

Installation

Use the Python package index (PyPI)

The package can be installed from the Python package index (PyPI) using pip.

pip install dompap

Download source from GitHub

Clone the repository from github at https://github.com/urpedersen/dompap.git (into a working directory of your choice), and add the package to your python path.

# Clone repository into some directory (replace [dir])
cd [dir]
git clone https://github.com/urpedersen/dompap.git

# Add to python path
export PYTHONPATH=$PYTHONPATH:$PWD/dompap

# Install requirements
pip install -r ./dompap/requirements.txt

Usage example

Python script

Below is an example of a 3D system of harmonic repulsive particles with the pair potential $$v(r) = (1 - r)^2$$ for $r<1$ and zero otherwise. The initial positions are set to a face-centered cubic (fcc) lattice, with five unit cells in each direction. The simulation is run for 100 steps (constant $NVT$ with Langevin thermostat), and the potential energy is printed every 10 steps.

from dompap import Simulation

# Initialize simulation object
sim = Simulation()

# Setup simulation
fcc_unit_cell = ([0.0, 0.0, 0.0], 
                 [0.5, 0.5, 0.0], 
                 [0.5, 0.0, 0.5], 
                 [0.0, 0.5, 0.5])
sim.set_positions(unit_cell_coordinates=fcc_unit_cell,
                  cells=(5, 5, 5), 
                  lattice_constants=(1.0, 1.0, 1.0))
sim.set_density(density=1.0)
sim.set_masses(masses=1.0)
sim.set_random_velocities(temperature=1.0)
sim.set_pair_potential(pair_potential_str='(1-r)**2', 
                       r_cut=1.0,
                       force_method='neighbor list', 
                       energy_method='neighbor list')
sim.set_pair_potential_parameters(sigma=1.0, epsilon=1.0)
sim.set_neighbor_list(skin=0.7, 
                      max_number_of_neighbors=128, 
                      method_str='double loop')
sim.set_integrator(time_step=0.01, 
                   target_temperature=1.0, 
                   temperature_damping_time=0.1)

# Run simulation
steps = 100
for step in range(steps):
    sim.step()
    if step % 10 == 0:
        print(f'Energy after {step} steps: {sim.get_potential_energy()}')

This simulation produces the output

Energy after 0 steps: 0.0
Energy after 10 steps: 3.1573922419447613
Energy after 20 steps: 16.330136084973663
Energy after 30 steps: 31.47341041787513
Energy after 40 steps: 43.913179017390576
Energy after 50 steps: 52.04197939534787
Energy after 60 steps: 57.968309542867964
Energy after 70 steps: 61.752276744879524
Energy after 80 steps: 67.33278804505039
Energy after 90 steps: 72.00507120397305

See examples for more examples of the capabilities of the dompap package.

Use dompap from the command line

The dompap package can also be used from the command line limited functionality. Type dompap --help for more information and to see the available options:

python3 -m dompap --help

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

dompap-0.0.5.tar.gz (855.9 kB view details)

Uploaded Source

Built Distribution

dompap-0.0.5-py3-none-any.whl (32.7 kB view details)

Uploaded Python 3

File details

Details for the file dompap-0.0.5.tar.gz.

File metadata

  • Download URL: dompap-0.0.5.tar.gz
  • Upload date:
  • Size: 855.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.11.2

File hashes

Hashes for dompap-0.0.5.tar.gz
Algorithm Hash digest
SHA256 24579526864c9cf27e7aabd0f31a41b07dcd02603523f24e878002cc515e36b9
MD5 ab954fe35029730dde2668a95339e48b
BLAKE2b-256 341c0d61ba0f323be1473c0464388c73fc772d7b9a45a474ba5511e665812789

See more details on using hashes here.

File details

Details for the file dompap-0.0.5-py3-none-any.whl.

File metadata

  • Download URL: dompap-0.0.5-py3-none-any.whl
  • Upload date:
  • Size: 32.7 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.11.2

File hashes

Hashes for dompap-0.0.5-py3-none-any.whl
Algorithm Hash digest
SHA256 4a28e0152075d8160a0a48d909f29e26fb07e51b25a0e0a2a618a8bbc78ef496
MD5 d0af250f37d77daa16036016a8ec147b
BLAKE2b-256 896de1eb649a61e2666e8f0545c33198daa26e8af15f4e7a2c91bc2ee0ac2750

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page