-
-
Notifications
You must be signed in to change notification settings - Fork 46.6k
/
Copy pathde_broglie
46 lines (37 loc) · 1.44 KB
/
de_broglie
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
PLANCK_CONSTANT_JS = 6.62607015e-34 # Planck's constant in Joule-seconds
PLANCK_CONSTANT_EVS = 4.135667696e-15 # Planck's constant in eV-seconds
def de_broglie_wavelength(momentum: float, in_ev: bool = False) -> str:
"""
Calculates the de Broglie wavelength of a particle using the given momentum.
Parameters:
momentum (float): Momentum of the particle.
in_ev (bool, optional): True if momentum is in eV·s.
If False, momentum is in kg·m/s.
Returns:
str: The calculated de Broglie wavelength of the particle in meters,
formatted in scientific notation.
Raises:
ValueError: If the momentum is zero or negative.
Usage example:
>>> de_broglie_wavelength(1e-24)
'6.62607015e-10'
>>> de_broglie_wavelength(1e-24, True)
'4.13566770e+09'
>>> de_broglie_wavelength(0)
Traceback (most recent call last):
...
ValueError: Momentum can't be zero or negative.
>>> de_broglie_wavelength(-1e-24)
Traceback (most recent call last):
...
ValueError: Momentum can't be zero or negative.
>>> de_broglie_wavelength(5e-20)
'1.32521403e-14'
"""
if momentum <= 0:
raise ValueError("Momentum can't be zero or negative.")
wavelength = PLANCK_CONSTANT_EVS / momentum if in_ev else PLANCK_CONSTANT_JS / momentum
return f"{wavelength:.8e}"
if __name__ == "__main__":
import doctest
doctest.testmod()