The stochastic self-consistent harmonic approximation (SSCHA) is a full computational python package that simulates thermodynamic and phononic properties of materials accounting for anharmonicity at a nonperturbative level, fully including quantum and thermal fluctuations.
If you are simulating transport or thermal properties of materials, phase diagrams, or phonon-related properties, then you need the SSCHA code. The method is based on a full-quantum variational method that optimizes the nuclear wave-functions (or density matrix at finite temperature) to minimize the free energy. In this way it can be used to calculate thermodynamic properties of strongly anharmonic systems and the conditions at which first- and second-order phase transitions occur, estimate the phonon frequencies expected experimentally, obtain spectral functions, phonon lifetimes, determine the structural changes imposed by quantum effects, the electron-phonon interaction in strongly anharmonic systems, and so on.
With the SSCHA code, you can easily obtain the following properties:
Phase diagrams (first and second-order phase transitions; Ferroelectric and Peierls transitions).
Thermodynamic properties at finite temperature (Entropy, thermal expansion, bulk modulus).
Vibrational spectroscopy (Raman & IR).
Thermal transport in presence of strong anharmonicity.
Electron-phonon coupling and superconductivity.
The method works even in presence of strong anharmonicity and with sizable quantum fluctuations of nuclei. Running each simulation is extremely easy.