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

get_EoRT_act calculates HoRT_act #187

Open
googhgoo opened this issue Jul 30, 2021 · 2 comments
Open

get_EoRT_act calculates HoRT_act #187

googhgoo opened this issue Jul 30, 2021 · 2 comments
Labels
bug Something isn't working

Comments

@googhgoo
Copy link
Member

Version of pMuTT
1.2.21

Describe the bug
Hello,
When we use get_EoRT_act to calculate activation energy, it actually calculates get_HoRT to calculate activation energy. This maybe due to the difference in the philosophy. When I use this function, I was more interested in the experimentally obtained value, which is the raw electronic energy terms + the zero point energy terms, but this function calculates enthalpy of activation instead, which is not necessarily technically correct. I think this need a clarification.

Another way to explain this is that, get_EoRT_act is the same as get_HoRT_act. These should not be the same values.

@googhgoo googhgoo added the bug Something isn't working label Jul 30, 2021
@jonlym
Copy link
Member

jonlym commented Aug 4, 2021

Hey Geun. Sorry for the delay. For some reason I stopped getting Github notifications.

I agree that I think it's a difference in philosophy. We wanted to use transition-state theory, which dictates the Ea is related to ΔH_TS and the A is related to ΔS_TS. See the exact formula we used here in the get_EoRT_act docstring. This has the added benefit of working with empirical objects that do not have electronic/vibrational properties.

We could add an option to get specifically the electronic + ZPE terms. Although not ideal, until a new version comes out you could calculate it using:
Ea = reaction.get_EoRT_state(state='ts', include_ZPE=True) - reaction.get_EoRT_state(state='reactants', include_ZPE=True))

@googhgoo
Copy link
Member Author

googhgoo commented Aug 5, 2021

Ea is related to ΔH_TS and the A is related to ΔS_TS.

Is this true? When we get Ea from experiment using Arrhenius plot, aren't we only getting E0+ZPE? When we get slope (-Ea/R), the T is completely out, so Ea obtained using Arrhenius plot is different than ΔH_TS.

So technically, A is related q_TS/q_IS, terms of which exclude E0+ZPE. I get that we are working with transition state theory most of the time, but I think it can be confusing for some people.

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

No branches or pull requests

2 participants