Because of the ability to recognize patterns behind complex phenomena, neural networks are starting to make their way into the natural and social sciences. In this work we will use neural networks as potentials with which we can calculate various properties of molecular crystals. More specifically, we will use potentials based on ANAKIN-ME methodology. On top of that, we added Van der Waals interaction with D4 model. Firstly, we optimized unit cells of molecular crystals with the help of ANI1x, ANI2x, ANI1x + D4, and ANI2x + D4 potentials, and calculated volumes which were then compared to the experimental data. After achieving reasonable agreement, we decided to use these potentials in order to calculate vibrational properties and Helmholtz free energy. Compound on which we have done the calculations is N’-2-propylidene-4-hydroxybenzohydrazide. This material is very special since it shows particularly rare thermosalient effect. In its crystal structure, it has 104 atoms making typical ab initio calculations very expensive. We proceeded to calculate phonon density of states and Helmholtz free energy. Although potentials do recognize some important features of this crystal like characteristic energy of vibrations of hydrogen atoms or phase transitions, we conclude that these potentials still cannot completely replace ab initio methods. Nevertheless, application of neural networks in solid state phyisics is a hot topic since they possibly represent extremely valuable tool for the field.