type to search

How to compute the principal moments of inertia for a set of molecules

Asked by , Edited by Alfons Nonell-Canals

I have different sets of molecules and I would like to compute the principal moment o inertia for the molecules of each set. I am looking for a way to do it easy and fast.

Usually, to manage molecules and compute different descriptors I program my own scripts with python and pybel but now, I have no idea how to do it.

or Cancel

1 answer


matteo floris [ Editor ]

Hi, there is a CDK implementation of principal moment of inertia descriptor:


This is a piece of code (any 3D input file):

MomentOfInertiaDescriptor descriptor = new MomentOfInertiaDescriptor();

IMolecule m = null;

IteratingMDLReader reader = new IteratingMDLReader(

new FileReader(new File( "MMsInc_STRUCTURES.sdf)), 

NoNotificationChemObjectBuilder.getInstance() );

while (reader.hasNext()) {

m = (IMolecule);

DoubleArrayResult retval = (DoubleArrayResult) descriptor.calculate(m).getValue();

System.out.println( m.getProperty("cdk:Title") + " " + retval);    


NN comments

Matteo, a short code example would have been nice…

matteo floris
-bla-ics Ok Egon, but I din’t find how to put html tags (pre or blockquote)….

or Cancel

Your answer

You need to join Blue Obelisk eXchange to complete this action, click here to do so.