Page 1 of 1

Errror reading charge in a .mol file

Posted: Sun Nov 20, 2022 8:17 pm
by jwk
WebMO v. 22 (demo site) and v. 19 (local) insert +3 charge on S when importing a .mol file.
Try this .mol file, which was exported from WebMO in mol/sdf format:
----------------------

WebMO

3 2 0 0 1 V2000
0.0000 0.6833 0.0000 S 0 1 0 0 0 0 0 0 0 0 0 0
-0.5726 -0.6896 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
1.7356 0.9075 0.0000 O 0 -1 0 0 0 0 0 0 0 0 0 0
1 2 2 0
1 3 1 0
M END
-----------------------
Of course this is easy to correct manually, but I would like to provide for beginning students a .mol file that has idealized geometry and reasonable formal charges - as a starting point for further calculations. Is there another way to do this in WebMO so we get the correct formal charges after importing?

Re: Errror reading charge in a .mol file

Posted: Mon Nov 21, 2022 5:04 pm
by schmidt
I can't import the snippet above, but it looks like SO2. When I export SO2 from WebMO as a MOL and then import, I get no charges, as expected. Can you elaborate?

Re: Errror reading charge in a .mol file

Posted: Mon Nov 21, 2022 5:57 pm
by jwk
(1) Build SO2 with 1 db and 1 sb. Add charges by right-clicking S (+1) and one O (-1). Then File, Export Molecule as mol/sdf format, creating a "molecule.mol" file in the download folder.
(2) Now start over: File, New. File, Import Molecule as mol/sdf, Choose File, find the above molecule.mol file, click OK.
Now display charges on the resulting molecule, showing +3 on S.

Step (1) would be what an instructor would do, wishing to distribute correctly-built-including-charges .mol files to students. With the .mol file (or text thereof) publicly available.
Step (2) would be what the student would do in order to import said molecule with charges.

Re: Errror reading charge in a .mol file

Posted: Mon Nov 21, 2022 6:19 pm
by jwk
More:
After assigning charges, do Cleanup, Comprehensive-Idealized.
The .mol text in the above window was modified by the Forum program so that it cannot be imported. To test my assertion, one would therefore need to rebuild SO2 per the directions given. Note that crucially the .mol file does record the formal charge assigned to each atom in the 6th field.

Re: Errror reading charge in a .mol file

Posted: Mon Nov 21, 2022 7:32 pm
by jwk
See the following web page for the formatted SO2 .mol test:

https://chem4.cns.uaf.edu/webmo1/so2-formatted-mol.htm

Re: Errror reading charge in a .mol file

Posted: Wed Nov 30, 2022 9:46 pm
by jwk
Avogadro 1.2 and 1.97 both put a +3 charge on S when opening the above .mol file. OpenBabel 3.1.1 converts this SO2.mol file into Sybyl-mol2 format as follows:
@<TRIPOS>MOLECULE
*****
3 2 0 0 0
SMALL
GASTEIGER

@<TRIPOS>ATOM
1 S -1.0184 1.3212 -23.9010 S.3 1 UNL1 0.1816
2 O 0.7228 1.4965 -23.9010 O.3 1 UNL1 -0.0423
3 O -1.6294 -0.0351 -23.9010 O.2 1 UNL1 -0.1393
@<TRIPOS>UNITY_ATOM_ATTR
1 1
charge 3
@<TRIPOS>BOND
1 2 1 1
2 1 3 2

That is, atom #1 (S) is given +3 charge.
These results suggest that both WebMO and Avogadro use an openbabel plugin to translate various molecular formats, rather than directly importing the .mol file contents. In the final analysis, perhaps it's best to let students assign formal charges, rather than rely on charges obtained in WebMO imports.

Re: Errror reading charge in a .mol file

Posted: Thu Dec 28, 2023 5:15 am
by jwk
A workaround tricks WebMO into assigning the correct formal charge on S the O=S-O resonance form. In the .mol file, enter 3 in place of 1. See the attached jpg. This also works for .mol files of iminium ions like CH2=NH2+ which have the same problem.