def minimum_image_displacement(crd_0, crd_1, box_size):
    """Find displacement between nearest periodic images of atom pair"""

    displacement = crd_0 - crd_1
    displacement[displacement < -box_size / 2] += box_size
    displacement[displacement > box_size / 2] -= box_size
    return displacement