Frequently Asked Questions¶
When Should I Change the Memory or Thread Count?¶
The memory and number of threads parameters can be very important for running the most efficient QM calculations. The defaults of 14GB and 8 CPUs were set to be conservative for our default methods and basis sets (HF-3c for geometry optimizations and B3LYP-D3MBJ/6-31G* for single point energies). for most drug like molecules (up to 30 heavy atoms). If you are going to run calculations with a non-default method or basis set, we recommend performing smaller scale benchmark calculations to make sure your resource settings are sufficient. How to Run Benchmark Floes for Cost Estimation
How do I check the metrics (memory, CPU usage) in a QM calculation?¶
By default, the memory and CPU usage metrics are turned on for Psi4 Property Calculation and Psi4 Geometry Optimization Cubes. That means you can see how much memory the cube actually used during a Floe by following these steps:
Go to the Floe tab, click on Jobs, then select your running or finished Floe.
Zoom in to the “Geometry Opt” or “Single-Point Energy” Cube in the Floe
Click on the ellipsis
Select Cube Metrics
A graph will appear showing how each cube metric changes over time.
You can deactivate a metric you do not want to see by clicking on its name in the top right.
When scrolling over the graph the value for each line appears, this is helpful since the units for each metric are different.
How Much Will My Floe Cost?¶
This is a difficult question to answer. The computational cost of QM calculations are going to depend on both the size of your molecule and the method and basis set you have chosen for your calculation. Before changing the method or basis set or moving to much larger molecules (>30 heavy atoms) it is best to run benchmark calculations. How to Run Benchmark Floes for Cost Estimation.
In our Floes, we have done our best to chose methods and basis sets that we believe have the best trade of for computational cost and accuracy. We have chosen HF-3c as a default method (with baked in minix basis set) for most of our calculations. The very small basis set makes calculations with HF-3c very affordable, and in our experience the geometries of optimized conformers are sufficient for most uses. You will see in Floes, such as Psi4 QM Conformer Ensemble and Psi4 QM Local Minima Search where we want to compare relative conformer energies we typically use a higher level of theory for the single point energy calculations.
What Method and Basis Set Should I use?¶
There is no universal right answer to this question. Before making changes, consider (1) if the default method and basis set are appropriate, (2) the type of calculation being performed, and (3) what elements and ions are included in your system.
The defaults chosen for each Floe have the goal of be affordable and reasonably accurate for most neutral drug like molecules. HF-3c – a corrected Hartree-Fock method from Sure and Grimme – is our default method (minix basis set) for most calculations ([Sure-2013]). We have found HF-3c to be very reliable for geometry optimizations. However, this basis set is quite small so if your molecule has atoms larger than Argon, you may want to consider a larger basis set. It is common practice to follow up geometry optimizations with single point energy calculations at a higher level of theory. Our default for these calculations is B3LYP-D3MBJ method with the 6-31G* basis set ([Grimme-2011]).
2. The type of calculation
Remember to think about the calculation being performed before changing the default method and basis set. If you are looking for a very accurate minimum geometry, then you may need a more sophisticated method and larger basis set. However, we have many Floes, such as Psi4 QM Conformer Ensemble, where all dihedral torsions are constrained during optimization. A more expensive method and basis set are unlikely to change the optimized structure in that case so the extra cost would be unjustified.
3. Heavy elements and ions
As mentioned above, the minix basis set (used to fit HF-3c) only really supports elements up to Argon. The 6-31G* basis set supports up to Krypton. Molecules with larger atoms, should use larger basis sets. In addition, most of our benchmarking has been performed with neutral molecules. If you are working with salts or other charged species then you may also need larger basis sets. Any time you change the defaults you should consider performing benchmark calculations to make sure the memory and thread count settings are sufficient and you know how expensive large scale jobs will cost.