Hello,
I’m attempting to run a system with rigid and non-rigid bodies.
I compress the box using NPT and then want to switch to isobaric but I get a “non-numeric pressure” when I switch to isobaric.
“”""
This runs with no problems
fix sphere_npt grp_sphere temp 1.5 0.1 10.0 iso 0.0 5.0 500.0 dilate all
fix plate_nvt grp_plate rigid/nvt molecule temp 1.5 0.1 10.0
run 200000
unfix sphere_npt
unfix plate_nvt
non-numeric pressure error occurs here, at first timestep
fix sphere_npt grp_sphere temp 0.1 0.1 10.0 iso 5.0 5.0 500.0 dilate all
fix plate_nvt grp_plate rigid/nvt molecule temp 0.1 0.1 10.0
run 100000
unfix sphere_npt
unfix plate_nvt
“”""
The pressure on the last timestep of the compression phase is 5.0 so I thought it would be fine once it switched to isobaric. Any ideas why things are going wrong?
If more information is needed, please let me know!
Thanks,
Stacey