code: def xyz2AC_vdW(atoms, xyz):
      # Get mol template
      mol = get_proto_mol(atoms)
      # Set coordinates
      conf = Chem.Conformer(mol.GetNumAtoms())
      for i in range(mol.GetNumAtoms()):
          conf.SetAtomPosition(
              i, (float(xyz[i][0]), float(xyz[i][1]), float(xyz[i][2]))
          )
      mol.AddConformer(conf)
      AC = get_AC(mol, xyz)
  
      return AC
  
imports: from rdkit import Chem
  import numpy
  import networkx
  import collections
  from collections import defaultdict
  import itertools
  import copy
  from scipy.spatial import distance_matrix
label_desc: Use image node_id 3 for calling a global function with example usage: xyz2AC_vdW(atoms, xyz) and returns: AC