I am performing a gaussian calculation of a cluster using ONIOM model. Here I am trying to use B3LYP method with explicit D3 correction. But I am facing a problem incorporating the D3 correction. The calculation execute normally when I use the method without incorporating the D3 correction. My route section without D3 correction reads as follows:

# opt=(maxcycle=300) freq ONIOM(B3LYP/aug-cc-pVTZ:hf/6-31g)

How to incorporate D3 correction. Any help in this regard is highly appreciated.

Thank you in advance