[AMBER] Triple bond

From: FyD via AMBER <amber.ambermd.org>
Date: Wed, 01 Jul 2026 10:01:15 +0200

Dear All,

I do not understand the energy values reported for '3-point' angles in
triple bonds by sander 2024 (automatic by REDServer Dev./PyRED)|pmemd
2024; done manually on the q4md-fft and MatriCS
https://www.matrics.u-picardie.fr/amber/ clusters):
For HF/6-31G(d) optimized ethyne H-C≡C-H (simple model), there are 2
angles HCC & CCH, and measured_angle_values ~ parm10 eq. value ~ 180
deg.; so the corresponding angle energy values ~ 0 kcal/mol:
sander|pmemd do not find '0'...

Similar 'results' are found for propyne and acetonitrile, and
discrepancies are observed between PyRED and sander for '3-point'
angles with triple bonds...

[xxxx.master P2875]# grep ANGLE dE_sander-vs-PyRED_SinglePointEnergyValues.log
    ANGLE sander vs PyRED single-point energy abs(differences) =
  .2000 [WARNING] i.e. C≡C_Ethyne
    ANGLE sander vs PyRED single-point energy abs(differences) =
  .2600 [WARNING] i.e. CC≡C_Propyne
    ANGLE sander vs PyRED single-point energy abs(differences) =
  .1600 [WARNING] i.e. CC≡N_Acetonitrile
    ANGLE sander vs PyRED single-point energy abs(differences) =
      0 [OK] i.e. C≈C_Toluene

Do I miss something?
Thank you
Best Francois

See data below:

WARNING: a job is kept _only_ 5 days on the public account...

See the home page of job P2875 (login=password='q4mdfft' for the
public account):
https://q4mdfft:q4mdfft.cluster.q4md-forcefieldtools.org/~ucpublic1/ADF2yADF10ADFNDAPwJLt0IXRSorhL/P2875.php
https://q4mdfft:q4mdfft.cluster.q4md-forcefieldtools.org/~ucpublic1/ADF2yADF10ADFNDAPwJLt0IXRSorhL/P2875/dE_sander-vs-PyRED_SinglePointEnergyValues.log

See energy values reported by PyRED:
https://q4mdfft:q4mdfft@cluster.q4md-forcefieldtools.org/~ucpublic1/ADF2yADF10ADFNDAPwJLt0IXRSorhL/P2875/Data-R.E.D.Server/Mol_MM/Energy_mm1-c1.txt
   versus energy values reported by sander|pmemd
https://cluster.q4md-forcefieldtools.org/~ucpublic1/ADF2yADF10ADFNDAPwJLt0IXRSorhL/P2875/Data-R.E.D.Server/Mol_MM/MM_sander_vs_pyred_SinglePoint/Mol_mm1-c1.log
    or
https://cluster.q4md-forcefieldtools.org/~ucpublic1/ADF2yADF10ADFNDAPwJLt0IXRSorhL/P2875/Data-R.E.D.Server/Mol_MM/MM_pmemd_vs_pyred_SinglePoint/Mol_mm1-c1.log
    NSTEP ENERGY RMS GMAX NAME NUMBER
       1 7.6586E+00 1.0429E+01 2.5501E+01 C1 1

  BOND = 0.2542 ANGLE = 0.2000 DIHED = 0.0000
  VDWAALS = 0.0000 EEL = 0.0000 HBOND = 0.0000
  1-4 VDW = -0.0055 1-4 EEL = 7.2098 RESTRAINT = 0.0000

  !!! ANGLE = 0.2000 !!!

[xxxx.master P2875]# cat Project.config
#_______Keyword(s) added in the Project.config file by the user before
the input submission procedure:
MOLECULE1-TITLE = C≡C_Ethyne
MOLECULE2-TITLE = CC≡C_Propyne
MOLECULE3-TITLE = CC≡N_Acetonitrile
MOLECULE4-TITLE = C≈C_Toluene

[xxxx.master P2875]# read_mol2_atom
Data-R.E.D.Server/Mol_MM/Mol_mm1-c1.mol2 (i.e. ethyne)
    1 C1 -0.592723 0.000000 0.000000 CZ 1 EYN -0.2932 0.0000 ****
    2 H1 -1.649675 0.000000 0.000000 HZ 1 EYN 0.2932 0.0000 ****
    3 C2 0.592723 0.000000 0.000000 CZ 1 EYN -0.2932 0.0000 ****
    4 H2 1.649675 0.000000 0.000000 HZ 1 EYN 0.2932 0.0000 ****

[xxxx.master P2875]# cat
Data-R.E.D.Server/Mol_MM/MM_sander_vs_pyred_SinglePoint/Mol_mm1-c1-in.rst7
(i.e. ethyne)
EYN
      4
   -0.5927230 0.0000000 0.0000000 -1.6496750 0.0000000 0.0000000
    0.5927230 0.0000000 0.0000000 1.6496750 0.0000000 0.0000000

    -> the cart.coordinates are identical; the rst7 is generated from
the mol2 using LEaP
    -> parmed reports what should be reported for 3-point angles

[xxxx.master P2875]# source /usr/local/amber24/amber.sh
[xxxx.paster P2875]# parmed MM_sander_vs_pyred_SinglePoint/Mol_mm1-c1-in.parm7
  [...]
ParmEd: a Parameter file Editor
Loaded Amber topology file
Data-R.E.D.Server/Mol_MM/MM_sander_vs_pyred_SinglePoint/Mol_mm1-c1.parm7
Reading input from STDIN...
> printangle
printangle command not recognized
> printangles
              Atom 1 Atom 2 Atom 3 Frc
Cnst Theta eq
       1 C1 ( CZ) 3 C2 ( CZ) 4 H2 ( HZ)
50.0000 180.0001
       2 H1 ( HZ) 1 C1 ( CZ) 3 C2 ( CZ)
50.0000 180.0001



_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jul 01 2026 - 01:30:03 PDT
Custom Search