-
Notifications
You must be signed in to change notification settings - Fork 141
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
Update atom enthalpies for H, C, N, O, F, Cl, Br, I with data from ATcT #416
base: main
Are you sure you want to change the base?
Conversation
e7ac70e
to
5719e06
Compare
5719e06
to
179e6c2
Compare
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.
Bill suggested we need to make sure we have the best possible numbers. atom_hf numbers are correct. The problem is atom_thermal. Give me couple of days to make sure we have the correct numbers.
'Li': 1.1, 'Be': 0.46, 'B': 0.29, 'C': 0.25, 'N': 1.04, 'O': 1.04, 'F': 1.05, 'Ne': 1.481, | ||
'Na': 1.54, 'Mg': 1.19, 'Al': 1.08, 'Si': 0.76, 'P': 1.28, 'S': 1.05, 'Cl': 1.1, | ||
'K': 1.481, 'Ca': 1.481, 'Ti': 1.802, 'Cu': 1.481, 'Zn': 1.481, 'Ge': 1.768, 'Br': 1.481, | ||
atom_thermal = {'H': 1.012, 'He': 1.481, |
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.
There is a problem with how H( 0 k) and H(298K) presented, each uses their own standard state. So it Looks like I have to email ATCT and request thermal correction relative to H(298 K)
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.
Aren't atom_thermal just ideal gas corrections? They shouldn't change, right?
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 think they depend on Cp (int Cp=0 to Cp 298)
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.
Let's make sure.
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.
Yes, they're the enthalpy increment of the elements (or 1/2 for diatomic reference states). Translational and rotational contributions are constant. If the vibrational frequencies of diatomics have changed in recent times, then the vibrational contribution could have changed, and if the excited state energy levels have changed in recent times, then the electronic contribution could have changed. My assumption would be that the frequencies and energy levels of diatomics are very well known, so these numbers should be accurate. I'll try to find recent frequencies and energy levels to calculate some of these numbers and see if my calculations agree.
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.
We do 2 for loops in [RMG-Py/arkane/encorr/corr.py](https://github.com/ReactionMechanismGenerator/RMG-Py/blob/74cee80a8c688b04dba265b2bb0b2a10724bee41/arkane/encorr/corr.py#L140) line 120 to 140. I Don't see a problem replacing with one for loop and use the J as units.
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.
Yes, it'll be exactly the same combining them into one number. What I'm trying to say is that it's better having them separate in the database because it makes for easier comparison with literature. For example, ATcT only lists atom_hf so a direct comparison to ATcT wouldn't be easy if we combine the numbers in the database.
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 agree with your comment that it would be easy to compare but the way we use/code we would never able to use that data. So why not combine?
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.
The way that Arkane separates all the terms (and then adds them back together) is the standard way to compute thermochemistry in literature, so I think we should keep it the way it is.
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.
By the way, I calculated the thermal corrections for all the gas phase diatomics, and my calculations agree with the values in the database.
179e6c2
to
b8fe881
Compare
b8fe881
to
795d984
Compare
795d984
to
3dd4b92
Compare
I'm going to get this PR going. I used updated data from ATcT. Some things I noticed:
@dranasinghe, I assigned you to this PR. None of the existing BACs have been corrected yet based on the updated data.
Let's discuss.
RMG-Py PR: ReactionMechanismGenerator/RMG-Py#1955