def read_xyz(file_path):
    """
    Function to read in xyz file and output as list of atoms and positions

    ARGS:
    file_path (str): path to xyz file

    RETURNS:
    list
    """

    with open(file_path, "r") as f:
        xyz_file = f.readlines()[2:]


    atom_list = []
    for line in xyz_file:
        element, pos = line.split()[0], [float(x) for x in line.split()[1:]]
        atom_list.append([element, pos])

    return atom_list