Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

crest xyz energies #9

Closed
joegair opened this issue Mar 16, 2023 · 4 comments
Closed

crest xyz energies #9

joegair opened this issue Mar 16, 2023 · 4 comments
Assignees
Labels
enhancement New feature or request

Comments

@joegair
Copy link

joegair commented Mar 16, 2023

It would be really useful to have a plot of energy vs conformer number pop up when a crest_conformers.xyz file is opened. The same structure as the plot that opens for geometry optimizations would be ideal (plotting energy vs conformer number).

The output would look something like this:

Conformer Number vs Energy

It would be useful to plot energy relative to the lowest energy structure rather than absolute energy.

Just like the plot for geometry optimizations, it would be great to have a second tab that just gives the conformer number and relative energy. I like to use this tab for scrolling through lists of geometries using the down arrow rather than clicking on each point in the plot.

Below is an example crest_conformers.xyz file. By default it contains geometries and energies, sorted from lowest energy structure to highest energy structure with energies in Hartree.

crest_conformers.txt

@joegair joegair added bug Something isn't working documentation Improvements or additions to documentation duplicate This issue or pull request already exists enhancement New feature or request question Further information is requested wontfix This will not be worked on labels Mar 16, 2023
@ajs99778
Copy link
Contributor

I've started working on this. I'd like to make it more general than just crest files, as I've seen other XYZ files with data in the comment line. Now when the plot opens, it shows any floating point number if finds in the comment lines. You'll have the choice to set the units of whatever the values are to Hartrees or kcal/mol. There's images of the plots below.

I think for plotting relative energies, I'll have it convert to kcal/mol. With Hartrees, the y-axis range looks jank (i.e. about the same as it currently looks). For what it's worth, when you hover over a point, it shows the energy in kcal/mol relative to the first point, which should be the lowest for crest.

If you want to test it out, you can install SEQCROW from GitHub. There's instructions for that here: https://github.com/QChASM/SEQCROW/tree/dev#installation

original
image

after choosing Hartrees
image

@ajs99778 ajs99778 removed bug Something isn't working documentation Improvements or additions to documentation duplicate This issue or pull request already exists question Further information is requested wontfix This will not be worked on labels Mar 17, 2023
@joegair
Copy link
Author

joegair commented Mar 17, 2023

This is fantastic! Generalizing to all xyz files is a great idea. Thank you!

@joegair
Copy link
Author

joegair commented Mar 21, 2023

Thanks so much for sharing this. I finally got around to trying it and it works exactly as I was hoping.

I think there is a minor typo in the installation instructions. You may want to change toolshed install /path/to/SEQCROW-1.6.2-none-any.whl, to toolshed install /path/to/SEQCROW-1.6.2-py3-none-any.whl

@ajs99778
Copy link
Contributor

Thanks for pointing out the typo, it's fixed now. I've been a bit busy this week, but I just put the update with this on the ChimeraX tool shed. I changed the unit buttons to a combobox so more options can fit.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants