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

xtbout.json not written for GFN-FF #745

Closed
nielskm opened this issue Jan 3, 2023 · 4 comments
Closed

xtbout.json not written for GFN-FF #745

nielskm opened this issue Jan 3, 2023 · 4 comments
Labels
easy to fix This issue has a low difficulty and little entry-barrier, great to get started with this project enhancement New feature or request

Comments

@nielskm
Copy link

nielskm commented Jan 3, 2023

The xtbout.json file that should be written when passing the --json flag is not generated when running GFN-FF calculations.

To Reproduce
Running xtb <some-molecule.xyz> --gfnff --json produces no JSON output file.

Additional context
xtb version: 6.5.1

@nielskm nielskm added the unconfirmed This report has not yet been confirmed by the developers label Jan 3, 2023
@awvwgk awvwgk added bug Something isn't working and removed unconfirmed This report has not yet been confirmed by the developers labels Jan 4, 2023
@MtoLStoN
Copy link
Member

MtoLStoN commented Jan 12, 2023

The "--json" flag is, at the moment, not used in conjunction with GFN-FF. We could use this flag to write all available information (replacing the need to use the --wrtopo flag with keywords). See #750.
Is that what you had in mind?

@nielskm
Copy link
Author

nielskm commented Jan 16, 2023

@MtoLStoN Thank you for taking the time to look at this. I was simply thinking if you could output a JSON containing total energy, etc. as done for the xTB methods.

@MtoLStoN MtoLStoN added enhancement New feature or request and removed bug Something isn't working labels Jan 17, 2023
@MtoLStoN
Copy link
Member

MtoLStoN commented Jan 17, 2023

I understand; I'll look into this. This should be relatively easy to do, if I find the time. If this is something you want faster, you are more than welcome to submit a patch.

The relevant call can be found here:

xtb/src/prog/main.F90

Lines 860 to 868 in 863dcdc

if (set%pr_json) then
select type(calc)
type is(TxTBCalculator)
call open_file(ich,'xtbout.json','w')
call main_json(ich, &
mol,chk%wfn,calc%basis,res,fres)
call close_file(ich)
end select
endif

And we would need an additional writer routine for the GFN-FF json file.

https://github.com/grimme-lab/xtb/blob/main/src/main/json.f90

This could probably be done with an explicit interface and a second writer routine just for GFN-FF.

@MtoLStoN MtoLStoN added the easy to fix This issue has a low difficulty and little entry-barrier, great to get started with this project label Jan 17, 2023
@MtoLStoN
Copy link
Member

Thanks to @Thomas3R, we have this implemented now!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
easy to fix This issue has a low difficulty and little entry-barrier, great to get started with this project enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

3 participants