Skip to content

Commit 3b0ccec

Browse files
authored
Small tweaks to MTS-RPC tutorial (#100)
* Small tweaks to MTS-RPC tutorial
1 parent d605cfc commit 3b0ccec

File tree

8 files changed

+6
-4
lines changed

8 files changed

+6
-4
lines changed

examples/pi-mts-rpc/data/h2o_md.xml

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@
44
<trajectory filename='pos' stride='100' format='xyz' cell_units='angstrom'> positions{angstrom} </trajectory>
55
<checkpoint filename='checkpoint' stride='1000' overwrite='True'/>
66
</output>
7-
<total_steps> 10000 </total_steps>
7+
<total_steps> 20000 </total_steps>
88
<ffsocket mode='unix' name='qtip4pf' pbc='false'>
99
<address>qtip4pf-md</address>
1010
</ffsocket>
@@ -23,7 +23,7 @@
2323
<dynamics mode='nvt'>
2424
<timestep units='femtosecond'> 0.5 </timestep>
2525
<thermostat mode="svr">
26-
<tau units="femtosecond"> 400 </tau>
26+
<tau units="femtosecond"> 500 </tau>
2727
</thermostat>
2828
</dynamics>
2929
</motion>

examples/pi-mts-rpc/data/h2o_mts.xml

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@
44
<trajectory filename='pos' stride='100' format='xyz' cell_units='angstrom'> positions{angstrom} </trajectory>
55
<checkpoint filename='checkpoint' stride='1000' overwrite='True'/>
66
</output>
7-
<total_steps> 2500 </total_steps>
7+
<total_steps> 5000 </total_steps>
88
<ffsocket mode='unix' name='qtip4pf' pbc='false'>
99
<address>qtip4pf-mts-full</address>
1010
</ffsocket>
@@ -32,7 +32,7 @@
3232
<timestep units='femtosecond'> 2.0 </timestep>
3333
<nmts>[1,4]</nmts>
3434
<thermostat mode="svr">
35-
<tau units="femtosecond"> 400 </tau>
35+
<tau units="femtosecond"> 500 </tau>
3636
</thermostat>
3737
</dynamics>
3838
</motion>

examples/pi-mts-rpc/mts-rpc.py

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -663,6 +663,8 @@ def autocorrelate(x, xbar=None, normalize=True):
663663
ax[1].plot(rpcmts_output["time"], rpcmts_output["kinetic_cv"], "m.", label="RPC-MTS")
664664
ax[1].set_xlabel("t / ps")
665665
ax[1].set_ylabel("K / eV")
666+
ax[0].set_xlim(0.0, 5.0)
667+
ax[1].set_xlim(0.0, 5.0)
666668
ax[0].legend()
667669

668670
# %%
3.14 KB
Loading
5.6 KB
Loading
-34.8 KB
Binary file not shown.
3.57 KB
Loading

examples/pi-mts-rpc/rpc-4.png

5.08 KB
Loading

0 commit comments

Comments
 (0)