Как исправить единицы СИ в этой программе? ⇐ Python

[Anonymous]

Я написал следующий сценарий для запуска MC-моделирования полимеров.
Скрипт производит нулевую полную энергию.
Что я такое делаете неправильно? Я подозреваю, что преобразование единиц измерения неверно. Однако я конкретно не знаю, что не так, поскольку я использовал наиболее подходящие единицы измерения и начальные значения.
Можете ли вы помочь мне найти ошибку?
< pre class="lang-py Prettyprint-override">import random
import math
import matplotlib.pyplot as plt
import datetime

# Global Parameters (Updated to SI Units)
polymer_chain_vec = []
N = 50 # Number of beads in the polymer chain
sim_steps_int = 100000 # Number of Monte Carlo steps
sigma = 1.0e-10 # Lennard-Jones sigma (m)
periodic_boundary_int = 1.0e-9 # Size of the periodic box (m)
cutoff_distance = 2.5 * sigma # Lennard-Jones cutoff distance (m)
max_atom_distance_float = 3.8e-10 # Maximum atom displacement (m)
epsilon = 1.0 # Lennard-Jones epsilon (J)
temperature_float = 500 # Temperature (K)
k_B = 1.380649e-23 # Boltzmann constant (J/K)
avogadro_number = 6.022e23 # Avogadro's number (mol^-1)

def write_txt(file_name_str, text_str):
"""Write text to a file."""
with open(file_name_str, 'a') as f:
f.writelines(text_str)

def apply_boundary_condition(point_pt):
"""Apply periodic boundary conditions to a point."""
x = point_pt[0] % periodic_boundary_int
y = point_pt[1] % periodic_boundary_int
return [x, y]

def get_distance(point_one_pt, point_two_pt):
"""Calculate the minimum image distance between two points."""
dx = abs(point_one_pt[0] - point_two_pt[0])
dy = abs(point_one_pt[1] - point_two_pt[1])

# Apply minimum image convention
dx = dx if dx

Подробнее здесь: https://stackoverflow.com/questions/792 ... is-program