calculate_G4 and cutoff functions

Issue #15 resolved
andrew_peterson repo owner created an issue

I think that in the Behler fingerprint scheme, calculate_G4 is incorrect. According to Artirth:

http://pubs.acs.org/doi/suppl/10.1021/nl5005674/suppl_file/nl5005674_si_001.pdf

each term should be multiplied by the product of three cutoff functions. Our code looks like it is multiplying by the sum of the three cutoff functions.

I noticed a discontinuity when plotting an angular fingerprint term versus bond length (abruptly goes to zero as the bond length exceeds the cutoff threshold), and think this explains it.

Comments (2)

  1. Log in to comment