[AMBER] MCPB.py Failing with OPC water?

From: Gustavo Seabra via AMBER <amber.ambermd.org>
Date: Thu, 13 Oct 2022 17:29:47 -0400

Hi,

I'm trying to build a system with an Fe2+ center, and using OPC water. The
manual states that, in AmberTools 22, MCPB.py should be compatible with OPC
water.

However, when trying to build the system with OPC water, I get
the following message:
<--
$ MCPB.py -i MCPB.in -s 1
[...]
Totally there are 105 atoms in the large model.
Totally there are 450 electrons in the large model.
Traceback (most recent call last):
  File "/opt/miniconda3/envs/AmberTools22/bin/MCPB.py", line 651, in
<module>
    gene_model_files(orpdbf, ionids, addres, addbpairs, gname, ff_choice,
  File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mcpb/gene_model_files.py",
line 1909, in gene_model_files
    build_large_model(mol, reslist, lmsresids, lmsresace, lmsresnme,
lmsresgly,
  File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mcpb/gene_model_files.py",
line 1600, in build_large_model
    IonLJParaDict = get_ionljparadict(watermodel)
  File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mol/element.py",
line 325, in get_ionljparadict
    IonLJParaDict.update(diljpara2)
UnboundLocalError: local variable 'diljpara2' referenced before assignment
-->

If I just replace the OPC model by TIP3P in the water_model line, it works
just fine.

--
Gustavo Seabra.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Oct 13 2022 - 15:00:02 PDT
Custom Search