From: SourceForge.net <no...@so...> - 2010-03-25 14:01:44
|
Feature Requests item #2975200, was opened at 2010-03-23 07:44 Message generated for change (Comment added) made by hansonr You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=379136&aid=2975200&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: Granger Hermitage (ghermitage) Assigned to: Nobody/Anonymous (nobody) Summary: Accept NBO PLOT format Initial Comment: Actually, there is a neater way to visualize Natural Bond Orbitals than I asked for in feature request 2974999. When PLOT (it has no subparameters) is specified in the NBO parameter list, all files needed for visualization are produced. The NBO parameter list ($NBO ... $END) appears in the electronic structure system input if NBO is embedded in the ESS, and appears in the .47 input file to GenNBO if NBO is being run standalone. $NBO AONBO=P $END is not necessary, for $NBO PLOT $END will output what is needed for visualization of all natural orbital types not just NBO in AO basis. Indeed, PNBO (prenormalized NBO) may be more useful for visualization of overlaps than NBO. PLOT causes the generation of these file types: molname.31 AO molname.32 PNAO molname.33 NAO molname.34 PNHO molname.35 NHO molname.36 PNBO molname.37 NBO molname.38 PNLMO molname.39 NLMO molname.40 MO molname.41 AO density matrix molname.46 Basis label file Regardless of what ESS NBO 5 is imbedded in or what ESS GenNBO 5 derives its input from the format of these files is unchanged, so this is a completely uniform input for natural orbital visualization. Also, these inputs make a visualization tool fully functional for all natural orbital types. Also, all these types of natural orbitals (including PNBO) are available in a machine oriented format rather than having to be picked out of various places in a somewhat more human oriented format of an ESS log or output file or nbo output file (.nbo). Also, the formats .31 .. .41, .46 could be expected to be more stable across NBO updates than that of log files that might see modification in the interests of human readability. The .31 file must be present for visualization as it contains such good things as atomic number. The different types of natural orbitals are useful depending on what is being investigated. The thing to do would be to supply a directoryname/molname and let all the .31-.41,.46 be loaded with a single command, then specify what basis you wanted to be working in at that moment. I will attach a zipped sample file of each file type, though the .33 file might have to be compressed with another method, say zipx, to get under the upload limit. I can imagine Jmol quite rapidly becoming the best tool to visualize natural orbitals. ---------------------------------------------------------------------- >Comment By: Bob Hanson (hansonr) Date: 2010-03-25 09:01 Message: for each new file you add, also add a full Gaussian/NBO output please. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-24 08:24 Message: good. That's what I needed. I'll experiment with appending to Gaussian output, since I have both the log file and the nbo file now. If you can generate more of these with the range of calculation I listed earlier, that would be great. A Gaussian file with associated nbo and .nn files would be perfect. I suggest doing a set of these, zipping them up into a single ZIP file and sending them to me via a service such as YouSendIt. Will get back to this later today. Bob ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-24 07:59 Message: Uploaded 14_a_b3lyp.nbo.zip which has AONBO=P in the NBO parameter list. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-24 07:39 Message: By the way, in case you hadn't noticed, it is really nice that we can drag the .XX files into Jmol and have instant access to the full range of orbitals there. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-24 07:37 Message: I think you want AONBO=P in the NBO parameter list. I am not really expecting much different than appears in the .31 file, but I will run it and upload the output. That sounds right. Key is that it is there so that I can compare how Jmol performs reading the Gaussian file (which I know is correct) with reading the GenNBO output. I could be putting coefficients or exponents in an incorrect order since I'm just guessing on the format of the .31 file. (If you have code for reading that, I'd love to see it!) ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-24 07:27 Message: I think you want AONBO=P in the NBO parameter list. I am not really expecting much different than appears in the .31 file, but I will run it and upload the output. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-24 07:22 Message: sorry, still no go. Gaussian 03 log with GFPrint on uploaded 14_a_gfprint.log.zip (for 14_a_b3lyp) does not contain information about NBOs. Expecting a line like: NBOs in the AO basis: AO 1 2 3 4 5 6 7 8 ---------- ------- ------- ------- ------- ------- ------- ------- ------- 1. C 1 (s) -0.0399 -0.0455 -0.0439 -0.0439 0.0062 -0.0060 0.0000 -0.0013 ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-24 07:22 Message: I am using Java 6 Update 18 here. The orbitals shown for .36, .37 files are simpler and more regular than the Natural Boind Orbitals I saw with GaussView 4.1.2 back when that worked here. Don't know if it a problem, might be the default contours are different and a particular level of contour shows the complexity. Could also be the extremely limited basis set I used to get a small and hence uploadable example. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-24 06:57 Message: working ok for me. I'll upload another version. ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 23:34 Message: Error reading file at line 525: NBO 198 String index out of range: 2009 for file D:/temp/14_a_b3lyp/14_a_b3lyp.35 type NBO The PLOT output files are those already uploaded to the site. The error doesn't occur with the .36 file. The .35 file is reported as being type NBO, but should be type NHO. ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 22:52 Message: Gaussian 03 log with GFPrint on uploaded 14_a_gfprint.log.zip (for 14_a_b3lyp) ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 22:20 Message: Unexpected end of file after line 2559: C 7- H89* C 7- H90* C 7- H78* for file D:/temp/1_a_mp2/1_a_mp2.36 type NBO -- load >> "D:temp\1_a_mp2\1_a_mp2.36" Same error for the .35 file. Uploading 1_a_mp2.31, 1_a_mp2.36, 1_a_mp2.46 now. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 21:52 Message: http://chemapps.stolaf.edu/jmol/docs/examples-11/Jmol-11.zip ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 21:51 Message: The standard default is just the number -- that's what we have here. I'll stick with that. I need the gfprint option in Gaussian so I can compare the bases. So I'm pretty sure this is working. Final tests, with Gaussian to compare, would be: a calculation with spherical D orbitals. a calculation involving chlorine. a calculation with cartesian D orbitals a calculation with F orbitals. a calculation with S and P but not SP atomic orbitals an open-shell calculation of a radical a calculation having both alpha and beta orbitals (and what you want to do about that) That should just about do it. You can run tests yourself using the version found in http://chemapps.stolaf.edu/jmol/docs/examples-11 Only thing missing there is I just got .31 files showing the AOs. Bob ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 19:14 Message: (a) I would just give the MOs a number as they are read in and use that number as a label. (b) Yes, give the PXXX the same names as XXX. For a 14_a.log basis functions listing, do you want an option switched on in the ESS parameter list or the NBO parameter list? ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 18:56 Message: Question: The .46 file has names for only AO, NAO, NHO, and NBO orbitals. Then (a) there are no names for MO orbitals, and (b) the names of PXXX are the same as XXX? Names do seem to match pretty well with depictions... Bob ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 18:54 Message: 14_a.log needs basis functions listed, please. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 18:52 Message: Jmol 11.9.36_dev can provisionally read .nbo files provided 0.31 file is present, 0.31 file for structure only, and 0.32-0.41 files provided 0.31 and 0.46 are present in the directory. See http://chemapps.stolaf.edu/jmol/docs/examples-11/Jmol-11.zip ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 17:00 Message: 14_a_b3lyp.nbo.zip (GenNBO output) and 14_a.log.zip (Gaussian 03 output) uploaded, 14_a_b3lyp.33.zip sent. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 13:12 Message: OK, Jmol is now reading the .nbo file and .31-.41 files. At least, provisionally. No orbitals yet, because I don't have any idea what the format of those files is. I'll need the .nbo output file for 14a, please. Also, if this comes from an ESS calc, can you please send me the ESS output file as well, showing the AO basis functions used? This will allow me to verify the file reading. Thanks Thanks. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 12:35 Message: OK. We could do that. Where is the description of the .31 file? Bob ---------------------------------------------------------------------- Comment By: Granger Hermitage (ghermitage) Date: 2010-03-23 10:52 Message: It would be okay to load one file in the range .32 to .40 at a time For files .32 to .40 as .yy : load "molname.yy" this load could get the molname.31 and molname.46 behind the scenes. A visualization session most often (perhaps always) compares and contrasts orbitals of the same type e.g. orbitals appearing in just one type in the range .32 .. .40 Files .32 to .40 are likely to be of the same format, so code for one should work for all. A visualization tool isn't aware of the difference in meaning of orbital type. ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 09:24 Message: More thoughts. What do you think of using the OUT file as the file that is indicated (or dragged) to read. Then, if the information such as AO Basis is missing, because this is NOT part of an ESS output, then we have Jmol automatically look for and load the .31 file. And if certain FILTER options are given (which we would agree upon), then it will also read another .NN file. Could be something like load "xxx.nbo" FILTER "PNBO" load "xxx.nbo" FILTER "PNLMO" That could work. The reader would know then to go grab the data in whichever file is appropriate. Bob ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 08:55 Message: I can't decompress the 33 file. (gzip would work) -- you can send that to ha...@st... directly if you want. We can look into this as a long-term project. It would be a pretty big job that would need quite a bit more information to create a totally new reader for these. While I appreciate the fact that it adds capabilities that the ESS programs don't have, I would have to think hard about how to implement this properly considering the complexity of having all those files to worry about. OK, there are two possibilities. We do also read Spartan directories, so the idea of reading a directory or reading a zip file with included files is possible, but it relies on Jmol first reading ALL the files and then parsing the concatenated script. It's really not ideal with these big files. In principle that could be adapted a bit. Jmol does not use file extensions for identifying files, so it's a bit problematic to think about doing this just from a list of file names. For example, you can drag and drop into Jmol, but that doesn't work if you want a specific SET of files. By the way, where are the energies stored? Or doesn't that apply in the case of NBOs? ---------------------------------------------------------------------- Comment By: Bob Hanson (hansonr) Date: 2010-03-23 07:57 Message: ok, I will look at that. Not sure that's exactly "easier"! ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=379136&aid=2975200&group_id=23629 |