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

unique conformers relative energies print on the screen #3

Closed
lizhiq16 opened this issue May 19, 2020 · 4 comments
Closed

unique conformers relative energies print on the screen #3

lizhiq16 opened this issue May 19, 2020 · 4 comments

Comments

@lizhiq16
Copy link

In the previous version of crest-2.9, the relative energy of each “unique conformers” will be output on the screen at the end of calculation. Now, this output is cancelled, and instead a "crest.energies" file is output. This function is what I really want, but I still want this energy to be displayed on the screen directly sometimes. I hope you can make adjustments.

Two days ago I asked some unrelated questions under another post, I apologize for that. But I still want to ask:
what's the function of automatic bond length constraint (-cbonds) in crest-2.10?

Best wishes!

@fabothch
Copy link

To answer the second part: the -cbonds option allows the user to constrain all bond lengths automatically by providing a force constant. This option has been implemented for being used together with the GFN-FF force field e.g. for metal organic complexes. This is used to avoid dissociation if the forces in the metadynamics simulation are very strong. This option is just for convenience, the user can supply any constraint manually if desired.

The first part has to be answered by @pprcht .

@pprcht
Copy link
Contributor

pprcht commented May 19, 2020

I surely can add some kind of verbose flag for that, but since this printout is (was) right at the end of the program anyways, it would probably be simpler to use cat crest.energies in the command line, just after execution of the program.

@awvwgk
Copy link
Member

awvwgk commented May 19, 2020

I have to admit, that I also liked this particular printout in the standard output, made it easy to quickly get a rough idea of the ensemble.

@lizhiq16
Copy link
Author

I often calculate a batch of molecules. Generally, I directly save the screen output separately as a file for further check, so I don’t want to save another more file (because there will be too many files). It's OK to keep current output, and frankly speaking, I prefer the current form output when dealing with a single task. I will add "crest.energies" to the end of the file I saved.
Thank you very much for your kindly reply.

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

No branches or pull requests

4 participants