qrotor 4.2.1__tar.gz → 4.2.2__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Potentially problematic release.
This version of qrotor might be problematic. Click here for more details.
- {qrotor-4.2.1 → qrotor-4.2.2}/PKG-INFO +1 -1
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/_version.py +1 -1
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/constants.py +33 -2
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor.egg-info/PKG-INFO +1 -1
- {qrotor-4.2.1 → qrotor-4.2.2}/LICENSE +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/README.md +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/__init__.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/plot.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/potential.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/rotate.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/solve.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/system.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor/systems.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor.egg-info/SOURCES.txt +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor.egg-info/dependency_links.txt +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor.egg-info/requires.txt +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/qrotor.egg-info/top_level.txt +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/setup.cfg +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/setup.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/__init__.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/test_constants.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/test_potential.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/test_rotate.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/test_solve.py +0 -0
- {qrotor-4.2.1 → qrotor-4.2.2}/tests/test_system.py +0 -0
|
@@ -53,7 +53,7 @@ r_CH = distance_CH * np.sin(np.deg2rad(angle_CH)) * 1e-10
|
|
|
53
53
|
r_NH = distance_NH * np.sin(np.deg2rad(angle_NH)) * 1e-10
|
|
54
54
|
"""Rotation radius of the amine group, in meters."""
|
|
55
55
|
|
|
56
|
-
#
|
|
56
|
+
# Inertias, SI units
|
|
57
57
|
I_CH3 = 3 * (periodictable.H.mass * amu_to_kg * r_CH**2)
|
|
58
58
|
"""Inertia of CH3, in kg·m^2."""
|
|
59
59
|
I_CD3 = 3 * (periodictable.D.mass * amu_to_kg * r_CH**2)
|
|
@@ -63,6 +63,7 @@ I_NH3 = 3 * (periodictable.H.mass * amu_to_kg * r_NH**2)
|
|
|
63
63
|
I_ND3 = 3 * (periodictable.D.mass * amu_to_kg * r_NH**2)
|
|
64
64
|
"""Inertia of ND3, in kg·m^2."""
|
|
65
65
|
|
|
66
|
+
# Inertias of the co-rotation, SI units
|
|
66
67
|
I_CH3NH3 = I_CH3 + I_NH3
|
|
67
68
|
"""Inertia of CH3NH3+, in kg·m^2."""
|
|
68
69
|
I_CD3NH3 = I_CD3 + I_NH3
|
|
@@ -72,7 +73,17 @@ I_CH3ND3 = I_CH3 + I_ND3
|
|
|
72
73
|
I_CD3ND3 = I_CD3 + I_ND3
|
|
73
74
|
"""Inertia of CD3ND3+, in kg·m^2."""
|
|
74
75
|
|
|
75
|
-
#
|
|
76
|
+
# Inertias of the dis-rotation, SI units
|
|
77
|
+
I_CH3NH3_dis = 1 / ((1/I_CH3) + (1/I_NH3))
|
|
78
|
+
"""Inertia of the disrotatory torsion of CH3NH3+, in kg·m^2."""
|
|
79
|
+
I_CD3NH3_dis = 1 / ((1/I_CD3) + (1/I_NH3))
|
|
80
|
+
"""Inertia of the disrotatory torsion of CD3NH3+, in kg·m^2."""
|
|
81
|
+
I_CH3ND3_dis = 1 / ((1/I_CH3) + (1/I_ND3))
|
|
82
|
+
"""Inertia of the disrotatory torsion of CH3ND3+, in kg·m^2."""
|
|
83
|
+
I_CD3ND3_dis = 1 / ((1/I_CD3) + (1/I_ND3))
|
|
84
|
+
"""Inertia of the disrotatory torsion of CD3ND3+, in kg·m^2."""
|
|
85
|
+
|
|
86
|
+
# Inertias, atomic units
|
|
76
87
|
I_CH3_amu = I_CH3 / (amu_to_kg * 1e-20)
|
|
77
88
|
"""Inertia of CH3, in amu·AA^2."""
|
|
78
89
|
I_CD3_amu = I_CD3 / (amu_to_kg * 1e-20)
|
|
@@ -82,6 +93,7 @@ I_NH3_amu = I_NH3 / (amu_to_kg * 1e-20)
|
|
|
82
93
|
I_ND3_amu = I_ND3 / (amu_to_kg * 1e-20)
|
|
83
94
|
"""Inertia of ND3, in amu·AA^2."""
|
|
84
95
|
|
|
96
|
+
# Inertias of the co-rotation, amu units
|
|
85
97
|
I_CH3NH3_amu = I_CH3_amu + I_NH3_amu
|
|
86
98
|
"""Inertia of CH3NH3+, in amu·AA^2."""
|
|
87
99
|
I_CD3NH3_amu = I_CD3_amu + I_NH3_amu
|
|
@@ -91,6 +103,16 @@ I_CH3ND3_amu = I_CH3_amu + I_ND3_amu
|
|
|
91
103
|
I_CD3ND3_amu = I_CD3_amu + I_ND3_amu
|
|
92
104
|
"""Inertia of CD3ND3+, in amu·AA^2."""
|
|
93
105
|
|
|
106
|
+
# Inertias of the dis-rotation, amu units
|
|
107
|
+
I_CH3NH3_dis_amu = I_CH3NH3_dis / (amu_to_kg * 1e-20)
|
|
108
|
+
"""Inertia of the disrotatory torsion of CH3NH3+, in amu·AA^2."""
|
|
109
|
+
I_CD3NH3_dis_amu = I_CD3NH3_dis / (amu_to_kg * 1e-20)
|
|
110
|
+
"""Inertia of the disrotatory torsion of CD3NH3+, in amu·AA^2."""
|
|
111
|
+
I_CH3ND3_dis_amu = I_CH3ND3_dis / (amu_to_kg * 1e-20)
|
|
112
|
+
"""Inertia of the disrotatory torsion of CH3ND3+, in amu·AA^2."""
|
|
113
|
+
I_CD3ND3_dis_amu = I_CD3ND3_dis / (amu_to_kg * 1e-20)
|
|
114
|
+
"""Inertia of the disrotatory torsion of CD3ND3+, in amu·AA^2."""
|
|
115
|
+
|
|
94
116
|
# Rotational energy
|
|
95
117
|
_hbar = const.physical_constants['reduced Planck constant'][0]
|
|
96
118
|
B_CH3 = ((_hbar**2) / (2 * I_CH3)) * (1000 / const.eV)
|
|
@@ -111,6 +133,15 @@ B_CH3ND3 = ((_hbar**2) / (2 * I_CH3ND3)) * (1000 / const.eV)
|
|
|
111
133
|
B_CD3ND3 = ((_hbar**2) / (2 * I_CD3ND3)) * (1000 / const.eV)
|
|
112
134
|
"""Kinetic rotational energy of CD3ND3+, in meV·s/kg·m^2."""
|
|
113
135
|
|
|
136
|
+
B_CH3NH3_dis = ((_hbar**2) / (2 * I_CH3NH3_dis)) * (1000 / const.eV)
|
|
137
|
+
"""Kinetic rotational energy of the disrotatory torsion of CH3NH3+, in meV·s/kg·m^2."""
|
|
138
|
+
B_CD3NH3_dis = ((_hbar**2) / (2 * I_CD3NH3_dis)) * (1000 / const.eV)
|
|
139
|
+
"""Kinetic rotational energy of the disrotatory torsion of CD3NH3+, in meV·s/kg·m^2."""
|
|
140
|
+
B_CH3ND3_dis = ((_hbar**2) / (2 * I_CH3ND3_dis)) * (1000 / const.eV)
|
|
141
|
+
"""Kinetic rotational energy of the disrotatory torsion of CH3ND3+, in meV·s/kg·m^2."""
|
|
142
|
+
B_CD3ND3_dis = ((_hbar**2) / (2 * I_CD3ND3_dis)) * (1000 / const.eV)
|
|
143
|
+
"""Kinetic rotational energy of the disrotatory torsion of CD3ND3+, in meV·s/kg·m^2."""
|
|
144
|
+
|
|
114
145
|
# Potential constants from titov2023 [C1, C2, C3, C4, C5]
|
|
115
146
|
constants_titov2023 = [
|
|
116
147
|
[2.7860, 0.0130,-1.5284,-0.0037,-1.2791], # ZIF-8
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|