@@ -1260,7 +1260,6 @@ def run_sp_job(self,
12601260 """
12611261 level = level or self .sp_level
12621262
1263- # determine_occ(xyz=self.xyz, charge=self.charge)
12641263 if level == self .opt_level and not self .composite_method \
12651264 and not (level .software == 'xtb' and self .species_dict [label ].is_ts ) \
12661265 and 'paths' in self .output [label ] and 'geo' in self .output [label ]['paths' ] \
@@ -1295,40 +1294,6 @@ def run_sp_job(self,
12951294 self .job_dict [label ]['sp' ] = dict ()
12961295 if self .composite_method :
12971296 raise SchedulerError (f'run_sp_job() was called for { label } which has a composite method level of theory' )
1298- if 'mrci' in level .method :
1299- if self .job_dict [label ]['sp' ]:
1300- # Parse orbital information from the CCSD job, then run MRCI
1301- job0 = None
1302- job_name_0 = 0
1303- for job_name , job in self .job_dict [label ]['sp' ].items ():
1304- if int (job_name .split ('_a' )[- 1 ]) > job_name_0 :
1305- job_name_0 = int (job_name .split ('_a' )[- 1 ])
1306- job0 = job
1307- with open (job0 .local_path_to_output_file , 'r' ) as f :
1308- lines = f .readlines ()
1309- core = val = 0
1310- for line in lines :
1311- if 'NUMBER OF CORE ORBITALS' in line :
1312- core = int (line .split ()[4 ])
1313- elif 'NUMBER OF VALENCE ORBITALS' in line :
1314- val = int (line .split ()[4 ])
1315- if val * core :
1316- break
1317- else :
1318- raise SchedulerError (f'Could not determine number of core and valence orbitals from CCSD '
1319- f'sp calculation for { label } ' )
1320- self .species_dict [label ].occ = val + core # the occupied orbitals are the core and valence orbitals
1321- self .run_job (label = label ,
1322- xyz = self .species_dict [label ].get_xyz (generate = False ),
1323- level_of_theory = 'ccsd/vdz' ,
1324- job_type = 'sp' )
1325- else :
1326- # MRCI was requested but no sp job ran for this species, run CCSD first
1327- logger .info (f'running a CCSD job for { label } before MRCI' )
1328- self .run_job (label = label ,
1329- xyz = self .species_dict [label ].get_xyz (generate = False ),
1330- level_of_theory = 'ccsd/vdz' ,
1331- job_type = 'sp' )
13321297 if self .job_types ['sp' ]:
13331298 if self .species_dict [label ].multi_species :
13341299 if self .output_multi_spc [self .species_dict [label ].multi_species ].get ('sp' , False ):
0 commit comments