9
$\begingroup$

I'm running a series of calculation for nonlinear chromophores. The general pipeline is geometry optimization (B3LYP/def2-TZVP) -> Frequencies -> Hyperpolarizability calculation (M06-2X (or CAM-B3LYP or MP2)/aug-cc-TZVP).

At the moment I plan to run Orca for optimization and frequencies because it is faster and GAMESS-US for hyperpolarizability, but I want to switch to the software that can do both tasks.

I can run anything that is free (at least for academia) and can be parallelized. Please help me find a suitable software, I could only think of NWChem/LSDalton (but that's also two separate codes).

$\endgroup$
5
  • 5
    $\begingroup$ Just out of curiosity: why is it a problem running two different (free) software? It doesn't seem very difficult to automate your workflow. $\endgroup$
    – Greg
    Commented Dec 5, 2020 at 18:16
  • $\begingroup$ It's mostly my perfectionism, it would be better to reformulate the question as "What is a good software for DFT and MP2 calculations of hyperpolarizability with various functionals?". AFAIK Orca can calculate only polarizability. $\endgroup$
    – Roman
    Commented Dec 6, 2020 at 6:20
  • 2
    $\begingroup$ I haven't used this program, so I wouldn't consider this a full answer, but Dalton should be able to do frequencies and hyperpolarizability. $\endgroup$
    – Tyberius
    Commented Dec 7, 2020 at 1:15
  • $\begingroup$ I've finally gave up and used Gaussian, since my university has a license. $\endgroup$
    – Roman
    Commented Jan 23, 2021 at 13:27
  • $\begingroup$ Maybe you can give a self-answer showing how you did it in GAUSSIAN, with minimalistic input files for the things that your questions asked for? $\endgroup$ Commented Feb 28, 2021 at 5:19

1 Answer 1

5
$\begingroup$

So I've used Gaussian as it is one of the few codes that can calculate hyperpolarizability for a lot of DFT functionals.

My input is very staightforward, hope it will help someone.

%NProcShared=24
%mem=192Gb
%Chk=optfreq.chk
#P B3LYP/Def2TZVP Opt Freq

 Optfreq

0 1
*Here goes initial geometry*

--Link1--
%NProcShared=24
%mem=192Gb
%OldChk=optfreq.chk
%Chk=polar.chk
#P M062X/aug-cc-pVDZ Guess=Read Geom=Checkpoint GFINPUT POP(Regular, NBO)  Polar

 Polar

0 1

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .