@@ -4039,7 +4039,7 @@ \subsubsection{Creating the system}
40394039\begin {verbatim }
40404040molecule CNT cnt.molecule
40414041create_atoms 0 random 5 8305 NULL overlap 3 &
4042- maxtry 500 mol CNT 7687
4042+ maxtry 500 mol CNT 7687
40434043\end {verbatim }
40444044}
40454045Let us use the \textit {minimize } command to reduce the energy of the system:
@@ -4072,35 +4072,41 @@ \subsubsection{Creating the system}
40724072 adiam n 0.3 adiam na 0.3
40734073\end {verbatim }
40744074}
4075- Then, let us perform a short equilibration of system using \textit {fix npt }:
4075+ Then, let us perform a short equilibration of system using \textit {fix npt }
4076+ with a isotropic couplage. To speed up the equilibration of the system, let us
4077+ impose a relatively large pressure of 1000\, atm for the first 10\, ps:
40764078{\normalsize
40774079\begin {verbatim }
40784080velocity all create 300 1829 dist gaussian &
4079- mom yes rot yes
4080- fix mynpt all npt temp 300 300 100 iso 1000 1000 1000
4081+ mom yes rot yes
4082+ fix mynpt all npt temp 300 300 100 &
4083+ iso 1000 1000 1000
40814084
40824085thermo 1000
40834086thermo_style custom step temp pe etotal press density vol
40844087
40854088run 10000
4086- unfix mynpt
4087-
4089+ \end {verbatim }
4090+ }
4091+ Then, for the following 20\, ps, let us continue the equilibration
4092+ with an imposed pressure of 1\, atm, and write the final state of the
4093+ system in a file named \textit {cnt-nylon-mix.data }:
4094+ {\normalsize
4095+ \begin {verbatim }
40884096fix mynpt all npt temp 300 300 100 iso 1 1 1000
40894097
40904098run 20000
40914099write_data cnt-nylon-mix.data
40924100\end {verbatim }
40934101}
4094- To speed up the equilibration of the system, a relatively large pressure of
4095- 1000\, atm is imposed for the first 10\, ps. Then for the following 20\, ps, the
4096- desired pressure of 1\, atm is imposed.
40974102
40984103\subsubsection {Atom maps and molecule templates }
40994104
4100- The REACTER protocol requires an \textit {atom map } to be provided by the user.
4101- This map contains information about the reaction mechanism, the atoms that are
4102- the initiators to the reaction, as well as distance cutoffs. Pre-reaction and
4103- post-reaction molecule templates must also be provided.
4105+ The REACTER protocol relies on classical force fields. For the reaction to
4106+ occur, an \textit {atom map } containing information about the reaction mechanism
4107+ must be provided by the user. The atom map specifies which atoms are
4108+ the initiators to the reaction and what are the distance cutoffs. In addition,
4109+ pre-reaction and post-reaction molecule templates must be provided.
41044110
41054111\section* {Author Contributions }
41064112% S.G. to update
0 commit comments