[AMBER] Electrostatic calculation

From: angad sharma via AMBER <amber.ambermd.org>
Date: Tue, 14 Feb 2023 19:33:38 +0530

Dear amber users,
I am trying to calculate the non-bonded interactions of Cu2+ ions with DNA
around the 5 angstrom of the DNA. In results i am getting positive
van-der-waal energy, i am not able to figure weather there is some mistakes
in my input file or is it something else.
Following are my input commands:-
---------------------------------------------------------------------------------------------------------------
echo "" > op1-ene_nb.txt

for i in {1..5000}
do
        echo "parm ../../dna-cu.prmtop
trajin ../tr2-5kfr.nc $i $i

reference ../tr2-5kfr.nc $i $i
strip !(:1-24<:5.0) outprefix dna-cu-op1 nobox
energy nonb_el_vdw .OP1,CU out nonbond_frm.txt nonbond elec vdw
go
quit

" > cppt.in
cpptraj -i cppt.in
tail -n 1 nonbond_frm.txt >> op1-ene_nb.txt

done

rm cppt.in nonbond_frm.txt
--------------------------------------------------------------------------------------------------------------------
Some results of this command is:-
                     vdw eletro total
       1 -0.1009 -638.8666 -638.9676
       1 -0.1069 -1101.7432 -1101.8501
       1 -0.1213 -1016.8659 -1016.9872
       1 -0.1194 -1225.8950 -1226.0144
       1 5.8338 -1315.7839 -1309.9501
       1 10.3802 -1231.8396 -1221.4594
       1 14.5967 -1336.6635 -1322.0668
       1 9.9898 -1046.7325 -1036.7427
       1 9.2984 -1408.1770 -1398.8786
       1 5.5465 -1274.6634 -1269.1169
       1 4.6642 -1552.3137 -1547.6495
       1 8.3172 -1353.3109 -1344.9936
       1 4.9239 -1540.5821 -1535.6582







*-------------------------------------------------*
*With Regards*
*Angad Sharma*

*Research Scholar*
*School of Computational & Integrative Sciences*

*Jawaharlal Nehru University*
*New Delhi 110067, (INDIA)*
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Feb 14 2023 - 06:30:05 PST
Custom Search