Feature Requests item #2791729, was opened at 2009-05-14 09:14
Message generated for change (Comment added) made by ldm001
You can respond by visiting: 
https://sourceforge.net/tracker/?func=detail&atid=379136&aid=2791729&group_id=23629

Please note that this message will contain a full copy of the comment thread,
including the initial issue submission, for this request,
not just the latest update.
Category: New IO Format
Group: None
Status: Open
Priority: 5
Private: No
Submitted By: Laurence Marks (ldm001)
Assigned to: Nobody/Anonymous (nobody)
Summary: Wien2k + vibrations

Initial Comment:
Can you please provide a reader for Wien2k DFT files www.wien2k.at (I can 
provide format info etc). Also, can you provide information on how to input 
vibrational modes; I am writing an interface for these so can produce a file in 
the "appropriate" format, if someone tells me what!

----------------------------------------------------------------------

>Comment By: Laurence Marks (ldm001)
Date: 2009-05-14 16:26

Message:
Web source:
http://www.wien2k.at/reg_user/textbooks/usersguide.pdf

I've attached Ti2O3.struct converted to a standard cif

Hunting a bit more, I think Angstroms have vanished from all the cases, so
please assume that everything is in BOHR (even if that key is not there)
and just ignore the units.

I'll put some additional cases in shortly (alas, it deletes the message if
not sent).




----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 16:14

Message:
Oh, just saw that latest comment ---

ANG like that, in caps? 

unit=ANG

or 

unit=ang

right now I see

unit=bohr

not 

unit=BOHR

OK, I'll try fixing the coordinate values for R. 

Some examples of the F, B, I, C$$ cells would be helpful. So you mean the
symmetry operations are not complete? Really?


----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 16:11

Message:
I could use a bit more information 

-- web source for that PDF guide? 
-- check of the Ti2O3 structure.

----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 15:09

Message:
OK, check out http://chemapps.stolaf.edu/jmol/docs/examples-11/new.htm

http://chemapps.stolaf.edu/jmol/docs/examples-11/Jmol-11.jar contains the
files. I think there is a problem with Ti2O3.struct -- it doesn't say BOHR,
but it looks like Bohr, and I'm not convinced the rhombahedral cell is
loading properly. I might need some help with that. But please check all of
them to see if they are loading properly. 

Bob Hanson


----------------------------------------------------------------------

Comment By: Laurence Marks (ldm001)
Date: 2009-05-14 13:27

Message:
Nobody was ldm001 (different computer, forgot to login).

----------------------------------------------------------------------

Comment By: Nobody/Anonymous (nobody)
Date: 2009-05-14 13:26

Message:
Yes, units=Bohr means atomic units or Bohr's, conversion a0=0.529177E0.

----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 13:11

Message:
When it says units=Bohr, that does mean this unit cell is in Bohr?

----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 13:10

Message:
Looks reasonable. The files you have provided should provide plenty of
examples. Jmol should be able to read the unit cell and symmetries and
derive all additional atoms from just the first atom in the set. This will
be a good test, of course, since the additional atom positions are also
given. I'll give it a try.



----------------------------------------------------------------------

Comment By: Laurence Marks (ldm001)
Date: 2009-05-14 10:23

Message:
I've attached some examples, as well as the relevant part of the
documentation. Fairly self evident and the Space Group numbers/names are
normally correct although it is more robust to use the symmetry.

Are all the atoms needed (as against the symmetry reduced set)? It is
simpler to use the reduced set although I can also do the full set. Also,
does Jmol understand about "F", "B", "I" cells etc to know that these
operations are implied or will it need more information? (For assorted
reasons A, B, or C centered cells are coded differently in Wien2k.)

----------------------------------------------------------------------

Comment By: Jonathan Gutow (gutow)
Date: 2009-05-14 10:15

Message:
As there is money being collected for distribution of Wien2k, I suggest
that the Wien2k project should put some programmer time into developing the
file reader if they want Jmol to read their files.

Jonathan

----------------------------------------------------------------------

Comment By: Bob Hanson (hansonr)
Date: 2009-05-14 09:47

Message:
please attach a few file examples. Vibrational modes are read along with
file data, so if your file format included vibrational mode information
then that could be automatically incorporated. If your vibrational
information must be separate, then I recommend the extended XYZ format.
This is simply:

# of atoms
comment line
[atom line]
[atom line]
[atom line]
[atom line]
...


where the atom line is simply:

element x y z vibx viby vibz

For example:

3
CoolMolecules Vibration Generator
O 0 0 0 0 0.005799922848714869 0
H 0.7698177590299584 -0.5480927092035551 0 0.06516975737819754
-0.04639938278971895 0
H -0.7698177590299585 -0.548092709203555 0 -0.06516975737819754
-0.04639938278971894 0


For multiple modes, you just string as many of those together as you
want:


3
CoolMolecules Vibration Generator mode 1
O 0 0 0 0 0.005799922848714869 0
H 0.7698177590299584 -0.5480927092035551 0 0.06516975737819754
-0.04639938278971895 0
H -0.7698177590299585 -0.548092709203555 0 -0.06516975737819754
-0.04639938278971894 0
3
CoolMolecules Vibration Generator mode 2
O 0 0 0 0 0.005799922848714869 0
H 0.7698177590299584 -0.5480927092035551 0 0.06516975737819754
-0.04639938278971895 0
H -0.7698177590299585 -0.548092709203555 0 -0.06516975737819754
-0.04639938278971894 0
3
CoolMolecules Vibration Generator mode 3
O 0 0 0 0 0.005799922848714869 0
H 0.7698177590299584 -0.5480927092035551 0 0.06516975737819754
-0.04639938278971895 0
H -0.7698177590299585 -0.548092709203555 0 -0.06516975737819754
-0.04639938278971894 0


And if you want to include charge, you can do that between z and vibx:

element x y z charge vibx viby vibz



----------------------------------------------------------------------

You can respond by visiting: 
https://sourceforge.net/tracker/?func=detail&atid=379136&aid=2791729&group_id=23629

------------------------------------------------------------------------------
Crystal Reports - New Free Runtime and 30 Day Trial
Check out the new simplified licensing option that enables 
unlimited royalty-free distribution of the report engine 
for externally facing server and web deployment. 
http://p.sf.net/sfu/businessobjects
_______________________________________________
Jmol-developers mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/jmol-developers

Reply via email to