def get_displacement_vector(molecule, start_atom, end_atom):
    """
    Get the displacement vector between `start_atom` and `end_atom`.

    """

    position1, position2 = (
        molecule.get_atomic_positions((start_atom, end_atom))
    )
    return position2 - position1