Compressing the System¶
How do I compress the system to a target density?
What is a packing fraction?
Show how to compute the packing fraction of a system.
Explain how how an Updater is an operation that modifies the system when its Trigger returns
Demonstrate using the QuickCompress updater to achieve a target packing fraction.
Demonstrate using the MoveSize tuner to adjust the trial move size.
import copy import math import hoomd
render function in the next (hidden) cell will render a snapshot using fresnel.
This is not intended as a full tutorial on fresnel - see the fresnel user documentation if you would like to learn more.
Self-assembly in hard particle systems typically occurs at a packing fraction above 0.5. The packing fraction is the ratio of the volume occupied by the particles to the volume of the periodic box.
So far, this tutorial as randomized a system of N octahedra in a box with a very low packing fraction and stored that in
random.gsd. Initialize a Simulation with this configuration and see what packing fraction it is at:
cpu = hoomd.device.CPU() sim = hoomd.Simulation(device=cpu, seed=20) sim.create_state_from_gsd(filename='random.gsd')
Compute the volume of the octahedron:
a = math.sqrt(2) / 2 V_particle = 1 / 3 * math.sqrt(2) * a**3
Compute the packing fraction:
initial_packing_fraction = (sim.state.N_particles * V_particle / sim.state.box.volume) print(initial_packing_fraction)
As you can see, this packing fraction is very low and the box volume needs to be significantly reduced to achieve a packing fraction above 0.5.
Use HPMC to move particles into non-overlapping configurations while you compress the system. Set up the HPMC integrator for the octahedron simulations:
mc = hoomd.hpmc.integrate.ConvexPolyhedron() mc.shape['octahedron'] = dict(vertices=[ (-0.5, 0, 0), (0.5, 0, 0), (0, -0.5, 0), (0, 0.5, 0), (0, 0, -0.5), (0, 0, 0.5), ]) sim.operations.integrator = mc
The QuickCompress updater¶
An Updater is a type of operation in HOOMD-blue that makes changes to the state. To use an Updater, first instantiate the object, assign a Trigger, and add it to the Simulation. Simulation will apply the Updater on time steps where the Trigger returns
True. The Periodic trigger returns
QuickCompress is an Updater that works with HPMC to quickly compress the box to a target volume. When triggered, QuickCompress reduces the box volume by a scale factor, while allowing slight overlaps between the particles. It then waits for the translation and rotation trial moves to remove these overlaps before it reduces the volume again. This process temporarily produces invalid system configurations, but is much quicker than a process that does not allow temporary overlaps.
Compute the final box size with a packing fraction above 0.5 and configure a QuickCompress to trigger every 10 time steps.
initial_box = sim.state.box final_box = hoomd.Box.from_box(initial_box) final_packing_fraction = 0.57 final_box.volume = sim.state.N_particles * V_particle / final_packing_fraction compress = hoomd.hpmc.update.QuickCompress(trigger=hoomd.trigger.Periodic(10), target_box=final_box)
Add the Updater to the Simulation:
The MoveSize tuner¶
A Tuner is another type of operation. Tuners make changes to other operations to improve performance. In HPMC, the translation and rotation trial move sizes have a huge performance impact. When the move size is too small it takes many time steps to make appreciable changes to the system. When the move size is too large very few moves are accepted and it again takes many time steps to make appreciable changes. The system makes the most progress at moderate move sizes and, in
most cases, the optimal acceptance ratio is 20%. The MoveSize tuner monitors the acceptance ratio and adjusts
a to achieve the target.
The optimal move size depends on the density of the system. QuickCompress changes the density rapidly during compression, so use the MoveSize tuner to adjust the move sizes regularly:
periodic = hoomd.trigger.Periodic(10) tune = hoomd.hpmc.tune.MoveSize.scale_solver(moves=['a', 'd'], target=0.2, trigger=periodic, max_translation_move=0.2, max_rotation_move=0.2) sim.operations.tuners.append(tune)
Run until complete¶
When the QuickCompress updater achieves the target box size and there are no overlaps between particles, the compression process is complete. The number of time steps needed to achieve this varies based on parameters. Check the
compress.complete property regularly and stop running the simulation when the compression completes:
while not compress.complete and sim.timestep < 1e6: sim.run(1000)
sim.timestep < 1e6 check ensures that this while loop will not waste resources in cases where the compression will never complete. The loop should complete before that point:
Check to see if the compression completed successfully:
if not compress.complete: raise RuntimeError("Compression failed to complete")
The MoveSize tuner should have adjusted the move sizes to relatively small values at the final density:
Now that the compression is complete, the particles are much closer together and not overlapping, but still arranged randomly.
Save the final configuration to a GSD file for use in the next stage of the simulation:
hoomd.write.GSD.write(state=sim.state, mode='xb', filename='compressed.gsd')