-
Notifications
You must be signed in to change notification settings - Fork 110
MOF workflow #1163
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
base: main
Are you sure you want to change the base?
MOF workflow #1163
Conversation
src/atomate2/common/jobs/mof.py
Outdated
This class wraps the zeo++ executable to calculate pore properties | ||
using given sorbate species. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Could you add the properties that one gets?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I just updated the information. There are many information you can get so one need to have a look at the zeo++ website
src/atomate2/common/jobs/mof.py
Outdated
|
||
output["is_mof"] = False | ||
if all( | ||
k in output["N2"] |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
why only N2?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I just made it more general, now you can input a dictionary
Happy to review once this is ready. |
You may wish to consider using the conda-installable fork of zeo++ that has been patched by EPFL: https://github.com/lsmo-epfl/zeopp-lsmo |
src/atomate2/common/jobs/mof.py
Outdated
output["is_mof"] = ( | ||
output["N2"]["PLD"] > 2.5 | ||
and output["N2"]["POAV_Volume_fraction"] > 0.3 | ||
and output["N2"]["POAV_A^3"] > output["N2"]["PONAV_A^3"] | ||
and output["N2"]["POAV_Volume_fraction"] | ||
> output["N2"]["PONAV_Volume_fraction"] | ||
) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I would suggest calling this is_porous
and not is_mof
. What if you want to run this on a zeolite?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I just made it more general, now you can input a dictionary. It is not necessary is_mof anymore!
Summary
Include a summary of major changes in bullet points:
Additional dependencies introduced (if any)
TODO (if any)
Add the full workflow for MOF optimization, filtering step, combination with MLFF (more broadly ase calculator), widom insertion Monte Carlo Simulation to compute Henry coefficients, averaged interaction energy, and heat of adsorption with MLFF. Create Documentation.