Can we optimize a molecular structure using Moller-Plesset MP4 method on GAUSSIAN software? If yes what keyword do we have to add?

  • 1
    $\begingroup$ Is this what you needed? $\endgroup$
    – PBH
    Oct 11, 2021 at 12:03
  • $\begingroup$ @PBH it is, so perhaps you could just write an answer with a sample input file. $\endgroup$ Oct 11, 2021 at 13:53
  • $\begingroup$ @PBH I believe it's very important to provide answers in a timely fashion (especially for new users!) if we want people to feel that it's worthwhile to continue using this new site. So I wrote an answer, but if you'd like to add anything or to provide an example input file, I'd highly encourage that, especially since we need to improve our answer-to-question ratio!. $\endgroup$ Oct 11, 2021 at 17:40
  • 4
    $\begingroup$ MP4 is quite outdated. While it used to be a recommendable method in the 1990s, today DFT methods (especially double-hybrid functionals) offer much more robust performance and much smaller computational resource demands than MP4. There is practically no reason to do any MP4 calculation nowadays, including but not limited to geometry optimization, unless you want to reproduce the results of an old paper. $\endgroup$
    – wzkchem5
    Oct 11, 2021 at 19:13
  • 1
    $\begingroup$ @NikeDattani, I was busy at the moment and didn't have sufficient time to write an answer. Also, I haven't used GAUSSIAN in a while (and am not an expert either) which is why I didn't try to answer the question as I didn't want to write an answer based on what I remember from some time back (which means that the answer could be wrong even if I wrote one). However, I believe that the link I provided should have been sufficient to get a start on solving the problem. $\endgroup$
    – PBH
    Oct 12, 2021 at 7:13

2 Answers 2


Gaussian can do all the way up to MP5 for single-point energies, but analytic gradients are only available for MP2, MP3, MP4(DQ) and MP4(SDQ), the latter including only single, double and quadruple substitutions.

MP4(SDTQ) which is like MP4(SDQ) but also with triple substitutions, and MP5, are listed in the manual as being able to do numerical gradients. Keep in mind that the Gaussian documentation already linked above says "MP4: Fourth-order MP theory correction [Raghavachari78], which defaults to full MP4 with single, double, triple and quadruple substitutions [Raghavachari78, Raghavachari80] (MP4(SDTQ))", meaning that MP4(SDTQ) and MP4 mean the same thing in Gaussian's terminology.

All of what I wrote above applies only to RHF and/or UHF reference states, since Gaussian only provides single-point energies for MP methods with an ROHF reference (you can see some discussion about that at the Chemistry SE here: Is it possible to perform ROMP2 numerical optimizations in Gaussian09?).

Beware that MP4 calculations are quite a lot slower than MP2 ones, which can make geometry optimization quite prohibitive depending on the size of your system, so you would probably want to use analytic gradients with RHF or UHF (instead of ROHF) and with MP4(DQ) or MP4(SDQ) rather than trying to do numericaly gradients in Gaussian with MP4 or numerical gradients manually on your own using ROMP4.

Furthermore, the MP series (MP2, MP3, MP4, MP5, etc.) is rarely used beyond MP2, largely because methods like CCSD can be far more accurate, and many DFT functionals may be more accurate while also being less computationally expensive (especially double-hybrid functionals which combine DFT with MP2-style calculations and can be astonishingly close to CCSD(T) in accuracy while being only a bit more expensive than MP2, see my answer to this: What are some recent developments in density functional theory?). So if you're about to launch 100s of geometry optimization calculations at the MP3 or MP4 level, I would encourage you to consider carefully the purpose of the project.

  • $\begingroup$ Thank you very much for your answer. From the links that you proposed I found one (somewhereville.com/2012/04/29/…) that talked about the same error that I had: ** ERROR IN INITNF. NUMBER OF VARIABLES ( 0) ** ** INCORRECT (SHOULD BE BETWEEN 1 AND 50) ** but when I tried the same input file that they proposed the calculation ended without the .chk file; is it normal? $\endgroup$
    – sarra6
    Oct 12, 2021 at 9:43
  • 2
    $\begingroup$ I solved the error thank you. The calculation ended noramally. $\endgroup$
    – sarra6
    Oct 12, 2021 at 10:03
  • 3
    $\begingroup$ Let's also not forget one of the biggest problems with the MP series: It's not even guaranteed to converge with higher perturbation orders. $\endgroup$
    – Antimon
    Oct 12, 2021 at 20:20
  • 1
    $\begingroup$ @Antimon thanks for the input! Sounds like we could benefit from having you around MMSE more often !!! $\endgroup$ Oct 12, 2021 at 20:41
  • 1
    $\begingroup$ I'll see what I can contribute! :) $\endgroup$
    – Antimon
    Oct 12, 2021 at 21:01

Complementing the @Nike answer:

Can we optimize a molecular structure using Moller-Plesset MP4 method on GAUSSIAN software?


If yes what keyword do we have to add?

You need to add the corresponding keywords like below:

# opt mp4(sdtq)/6-31++g(d,p) geom=connectivity

From Gaussian site:

This keyword requests that a geometry optimization be performed. The geometry will be adjusted until a stationary point on the potential surface is found. Analytic gradients will be used if available. For the Hartree-Fock, CIS, MP2, MP3, MP4(SDQ), CID, CISD, CCD, CCSD, QCISD, BD, CASSCF, and all DFT and semi-empirical methods, the default algorithm for both minimizations (optimizations to a local minimum) and optimizations to transition states and higher-order saddle points is the Berny algorithm using GEDIIS [Li06] in redundant internal coordinates [Pulay79, Fogarasi92, Pulay92, Baker93, Peng93, Peng96] (corresponding to the Redundant option).


The MPn method keywords request a Hartree-Fock calculation (by default, RHF for singlets, UHF for higher multiplicities) followed by a Møller-Plesset correlation energy correction (...) MP4: Fourth-order MP theory correction, which defaults to full MP4 with single, double, triple and quadruple substitutions (MP4(SDTQ)).

6-31++g(d,p): Is the chosen basis set. You need to select the basis set that better represent your system.


The Geom keyword specifies the source of the molecule specification input, options related to coordinate definitions, and geometry-related output. By default, it is read from the input stream, as described previously. Geom may be used to specify an alternate input source. It also controls what geometry-related information is printed and use of internal consistency checks on the Z-matrix. The Geom keyword is not meaningful without at least one item selection option.

  • $\begingroup$ Thank you very much for your answer I tried the same opt keyword but the calculation ended with the following error: ** ERROR IN INITNF. NUMBER OF VARIABLES ( 0) ** ** INCORRECT (SHOULD BE BETWEEN 1 AND 50) **. Do you have any suggestions to solve it? $\endgroup$
    – sarra6
    Oct 12, 2021 at 9:42
  • $\begingroup$ The above comment may now be irrelevant (see the comments on my answer). $\endgroup$ Oct 12, 2021 at 14:30
  • 3
    $\begingroup$ Concerning the basis set: If you're already going to the effort of an MP4 calculation with all substitutions, I'd suggest to choose a better basis set than a double-zeta Pople one to match that performance. But as mentioned in Nike's answer, MP4 itself is a somewhat questionable choice by today's standards. $\endgroup$
    – Antimon
    Oct 12, 2021 at 20:18

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.