Skip to content

Add box resizing to MD integrators #234

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
clpetix wants to merge 12 commits into
base: main
Choose a base branch
from
Draft

Add box resizing to MD integrators #234

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
28696cc
Add volume scaling for Langevin ynamics.
clpetix Feb 7, 2024
1a846ab
Update implementation to match HOOMD required syntax.
clpetix Feb 7, 2024
52a7edb
Update period to prevent triggering repeatedly.
clpetix Feb 7, 2024
b2a16d9
Check if barostat is set before trying to scale box.
clpetix Feb 10, 2024
8576067
Add barostat to LAMMPS and Dilute simulation methods.
clpetix Feb 12, 2024
ac0d4c3
Remove default argument from backend operation.
clpetix Feb 12, 2024
798586d
Set freud version less than 3
clpetix Apr 6, 2024
b87c3c3
Require freud version 2.
clpetix Apr 6, 2024
28f55ce
Revert freud change.
clpetix Apr 6, 2024
4636ea9
Merge branch 'main' into feature/box-resize
clpetix May 3, 2024
4d0567c
Merge branch 'main' into feature/box-resize
clpetix May 3, 2024
d4b929a
Merge branch 'main' into feature/box-resize
mphoward Jun 27, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
60 changes: 59 additions & 1 deletion src/relentless/simulate/hoomd.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -656,9 +656,10 @@ class RunLangevinDynamics(_Integrator):

"""

def __init__(self, steps, timestep, T, friction, seed, analyzers):
def __init__(self, steps, timestep, T, friction, seed, analyzers, barostat=None):
super().__init__(steps, timestep, analyzers)
self.T = T
self.barostat = barostat
self.friction = friction
self.seed = seed

Expand Down Expand Up @@ -709,6 +710,63 @@ def _call_v2(self, sim):
for analyzer in self.analyzers:
analyzer.pre_run(sim, self)

if self.barostat is not None:
if isinstance(self.barostat, extent.Extent):
period = None
if sim.dimension == 2:
Lx, Ly, xy = self.barostat.as_array("HOOMD")
Lz, xz, yz = None, None, None
else:
Lx, Ly, Lz, xy, xz, yz = self.barostat.as_array("HOOMD")
elif all(isinstance(n, extent.Extent) for n in self.barostat):
period = 1
if sim.dimension == 2:
Lx1, Ly1, xy1 = self.barostat[0].as_array("HOOMD")
Lz1, xz1, yz1 = None, None, None

Lx2, Ly2, xy2 = self.barostat[1].as_array("HOOMD")
Lz2, xz2, yz2 = None, None, None

Lz, xz, yz = None, None, None
Lx = hoomd.variant.linear_interp(
[(0, Lx1), (self.steps - 1, Lx2)]
)
Ly = hoomd.variant.linear_interp(
[(0, Ly1), (self.steps - 1, Ly2)]
)
xy = hoomd.variant.linear_interp(
[(0, xy1), (self.steps - 1, xy2)]
)
else:
Lx1, Ly1, Lz1, xy1, xz1, yz1 = self.barostat[0].as_array(
"HOOMD"
)
Lx2, Ly2, Lz2, xy2, xz2, yz2 = self.barostat[1].as_array(
"HOOMD"
)

Lx = hoomd.variant.linear_interp(
[(0, Lx1), (self.steps - 1, Lx2)]
)
Ly = hoomd.variant.linear_interp(
[(0, Ly1), (self.steps - 1, Ly2)]
)
Lz = hoomd.variant.linear_interp(
[(0, Lz1), (self.steps - 1, Lz2)]
)
xy = hoomd.variant.linear_interp(
[(0, xy1), (self.steps - 1, xy2)]
)
xz = hoomd.variant.linear_interp(
[(0, xz1), (self.steps - 1, xz2)]
)
yz = hoomd.variant.linear_interp(
[(0, yz1), (self.steps - 1, yz2)]
)

hoomd.update.box_resize(
Lx=Lx, Ly=Ly, Lz=Lz, xy=xy, xz=xz, yz=yz, period=period
)
hoomd.run(self.steps)

for analyzer in self.analyzers:
Expand Down
6 changes: 5 additions & 1 deletion src/relentless/simulate/md.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -118,9 +118,12 @@ class RunLangevinDynamics(_Integrator):

"""

def __init__(self, steps, timestep, T, friction, seed, analyzers=None):
def __init__(
self, steps, timestep, T, friction, seed, barostat=None, analyzers=None
):
super().__init__(steps, timestep, analyzers)
self.T = T
self.barostat = barostat
self.friction = friction
self.seed = seed

Expand All @@ -130,6 +133,7 @@ def _make_delegate(self, sim):
steps=self.steps,
timestep=self.timestep,
T=self.T,
barostat=self.barostat,
friction=self.friction,
seed=self.seed,
analyzers=self.analyzers,
Expand Down