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

hess and ohess don't support diatomic molecules? #1052

Open
longdragoon opened this issue Jun 16, 2024 · 2 comments · May be fixed by #1077
Open

hess and ohess don't support diatomic molecules? #1052

longdragoon opened this issue Jun 16, 2024 · 2 comments · May be fixed by #1077
Labels
bug Something isn't working

Comments

@longdragoon
Copy link

longdragoon commented Jun 16, 2024

Describe the bug
Hi I'm constructing a small molecule library including the H2 and HCl. The hess (or ohess, as well as the thermo submodule) stops working for these diatomic molecules. xtb can output the hessian and the vibspectrum file correctly but does not give the thermodynamic correction.

To Reproduce
xtb installed with conda from conda-forge
run with "xtb --gfn 1 xtbopt.xyz --hess --verbose >> xtbopt.out 2>> xtbopt.err"
h2.zip
stdout stopped after detecting this molecule has an Dinfh symmetry. stderr does not have any useful information.
all input and output files provided

Expected behaviour
Please support the thermodynamic corrections for linear molecules

@longdragoon longdragoon added the unconfirmed This report has not yet been confirmed by the developers label Jun 16, 2024
@haneug haneug added bug Something isn't working and removed unconfirmed This report has not yet been confirmed by the developers labels Jun 20, 2024
@haneug
Copy link
Member

haneug commented Jun 20, 2024

I could not reproduce this with the downloadable binary but only with the conda version. I submitted a PR that should hopefully resolve this issue but I really do not know much about the conda build.

@longdragoon
Copy link
Author

I could not reproduce this with the downloadable binary but only with the conda version. I submitted a PR that should hopefully resolve this issue but I really do not know much about the conda build.

Yes That's what I noticed, too.
Additionally, I think the two variables magic4 and magic5 (also conv4 and conv5) for calculating the translational entropy (s_tr) is not very correct in xtb/src/thermo.f90

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants