qrotor 4.2.1__py3-none-any.whl → 4.2.2__py3-none-any.whl

qrotor/_version.py

@@ -11,5 +11,5 @@ https://semver.org/
 ---
 """
 
-__version__ = "v4.2.1"
+__version__ = "v4.2.2"
 

qrotor/constants.py

@@ -53,7 +53,7 @@ r_CH = distance_CH * np.sin(np.deg2rad(angle_CH)) * 1e-10
 r_NH = distance_NH * np.sin(np.deg2rad(angle_NH)) * 1e-10
 """Rotation radius of the amine group, in meters."""
 
-# Inertia, SI units
+# Inertias, SI units
 I_CH3 = 3 * (periodictable.H.mass * amu_to_kg * r_CH**2)
 """Inertia of CH3, in kg·m^2."""
 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)
 I_ND3 = 3 * (periodictable.D.mass * amu_to_kg * r_NH**2)
 """Inertia of ND3, in kg·m^2."""
 
+# Inertias of the co-rotation, SI units
 I_CH3NH3 = I_CH3 + I_NH3
 """Inertia of CH3NH3+, in kg·m^2."""
 I_CD3NH3 = I_CD3 + I_NH3
@@ -72,7 +73,17 @@ I_CH3ND3 = I_CH3 + I_ND3
 I_CD3ND3 = I_CD3 + I_ND3
 """Inertia of CD3ND3+, in kg·m^2."""
 
-# Inertia, atomic units
+# Inertias of the dis-rotation, SI units
+I_CH3NH3_dis = 1 / ((1/I_CH3) + (1/I_NH3))
+"""Inertia of the disrotatory torsion of CH3NH3+, in kg·m^2."""
+I_CD3NH3_dis = 1 / ((1/I_CD3) + (1/I_NH3))
+"""Inertia of the disrotatory torsion of CD3NH3+, in kg·m^2."""
+I_CH3ND3_dis = 1 / ((1/I_CH3) + (1/I_ND3))
+"""Inertia of the disrotatory torsion of CH3ND3+, in kg·m^2."""
+I_CD3ND3_dis = 1 / ((1/I_CD3) + (1/I_ND3))
+"""Inertia of the disrotatory torsion of CD3ND3+, in kg·m^2."""
+
+# Inertias, atomic units
 I_CH3_amu = I_CH3 / (amu_to_kg * 1e-20)
 """Inertia of CH3, in amu·AA^2."""
 I_CD3_amu = I_CD3 / (amu_to_kg * 1e-20)
@@ -82,6 +93,7 @@ I_NH3_amu = I_NH3 / (amu_to_kg * 1e-20)
 I_ND3_amu = I_ND3 / (amu_to_kg * 1e-20)
 """Inertia of ND3, in amu·AA^2."""
 
+# Inertias of the co-rotation, amu units
 I_CH3NH3_amu = I_CH3_amu + I_NH3_amu
 """Inertia of CH3NH3+, in amu·AA^2."""
 I_CD3NH3_amu = I_CD3_amu + I_NH3_amu
@@ -91,6 +103,16 @@ I_CH3ND3_amu = I_CH3_amu + I_ND3_amu
 I_CD3ND3_amu = I_CD3_amu + I_ND3_amu
 """Inertia of CD3ND3+, in amu·AA^2."""
 
+# Inertias of the dis-rotation, amu units
+I_CH3NH3_dis_amu = I_CH3NH3_dis / (amu_to_kg * 1e-20)
+"""Inertia of the disrotatory torsion of CH3NH3+, in amu·AA^2."""
+I_CD3NH3_dis_amu = I_CD3NH3_dis / (amu_to_kg * 1e-20)
+"""Inertia of the disrotatory torsion of CD3NH3+, in amu·AA^2."""
+I_CH3ND3_dis_amu = I_CH3ND3_dis / (amu_to_kg * 1e-20)
+"""Inertia of the disrotatory torsion of CH3ND3+, in amu·AA^2."""
+I_CD3ND3_dis_amu = I_CD3ND3_dis / (amu_to_kg * 1e-20)
+"""Inertia of the disrotatory torsion of CD3ND3+, in amu·AA^2."""
+
 # Rotational energy
 _hbar = const.physical_constants['reduced Planck constant'][0]
 B_CH3 = ((_hbar**2) / (2 * I_CH3)) * (1000 / const.eV)
@@ -111,6 +133,15 @@ B_CH3ND3 = ((_hbar**2) / (2 * I_CH3ND3)) * (1000 / const.eV)
 B_CD3ND3 = ((_hbar**2) / (2 * I_CD3ND3)) * (1000 / const.eV)
 """Kinetic rotational energy of CD3ND3+, in meV·s/kg·m^2."""
 
+B_CH3NH3_dis = ((_hbar**2) / (2 * I_CH3NH3_dis)) * (1000 / const.eV)
+"""Kinetic rotational energy of the disrotatory torsion of CH3NH3+, in meV·s/kg·m^2."""
+B_CD3NH3_dis = ((_hbar**2) / (2 * I_CD3NH3_dis)) * (1000 / const.eV)
+"""Kinetic rotational energy of the disrotatory torsion of CD3NH3+, in meV·s/kg·m^2."""
+B_CH3ND3_dis = ((_hbar**2) / (2 * I_CH3ND3_dis)) * (1000 / const.eV)
+"""Kinetic rotational energy of the disrotatory torsion of CH3ND3+, in meV·s/kg·m^2."""
+B_CD3ND3_dis = ((_hbar**2) / (2 * I_CD3ND3_dis)) * (1000 / const.eV)
+"""Kinetic rotational energy of the disrotatory torsion of CD3ND3+, in meV·s/kg·m^2."""
+
 # Potential constants from titov2023 [C1, C2, C3, C4, C5]
 constants_titov2023 = [
     [2.7860, 0.0130,-1.5284,-0.0037,-1.2791],  # ZIF-8

{qrotor-4.2.1.dist-info → qrotor-4.2.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: qrotor
-Version: 4.2.1
+Version: 4.2.2
 Summary: QRotor
 Author: Pablo Gila-Herranz
 Author-email: pgila001@ikasle.ehu.eus

{qrotor-4.2.1.dist-info → qrotor-4.2.2.dist-info}/RECORD

@@ -1,20 +1,20 @@
 qrotor/__init__.py,sha256=rG2dH4QjsVUOMBhFnv5gXs3QnrUg7fywd5pIDmMBXcQ,246
-qrotor/_version.py,sha256=YSHr6C919Wy5fXxFzjIN6vSskUpecEZNWKepHCIAsy4,198
-qrotor/constants.py,sha256=tX_mM7WmP2wN4WYaL551jPcI-0Ca5SNOHATUxiUJQCQ,4934
+qrotor/_version.py,sha256=xJziIF1hUqXZ6UNnYfvk3UOaM4YT6tHCKcJ7dYAuyRQ,198
+qrotor/constants.py,sha256=YRri136Zg5dqqDS-180sxOB5ytBJW74msUbHYbrfiMc,6655
 qrotor/plot.py,sha256=loyd-1sXvfD8_PyHCThE3VQpRC0qs1iONGY4LXjJyZ4,14086
 qrotor/potential.py,sha256=2HjSDVJWSvwplyH5pq0MMPoQsLIBJbbTwkr350u949I,18491
 qrotor/rotate.py,sha256=Wje9Q9SFhDvizz58MzNGBwsmgV-3wN9z2SnUNTIXzeg,8107
 qrotor/solve.py,sha256=YkOR1SJlpk41PCNEhslv6X3wV1TWMNztT78qX3Pngf0,10722
 qrotor/system.py,sha256=ahYurNUmVOV7B6aZSe7rhcruagj5rW9UClqG-H1vVvY,11454
 qrotor/systems.py,sha256=Hcx0QvMWpaPMfC6HWpkZPPWDyHk9rxWKdAxWNnD2NMg,8184
-qrotor-4.2.1.dist-info/licenses/LICENSE,sha256=hIahDEOTzuHCU5J2nd07LWwkLW7Hko4UFO__ffsvB-8,34523
+qrotor-4.2.2.dist-info/licenses/LICENSE,sha256=hIahDEOTzuHCU5J2nd07LWwkLW7Hko4UFO__ffsvB-8,34523
 tests/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 tests/test_constants.py,sha256=YHKkPyZlzjchxxzON_VSNsQdKnpkknsFVoIA6TcUk70,399
 tests/test_potential.py,sha256=_Vq9t9Xm59kNbyYwXlRnvKcxwL7vntD2j14W2aUtF6I,1302
 tests/test_rotate.py,sha256=2On2d1E82hdisFC5DXpaqqYNnteX7ZP3PAnGa_oGm2M,1896
 tests/test_solve.py,sha256=tEjLUZC7oe6LCQD5b2xf2aaK9lu-zI4lzuPXOGR2GAs,861
 tests/test_system.py,sha256=36d-8AdoJdzq0O9_O3s8wwBPGa-M7A86YiHqhhAsCZ8,742
-qrotor-4.2.1.dist-info/METADATA,sha256=pkuA9iv1jSYiOUhxf9aYdRlKNHHyaLwJaiuWZACJyg0,9649
-qrotor-4.2.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-qrotor-4.2.1.dist-info/top_level.txt,sha256=mLnYs07-amqX4TqbDV2_XvgdpHfgrYmzmYb7dwoh6EQ,13
-qrotor-4.2.1.dist-info/RECORD,,
+qrotor-4.2.2.dist-info/METADATA,sha256=q11ZEF8wbyCfEJ36y_1KkO4Sxge0KnSMOSN78Ha4uKU,9649
+qrotor-4.2.2.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+qrotor-4.2.2.dist-info/top_level.txt,sha256=mLnYs07-amqX4TqbDV2_XvgdpHfgrYmzmYb7dwoh6EQ,13
+qrotor-4.2.2.dist-info/RECORD,,