From: linefinc (linefinc_at_gmail.com)
Date: Tue Jun 20 2006 - 11:02:28 CDT
TNX TNX TNX
I have use a short simulation with 100 steps and after the usal
simulation and all work!
2006/6/20, Giovanni Bellesia <gbellesia_at_chem.ucsb.edu>:
> Hi there,
> I quickly went through your files.
> As far as I see your molecules are perfectly linear and this may create
> problems in the calculation of the potential energy (bond angle potential).
>
> As a matter of fact, if you put to zero the value of Ktheta in the
> parameters file you will see that your simulation will start correctly.
> In conclusion, a quick solution if you want to use these initial
> coordinates might be to run a short minimization (100 steps or less)
> using Ktheta=0 then you can use the newly-obtained coordinates to start
> your simulation putting back Ktheta to its right value.
>
> Giovanni
>
> > I work as usual on my CO2 molecule bat i don't find the problem in my
> > simulation.
> > I have made too mach experiment bat i don't resolve the problem.
> >
> > I wont made a CO2 solvent box, it is cubic and with 36A of side. In
> > this box there's 300 molecules.
> > I use packmol to generate it. To avoid problem with periodic boundary
> > conditions i generate a box 4A smaller.
> > The box start in {2.0 2.0 2.0} and end at {34.0 34.0 34.0}, the
> > tolerance is 4.5A between any molecule.
> >
> > I generate correctly the psf. On this solvent box, I measure bounding
> > box that correspond with the input data of packmol.
> >
> >> Main< (Desktop) 51 % measure minmax [atomselect top all]
> > {2.0 2.0 2.0} {34.0 34.0 34.0}
> >
> > now i run the simulation on windows and namd plot out (below).
> > I run the simulation on linux end plot a different error!
> > Possible !?!?!?
> >
> > I don't understend where is located the error.
> >
> > TNX
> >
> > Sassaro
> >
> >
> >
> > ############### WINDOWS ###############
> >
> >
> >
> > Charm++: standalone mode (not using charmrun)
> > Info: NAMD 2.6b1 for Win32-i686
> > Info:
> > Info: Please visit http://www.ks.uiuc.edu/Research/namd/
> > Info: and send feedback or bug reports to namd_at_ks.uiuc.edu
> > Info:
> > Info: Please cite Kale et al., J. Comp. Phys. 151:283-312 (1999)
> > Info: in all publications reporting results obtained with NAMD.
> > Info:
> > Info: Based on Charm++/Converse 50900 for net-win32-smp
> > Info: Built Fri Jul 29 12:33:47 2005 by administrator on malta
> > Info: Running on 1 processors.
> > Info: 0 kB of memory in use.
> > Info: Changed directory to C:\Documents and Settings\Matteo\Desktop\co2
> > Info: Configuration file is min_solv.conf
> > TCL: Suspending until startup complete.
> > Info: SIMULATION PARAMETERS:
> > Info: TIMESTEP 1
> > Info: NUMBER OF STEPS 0
> > Info: STEPS PER CYCLE 20
> > Info: PERIODIC CELL BASIS 1 36 0 0
> > Info: PERIODIC CELL BASIS 2 0 36 0
> > Info: PERIODIC CELL BASIS 3 0 0 36
> > Info: PERIODIC CELL CENTER 0 0 0
> > Info: WRAPPING ALL CLUSTERS AROUND PERIODIC BOUNDARIES ON OUTPUT.
> > Info: LOAD BALANCE STRATEGY Other
> > Info: LDB PERIOD 4000 steps
> > Info: FIRST LDB TIMESTEP 100
> > Info: LDB BACKGROUND SCALING 1
> > Info: HOM BACKGROUND SCALING 1
> > Info: MAX SELF PARTITIONS 50
> > Info: MAX PAIR PARTITIONS 20
> > Info: SELF PARTITION ATOMS 125
> > Info: PAIR PARTITION ATOMS 200
> > Info: PAIR2 PARTITION ATOMS 400
> > Info: MIN ATOMS PER PATCH 100
> > Info: INITIAL TEMPERATURE 300
> > Info: CENTER OF MASS MOVING? NO
> > Info: DIELECTRIC 1
> > Info: EXCLUDE SCALED ONE-FOUR
> > Info: 1-4 SCALE FACTOR 1
> > Info: DCD FILENAME sim_out.dcd
> > Info: DCD FREQUENCY 10
> > Info: DCD FIRST STEP 10
> > Info: DCD FILE WILL CONTAIN UNIT CELL DATA
> > Info: NO EXTENDED SYSTEM TRAJECTORY OUTPUT
> > Info: NO VELOCITY DCD OUTPUT
> > Info: OUTPUT FILENAME sim_out
> > Info: BINARY OUTPUT FILES WILL BE USED
> > Info: NO RESTART FILE
> > Info: SWITCHING ACTIVE
> > Info: SWITCHING ON 10
> > Info: SWITCHING OFF 12
> > Info: PAIRLIST DISTANCE 14.5
> > Info: PAIRLIST SHRINK RATE 0.01
> > Info: PAIRLIST GROW RATE 0.01
> > Info: PAIRLIST TRIGGER 0.3
> > Info: PAIRLISTS PER CYCLE 2
> > Info: PAIRLISTS ENABLED
> > Info: MARGIN 0
> > Info: HYDROGEN GROUP CUTOFF 2.5
> > Info: PATCH DIMENSION 17
> > Info: ENERGY OUTPUT STEPS 1000
> > Info: USING VERLET I (r-RESPA) MTS SCHEME.
> > Info: C1 SPLITTING OF LONG RANGE ELECTROSTATICS
> > Info: PLACING ATOMS IN PATCHES BY HYDROGEN GROUPS
> > Info: NONBONDED FORCES EVALUATED EVERY 2 STEPS
> > Info: RANDOM NUMBER SEED 1150733412
> > Info: USE HYDROGEN BONDS? NO
> > Info: COORDINATE PDB mixture.pdb
> > Info: STRUCTURE FILE mixture.psf
> > Info: PARAMETER file: CHARMM format!
> > Info: PARAMETERS par_all27_prot_lipid.prm
> > Info: USING ARITHMETIC MEAN TO COMBINE L-J SIGMA PARAMETERS
> > Info: SUMMARY OF PARAMETERS:
> > Info: 1 BONDS
> > Info: 1 ANGLES
> > Info: 0 DIHEDRAL
> > Info: 0 IMPROPER
> > Info: 2 VDW
> > Info: 0 VDW_PAIRS
> > Info: ****************************
> > Info: STRUCTURE SUMMARY:
> > Info: 900 ATOMS
> > Info: 600 BONDS
> > Info: 300 ANGLES
> > Info: 0 DIHEDRALS
> > Info: 0 IMPROPERS
> > Info: 0 EXCLUSIONS
> > Info: 2697 DEGREES OF FREEDOM
> > Info: 900 HYDROGEN GROUPS
> > Info: TOTAL MASS = 13202.7 amu
> > Info: TOTAL CHARGE = 0 e
> > Info: *****************************
> > Info: Entering startup phase 0 with 0 kB of memory in use.
> > Info: Entering startup phase 1 with 0 kB of memory in use.
> > Info: Entering startup phase 2 with 0 kB of memory in use.
> > Info: Entering startup phase 3 with 0 kB of memory in use.
> > Info: PATCH GRID IS 2 (PERIODIC) BY 2 (PERIODIC) BY 2 (PERIODIC)
> > Info: REMOVING COM VELOCITY -0.0452495 -0.157624 -0.0325843
> > Info: LARGEST PATCH (6) HAS 117 ATOMS
> > Info: Entering startup phase 4 with 0 kB of memory in use.
> > Info: Entering startup phase 5 with 0 kB of memory in use.
> > Info: Entering startup phase 6 with 0 kB of memory in use.
> > Measuring processor speeds... Done.
> > Info: Entering startup phase 7 with 0 kB of memory in use.
> > Info: NONBONDED TABLE R-SQUARED SPACING: 0.0625
> > Info: NONBONDED TABLE SIZE: 769 POINTS
> > Info: ABSOLUTE IMPRECISION IN FAST TABLE ENERGY: 3.38813e-021 AT 11.8295
> > Info: RELATIVE IMPRECISION IN FAST TABLE ENERGY: 5.45773e-014 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN FAST TABLE FORCE: 9.09495e-013 AT 0.03125
> > Info: RELATIVE IMPRECISION IN FAST TABLE FORCE: 1.6856e-015 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN VDWA TABLE ENERGY: 1.75 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWA TABLE ENERGY: 6.75473e-014 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN VDWA TABLE FORCE: 19968 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWA TABLE FORCE: 6.16499e-015 AT 0.0441942
> > Info: ABSOLUTE IMPRECISION IN VDWB TABLE ENERGY: 4.1359e-025 AT 11.8295
> > Info: RELATIVE IMPRECISION IN VDWB TABLE ENERGY: 7.56853e-015 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN VDWB TABLE FORCE: 3.87741e-026 AT 11.9138
> > Info: RELATIVE IMPRECISION IN VDWB TABLE FORCE: 5.97409e-016 AT 11.9138
> > Info: Entering startup phase 8 with 0 kB of memory in use.
> > Info: Finished startup with 0 kB of memory in use.
> > TCL: Minimizing for 1000 steps
> > ETITLE: TS BOND ANGLE DIHED
> > IMPRP ELECT VDW BOUNDARY MISC
> > KINETIC TOTAL TEMP TOTAL2
> > TOTAL3 TEMPAVG PRESSURE GPRESSURE
> > VOLUME PRESSAVG GPRESSAVG
> >
> > ENERGY: 0 0.0607 -1.#IND 0.0000
> > 0.0000 -11.8639 -247.8474 0.0000
> > 0.0000 0.0000 -1.#IND 0.0000
> > -1.#IND -1.#IND 0.0000 -1.#IND
> > -1.#IND 46656.0000 -1.#IND -1.#IND
> >
> > INITIAL STEP: 1e-006
> > GRADIENT TOLERANCE: -1.#IND
> >
> > ############### LINUX ###############
> > Charm++: standalone mode (not using charmrun)
> > Info: NAMD 2.6b1 for Linux-i686
> > Info:
> > Info: Please visit http://www.ks.uiuc.edu/Research/namd/
> > Info: and send feedback or bug reports to namd_at_ks.uiuc.edu
> > Info:
> > Info: Please cite Kale et al., J. Comp. Phys. 151:283-312 (1999)
> > Info: in all publications reporting results obtained with NAMD.
> > Info:
> > Info: Based on Charm++/Converse 50900 for net-linux-iccstatic
> > Info: Built Fri Jul 29 11:30:30 CDT 2005 by jim on kyoto.ks.uiuc.edu
> > Info: Sending usage information to NAMD developers via UDP. Sent data
> > is:
> > Info: 1 NAMD 2.6b1 Linux-i686 1 delta1 uno
> > Info: Running on 1 processors.
> > Info: 7723 kB of memory in use.
> > Info: Configuration file is min_solv.conf
> > TCL: Suspending until startup complete.
> > Info: SIMULATION PARAMETERS:
> > Info: TIMESTEP 1
> > Info: NUMBER OF STEPS 0
> > Info: STEPS PER CYCLE 20
> > Info: PERIODIC CELL BASIS 1 36 0 0
> > Info: PERIODIC CELL BASIS 2 0 36 0
> > Info: PERIODIC CELL BASIS 3 0 0 36
> > Info: PERIODIC CELL CENTER 0 0 0
> > Info: WRAPPING ALL CLUSTERS AROUND PERIODIC BOUNDARIES ON OUTPUT.
> > Info: LOAD BALANCE STRATEGY Other
> > Info: LDB PERIOD 4000 steps
> > Info: FIRST LDB TIMESTEP 100
> > Info: LDB BACKGROUND SCALING 1
> > Info: HOM BACKGROUND SCALING 1
> > Info: MAX SELF PARTITIONS 50
> > Info: MAX PAIR PARTITIONS 20
> > Info: SELF PARTITION ATOMS 125
> > Info: PAIR PARTITION ATOMS 200
> > Info: PAIR2 PARTITION ATOMS 400
> > Info: MIN ATOMS PER PATCH 100
> > Info: INITIAL TEMPERATURE 300
> > Info: CENTER OF MASS MOVING? NO
> > Info: DIELECTRIC 1
> > Info: EXCLUDE SCALED ONE-FOUR
> > Info: 1-4 SCALE FACTOR 1
> > Info: DCD FILENAME sim_out.dcd
> > Info: DCD FREQUENCY 10
> > Info: DCD FIRST STEP 10
> > Info: DCD FILE WILL CONTAIN UNIT CELL DATA
> > Info: NO EXTENDED SYSTEM TRAJECTORY OUTPUT
> > Info: NO VELOCITY DCD OUTPUT
> > Info: OUTPUT FILENAME sim_out
> > Info: BINARY OUTPUT FILES WILL BE USED
> > Info: NO RESTART FILE
> > Info: SWITCHING ACTIVE
> > Info: SWITCHING ON 10
> > Info: SWITCHING OFF 12
> > Info: PAIRLIST DISTANCE 14.5
> > Info: PAIRLIST SHRINK RATE 0.01
> > Info: PAIRLIST GROW RATE 0.01
> > Info: PAIRLIST TRIGGER 0.3
> > Info: PAIRLISTS PER CYCLE 2
> > Info: PAIRLISTS ENABLED
> > Info: MARGIN 0
> > Info: HYDROGEN GROUP CUTOFF 2.5
> > Info: PATCH DIMENSION 17
> > Info: ENERGY OUTPUT STEPS 1000
> > Info: USING VERLET I (r-RESPA) MTS SCHEME.
> > Info: C1 SPLITTING OF LONG RANGE ELECTROSTATICS
> > Info: PLACING ATOMS IN PATCHES BY HYDROGEN GROUPS
> > Info: NONBONDED FORCES EVALUATED EVERY 2 STEPS
> > Info: RANDOM NUMBER SEED 1150733713
> > Info: USE HYDROGEN BONDS? NO
> > Info: COORDINATE PDB mixture.pdb
> > Info: STRUCTURE FILE mixture.psf
> > Info: PARAMETER file: CHARMM format!
> > Info: PARAMETERS par_all27_prot_lipid.prm
> > Info: USING ARITHMETIC MEAN TO COMBINE L-J SIGMA PARAMETERS
> > Info: SUMMARY OF PARAMETERS:
> > Info: 1 BONDS
> > Info: 1 ANGLES
> > Info: 0 DIHEDRAL
> > Info: 0 IMPROPER
> > Info: 2 VDW
> > Info: 0 VDW_PAIRS
> > Info: ****************************
> > Info: STRUCTURE SUMMARY:
> > Info: 900 ATOMS
> > Info: 600 BONDS
> > Info: 300 ANGLES
> > Info: 0 DIHEDRALS
> > Info: 0 IMPROPERS
> > Info: 0 EXCLUSIONS
> > Info: 2697 DEGREES OF FREEDOM
> > Info: 900 HYDROGEN GROUPS
> > Info: TOTAL MASS = 13202.7 amu
> > Info: TOTAL CHARGE = 0 e
> > Info: *****************************
> > Info: Entering startup phase 0 with 7943 kB of memory in use.
> > Info: Entering startup phase 1 with 7943 kB of memory in use.
> > Info: Entering startup phase 2 with 8047 kB of memory in use.
> > Info: Entering startup phase 3 with 8054 kB of memory in use.
> > Info: PATCH GRID IS 2 (PERIODIC) BY 2 (PERIODIC) BY 2 (PERIODIC)
> > Info: REMOVING COM VELOCITY -0.289994 -0.0983423 -0.125683
> > Info: LARGEST PATCH (6) HAS 117 ATOMS
> > Info: Entering startup phase 4 with 8195 kB of memory in use.
> > Info: Entering startup phase 5 with 8195 kB of memory in use.
> > Info: Entering startup phase 6 with 8195 kB of memory in use.
> > Measuring processor speeds... Done.
> > Info: Entering startup phase 7 with 8196 kB of memory in use.
> > Info: NONBONDED TABLE R-SQUARED SPACING: 0.0625
> > Info: NONBONDED TABLE SIZE: 769 POINTS
> > Info: ABSOLUTE IMPRECISION IN FAST TABLE ENERGY: 3.38813e-21 AT 12.0805
> > Info: RELATIVE IMPRECISION IN FAST TABLE ENERGY: 1.04308e-13 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN FAST TABLE FORCE: 1.69407e-21 AT 11.9974
> > Info: RELATIVE IMPRECISION IN FAST TABLE FORCE: 1.09878e-16 AT 11.9974
> > Info: ABSOLUTE IMPRECISION IN VDWA TABLE ENERGY: 2.75 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWA TABLE ENERGY: 6.72641e-14 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN VDWA TABLE FORCE: 6656 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWA TABLE FORCE: 2.055e-15 AT 0.0441942
> > Info: ABSOLUTE IMPRECISION IN VDWB TABLE ENERGY: 7.45058e-09 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWB TABLE ENERGY: 7.56853e-15 AT 11.9138
> > Info: ABSOLUTE IMPRECISION IN VDWB TABLE FORCE: 7.62939e-06 AT 0.0441942
> > Info: RELATIVE IMPRECISION IN VDWB TABLE FORCE: 7.96546e-16 AT 11.9138
> > Info: Entering startup phase 8 with 8713 kB of memory in use.
> > Info: Finished startup with 9867 kB of memory in use.
> > TCL: Minimizing for 1000 steps
> > ETITLE: TS BOND ANGLE DIHED
> > IMPRP ELECT VDW BOUNDARY MISC
> > KINETIC TOTAL TEMP TOTAL2
> > TOTAL3 TEMPAVG PRESSURE GPRESSURE
> > VOLUME PRESSAVG GPRESSAVG
> >
> > ENERGY: 0 0.0607 0.0056 0.0000
> > 0.0000 -11.8639 -247.8474 0.0000
> > 0.0000 0.0000 -259.6450 0.0000
> > -259.6450 -259.6450 0.0000 99999999.9999
> > 99999999.9999 46656.0000 99999999.9999 99999999.9999
> >
> > INITIAL STEP: 1e-06
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > RESTARTING CONJUGATE GRADIENT ALGORITHM
> > INITIAL STEP: 2.5e-07
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > RESTARTING CONJUGATE GRADIENT ALGORITHM
> > INITIAL STEP: 6.25e-08
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 1.5625e-08
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 3.90625e-09
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 9.76562e-10
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 2.44141e-10
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 6.10352e-11
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 1.52588e-11
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 3.8147e-12
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 9.53674e-13
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 2.38419e-13
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 5.96046e-14
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 1.49012e-14
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 3.72529e-15
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 9.31323e-16
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 2.32831e-16
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 5.82077e-17
> > GRADIENT TOLERANCE: nan
> > BRACKET: 0 0 nan nan nan
> > NEW SEARCH DIRECTION
> > INITIAL STEP: 1.45519e-17
>
>
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:42:14 CST