TiO2 MD simulation


Forum Vet
(Posted for Saeed that is having trouble using this website)

I am new with both the NWCHEM and molecular dynamics
calculations. I want to simulate TiO2 in the crystalline state. I have
written .par file for Ti and O atoms using literature data which just
include electrostatic and van der Walls parameters and not the bond or
angle parameters. This is because I found that in this ionic system we
do not need to define parameters for bond, angle and etc. I have
written .pdb file using crystalline structure data. I have also
included just the name of my single segment in a .seq file to avoid
automatic inclusion of solvent by the software. The .frg, .sgm, .rst
and .top files were generated using prepare task. Of course, in the
.top and .rst file generation, the bond, angle and etc. definitions
were deleted from .sgm file as it was urged before. Unfortunately when
I want to run a classical molecular dynamics simulation using these
.rst and .top files, I get that "EOF was encountered in the topology
file". Of course, out file just show error that "* 0: Error in
number of solute fractions 0". Hereby I would like to help me to
overcome this problem. All files are attached too. Please let me know
if you need more information.

Saeed

Clicked A Few Times
First, I would strongly discourage using NWChem classical MD module to simulate TiO2 crystalline state, especially if you are new to this area. MD module is oriented more towards soft mater systems. You would be much better off looking into plane wave module.

Some of the problems you see likely coming from the fragment file. You should remove all the connectivity fields from there (the ones that have two numbers, e.g. "1 2"). You also should put charges into the fragment file. Your setup also has no information about the periodic nature of the system (e.g. cell size).

Marat

Clicked A Few Times
Request
Dear Marat

Thanks a lot. We want to use NWCHEM for classical MD simulation of organic molecules' adsorption from solution on the TiO2 surface. At first we want to be able to use NWCHEM for classical MD simulation. We have removed all the connectivity fields from fragment file. unfortunately we did not find that how add cell size ( Is it necessary, if we do not want to apply periodicity?). Following are used and created files.

Best Regards
Saeed

Clicked A Few Times
Request
Dear Marat

Thanks a lot. We want to use NWCHEM for classical MD simulation of organic molecules' adsorption from solution on the TiO2 surface. At first we want to be able to use NWCHEM for classical MD simulation. We have removed all the connectivity fields from fragment file. unfortunately we did not find that how add cell size ( Is it necessary, if we do not want to apply periodicity?). Following are used and created files.

Best Regards
Saeed

pdb file;(TiO2.pdb)

ATOM 1 Ti1 TIO 1 5.164 0.702 -0.905
ATOM 2 O1 TIO 1 7.279 0.372 -0.689
ATOM 3 O2 TIO 1 3.184 0.430 -1.379
ATOM 4 Ti2 TIO 1 1.178 0.699 -0.892
ATOM 5 O3 TIO 1 1.158 2.357 0.210
ATOM 6 O4 TIO 1 1.033 1.887 -2.360
ATOM 7 Ti3 TIO 1 1.117 3.460 -1.242
ATOM 8 O5 TIO 1 5.193 2.406 0.075
ATOM 9 O6 TIO 1 5.325 1.847 -2.470
ATOM 10 Ti4 TIO 1 5.242 3.430 -1.446
ATOM 11 O7 TIO 1 3.149 3.820 -1.537
ATOM 12 O8 TIO 1 7.306 3.822 -1.502
END

fragment file;(TIO.frg)

# This is an automatically generated fragment file
  1. Atom types and connectivity were derived from coordinates
  2. Atomic partial charges are crude estimates
  3. 04/30/14 18:11:29

$TIO
  12    1    1    0
TIO
   1Ti1   Ti        0    0    0    1    1   +2.196000    0.000000
2 O1 Oi 0 0 0 1 1 -1.098000 0.000000
3 O2 Oi 0 0 0 1 1 -1.098000 0.000000
4Ti2 Ti 0 0 0 1 1 +2.196000 0.000000
5 O3 Oi 0 0 0 1 1 -1.098000 0.000000
6 O4 Oi 0 0 0 1 1 -1.098000 0.000000
7Ti3 Ti 0 0 0 1 1 +2.196000 0.000000
8 O5 Oi 0 0 0 1 1 -1.098000 0.000000
9 O6 Oi 0 0 0 1 1 -1.098000 0.000000
10Ti4 Ti 0 0 0 1 1 +2.196000 0.000000
11 O7 Oi 0 0 0 1 1 -1.098000 0.000000
12 O8 Oi 0 0 0 1 1 -1.098000 0.000000


segment file; (TIO.sgm)

  1. This is an automatically generated segment file
  2. 04/13/14 13:20:08
   4.600000
12 0 0 0 0 0 1 1
0.000000
1Ti1 0 0 0 1 1
Ti 0.000000 0.000000
2 O1 0 0 0 1 1
Oi 0.000000 0.000000
3 O2 0 0 0 1 1
Oi 0.000000 0.000000
4Ti2 0 0 0 1 1
Ti 0.000000 0.000000
5 O3 0 0 0 1 1
Oi 0.000000 0.000000
6 O4 0 0 0 1 1
Oi 0.000000 0.000000
7Ti3 0 0 0 1 1
Ti 0.000000 0.000000
8 O5 0 0 0 1 1
Oi 0.000000 0.000000
9 O6 0 0 0 1 1
Oi 0.000000 0.000000
10Ti4 0 0 0 1 1
Ti 0.000000 0.000000
11 O7 0 0 0 1 1
Oi 0.000000 0.000000
12 O8 0 0 0 1 1
Oi 0.000000 0.000000

Clicked A Few Times
Dear Marat

Thanks a lot. We want to use NWCHEM for classical MD simulation of organic molecules' adsorption from solution on the TiO2 surface. At first we want to be able to use NWCHEM for classical MD simulation. We have removed all the connectivity fields from fragment file. unfortunately we did not find that how add cell size ( Is it necessary, if we do not want to apply periodicity?). Following are used and created files.

Best Regards
Saeed

parameter file; (amber.par)

This is the AMBER96 user defined parameter file for NWChem 3.2 and ARGOS 7.0
Electrostatic 1-4 scaling factor 0.833333
Relative dielectric constant 1.000000
Parameters epsilon R*
Atoms
Ti 47.86700 2.58660E+00 2.20000E-01 1 1111111111
            22 1.29330E+00 2.20006E-01                              TPS000106 JCP 85, 6720-6727 (1986)
Oi 15.99900 1.41810E+00 3.24500E-01 1 1111111111
             8 0.70905E+00 3.24500E-01                              TPS000106 JCC 12, 1125-1128 (1991)
Cross
Bonds
Angles
Proper dihedrals
Improper dihedrals
Atom types
Ti 22 0 0 0 0 0 0 0 0 0
                                0  0      0      0      0
0 0 0 0 0

Oi 8 0 0 0 0 0 0 0 0 0
                                0  0      0      0      0
0 0 0 0 0

End


sequence file; (TiO2.seq)

   1TIO   
stop



prepare nw file; (TiO2.nw)

echo
memory noverify heap 1 mb stack 32 mb global 8 mb

start TiO2

prepare
amber
source TiO2.pdb
new_top new_rst
end
task prepare

Clicked A Few Times
a
aa


Forum >> NWChem's corner >> QMMM