Skip to content

Commit f7fce02

Browse files
committed
Send the species multiplicity to Arkane's input file
to allow it to process (and differentiate between) singlet and triplet forms of the same structure, see RMG-Py PR #2854
1 parent 79c2aca commit f7fce02

File tree

1 file changed

+4
-2
lines changed

1 file changed

+4
-2
lines changed

arc/statmech/arkane.py

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -50,7 +50,7 @@
5050
% for spc in species_list:
5151
% if spc['smiles']:
5252
species('${spc['label']}', '${spc['path']}'${spc['pdep_data'] if 'pdep_data' in spc else ''},
53-
structure=SMILES('${spc['smiles']}'))
53+
structure=SMILES('${spc['smiles']}'), spinMultiplicity=${spc['multiplicity']})
5454
% else:
5555
species('${spc['label']}', '${spc['path']}'${spc['pdep_data'] if 'pdep_data' in spc else ''})
5656
% endif
@@ -325,7 +325,9 @@ def render_arkane_input_template(self,
325325
if e0_only or spc.compute_thermo:
326326
species_list.append({'label': spc.label,
327327
'path': spc.yml_path or os.path.join(statmech_dir, 'species', f'{spc.label}.py'),
328-
'smiles': spc.mol.copy(deep=True).to_smiles() if not spc.is_ts else ''})
328+
'smiles': spc.mol.copy(deep=True).to_smiles() if not spc.is_ts else '',
329+
'multiplicity': spc.multiplicity,
330+
})
329331
ts_list = [{'label': rxn.ts_species.label,
330332
'path': rxn.ts_species.yml_path or os.path.join(statmech_dir, 'TSs', f'{rxn.ts_species.label}.py')}
331333
for rxn in self.reactions] if self.reactions else list()

0 commit comments

Comments
 (0)