Skip to content

Quick Fix to Bondi radii #25

@dunphysd

Description

@dunphysd

Describe the bug
Current buried volume calculations include H-atoms by default. However, the Bondi radii for H-atoms (1.20 A) is not the radii used on popular buried volume tools such as SambVca. Instead, a radius of 1.10 A is used as "carbon bonded to hydrogen". Given the widespread use of SambVca as a tool for buried volume calculations, this discrepancy may cause confusion when different Vbur volumes are produced with SEQCROW.

On a separate note, I am not aware of any way to turn off H-atom inclusion, this would also be a good feature to add given they are frequently omitted in literature reported buried volumes.

To Reproduce
Steps to reproduce the behavior:

  1. Calculate buried volume in SEQCROW as you would normally, making sure there is some H-atom contribution to the buried volume
  2. Calculate buried volume using SambVca on the same .xyz file, making sure to include H-atoms in the Vbur calculation
  3. Compare values. SEQCROW will consistently produce a larger Vbur value due to the larger H-atom radius
    Expected behavior
    The buried volume produced by SEQCROW should match that of SambVca

Metadata

Metadata

Assignees

Labels

bugSomething isn't workingdocumentationImprovements or additions to documentationduplicateThis issue or pull request already existsenhancementNew feature or requestquestionFurther information is requestedwontfixThis will not be worked on

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions