We should add the ability to resize the box to MD integration methods. We can use the barostat argument to do this:

• If barostat is an Extent, resize the box immediately at the start of the run.
• If barostat is a tuple(Extent, Extent), do linear interpolation between the two box sizes from the beginning to the end of the run.