Errror reading charge in a .mol file

Post Reply
jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Errror reading charge in a .mol file

Post 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?

schmidt
Posts: 82
Joined: Sat May 30, 2020 3:00 pm
Full Name: JR Schmidt
Organization: WebMO, LLC

Re: Errror reading charge in a .mol file

Post 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?

jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Re: Errror reading charge in a .mol file

Post 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.

jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Re: Errror reading charge in a .mol file

Post 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.

jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Re: Errror reading charge in a .mol file

Post by jwk »

See the following web page for the formatted SO2 .mol test:

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

jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Re: Errror reading charge in a .mol file

Post 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.

jwk
Posts: 72
Joined: Sat Jun 20, 2020 3:37 am
Full Name: John Keller
Organization: University of Alaska Fairbanks
Subdiscipline: Chemistry

Re: Errror reading charge in a .mol file

Post 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.
Attachments
SO2-molfile.jpg
SO2-molfile.jpg (43.75 KiB) Viewed 1705 times

Post Reply