This paper reports on the development of a collection of dynamical models of one-dimensional peak profile functions used to fit dynamic structure factors S (Q, ħω) of Quasi Elastic Neutron Scattering (QENS) data. The objective of this development is to create a maintainable and interoperable Python library with models reusable in other projects related to the analysis of data from Quasi Elastic Neutron Scattering experiments. The ambition is that the library also will serve as a platform where scientists can make their models available for others. We illustrate how the library can be used by newcomers to the field as well as by experts via different examples. These examples, provided as Jupyter notebooks, show how the QENS models can be integrated in the whole QENS data processing pipeline.