code: def read_structure(filepath, outpath="."):
      r"""Read the structure of the polyatomic system from a file and returns
      a list with the symbols of the atoms in the molecule and a 1D array
      with their positions :math:`[x_1, y_1, z_1, x_2, y_2, z_2, \dots]` in
      atomic units (Bohr radius = 1).
  
      The atomic coordinates in the file must be in Angstroms.
      The `xyz <https://en.wikipedia.org/wiki/XYZ_file_format>`_ format is supported. Additionally,
      the new file ``structure.xyz``, containing the input geometry, is created in a directory with
      path given by ``outpath``.
  
      Args:
          filepath (str): name of the molecular structure file in the working directory
              or the absolute path to the file if it is located in a different folder
          outpath (str): path to the output directory
  
      Returns:
          tuple[list, array]: symbols of the atoms in the molecule and a 1D array with their
          positions in atomic units.
  
      **Example**
  
      >>> symbols, coordinates = read_structure('h2.xyz')
      >>> print(symbols, coordinates)
      ['H', 'H'] [0.    0.   -0.66140414    0.    0.    0.66140414]
      """
      file_in = filepath.strip()
      file_out = os.path.join(outpath, "structure.xyz")
  
      copyfile(file_in, file_out)
  
      symbols = []
      coordinates = []
      with open(file_out, encoding="utf-8") as f:
          for line in f.readlines()[2:]:
              symbol, x, y, z = line.split()
              symbols.append(symbol)
              coordinates.append(float(x))
              coordinates.append(float(y))
              coordinates.append(float(z))
  
      return symbols, np.array(coordinates) / bohr_angs
  
file_path: pennylane/pennylane/qchem/structure.py
imports: import os
  import re
  from shutil import copyfile
  import numpy
label_desc: Use image node_id 1 for calling a global function with example usage: read_structure(filepath, outpath) and returns: symbols,  unknown