def lennard_jones(r, sig, eps):
    """
    Calculate Lennard Jones potential for given distance, sigma, and epsilon values.
    Energy unit: (kB)
    """
    return 4 * eps * ((sig / r)**12 - (sig / r)**6)