“…For the exchange-correlation functional, we employ B3LYP [31] and CAM-B3LYP [32] which have been integrated into the software as well. The former has provided good prediction in our previous studies on molecular orbital interaction problems [33], [34], [35], [36], excitation in solvent problems [37], [38], and reaction path coordinate to calculate the tunneling probability [39], [40], [41]. Meanwhile, the latter is needed to consider the correction of the long-range interaction in states 1a , 1b , and 1c .…”