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

Not identifing H_abstraction family #606

Open
MichalKesl opened this issue Mar 6, 2023 · 4 comments
Open

Not identifing H_abstraction family #606

MichalKesl opened this issue Mar 6, 2023 · 4 comments
Assignees

Comments

@MichalKesl
Copy link
Contributor

Describe the bug
I tried to calculate the rate of the H_abstraction reaction : H2N2_116 + ammonia_5 <=> NH2_3 + H3N2_4, without including TS guess. I got this reaction from RMG output that identified it as H_abstraction. While running ARC the reaction wasn't identified as H_abstraction as expected.

How to reproduce
Here is some more examples that ARC didn't identified to be H_abstraction but they are.


job_types:
  rotors: false

allow_nonisomorphic_2d: True  

compute_thermo: false
freq_scale_factor: 1

species:
  - label: 'H4N2_1'
    smiles: 'NN'
  - label: 'H2_2'
    smiles: '[H][H]'
  - label: 'NH2_3'
    smiles: '[NH2]'
  - label: 'H3N2_4'
    smiles: '[NH]N'
  - label: 'ammonia_5'
    smiles: 'N'  
  - label: 'HN2_6'
    smiles: '[N]=N' 
  - label: 'N2H2_19'
    smiles: 'N=N'
  - label: 'H_8'
    smiles: '[H]'    
  - label: 'H2N2_116'
    smiles: '[NH][NH]'
  - label: 'H2N2_113'
    smiles: '[N]N'

reactions:
  - label: H2N2_116 + ammonia_5 <=> NH2_3 + H3N2_4
  - label: H2N2_116 + H4N2_1 <=> H3N2_4 + H3N2_4
  - label: N2H2_19 + H2N2_116 <=> HN2_6 + H3N2_4
  - label: H_8 + H3N2_4 <=> H2_2 + H2N2_116

ARC.log output:

Setting multiplicity of reaction NH2_3 + H3N2_4 <=> H2N2_116 + ammonia_5 to 3
Considering reaction: H2N2_116 + ammonia_5 <=> NH2_3 + H3N2_4



Setting multiplicity of reaction H3N2_4 + H3N2_4 <=> H2N2_116 + H4N2_1 to 3
Considering reaction: H2N2_116 + H4N2_1 <=> H3N2_4 + H3N2_4



Setting multiplicity of reaction HN2_6 + H3N2_4 <=> N2H2_19 + H2N2_116 to 3
Considering reaction: N2H2_19 + H2N2_116 <=> HN2_6 + H3N2_4



Setting multiplicity of reaction H2_2 + H2N2_116 <=> H_8 + H3N2_4 to 3
Considering reaction: H_8 + H3N2_4 <=> H2_2 + H2N2_116

Here is a graphic image of the species of the first reaction:
image

@alongd
Copy link
Member

alongd commented Mar 12, 2023

@MichalKesl, if you search the RMG website for this reaction, it doesn't find it either (which implies it thinks it's not H_Abstraction)

@alongd
Copy link
Member

alongd commented Mar 12, 2023

But in the reverse direction it does find it... Strange.

@kfir4444
Copy link
Collaborator

@Laxzal pointed out that these reactions are marked forbidden in the RMG-database, maybe that's why it can't find them?

@MichalKesl
Copy link
Contributor Author

The reactions are marked forbidden in the RMG-database

@alongd alongd reopened this Jun 11, 2023
@alongd alongd self-assigned this Dec 4, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants