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?
Errror reading charge in a .mol file
Re: Errror reading charge in a .mol file
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?
-
- Posts: 87
- 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
(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.
(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.
-
- Posts: 87
- 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
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.
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.
-
- Posts: 87
- 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
See the following web page for the formatted SO2 .mol test:
https://chem4.cns.uaf.edu/webmo1/so2-formatted-mol.htm
https://chem4.cns.uaf.edu/webmo1/so2-formatted-mol.htm
-
- Posts: 87
- 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
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.
@<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.
-
- Posts: 87
- 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
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 (43.75 KiB) Viewed 43055 times