mirror of
https://github.com/TheAlgorithms/Python.git
synced 2024-11-27 15:01:08 +00:00
Compare commits
4 Commits
3aaf6dc088
...
f740be8ac8
Author | SHA1 | Date | |
---|---|---|---|
|
f740be8ac8 | ||
|
f3f32ae3ca | ||
|
4200aa200b | ||
|
4fbdd8e171 |
|
@ -16,7 +16,7 @@ repos:
|
||||||
- id: auto-walrus
|
- id: auto-walrus
|
||||||
|
|
||||||
- repo: https://github.com/astral-sh/ruff-pre-commit
|
- repo: https://github.com/astral-sh/ruff-pre-commit
|
||||||
rev: v0.7.3
|
rev: v0.7.4
|
||||||
hooks:
|
hooks:
|
||||||
- id: ruff
|
- id: ruff
|
||||||
- id: ruff-format
|
- id: ruff-format
|
||||||
|
|
48
physics/translation_partition
Normal file
48
physics/translation_partition
Normal file
|
@ -0,0 +1,48 @@
|
||||||
|
import math
|
||||||
|
|
||||||
|
def translation_partition_function(mass: float, temperature: float,
|
||||||
|
volume: float) -> float:
|
||||||
|
"""
|
||||||
|
Calculates the translational partition
|
||||||
|
function using the formula:
|
||||||
|
q_trans = ((2 * pi * m * k_B * T)^(3/2) * V) / h^3
|
||||||
|
|
||||||
|
mass: Mass of the molecule in kg
|
||||||
|
temperature: Temperature in Kelvin
|
||||||
|
volume: Volume in cubic meters
|
||||||
|
|
||||||
|
>>> round(translation_partition_function(2e-26, 300, 1e-3), 4)
|
||||||
|
4.081816847078438e+28
|
||||||
|
>>> round(translation_partition_function(4e-26, 300, 1e-3), 4)
|
||||||
|
1.1545121488522627e+29
|
||||||
|
>>> round(translation_partition_function(-4e-26, 300, 1e-3), 4)
|
||||||
|
Traceback (most recent call last):
|
||||||
|
...
|
||||||
|
ValueError: Mass must be positive
|
||||||
|
>>> round(translation_partition_function(4e-26, 0, 1e-3), 4)
|
||||||
|
Traceback (most recent call last):
|
||||||
|
...
|
||||||
|
ValueError: Temperature must be positive
|
||||||
|
>>> round(translation_partition_function(4e-26, 300, 0), 4)
|
||||||
|
Traceback (most recent call last):
|
||||||
|
...
|
||||||
|
ValueError: Volume must be positive
|
||||||
|
"""
|
||||||
|
|
||||||
|
if mass <= 0:
|
||||||
|
raise ValueError("Mass must be positive")
|
||||||
|
if temperature <= 0:
|
||||||
|
raise ValueError("Temperature must be positive")
|
||||||
|
if volume <= 0:
|
||||||
|
raise ValueError("Volume must be positive")
|
||||||
|
|
||||||
|
h = 6.62607015e-34 # Planck's constant in J·s
|
||||||
|
k_B = 1.380649e-23 # Boltzmann constant in J/K
|
||||||
|
|
||||||
|
prefactor = (2 * math.pi * mass * k_B * temperature) ** (3 / 2)
|
||||||
|
denominator = h ** 3
|
||||||
|
return (prefactor * volume) / denominator
|
||||||
|
|
||||||
|
if __name__ == "__main__":
|
||||||
|
import doctest
|
||||||
|
doctest.testmod(name="translation_partition_function")
|
Loading…
Reference in New Issue
Block a user