Skip to content

Commit

Permalink
ZXY2
Browse files Browse the repository at this point in the history
  • Loading branch information
Trovemaster committed Jun 7, 2023
1 parent fce2771 commit 90ffd0b
Show file tree
Hide file tree
Showing 2 changed files with 51 additions and 0 deletions.
48 changes: 48 additions & 0 deletions docs/source/potentials.rst
Original file line number Diff line number Diff line change
Expand Up @@ -321,6 +321,54 @@ The vibrational coordinates are
\xi_3 &= \alpha-\alpha_{\rm eq},
\end{split}
where :math:`r_1 = r_1^{\rm eq}` and :math:`r_2 = r_2^{\rm eq}` are the internal stretching coordinates and :math:`\alpha` is the interbond angle, and the equilibrium parameters are :math:`r_1^{\rm eq}`, :math:`r_2^{\rm eq}` and :math:`\alpha_{\rm eq}`. Note that the exponent :math:`k` associated with the bending coordinate :math:`\xi_3` assumes only even values because of the symmetry of the XYZ molecule.


The input file example is `CaOH_Koput_step1.inp <https://raw.githubusercontent.com/Trovemaster/TROVE/develop/docs/source/input/CaOH_Koput_step1.inp>`_ where this PES is used.



ZXY2 type
---------


``poten_zxy2_morse_cos``
^^^^^^^^^^^^^^^^^^^^^^^^


The PEF is expressed analytically as (see [15AlOvPo]_, [23MeOwTe]_)

.. math::
V = \sum_{ijklmn} a_{ijklmn} \xi_{1}^{i} \xi_{2}^{j} \xi_{3}^{k} \xi_{4}^{l} \xi_{5}^{m} \xi_6^{n},
The vibrational coordinates are

.. math::
\begin{split}
\xi_i &= 1 - \exp[-b_i(r_i-r_i^{\rm e})],\quad i={\rm 1, 2, 3}, \\
\xi_4 &= \alpha_1-\alpha_{\rm e}, \\
\xi_5 &= \alpha_2-\alpha_{\rm e}, \\
\xi_6 &= 1+\cos\tau.
\end{split}
where :math:`a_{ijklmn}:math:` are the expansion parameters. Here, :math:`r_1:math:`, :math:`r_2:math:` and :math:`r_3:math:` are the bond lengths, :math:`\alpha_{1}:math:` and :math:`\alpha_{2}:math:` are the bond angles, :math:`\tau:math:` is the dihedral angle between two bond planes, and :math:`b_i:math:` is a Morse-oscillator parameter.














where :math:`r_1 = r_1^{\rm eq}` and :math:`r_2 = r_2^{\rm eq}` are the internal stretching coordinates and :math:`\alpha` is the interbond angle, and the equilibrium parameters are :math:`r_1^{\rm eq}`, :math:`r_2^{\rm eq}` and :math:`\alpha_{\rm eq}`. Note that the exponent :math:`k` associated with the bending coordinate :math:`\xi_3` assumes only even values because of the symmetry of the XYZ molecule.


Expand Down
3 changes: 3 additions & 0 deletions docs/source/references.rst
Original file line number Diff line number Diff line change
Expand Up @@ -68,3 +68,6 @@ References
.. [17HuScFr] X. Huang, D.W. Schwenke, R.S. Freedman, Т.J. Lee, J. Quant. Spectrosc. Radiat. Transf., 203, 224 (2017), Ames-2016 line lists for 13 isotopologues of CO\ :sub:`2`\ : Updates, consistency, and remaining issues.
.. [21ClYu] V.H. J. Clark and S.N. Yurchenko. Phys. Chem. Chem. Phys., 23, 11990 (2021), Modelling the non-local thermodynamic equilibrium spectra of silylene (SiH\ :sub:`2`\ ).
.. [23MeOwTe] T. Mellor, A. Owens, J. Tennyson, and S.N. Yurchenko. MNRAS, 520, 1997 (2023), ExoMol line lists - XLVIII. High-temperature line list of thioformaldehyde (H\ :sub:`2`\ CS).

0 comments on commit 90ffd0b

Please sign in to comment.