Determining excited state processes for small nanoclusters, specifically gold, aids in our ability to fine-tune luminescent materials and optical devices. Using TDDFT and TDDFT + TB, we present a detailed theoretical explanation for the dual emission peaks displayed in Au₁₄Cd(S-Adm)₁₂ (Adm = adamantane). As dual emission is relatively rare, we decipher whether the mechanism originates from two different excited states or from two different minima on the same excited state surface. This unique mechanism, which proposes that the dual emission results from two minima on the first excited state, stems from geometrical changes in the bi-tetrahedron core during the emission process.

Understanding photoluminescent mechanisms has become essential for photocatalytic, biological, and electronic applications. Unfortunately, analyzing excited state potential energy surfaces (PESs) in large systems is computationally expensive, and hence limited with electronic structure methods such as time-dependent density functional theory (TDDFT). Inspired by the sTDDFT and sTDA methods, time-dependent density functional theory plus tight binding (TDDFT + TB) has been shown to reproduce linear response TDDFT results much faster than TDDFT, particularly in large nanoparticles. For photochemical processes, however, methods must go beyond the calculation of excitation energies. Herein, this work outlines an analytical approach to obtain the derivative of the vertical excitation energy in TDDFT + TB for more efficient excited state PES exploration. The gradient derivation is based on the Z vector method, which utilizes an auxiliary Lagrangian to characterize the excitation energy. The gradient is obtained when the derivatives of the Fock matrix, the coupling matrix, and the overlap matrix are all plugged into the auxiliary Lagrangian, and the Lagrange multipliers are solved. This article outlines the derivation of the analytical gradient, discusses the implementation in Amsterdam Modeling Suite, and provides proof of concept by analyzing the emission energy and optimized excited state geometry calculated by TDDFT and TDDFT + TB for small organic molecules and noble metal nanoclusters.