|
From: Rick M. <rm...@sa...> - 2006-07-12 20:32:23
|
cclib looks cool.
Should certainly be possible to make a more robust interface.
Currently the HF and DFT codes have the following interface:
>>> h2 = Molecule('H2',atomlist = [(1,(0,0,-0.7)),(1,(0,0,0.7))])
>>> energy, orbital_energy, orbitals = dft(h2)
If you also have the basis set, you can compute both orbital
amplitudes and densities at points in real space.
Let me know if I can help in any way. I've written code along these
lines before, and can probably dig up something close to what you need.
Rick
On Jul 12, 2006, at 5:19 AM, Noel O'Boyle wrote:
> Dear Rick,
>
> I am one of the developers of cclib (http://cclib.sf.net) which
> provides
> an interface to the results of comp chem calcs from various
> proprietary
> comp chem packages. This is currently GPLed, but we may be moving to a
> more liberal license.
>
> A further goal of cclib is to provide some algorithms that allow
> analysis of the results of comp chem calcs in a platform and package
> independent way. Typical examples of some useful algorithms that
> can be
> carried out after calculations are those that simply involve
> partitioning of the electron density, e.g. Hirshfeld charge analysis,
> AIM, NBO and so on.
>
> I have already coded a small bridge to PyQuante as part of cclib. I am
> interested in expanding this to allow 'recreation' and manipulation of
> the electron density given the basis set details extracted from log
> files. Could you give me any idea whether I will be able to use
> PyQuante
> to accomplish some of these goals?
|