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