Tuning the range separation in LC-wPBE for organic electronics
Optimizing TADF emission
Long-range corrected density functional theory can ameliorate some of the known problems of DFT dealing with charge transfer states (e.g. push-pull systems). These systems are of great interest for organic electronics, and most recently for the so-called thermally activated delayed fluorescence (TADF) where triplet excitons in OLEDs slowly convert to singlet excitons which can fluoresce. To maximize luminescence yield, usually the singlet-triplet gap of is minimized, while maximizing spin-orbit coupling matrix elements.
The optimal range separation parameter gamma yields a HOMO energy corresponding to the ionization potential (IP). The optimal gamma will differ from molecule to molecule. Scanning over a number of gamma’s and then further refining is a typical task for which python scripting with PLAMS can make your life easier.
Tuning gamma for LC-wPBE
To get you started with tuning range-separation hybrids (RSHs) for your own systems of interest, we provide here an example script to tune LC-wPBE for the push-pull system p-Nitroaniline. In ADF2017 only RSHs with Yukawa switching, available in XCFun, are amenable to changing gamma. So we will work with LCY-wPBE, which differs from the LC-wPBE with error function switching. Yukawa switching actually has higher accuracy for atomization energies and Rydberg excitations, as shown by Akinaga and Ten-no.
Run the PLAMS script from command line, e.g. $ADFBIN/plams gammascan.py > gammascan.out
Tip: Windows and Mac users can use the pre-configured command line as explained here.
The script first defines the target to optimize. It minimizes simultaneously the difference between HOMO and IP for the neutral and the anion system, the latter being a proxy for optimizing the LUMO as close as possible to the fundamental gap, following the suggestion of Kronik and coworkers.
Then the job type is defined and the molecule read from the Molecule database shipped with the ADF Modeling Suite. We start with basic ADF settings and a large stride in gamma-space, followed by a further refinement around the found minimum with a better basis set and higher numerical accuracy.
The first optimum value for 0.3-0.4-0.5 is found to be 0.4:
gamma J 0.3000 0.01767257 0.4000 0.00307860 0.5000 0.01202306
In the last part we switch to a TZP basis and run a finer grid around the first minimum:
gamma J 0.3500 0.01303417 0.3700 0.00987240 0.3900 0.00691411 0.4100 0.00420414 0.4300 0.00207710 0.4500 0.00232538 Refined optimal gamma value: 0.4300
Make sure to reach out to our support team if you have any questions or problems about scripting for tuning RSHs or for similar tasks which could benefit from automation!
You may also be interested in checking out our advanced tutorial on optimizing TADF emitters including singlet-triplet gap and spin-orbit coupling calculations.
