is a stochastic force of magnitude (see \[app:1\]) $$\bm{F}_{i,T} =\Lambda_{i,j}\sqrt{\frac{2Bk_{B}TL_{i}}{\Delta t_{j}}}\bm{s}_{i},$$ where $\Lambda_{i,j}$ is a normal distributed Gaussian random number; $k_{B}$ is Boltzmann’s constant; $T$ is absolute temperature; $\Delta t_{j}$ is the increment of time at simulation step $j$, and the unit vector $\bm{s}_{i}$ denotes the effective glide direction of node $i$.

   Temperature(K)   Drag coefficient($\mu$Pa$\bm{\cdot}$s)
  ---------------- ----------------------------------------
        300                           30
        600                           60
        900                           90
        1200                         120

  : Temperature dependence of drag coefficient, after [@Hiratani2002_JEMT].

\[table:1\]

[|c|c|c|c|c|&gt;p[1.2cm]{}|c|c|c|]{} & & & &\
& (i) &(ii) &(i) &(ii) &(i) &(ii) &(i) &(ii)\
300 & 299 & 299 &52.7 &16.6 &52.3 &17.2 & $\sim$31.6&$\sim$36.5\
600 & 596 & 596 &74.5 &23.5 &72.9 &23.6 & $\sim$21.2& $\sim$16.7\
900 & 894 & 896 &91.2 &28.8 &92.5 &29.8 & $\sim$13.1& $\sim$15.9\
1200 & 1189& 1189&105.3 &33.3 &105.7 & 32.4& $\sim$9.0& $\sim$8.2\

\[table:2\]

Validation: thermal equilibrium properties {#sec:3.2}
------------------------------------------

To validate our model implementation, we consider an initially straight edge dislocation which we assume to run in $x$ direction and move in the $xy$ plane under the influence of thermal forces. After sufficient simulation time, once thermal equilibrium has been reached, the mean kinetic energy for each degree of freedom should be $k_{B}T/2$ and the distribution of nodal velocities 