14
$\begingroup$

Gaussian 16 claims to have an external tool to calculate sigma-profile COSMO/RS (Klamt) [1]. In the support website, there is a SCRF flag for COSMORS [2].

How do I calculate sigma-profile using Gaussian 16 and COSMO/RS? Are there any examples available?

Improve this question
$\endgroup$
4
  • 6
    $\begingroup$ SCRF(COSMORS) (gaussian.com/scrf/?tabid=2) will produce some kind of interface file. Then the manual of Cosmotherm details how to use these files (pdf). Related also: ccl.net/chemistry/resources/messages/2006/02/21.005-dir/… $\endgroup$ Commented May 12, 2020 at 12:19
  • $\begingroup$ Hi @Martin-マーチン can you expand that into an answer? $\endgroup$ Commented May 18, 2020 at 1:18
  • $\begingroup$ I won't, I have too little experience, and I cannot test this. If anyone wants to use (parts of or all of) the information above to turn that into an answer, I'd welcome that very much. $\endgroup$ Commented May 19, 2020 at 9:22
  • 1
    $\begingroup$ Does Guassians implementation still have the typo en.wikipedia.org/wiki/COSMO_solvation_model? making it correct for ions but not neutrals? $\endgroup$ Commented Dec 6, 2021 at 18:02

1 Answer 1

Reset to default
6
$\begingroup$

Gaussian itself doesn't have any way of computing sigma-profiles on its own. As you noticed, it has the scrf=COSMORS keyword, but this just produces a file that is used by a dedicated program (e.g. COSMO-RS, COSMOtherm).

On the Gaussian side, it really is as straightforward as using that keyword. Here is a minimal example: 

#p B88-VWN-P86/TZVP scrf=cosmors peroxide 0 1 O -0.22216884 -0.6214829 -0.08709297 H 0.09625091 -1.07514452 0.69674375 O 0.22216884 0.6214829 -0.08709297 H -0.09625091 1.07514452 0.69674375 h2o2.cosmo

This produces the necessary .cosmo file. Note the choice of functional and basis set is important, as COSMOtherm only ships parameter sets for certain combinations. In this case, you need to be sure the geometry is converged at this level of theory to use the BP_TZVP_C30_1701.ctd parameter file. If you can't converge it at that level of theory, you can use the SVP basis set and a different parameter file (BP_SVP_AM1_C30_1701.ctd).

I don't have access to either of the dedicated COSMO programs, but at least the input for COSMOtherm seems pretty straightforward. In an input file (name.inp), you would at minimum write: 

ctd=name.ctd #Identifier f=name.cosmo

This would be run using cosmotherm name.inp from the commandline.

Improve this answer
$\endgroup$
0

You must log in to answer this question.

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.