7
$\begingroup$

In the recent publication:

"Density Functional Theory Calculation on the Structural, Electronic, and Optical Properties of Fluorene-Based Azo Compounds",

the author claimed to have used Gaussian16 on an Intel core i3-6006U CPU@2.00 GHz, 1.99 GHz. In my understanding of such basis sets, it is difficult to do this.

Is it possible to do TD-DFT calculations with 6-31+G(d,p) on a laptop?

$\endgroup$
  • 3
    $\begingroup$ The most important piece of information, the main memory, is missing. However, it would be possible to run a TDDFT calculation disk-based (don't know what Gaussian implements), as long as you can fit the matrix for the eigenvalue problem into memory. $\endgroup$ – TAR86 Aug 3 at 15:29
  • 2
    $\begingroup$ +1. But the size of the basis set is less important than the number of atoms. 6-31+G(d,p) is not so big if you're doing TD-DFT on the Li$_2$ molecule. The issue here is that there's 46 atoms in their largest XYZ file (see this PDF). $\endgroup$ – Nike Dattani Aug 3 at 15:50
  • 1
    $\begingroup$ @TAR86 Max memory supported by that CPU is 32 GB. Now, it seems more likely a laptop with a fairly low end CPU would have say 16 GB or less, but at least there is a strict constraint here. $\endgroup$ – Anyon Aug 3 at 16:41
11
$\begingroup$

Basis set name versus number of total orbitals

I would like to first address a part of the question that appears to be a misconception about the use of a 6-31+G(d,p) basis set, since you wrote:

"In my understanding of such basis sets, it is difficult to do this."

6-31+G(d,p) is not a "big" or "small" basis set, unless we're dealing with a big or small number of atoms (and aiming for a big or small amount of accuracy in our solution to the electronic Schrödinger equation). 6-31+G(d,p) is a fairly small basis set if we are dealing with a small molecule and only aiming to do TD-DFT. As currently written, the question only indicates that the paper was studying "azo" compounds, so the question alone does not tell us if the basis set is manageable or not manageable on a laptop computer. The question mentioned a journal publication, and we have to look into that paper to see what specific "azo" molecules are being discussed, before we can assess anything.

The calculation mentioned in the paper is possible on a laptop

The largest XYZ file in the articles Supplementary PDF is for a molecule with 46 atoms: bis(9H-fluoren-2-yl)diazene (which they call AzoFL). I just ran for you a CCSD(T)/6-31+G(d,p) calculation in MOLPRO, which seems to show that 6-31+G(d,p) seems manageable for a 46-atom molecule. I am aware the authors of the paper mentioned in the question, used Gaussian16, but I do not have access to Gaussian, so I used MOLPRO, which I do not think will affect whether or not the calculation is feasible on a laptop.

Since CSCD(T) is much more expensive than B3-LYP/TD-DFT (which is what is used in the paper), I conclude that B3-LYP/TD-DFT is feasible with the same or even smaller resources. The CCSD(T) calculation is still running, but the fact that it is running, is enough to know that it is possible (if the calculation is not possible due to lack of disk space or RAM, MOLPRO would have told us earlier in the output file).

What I have finished though is HF and MP2, and neither of them (nor the integrals) took more than 1000 seconds:

 PROGRAMS   *        TOTAL      RMP2       RHF       INT
 CPU TIMES  *      1872.67    488.34    832.73    551.49
 REAL TIME  *      3816.76 SEC
 DISK USED  *        84.74 GB

The disk space needed was only 85GB, which most laptops do have available. Unfortunately MOLPRO did not print the "peak RAM usage", and as user Anyon has pointed out, the CPU you mention in your question has a maximum capacity of 32GB of RAM, so this is the main limitation. I can re-run the calculation to only use 32GB of RAM and see what happens, but when the total disk usage is 85GB, it essentially means that a machine with 32GB of RAM will be enough, even if it means some of the disk has to act as SWAP space. I would not recommend doing this very often though, because if you're reading and writing 85GB too often, the lifetime of your disk might be affected.

Concluding remarks

  • It is possible to do a 46-atom TD-DFT-B3LYP/6-31+G(d,p) calculation on a laptop with an Intel Core i3-6006U CPU, but it requires 85GB of disk space, and I wouldn't recommend doing such calculations often if you care about the lifetime of your disk.
  • The name of the basis set, for example: 6-31+G(d,p), is not what matters in determining whether or not TD-DFT-B3LYP is possible, as much as the number of atoms, and more importantly the total number of orbitals.
| cite | improve this answer | |
$\endgroup$
  • 4
    $\begingroup$ @nike-datani I think even ~16-32 GB would be sufficient. G16 can optimize the disk space usage and if we solve for lower number of roots (~5), I think it is feasible, though one would have to save the restart *.wfn periodically if in case the hard drive fails! $\endgroup$ – mykd Aug 3 at 21:44

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.