Skip to content

Commit 5fad50b

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 da68824 commit 5fad50b

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)