“…For the sake of numerical convenience, a Lorentzian function of the form normalΓ α ν ( ω ) = normalΓ α ν W α 2 false( ω − normalΩ α false) 2 + W α 2 is adopted, where Γ αν is the effective hybridization strength between the νth d orbital on the Ni ion and the s bands in the αth reservoir and Ω α and W α are the band center and bandwidth of the αth reservoir, respectively. The involved energetic parameters (ϵ ν , U ν , Γ αν , W α , and Ω α ) are extracted from the DFT calculations ,,− (see section S3 of the Supporting Information for details). In particular, the magnitudes of s–d hybridization strength Γ αν are obtained by analyzing the Kohn–Sham Green’s functions. , The parametrized AIM is then solved by using the Fermionic HEOM method, which is capable of addressing the impurity–reservoir hybridization, electron transport, multielectron cotunneling, and strong correlation effect in a nonperturbative manner. ,, …”