From: Vermaas, Josh (
Date: Wed Apr 14 2021 - 18:35:05 CDT

What is the output of the log? Everything looks plausible in the inputs (although isn’t the ethanol box 82x82x82?).

What I would try is to do this once through the tkconsole and compare with the GUI output:

package require solvate
solvate input.psf input.pdb -o solvated -b 1.4 -minmax [list [list -40 -40 -40] [list 60 60 60]] -spsf etohbox.psf -spdb etohbox.pdb -stop top_all36_cgenff.rtf -ws 82 -ks “name C1”

This is what I used, so it should totally work.


From: Raisa Carmen Andeme Ela <>
Date: Wednesday, April 14, 2021 at 4:55 PM
To: "" <>, "Vermaas, Josh" <>
Subject: VMD: ethanol solvation box

Hello VMD Family,

I have been trying to create a solvation box around a macromolecule (lignin); I am utilizing pre-equilibrated boxes created by Vermaas et al (10.1021/acssuschemeng.0c07156<;!!HXCxUKc!lWQkTBLlZyHd1SwX6Q8WIiQqvQHK_6j1PxJDO4KE7W56QEtYvmt-r_DzbeHDPX0$>).

In VMD, I am inserting the following commands (please see image), but I am only getting carbon atoms making up the solvent box, and not full ethanol molecules (please see image attached); I have tried key selections "name C1, name O1, name HO1" and I still get this.

Can anyone please assist me?

Thank you so much in advance!