Simulating evaporative cooling of ultracold polar molecules.

Major progress has been made over the past few years in cooling polar molecules to ultracold temperatures. Common to these advances is evaporative cooling, where the gas distribution has its hot tail systematically truncated while elastic collisions work to maintain thermal equilibrium. However, the large 2 and 3-body loss mechanisms in kinetic molecular samples still pose a major hurdle for many groups in achieving stable quantum degenerate samples at their desired regimes. In light of this, we present a numerical tool for simulating evaporation in polar molecules. Our technique employs a Monte Carlo algorithm that allows for the inclusion of many-body quantum statistics. We also develop an efficient technique to include the accurate anisotropic energy dependent scattering cross sections for various regimes of operation, employing the use of Gaussian process interpolation. Our tools showcase good agreement with experimental data from the Max-Planck-Institute for Quantum Optics, paving the way for explorations of optimal evaporation schemes.