Skip to content

Commit

Permalink
Allow users to set a specific MRCI procedure
Browse files Browse the repository at this point in the history
e.g.: MP2_RS3_MRCI-F12
  • Loading branch information
alongd committed Aug 1, 2024
1 parent a7c1e1a commit e3fd2fc
Showing 1 changed file with 13 additions and 3 deletions.
16 changes: 13 additions & 3 deletions arc/job/adapters/molpro.py
Original file line number Diff line number Diff line change
Expand Up @@ -46,8 +46,8 @@
${cabs}
int;
{hf;${shift}
maxit,1000;
wf,spin=${spin},charge=${charge};}
maxit,1000;
wf,spin=${spin},charge=${charge};}
${restricted}${method};
Expand Down Expand Up @@ -252,13 +252,23 @@ def write_input_file(self) -> None:
if 'mrci' in self.level.method:
input_dict['orbitals'] = '\ngprint,orbitals;\n'
input_dict['restricted'] = ''
input_dict['method'] = f"""{{casscf;
if '_' in self.level.method:
methods = self.level.method.split('_')
input_dict['method'] = ''
for method in methods:
input_dict['method'] += f"""{{{method.lower()};
maxit,999;
wf,spin={input_dict['spin']},charge={input_dict['charge']};}}
"""
else:
input_dict['method'] = f"""{{casscf;
maxit,999;
wf,spin={input_dict['spin']},charge={input_dict['charge']};}}
{{mrci{"-f12" if "f12" in self.level.method.lower() else ""};
maxit,999;
wf,spin={input_dict['spin']},charge={input_dict['charge']};}}"""

input_dict['block'] += '\n\nE_mrci=energy;\nE_mrci_Davidson=energd;\n\ntable,E_mrci,E_mrci_Davidson;'

input_dict = update_input_dict_with_args(args=self.args, input_dict=input_dict)
Expand Down

0 comments on commit e3fd2fc

Please sign in to comment.