# Set up elements and positions using xyz format h2o = MoleculeConfiguration(xyz_format= """3 O 0.000 0.000 0.0 H 0.757 0.586 0.0 H -0.757 0.586 0.0""")