Full-colour 3D printing of molecular models and orbitals (wavefunctions).

We have been experimenting with full-colour 3D printing of molecular objects. I thought I might here share some of our observations. Firstly, I list the software used:

  1. Crystal structures as sources of ball&stick models (e.g. the CCDC database).
  2. Gaussian style cube files for sources of wavefunctions.
  3. There are plenty of other sources of course, in fact any file type which is supported by ..
  4. Jmol 13.3 (the app, not the applet) to render the above.
  5. And www.shapeways.com who provide the full colour 3D printing service.

My first attempt was to 3D print a molecular orbital; in fact the one shown in this post. To be printable, a 3D object must be fully connected, in other words not contain any disconnected components. To ensure this is true for a molecular orbital, one has to select a very low isosurface threshold; 0.005 au in this case. The orbital is computed as a Gaussian cube, and converted to a .wrl file as follows

  1. Start the Jmol.app.
  2. In the console, enter e.g. 
    isosurface sign color yellow green cutoff 0.005 "144.cub"
    (where 144.cub is the Gaussian cube file). Then
  3. write model.wrl
    to write out a VRML (.wrl) file. 
  4. Create an account at Shapeways, go to “make your own” and upload the .wrl file. After about five minutes, it should appear under “my models“, looking something like below:
  5. Select “full color sandstone” as the material and see what the quoted price is. For the above it was €3.25. It is actually quite tricky to anticipate what physical size this represents, but upon printing and receiving it by post, it turned out to be about 2 cm in size! I should point out that it takes a day or so for  Shapeways to confirm that the model has been sent to manufacturing (it is perfectly possible that it may be rejected at this stage because it cannot be manufactured for whatever reason).
  6. The thing now is to get the model to scale. To do this, a scale-factor has to be applied to the VRML file, which is easier said than done. Fortunately (thanks Bob Hanson!) this can be applied at the Jmol stage, and after step 2 above in the Jmol console, apply the command:
    set exportScale 8.0
    to scale up by 8-fold, before writing out the .wrl file. You know it has worked when the size is given as cm: 13.228 w x 13.288 d x 2.786 h and the price goes up to €69. I will assume (I have not printed at this scale yet) that it is less than 13 cm across, but something big enough to hand around a group of people. I should mention that we tried a variety of ways of scaling the model using programs such as Blender, but they all proved unreliable. You can of course also set the colours in the command above.
  7. It is possible to set the scale larger than 8, but say at 10-fold the materials available to you reduce significantly. Full-colour sandstone is no longer available, and you can only print in a single colour (the materials by the way include stainless steel, for €5700.71). 
  8. The resulting model is made from gypsum with a polymer binder. I have not had the courage yet to test how robust it is by breaking one of the models; I suspect it might be fairly fragile. Time will tell.
  9. I also tried a ball&stick. The molecule selected has been described in another post and is saved by converting the .cif file to .pdb (using in my case the program Mercury). This helps normalise the bond types, producing:
  10. In this case, converting the crystal coordinates via Jmol to .wrl required a scaling factor of 10.0 in this instance to produce a model costed as €24.53. It is very much a matter of adjusting the exportscale against the resulting printing price to get the desired result (by trial and error, but perhaps there is a more systematic way of doing this?). One of my colleagues (Paul) has had this sort of thing printed to about 25cm in size for about this price, and the model seems reasonably robust to physical handling. It may also benefit by fine-tuning of the bond radius; thicker would clearly be stronger.

If anyone reading this post has their own experience of 3D colour printing of chemical models, do please post comments here. And I dare say that in a few years time, students will simply press the “3D print” button on the tablet they are using to view lecture notes to get a copy. Mind you, I am somewhat ambivalent about such a process, having spent the last twenty years trying to discourage students from using the “2D print” button on their computer. Will I eventually come to adopt the same attitude to 3D print (if you click on the MO image above, you will get a virtual 3D model instead of a physical one). 


The width of a bond can be set using the Jmol command: set bondRadiusMilliAngstroms 300


Acknowledgments

This post has been cross-posted in PDF format at Authorea.

Henry Rzepa

Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.

View Comments

  • Afternoon all,

    Just wanted to contribute a little bit of info about our institution, as we work heavily with molecular models using 3D printing technology and might be a useful resource for those reading this article.

    At the MSOE Center for BioMolecular Modeling (http://cbm.msoe.edu/) we use Jmol, Pymol, VMD, Chimera, and other molecular visualization programs to create full-color 3D printed physical models of protein and molecular structures. These materials go directly into high school and undergraduate classrooms for educational purposes - including through our grant-funded "lending library", where educators can borrow our models free of charge.

    We have also created physical models using X-ray and MRI data, including hips, brains and other macroscopic objects.

    Happy modeling everyone!

    -Mark Hoelzer
    MSOE Center for BioMolecular Modeling

Recent Posts

Detecting anomeric effects in tetrahedral boron bearing four oxygen substituents.

In an earlier post, I discussed a phenomenon known as the "anomeric effect" exhibited by…

1 week ago

Internet Archeology: reviving a 2001 article published in the Internet Journal of Chemistry.

In the mid to late 1990s as the Web developed, it was becoming more obvious…

2 months ago

Detecting anomeric effects in tetrahedral carbon bearing four oxygen substituents.

I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…

2 months ago

Data Citation – a snapshot of the chemical landscape.

The recent release of the DataCite Data Citation corpus, which has the stated aim of…

2 months ago

Mechanistic templates computed for the Grubbs alkene-metathesis reaction.

Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…

3 months ago

3D Molecular model visualisation: 3 Million atoms +

In the late 1980s, as I recollected here the equipment needed for real time molecular…

3 months ago