Page 6 of 6 < 1 2 3 4 5 6
Topic Options
#24429 - 05/26/10 03:08 PM Re: generating parameters for pyruvate [Re: RD2015]
Kenno Offline
Forum Member

Registered: 12/19/05
Posts: 1535
Loc: Baltimore, MD
Don't get blinded by the averages on the last line, they are the reason why (at least some) automatic fitting programs perform so badly.
  • The interaction energies with the CH3 group got better, but they're totally irrelevant because they're repulsive.
  • Interaction 1) got substantially worse; I'd try to fix that up.
  • Interaction 5) got substantially better, at the cost of interaction 6), which got a little bit worse. Unfortunately, out of these two interactions, interaction 6) is more important (lower AI energy). Yet, improving interaction 6) should not lead to a worsening of interaction 3). The carboxylate group is a bit messy and I'm not sure you can improve it, but you can try. I'm not a psychic.

When you're done with that, yes, it's perfectly possible that the final unconstrained values are better than the constrained ones.

Top
#24495 - 06/07/10 11:48 AM Re: generating parameters for pyruvate [Re: Kenno]
RD2015 Offline
Forum Member

Registered: 01/29/10
Posts: 85
Sorry I haven't posted my results in a while. I was on break.

I adjusted the charges and replotted the scan.

Code:
PYRU WATER INTERACTION
QM DIPOLE (DEBYE)
MP2/6-31G*: X= 4.6182 Y= 1.9920 Z= 0.6033 TOT= 5.0655
EMPIRICAL DIPOLE: 6.17438 2.37572 0.68391 6.65092
1) O2...HOH, 180. DEGREES
A.I. -10.21 1.94
EMP. -9.98738 -11.5515 1.56408 1.7
ENE DIFF: 0.22262 DIST DIFF: -0.24
2) O2...HOH LINEAR, HOH: 90. DEGREES
A.I. -8.86 1.97
EMP. -10.2175 -11.673 1.45541 1.66
ENE DIFF: -1.3575 DIST DIFF: -0.31
3) O11...HOH, 180. DEGREES
A.I. -13.45 1.87
EMP. -11.6132 -13.1247 1.51148 1.72
ENE DIFF: 1.8368 DIST DIFF: -0.15
4) O11...HOH LINEAR, HOH: 90. DEGREES
A.I. -12.97 1.84
EMP. -13.9382 -16.0165 2.07834 1.68
ENE DIFF: -0.9682 DIST DIFF: -0.16
5) O12...HOH, 0. DEGREES
A.I. -11.03 2.01
EMP. -9.98958 -11.3142 1.32461 1.75
ENE DIFF: 1.04042 DIST DIFF: -0.26
6) O12...HOH LINEAR, HOH: 90. DEGREES
A.I. -12.65 1.86
EMP. -14.0225 -16.1759 2.1534 1.68
ENE DIFF: -1.3725 DIST DIFF: -0.18
7) C3H..OHH 90.0 DEGREES
A.I. 0 4.50
EMP. 2.22844 2.29245 -6.401592E-02 4.49
ENE DIFF: 2.22844 DIST DIFF: -1E-02
8) C3H..OHH 90.0 DEGREES
A.I. 0 4.50
EMP. 2.67073 2.99919 -0.328457 2.8
ENE DIFF: 2.67073 DIST DIFF: -1.7
9) C3H..OHH 90.0 DEGREES
A.I. 0 4.50
EMP. 3.03498 3.12064 -8.56574E-02 4.49
ENE DIFF: 3.03498 DIST DIFF: -1E-02
AVEDIFF, RMSDIFF, AVERAGE ABSOLUTE ERROR
0.815088 1.64993 1.63691


Attachments
CCCH.png

OCCC.png



Top
#24498 - 06/07/10 02:16 PM Re: generating parameters for pyruvate [Re: RD2015]
Kenno Offline
Forum Member

Registered: 12/19/05
Posts: 1535
Loc: Baltimore, MD
Hmm, the charges are not convincingly improved. This seems to be a bit of a tough case; the ranking of interactions 1-6 is significantly different in the QM and the MM. Maybe we should go for a compromise after all, like getting the average difference of 1 and 2 as close to 0 as possible, and same for the average difference of 3,4,5, and 6.

The scans are quite nice.

Top
#34581 - 11/20/14 09:34 PM Re: generating parameters for pyruvate [Re: Kenno]
yxb Offline
Forum Member

Registered: 11/12/14
Posts: 1
Hi Kenno,I'm generating parameters for my molecular following the "CGenFF tutorial" (file:///E:/documents/ug-tutorial/Charmm/cgenff-tutorial/cgenff_tutorial_old/supp/tutorial.html) and your guidance.

I have finished the part "Target Part".When I do the next "Initial charge optimization using the QM geometry",I'm confused with how to get the file water_constr.inp. Is there any generator to get it? Or should I write it manually? If it's the latter case, how can I prepare "the IC table for surface of interest"? Is there any principles to write this part?

Sorry for so many questions,maybe some are silly. frown I'm new for CHARMM and I'm really trapped here. Looking forward for your reply.

Top
#37644 - 08/29/19 05:48 AM Re: generating parameters for pyruvate [Re: RD2015]
thapa Offline
Forum Member

Registered: 08/08/19
Posts: 1
Can you please explain how the dummy atoms are chosen to be placed in the script ?

In this folder :pyru_water_gauss.tgz


Edited by thapa (08/29/19 05:51 AM)

Top
Page 6 of 6 < 1 2 3 4 5 6

Moderator:  alex, chmgr, John Legato