@@ -2930,12 +2930,13 @@ \subsubsection{Deform the structure}
29302930the \flecmd {relax.data} file. Copy the following lines into \flecmd {deform.lmp}:
29312931
29322932Next, let us use \lmpcmd {fix nvt} instead of \lmpcmd {fix npt} to apply a
2933- Nosé-Hoover thermostat without any barostat:
2933+ Nosé-Hoover thermostat without a barostat:
29342934\ begin{lstlisting}
29352935fix mynvt all nvt temp 300.0 300.0 100
29362936timestep 0.5
29372937\end {lstlisting }
2938- Here, no barostat is used because the box volume will be imposed by the \lmpcmd {fix deform}.
2938+ Here, a barostat is not used because the change in the box volume will be imposed
2939+ by the \lmpcmd {fix deform}.
29392940
29402941\begin {figure }
29412942\includegraphics [width=\linewidth ]{SIO-deformed-charge}
@@ -2946,60 +2947,58 @@ \subsubsection{Deform the structure}
29462947\label {fig:SIO-deformed-charge }
29472948\end {figure }
29482949
2949- Let us run for 5000 steps without deformation, then apply the \textit {fix deform }
2950- for elongating progressively the box along \textit { x } during 25000 steps. Add
2951- the following line to \flecmd {input .lmp}:
2950+ Let us run for 5000 steps without deformation, then apply the \lmpcmd {fix deform}
2951+ to progressively elongate the box along the $ x $ -axis during 25000 steps. Add
2952+ the following line to \flecmd {deform .lmp}:
29522953\ begin{lstlisting}
29532954run 5000
29542955
29552956fix mydef all deform 1 x erate 5e-5
29562957
29572958run 25000
29582959
2959- write_data silica-deformed .data
2960+ write_data deform .data
29602961\end {lstlisting }
2961- Run the \flecmd {input .lmp} file using LAMMPS. During the deformation, the charge
2962- values progressively evolve until the structure eventually breaks down. After the
2962+ Run the \lmpcmd {deform .lmp} file using LAMMPS. During the deformation, the charge
2963+ values progressively evolve until the structure eventually breaks down. After the
29632964structure breaks down, the charges equilibrate near new average values that differ
2964- from the starting averages (Fig.\, \ref {fig:SIO-deformed-charge }). The difference
2965+ from the initial averages (Fig.\, \ref {fig:SIO-deformed-charge }). The difference
29652966between the initial and the final charges can be explained by the presence of
2966- defects as well as new solid/vacuum interfaces, and the fact that surface atoms
2967- typically have different charges compared to bulk atoms. There is also a strong
2968- increase in temperature during the rupture of the material (Fig.\, \ref {fig:SIO-temperature }).
2969- At the end of the deformation, one can visualize the broken material using VMD.
2970- Notice the different charge values of the atoms located near the vacuum interfaces,
2971- compared to the atoms located in the bulk of the material (Fig.\, \ref {fig:SIO-deformed }).
2972-
2973- \begin {figure }
2974- \includegraphics [width=\linewidth ]{SIO-temperature}
2975- \caption {Evolution of the temperature $ T$ of the silica system over time $ t$ .
2976- The material ruptures near $ t = 10 ~\text {ps}$ .}
2977- \label {fig:SIO-temperature }
2978- \end {figure }
2967+ defects, as well as new solid/vacuum interfaces, and the fact that surface atoms
2968+ typically have different charges compared to bulk atoms (Fig.\, \ref {fig:SIO-deformed }).
2969+ There is also a sharp increase in temperature during the rupture of
2970+ the material (Fig.\, \ref {fig:SIO-temperature }).
29792971
29802972\begin {figure }
29812973\includegraphics [width=\linewidth ]{SIO-deformed}
29822974\caption {Amorphous silicon oxide after deformation. The atoms are colored by their
2983- charges. Dangling oxygen groups appear in greenish, bulk Si atoms with a charge of
2975+ charges. Dangling oxygen groups appear in greenish, bulk Si atoms with a charge of
29842976about $ 1.8 ~\text {e}$ appear in red/orange, and bulk O atoms with a charge of
29852977about $ -0.9 ~ \text {e}$ appear in blue.}
29862978\label {fig:SIO-deformed }
29872979\end {figure }
29882980
2989- One can have a look at the charge distribution after deformation, as well as during
2990- the deformation (Fig.\, \ref {fig:SIO-distribution-bis }). As expected, the final
2991- charge distribution slightly differs from the previously calculated. If
2992- no new species were formed during the simulation, the \textit {species.log } file
2993- should resemble:
2981+ \begin {figure }
2982+ \includegraphics [width=\linewidth ]{SIO-temperature}
2983+ \caption {Evolution of the temperature $ T$ of the silica system over time $ t$ .
2984+ The material ruptures near $ t = 10 ~\text {ps}$ .}
2985+ \label {fig:SIO-temperature }
2986+ \end {figure }
2987+
2988+ You can examine the charge distribution after deformation, as well as during
2989+ deformation (Fig.\, \ref {fig:SIO-distribution-bis }). As expected, the final
2990+ charge distribution slightly differs from the previously calculated one. If
2991+ no new species were formed during the simulation, the \flecmd {species.log} file
2992+ should look like this:
29942993\ begin{lstlisting}
29952994# Timestep No_Moles No_Specs O384Si192
29962995 5 1 1 1
29972996(...)
29982997# Timestep No_Moles No_Specs O384Si192
29992998 30000 1 1 1
30002999\end {lstlisting }
3001- Sometimes, $ \text {O}_2 $ molecules are formed during the deformation. If this is
3002- the case, a new column \textit {O2 } appears in the \textit {species.log } file.
3000+ Sometimes, $ \text {O}_2 $ molecules are formed during the deformation. If this occurs,
3001+ a new column \lmpcmd {O2} appears in the \flecmd {species.log} file.
30033002
30043003\begin {figure }
30053004\includegraphics [width=\linewidth ]{SIO-distribution-bis}
@@ -3010,14 +3009,13 @@ \subsubsection{Deform the structure}
30103009\end {figure }
30113010
30123011\subsubsection {Decorate the surface }
3013- In ambient conditions, some of the surface $ \text {SiO}_2 $ atoms are chemically
3012+
3013+ Under ambient conditions, some of the surface $ \text {SiO}_2 $ atoms become chemically
30143014passivated by forming covalent bonds with hydrogen (H) atom \cite {sulpizi2012silica }.
30153015Let us add hydrogen atoms randomly to the cracked silica and observe how the
3016- system evolves. Next to \flrcmd {RelaxSilica/} and \flrcmd {Deform/}, create a folder,
3017- and call it \flrcmd {Decorate/}. Then, let us modify the previously generated data
3018- file \flecmd {silica-deformed.data} and make space for a third atom type.
3019- Copy \flecmd {silica-deformed.data} from the \flrcmd {Deform/} folder, and modify
3020- the first lines as follows:
3016+ system evolves. To do so, we first need to modify the previously generated data
3017+ file \flecmd {deform.data.data} and make space for a third atom type.
3018+ Modify the first lines of \flecmd {deform.data} as follows:
30213019\ begin{lstlisting}
30223020576 atoms
302330213 atom types
@@ -3032,6 +3030,7 @@ \subsubsection{Decorate the surface}
30323030
30333031(...)
30343032\end {lstlisting }
3033+
30353034Create a file called \flecmd {input.lmp} into the \flrcmd {Decorate/} folder, and
30363035copy the following lines into it:
30373036\ begin{lstlisting}
0 commit comments