Photoluminescence of ligand-protected gold nanoclusters: progress in experimental and theoretical studies.

Photoluminescence of ligand-protected gold nanoclusters: progress in experimental and theoretical studies - PMC Skip to main content An official website of the United States government Here's how you know Here's how you know Official websites use .gov A .gov website belongs to an official government organization in the United States. Secure .gov websites use HTTPS A lock ( Lock Locked padlock icon ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites. Search Log in Dashboard Publications Account settings Log out Search… Search NCBI Primary site navigation Search Logged in as: Dashboard Publications Account settings Log in Search PMC Full-Text Archive Search in PMC Journal List User Guide PERMALINK Copy As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with, the contents by NLM or the National Institutes of Health. Learn more: PMC Disclaimer | PMC Copyright Notice Chem Sci . 2026 Mar 10;17(14):6907–6923. doi: 10.1039/d6sc00001k Search in PMC Search in PubMed View in NLM Catalog Add to search Photoluminescence of ligand-protected gold nanoclusters: progress in experimental and theoretical studies Kang Li Kang Li a Department of Chemistry, Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, Xiangtan University, Hunan Province 411105, China, Email: [email protected], Email: [email protected] Find articles by Kang Li a , Pu Wang Pu Wang a Department of Chemistry, Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, Xiangtan University, Hunan Province 411105, China, Email: [email protected], Email: [email protected] Find articles by Pu Wang a, ✉ , Yong Pei Yong Pei a Department of Chemistry, Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, Xiangtan University, Hunan Province 411105, China, Email: [email protected], Email: [email protected] b College of Chemistry & Chemical Engineering, Hunan Normal University, Changsha, Hunan 410081, China Find articles by Yong Pei a, b, ✉ Author information Article notes Copyright and License information a Department of Chemistry, Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, Xiangtan University, Hunan Province 411105, China, Email: [email protected], Email: [email protected] b College of Chemistry & Chemical Engineering, Hunan Normal University, Changsha, Hunan 410081, China ✉ Corresponding author. Received 2026 Jan 1; Accepted 2026 Feb 26; Collection date 2026 Apr 16. This journal is © The Royal Society of Chemistry This article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence . You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes. PMC Copyright notice PMCID: PMC12973312  PMID: 41816613 Abstract Ligand-protected gold nanoclusters, with their small sizes and atomically precise structures, exhibit unique photophysical properties that position them as promising candidates for applications in bioimaging, sensing, and optoelectronics. However, except for a few ligand-protected gold nanoclusters, their photoluminescence quantum yield (PLQY) is still relatively low, and the underlying emission mechanisms are not yet fully elucidated. This review highlights recent experimental and theoretical advances aimed at enhancing and elucidating the photoluminescence (PL) mechanisms of ligand-protected gold nanoclusters. Experimentally, strategies such as heterometallic doping, ligand engineering, and structural rigidification have been employed to promote radiative transitions and suppress nonradiative decay, resulting in substantial improvements in PLQY. Theoretically, methods including density functional theory (DFT), time-dependent density functional theory (TDDFT), and nonadiabatic molecular dynamics (NA-MD) have provided important insights into the PL origin, emission pathways, and excited-state dynamics of ligand-protected gold nanoclusters. These advances will deepen our understanding of structure–property relationships in ligand-protected gold nanoclusters and pave the way for the rational design of highly emissive gold nanoclusters. Experimental enhancement and theoretical mechanisms of photoluminescence in ligand-protected Au nanoclusters. 1. Introduction Owing to their ultrasmall size (1–3 nm) and atomically precise structures, ligand-protected gold nanoclusters exhibit unique physical and chemical properties that distinguish them from traditional gold nanoparticles (>3 nm, characterized primarily by surface plasmon resonance). In particular, they display rich photoluminescence (PL) behaviors across the visible and near-infrared (NIR) regions, which has established them as a research focus in the field of advanced functional nanomaterials. This atomic-level precision is specifically reflected in the well-defined “metal core—ligand shell” configuration of ligand-protected gold nanoclusters, which not only provides a structural basis for regulating PL performance but also renders them an ideal model for elucidating the correlation between structure and PL properties. Notably, the PL properties of ligand-protected gold nanoclusters are highly aligned with the demands of biological applications. Their excellent biocompatibility, relatively low toxicity, and persistent NIR emission capable of penetrating biological tissues render them ideal fluorescent probes for bioimaging, disease diagnostics, and therapeutic applications. 1–7 However, research on the PL of ligand-protected gold nanoclusters still faces two major challenges. First, from an experimental perspective, although systems with high photoluminescence quantum yield (PLQY) have been synthesized in recent years ( e.g. , Ir@Au 12 (PLQY = 87%) and Au 16 Cu 6 (PLQY ≈ 100%)), 8–13 the PLQY of most gold nanoclusters is still relatively low, particularly in the NIR-II region, which fails to meet the signal intensity requirements for practical applications. Second, from a mechanistic perspective, the complex electronic dynamics induced by core–ligand interactions in these gold nanoclusters make it challenging to experimentally determine the origin of PL (core-dominated vs. surface-ligand-involved). Furthermore, the strong intrinsic spin–orbit coupling (SOC) of Au atoms efficiently promotes intersystem crossing (ISC) and reverse intersystem crossing (RISC) between singlet and triplet states, 14–16 leading to fluorescence, phosphorescence, thermally activated delayed fluorescence (TADF), and the coexistence of multiple emission phenomena. 17–23 The overlapping emission peaks arising from the small singlet–triplet state energy gap (Δ E ST ) further hinder the differentiation of emission types (phosphorescence or TADF). Recently, experimental studies have made substantial progress in exploring the PL mechanism of ligand-protected gold nanoclusters via various characterization techniques, such as steady-state and time-resolved spectroscopy. 18,24,25 These efforts have provided clues to the origin of PL but cannot fully reveal the microscopic mechanisms underlying excited-state dynamics and nonradiative relaxation. First-principles theoretical calculations have become a key tool to address these bottlenecks. Specifically, they can directly overcome the limitations of experimental characterization and provide atomic-level insights for the direct exploration of electronic structures, excited-state relaxation dynamics, and energy/charge transfer processes. With the ongoing advancement of computational methods and the increasing computational efficiency and computing resources, theoretical calculations have become an indispensable partner in experiments, offering theoretical references for the modulation of PL properties in these gold nanoclusters. In this review, we summarize the latest experimental advances in enhancing the PL intensity of ligand-protected gold nanoclusters, as well as theoretical studies focused on exploring their excited-state properties. First, we summarize the experimental strategies employed to boost PL efficiency in recent years, including enhancing radiative decay ( e.g. , heterometal doping) and suppressing nonradiative pathways ( e.g. , ligand rigidification). Subsequently, we introduce the key computational methods utilized to study the excited-state properties of these gold nanoclusters, including time-dependent density functional theory (TDDFT) and nonadiabatic molecular dynamics (NA-MD), and analyze the applicability and accuracy boundaries of each method. Finally, we discuss theoretical investigations into the excited-state properties of these gold nanoclusters, including the origins of photoluminescence, SOC-driven emission pathways, and excited-state nonradiative relaxation dynamics, aiming to provide theoretical references for the precise design of ligand-protected gold nanoclusters featuring high PL performance. 2. Experimental advances in photoluminescence enhancement Research on the photoluminescence of ligand-protected gold nanoclusters faces two primary challenges, including low PLQY and unclear photoluminescence mechanisms, which severely constrain the rational design of high-performance gold nanoclusters. Over the past decade, substantial experimental efforts have been devoted to addressing these issues, resulting in a range of effective strategies for PLQY enhancement, including size control, 26–28 heterometallic doping, 8,29–31 tailored ligand engineering, 9,32–34 and aggregation-induced emission (AIE). 35–37 The availability of ligand-protected gold nanoclusters with atomically precise structures and systematic spectroscopic data has also provided an essential experimental foundation for theoretical studies. Such well-defined cluster systems serve as benchmark models for validating computational approaches. For example, in TDDFT calculations, the assessment and selection of exchange–correlation functionals still rely largely on comparisons with experimental optical spectra (see Section 3). Nevertheless, most existing reviews emphasize synthetic strategies and potential applications of gold nanoclusters with high PLQY, while discussions on how different experimental strategies enhance PLQY via photophysical modulation are fragmented across the literature. A unified photophysical framework linking experimental strategies to the key factors that determine PLQY is still lacking. In this context, the present section provides a systematic overview of recent experimental advances in improving the PLQY of ligand-protected gold nanoclusters from a photophysical perspective. By focusing on the key physical parameters governing PLQY, this section highlights major experimental breakthroughs, outlines unresolved challenges, and underscores the indispensable role of computational methods in elucidating PLQY enhancement mechanisms and guiding rational experimental design. 2.1. Photophysical basis for PLQY enhancement According to the PLQY expression, Φ = k r /( k r + k nr ), the decay of the excited states in ligand-protected gold clusters proceeds primarily through two competing pathways, i.e. , radiative and nonradiative transitions (as illustrated in Fig. 1 ). The relative magnitudes of the corresponding rate constants ultimately determine the PLQY of the system. Upon photon absorption, an electron is promoted from an occupied orbital to an unoccupied orbital with higher energy, driving the system from the ground state to an excited state. In accordance with Kasha's rule, excitation to higher-lying excited states is typically followed by rapid nonradiative relaxation processes, including internal conversion (IC), ISC, and vibrational relaxation, leading to population of the lowest excited state (S 1 or T 1 ), from which radiative decay to the ground state mainly occurs. In the radiative process, an excited-state molecule releases energy and returns to the ground state through photon emission. The efficiency of this process is quantified using Einstein's spontaneous emission rate formula: 1 where ω = E em /ℏ is the angular frequency of the emitted photon, e is the elementary charge, ε 0 is the vacuum permittivity, ℏ is the reduced Planck constant, c is the speed of light in vacuum, and µ is the transition dipole moment. Fig. 1. Schematic illustration of PL processes in ligand-protected gold nanoclusters (left) and the key factors influencing PLQY of ligand-protected gold nanoclusters (right). Open in a new tab Within the framework of Einstein's theory of spontaneous emission, both a larger emission energy E em and a larger transition dipole moment µ favor an increased radiative decay rate. Meanwhile, excited states can also return to the ground state through nonradiative pathways. The corresponding nonradiative rate constant k nr is commonly governed by the energy gap law and decreases exponentially with increasing excited–ground state energy gap, which is often correlated with a higher emission energy E em . Consequently, increasing E em not only enhances the radiative transition rate but also suppresses nonradiative decay. Because the equilibrium geometries of the ground and excited states differ, photoexcitation is followed by structural relaxation on the excited-state potential energy surface. This relaxation involves an energetic cost associated with geometric reorganization, which is commonly described by the reorganization energy λ . This relaxation leads to an emission energy E em that is lower than the vertical excitation energy E ex defined at the ground-state equilibrium geometry. Therefore, for a given electronic level structure, suppressing excited-state structural relaxation can reduce energy loss and increase E em . Beyond indirect modulation via emission energy, direct control of nonradiative decay channels is also critical. In particular, weakening nonadiabatic coupling between electronic states and reducing electron–vibrational coupling strength can effectively suppress nonradiative processes, thereby lowering k nr . In summary, enhancing the PLQY of ligand-protected gold nanoclusters can be achieved by engineering their electronic structures to increase the emission energy and transition dipole moment, while simultaneously suppressing excited-state structural relaxation and nonradiative decay through ligand engineering or structural rigidification. Such strategies enable larger radiative rates and smaller nonradiative rates, ultimately leading to improved photoluminescence efficiency. On this basis, the following sections will systematically review recent experimental efforts aimed at enhancing the PLQY of ligand-protected gold nanoclusters, with particular emphasis on strategies that promote radiative transitions and suppress nonradiative processes. 2.2. Promoting the efficiency of the radiative process In the radiative process, an excited-state molecule releases energy and returns to the ground state through photon emission. According to Einstein's theory of spontaneous emission, the radiative rate is mainly governed by two key factors: the emission energy E em and the transition dipole moment µ . Accordingly, strategies for enhancing the radiative rate can be broadly classified into two aspects: (i) increasing the emission energy of the lowest excited state through electronic structure modulation, and (ii) enhancing the transition dipole moment to accelerate radiative decay. 2.2.1. Increasing the emission energy Relative to the ground state, the excitation energy E ex , defined as the energy difference between the excited state and ground state, is determined by the energy gap between occupied and unoccupied orbitals and the electron–hole interaction following excitation. When the lowest excited state (S 1 or T 1 ) is dominated by a single HOMO → LUMO electronic transition, the excitation energy E ex shows a strong correlation with the HOMO–LUMO gap E g . In such cases, enlarging E g through electronic structure modulation is expected to increase the excitation energy of the lowest excited state and increase the emission energy E em . Recently, Buschmann et al. summarized the correlations among the HOMO–LUMO gap E g , the emission energy E em , and the PLQY for a series of Au 13 core clusters protected by different ligands. 38 Their results revealed that Au 13 clusters with larger E g simultaneously exhibit higher emission energies and higher PLQY. These observations provide empirical support for the design strategy of enhancing the emission energy via E g engineering, thereby promoting radiative decay while suppressing nonradiative relaxation. Previous studies have demonstrated that the HOMO–LUMO gap E g in ligand-protected gold nanoclusters can be effectively tuned through various experimental strategies, including cluster size control, 39–41 ligand modification, 33,42–44 and metal doping. 8,45,46 Among these strategies, metal doping is an effective strategy to enhance the luminescence efficiency by enlarging the E g , and it has been widely applied in experiments. Takano and co-workers showed that replacing the central atom of the Au 13 core with elements lighter than Au (in the order of Au > Pt > Ir in atomic weight) can widen the E g ( Fig. 2A and B ). 47 Using this principle, they successfully synthesized a series of heterometal-doped clusters, including [RuAu 12 (dppm) 6 ] 2+ , [RhAu 12 (dppm) 6 ] 3+ , and [IrAu 12 (dppm) 6 ] 3+ (dppm = 1,2-bis(diphenylphosphino)methane). Notably, [RuAu 12 (dppm) 6 ] 2+ exhibits intense room-temperature phosphorescence in degassed CH 2 Cl 2 solution, with a PLQY of 37%. 8 Using a similar doping strategy, a series of clusters with the general formula [MAu 12 (dppe) 5 Cl 2 ] (M = Au, Pd, Pt, Rh, Ir; dppe = Ph 2 PCH 2 CH 2 PPh 2 ) were also successfully prepared. 30 Experimental results showed that heterometal doping substantially enlarges the E g of the clusters ( Fig. 2B ), thereby accelerating radiative transitions while suppressing nonradiative decay pathways. As a result, the PLQY of these systems is significantly enhanced ( Fig. 2C ). Importantly, such doped cluster systems with well-defined atomic structures and systematic spectroscopic data not only validate the physical mechanism of enhancing radiative efficiency through HOMO–LUMO gap engineering, but also provide ideal experimental models and reliable benchmarks for subsequent excited-state electronic-structure calculations using TDDFT and methodological assessment. Fig. 2. (A) Two-step jellium potentials for doping lower (left) and higher (right) valent atoms into Au 13 . Reproduced with permission. 47 Copyright 2021, American Chemical Society. (B) Plot of the HOMO–LUMO gaps of M@Au 12 (dppe) 5 Cl 2 (M = Ir, Pt, Au, Rh, Pd). (C) Plot of the PLQYs of M@Au 12 (dppe) 5 Cl 2 (MAu 12 ) and M@Au 12 (dppm) 6 (MAu 12 -m) versus the energy of their PL peaks. Reproduced with permission. 30 Copyright 2022, Wiley-VCH. Open in a new tab 2.2.2. Increasing the transition dipole moment In addition to the emission energy E em , the transition dipole moment µ represents another key parameter governing the radiative transition rate. For a given emission energy, excited states with larger transition dipole moments generally exhibit a higher radiative transition rate. Recent studies on the Au 38 S 2 (S-Adm) 20 cluster (S-Adm = adamantanethiolate) have shown that this system displays a relatively high PLQY of 15%. 48 The enhanced radiative efficiency was attributed to symmetry-allowed transitions between frontier orbitals ( Fig. 3A ), which lead to an increase in the transition dipole moment. 48,49 Beyond orbital symmetry, the cluster's overall geometry, specifically its aspect ratio, exerts notable influence on the transition dipole moment. Luo and co-workers reported that a series of rod-shaped Au clusters (Au 42 , Au 60 , Au 78 , Au 96 , and Au 114 ) with hexagonal close-packed (HCP) cores exhibit unusually intense absorption peaks in their UV-vis spectra 40 ( Fig. 3B ). Theoretical calculations revealed that these absorption features originate from electronic transitions with large transition dipole moments ( Fig. 3C ). Among them, the Au 42 (PET) 32 cluster (PET = 2-phenylethanethiolate) shows a relatively high PLQY (11.9%), 50 and subsequent theoretical studies by Xie et al. confirmed that its large transition dipole moment significantly enhances the radiative decay channel. 51 Similar conclusions have also been supported by computational investigations of Au 60 clusters. 52 Similarly, in a theoretical study on the optical absorption properties of the Au 76 (SR) 44 cluster with a face-centered cubic (FCC) core, Ma et al. found that its larger aspect ratio leads to a stronger longitudinal transition dipole moment 53 ( Fig. 3D ). Although existing studies across different cluster systems have revealed several empirical correlations between geometric structure, orbital symmetry, and transition dipole moments, a comprehensive understanding of the relationship between the structure of a specific cluster and its transition dipole moment remains to be established. Fig. 3. (A) Isosurface and irreducible representation of the frontier orbitals of Au 38 S 2 (S-Adm) 20 and Au 30 (S-Adm) 18 . Reproduced with permission. 49 Copyright 2024, American Chemical Society. (B) Optical absorption spectra of Au 42 , Au 60 , Au 78 , Au 96 , and Au 114 ; (C) the HOMO and LUMO orbitals and electric transition dipole moment for transitions from the ground state to the excited states of Au 42 (SCH 3 ) 32 , Au 60 (SH) 44 , Au 78 (SH) 56 , and Au 96 (SH) 68 . Reproduced with permission. 40 Copyright 2024, National Academy of Sciences. (D) Transversal and longitudinal dipole moment transitions of Au 76 (SR) 44 . Reproduced with permission. 53 Copyright 2016, American Chemical Society. Open in a new tab The HOMO–LUMO gap E g can be indirectly obtained from experimental techniques such as optical absorption spectroscopy and can be effectively tuned through modulation of cluster size, choice of heterometal dopants, or ligand electronic effects. In contrast, the intensity of an apparent absorption peak in the experimental spectrum may arise from the superposition of multiple excited-state transitions. As a result, it is experimentally challenging to directly correlate the intensity of the longest-wavelength absorption peak with the transition dipole moment of the lowest excited state. Fundamentally, the transition dipole moment is determined by the spatial distribution and symmetry characteristics of the orbitals involved in the electronic transition, which generally require excited-state electronic structure calculations for detailed analysis. This issue becomes particularly important in the context of applications such as bioimaging, where NIR emission is desired and blue-shifting the emission ( i.e. , increasing the emission energy) is no longer an option. Under such conditions, enhancing the transition dipole moment represents a more viable approach for improving radiative efficiency. Therefore, establishing more general and design-oriented relationships between the cluster structure and transition dipole moment through theoretical calculations is essential for guiding the rational design of highly emissive gold clusters. 2.3. Suppressing nonradiative processes In addition to radiative transitions, excited-state ligand-protected gold nanoclusters can also undergo nonradiative deactivation pathways, which diminish the PLQY. Consequently, beyond enhancing radiative efficiency, suppressing nonradiative decay channels is another critical strategy for improving the PL performance of gold clusters. In fact, most reported ligand-protected gold nanoclusters still exhibit relatively low PLQY, indicating that nonradiative channels generally dominate the excited-state decay processes. Structurally, ligand-protected gold nanoclusters adopt a core–shell structure consisting of a metallic core surrounded by an organic ligand shell. 54 In such systems, various vibrational modes introduced by both the metal core and the ligands can couple with electronic excited states, providing efficient pathways for nonradiative energy dissipation. Recent studies have emphasized strategies like ligand shell rigidification and suppression of surface vibrations to mitigate nonradiative relaxation pathways and consequently enhance PLQY. 2.3.1. Ligand rigidification In ligand-protected gold nanoclusters, increasing the overall rigidity of the ligand shell represents an important strategy for suppressing nonradiative decay. To date, gold clusters protected by various types of ligands have been reported, including thiolates, 55–57 phosphines, 58 alkynyls, 59 and N-heterocyclic carbenes (NHCs). 60,61 By comparing Au 13 clusters protected by different ligand types, Buschmann et al. showed that Au 13 clusters protected by X-type ligands ( i.e. , thiolates and alkynyls) typically undergo larger structural relaxation in the excited state, whereas Au 13 clusters protected by L-type ligands ( i.e. , phosphines and NHCs) exhibit much smaller excited-state geometric changes. 38 In recent years, a series of NHC-protected gold clusters with exceptionally high PLQY have been reported. Luo et al. synthesized chiral NHC-protected Au 13 clusters with PLQY as high as 61% (ref. 62 ) ( Fig. 4A ). By combining ligand-exchange experiments with theoretical simulations, they revealed that noncovalent inter-ligand interactions, including C–H⋯N, C–H⋯π, and π⋯π interactions, play a crucial role in enhancing the overall rigidity of the ligand shell and suppressing excited-state structural relaxation. Subsequently, Kulkarni et al. reported bidentate NHC-protected Au 13 clusters with PLQY reaching 62% (ref. 10 ) ( Fig. 4B ), attributing the enhanced performance to increased ligand rigidity and strong interligand interactions. For diphosphine-protected gold nanoclusters, Mutoh and co-workers combined Ir doping with tailor-made diphosphine ligands to synthesize [IrAu 12 (bbpe) 6 ] 3+ (bbpe = bis[benzo[ b ]phosphindole]ethane), which exhibits an exceptionally high PLQY of 87%. 9 Similar ligand-engineering strategies have also been successfully applied to thiolate-protected gold nanoclusters. Enhancing the π⋯π interactions by reducing the methyl group at the end of the thiol ligand, Wang et al. improved the PLQY of the Au 52 (SR) 32 cluster from 3.8% to 18.3% (ref. 34 ) ( Fig. 4C ). And this approach was subsequently applied to the Au 36 (SR) 24 cluster, achieving a similar PLQY enhancement 63 ( Fig. 4D ). Fig. 4. (A) The chiral NHC-stabilized Au 13 cluster with high PLQY (61%). Reproduced with permission. 62 Copyright 2023, Wiley-VCH. (B) The bis-NHC-protected Au 13 clusters with π-extended aromatic systems and high PLQY. Reproduced with permission. 10 Copyright 2025, American Chemical Society. (C) The PLQY of Au 52 (SR) 32 with four different aromatic thiolate ligands. Reproduced with permission. 34 Copyright 2023, American Chemical Society. (D) Emission spectra of Au 36 (SR) 24 with four different aromatic thiolate ligands. Reproduced with permission. 63 Copyright 2024, American Chemical Society. Open in a new tab Collectively, these studies indicate that enhancing interligand interactions to rigidify the ligand shell is an effective strategy for suppressing excited-state structural relaxation and thereby reducing nonradiative decay. However, most ligand-modification approaches reported to date have focused on hydrophobic ligands. Meanwhile, hydrophilic ligands are essential for ligand-protected gold nanoclusters aimed at biological applications, 4,18,64,65 and the surrounding aqueous environment plays a decisive role in governing their PL properties. From a theoretical perspective, a faithful description of these aqueous environment effects requires going beyond implicit-solvent treatments of bulk polarity and instead adopting explicit, dynamic multiscale simulations ( e.g. , QM/MM), coupled with nonadiabatic excited-state dynamics to capture solvent-coupled structural relaxation as well as charge- and energy-transfer processes. 66–72 Notably, such multiscale and nonadiabatic approaches have been less frequently applied to ligand-protected gold nanoclusters and merit further development. 2.3.2. Suppression of staple motions Unlike phosphine and NHC ligands, which bind directly to the metallic core through anchoring atoms (P or C), thiolate ligands in thiolate-protected gold nanoclusters can form characteristic Au–S staple motifs involving a portion of surface gold atoms. The vibrational modes associated with these staple motifs can couple with electronic excited states, thereby facilitating nonradiative energy dissipation. Modulating the Au–S staple structures to restrict their vibrational relaxation has therefore been demonstrated as an effective strategy for enhancing the PLQY of thiolate-protected gold nanoclusters. 73 In 2021, Li et al. reported that the Au 38 S 2 (S-Adm) 20 cluster exhibited a relatively high PLQY of 15%. 48 They attributed this PL enhancement to the rigidifying effect of surface µ 3 -S atoms and eight bridging SR ligands ( Fig. 5A ), which significantly suppresses energy dissipation caused by Au–S vibration. Similar effects have been further validated in isomeric systems. Chen and Xia et al. identified two isomers of Au 28 (SR) 20 (Au 28i and Au 28ii ) with identical gold cores but different surface Au–S staple motif structures, leading to significant differences in the PL intensity. 74,75 In particular, the long-lived emission component of Au 28ii is significantly enhanced, indicating that the interlocked short motifs effectively inhibit nonradiative decay. Subsequently, Dong et al. synthesized the Au 28 (SCH 2 Ph- t Bu) 22 cluster, whose compact and rigid shell layer effectively reduced nonradiative energy loss and enhanced the fluorescence intensity. 76 Fig. 5. (A) Atomic structure of Au 38 S 2 (S-Adm) 20 determined from single-crystal diffraction. Reproduced with permission. 48 Copyright 2021, American Chemical Society. (B) PLQY of Au 42 dissolved in deaerated 2-MeTHF containing DMBA ( N , N -dimethylbenzamide) with different concentrations. (C) PL spectra of Au 42 dissolved in C 2 Cl 4 and embedded in a polystyrene (PS) film. Reproduced with permission. 12 Copyright 2024, American Chemical Society. (D) Schematic diagram of intramolecular cross-linking of ligands by the formation of bis-Schiff linkages on the Au 22 (SG) 18 surface. (E) PL spectra of Au 22 (SG) 18 and PDA-Au 22 (SG) 18 NCs. Reproduced with permission. 77 Copyright 2022, Springer Nature. (F) Comparison of the emission spectra of Au 16 Cu 6 and Au 22 in CH 2 Cl 2 . Reproduced with permission. 11 Copyright 2024, American Association for the Advancement of Science. Open in a new tab Au–S staple motion can also be constrained through external molecular environments. By introducing amide molecules to further restrict staple vibrations in Au 42 (PET) 32 clusters, the room-temperature PLQY was increased to 50% and even reached 75% when embedded in a polymer matrix containing amide molecules 12 ( Fig. 5B and C ). In addition to environmental confinement, surface modification engineering can achieve similar effects. Deng et al. employed bis-Schiff base cross-links to restrict the movement of Au–S staples in individual Au 22 (SG) 18 clusters (SG = glutathione), achieving a PLQY of 48% (ref. 77 ) ( Fig. 5D and E ). Similar to thiolate-protected systems, alkynyl-protected gold nanoclusters also feature Au–C Created by potrace 1.16, written by Peter Selinger 2001-2019 CR staple motifs. Recently, metal doping within staple motifs provides another effective means of suppressing nonradiative decay in alkynyl-protected gold nanoclusters. Shi et al. doped six Cu atoms into the surface staple units of Au 22 ( t BuPhC Created by potrace 1.16, written by Peter Selinger 2001-2019 C) 18 , yielding a more compact Au 16 Cu 6 ( t BuPhC Created by potrace 1.16, written by Peter Selinger 2001-2019 C) 18 cluster with a near-unity PLQY approaching 100% (ref. 11 ) ( Fig. 5F ). From a theoretical perspective, suppressing staple vibrations can be viewed as mitigating excited-state structural relaxation and vibronic coupling at the surface, thereby reducing nonradiative decay. In summary, experimental studies on the photoluminescence of ligand-protected gold nanoclusters in recent years have established a clear and reproducible set of structure–property relationships. Although the specific experimental strategies differ in their practical implementations, their underlying mechanisms can be consistently rationalized in terms of effective modulation of key parameters governing excited-state decay processes. On the one hand, radiative decay rates are enhanced by increasing the emission energy and the transition dipole moment. On the other hand, nonradiative decay is efficiently suppressed by strengthening the rigidity of ligands and surface motifs and by reducing the coupling between vibrational modes and electronic excited states. It should be emphasized that these experimental observations cannot be fully understood solely on the basis of static structural parameters or ground-state electronic structures. Instead, they critically depend on excited-state electronic structure characteristics, including energy levels, coupling parameters, and the dynamical evolution of excited-state decay processes. Therefore, to achieve a systematic mechanistic understanding and predictive capability for the PLQY of ligand-protected gold nanoclusters, it is essential to introduce theoretical and computational frameworks centered on excited states. Beyond identifying the key factors governing PLQY, theoretical methods are also indispensable for determining emission types, elucidating emission origins, and describing the dynamics of nonradiative relaxation processes. This includes not only calculations of excited-state energies and transition properties based on TDDFT, but also quantitative descriptions of SOC effects, ISC, and vibrational coupling accompanied by structural relaxation. On this basis, the following section will focus on the crucial role of excited-state theoretical approaches in addressing these issues. 3. Theoretical tools for understanding photoluminescence Although extensive experimental efforts have been devoted to enhancing the PLQY of ligand-protected gold nanoclusters and to elucidating their excited-state properties and emission pathways, 18,19,24,31,41 there is still a lack of a comprehensive understanding of their PL mechanisms, which hinders the rational design of clusters with specific emission wavelengths and lifetimes. Specifically, experimental characterization faces three limitations. Firstly, the origin of photoluminescence remains ambiguous due to the inherent core–shell structure of ligand-protected gold nanoclusters. Experimental researchers have proposed several competing explanations, including core-dominated emission, ligand participation, ligand-to-metal charge transfer (LMCT), and metal-to-ligand charge transfer (MLCT). These hypotheses remain under debate because the electronic distribution and transfer processes in the excited states of the clusters cannot be directly observed. Secondly, the strong SOC in gold atoms often results in the coexistence of multiple emission types, including fluorescence, phosphorescence, and TADF. However, the experimentally observed fluorescence and phosphorescence emission bands usually overlap owing to their broad bandwidths and small energy separations, making it challenging to unambiguously distinguish and assign the emission types. Finally, although ultrafast spectroscopy provides powerful means to probe the timescales of excited-state relaxation in these clusters, capturing the atomic-scale details of processes such as structural relaxation ( e.g. , core distortion and ligand vibration) and charge or energy transfer remains highly challenging. Theoretical calculations can complement these experimental limitations and serve as indispensable tools for elucidating PL mechanisms. In light of the challenges mentioned above, this section systematically reviews two major computational approaches employed to investigate the PL mechanisms of ligand-protected gold nanoclusters and clarifies their relevance to the unresolved experimental issues (as illustrated in Fig. 6 ). The first category involves electronic structure methods, including density functional theory (DFT), TDDFT, and multiconfigurational methods. By selecting an appropriate theoretical methods, theoretical calculations can reproduce experimental absorption and emission spectra, thereby providing a reliable foundation for subsequent rate-constant evaluations and excited-state property analyses. In particular, analytical tools such as natural transition orbital (NTO) and electron–hole analysis 78 enable direct visualization of charge distribution and transfer characteristics in excited states, helping to elucidate the origin of photoluminescence. The second category focuses on the simulation of nonradiative processes, encompassing static rate-constant calculations based on Fermi's golden rule and NA-MD methods that explicitly account for electron–nuclear coupling. By computing and comparing rate constants associated with different decay pathways, one can distinguish emission types, estimate PLQY, and unravel the underlying mechanisms. Moreover, NA-MD simulations provide time-resolved insights into electronic relaxation dynamics, revealing transition times and pathways among excited states. Together, these computational strategies establish a robust methodological framework for in-depth mechanistic analysis of photoluminescence and for the rational design of ligand-protected gold nanoclusters with enhanced PL performance. Fig. 6. Computational approaches for investigating PL mechanisms in ligand-protected gold nanoclusters. Open in a new tab 3.1. Electronic structure methods The selection of electronic structure methods is pivotal for accurately modeling the electronic configuration and deriving dependable photophysical characteristics, serving as a fundamental prerequisite for valid theoretical interpretations. DFT and TDDFT are widely used to analyze the electronic structure and excited-state properties of ligand-protected gold nanoclusters, owing to their favorable balance between computational efficiency and accuracy. 43,55,79,80 However, the reliability of theoretical predictions, such as UV-vis absorption spectra, emission energies, and excited-state lifetimes, strongly depends on the choice of the exchange–correlation functional. 80–83 Generalized gradient approximation (GGA) functionals, such as BP86 (ref. 84 ) and PBE, 85 are frequently selected in computational studies of ligand-protected gold nanoclusters due to their computational efficiency and acceptable accuracy. However, they tend to underestimate excitation energies, thus limiting their applicability in excited-state properties. Recently, in the computational investigation of X@Cu 14 (X = Au, Cl) clusters, Luo et al. assessed the performance of three GGA/meta-GGA functionals (BP86, PBE, and TPSS 86,87 ) and nine hybrid functionals (PBE0, 88 B3LYP, 89,90 revPBE38, 91 TPSSh, 92 TPSS0, 93 M06, M06-2X, 94 CAM-B3LYP, 95 and ωB97X-D4 (ref. 96 )) for the calculation of fluorescence and phosphorescence emission energies. 97 For Au@Cu 14 , hybrid functionals with high Hartree–Fock exchange ( e.g. , revPBE38) most accurately replicate the experimental T 1 → S 0 phosphorescence emission energy ( Fig. 7 ), while for Cl@Cu 14 , GGA/meta-GGA functionals (BP86, PBE, and TPSS) align more closely with the experiment. Notably, the calculated Δ E (S 1 –T 1 ) increased with higher HF exchange fractions, resulting in a growing deviation between the calculated fluorescence and phosphorescence emission properties and experimental data. Additionally, the oscillator strengths of fluorescence radiation show an increasing trend with rising HF exchange percentage in the functional. For phosphorescence radiation, except for M06-2X, other functionals showed minimal influence on the oscillator strengths. Qian and co-workers evaluated the performance of various functionals for ligand-protected noble-metal clusters, finding that the proportion of HF exchange predominantly influences simulated spectral features, regardless of the hybrid functional type. 98 Benchmarking against experimental spectra indicates that functionals incorporating approximately 10–20% HF exchange yield the most accurate spectral predictions. Furthermore, She et al. examined the fluorescence mechanisms and doping effects in heteroatom-doped M@Au 12 (dppe) 6 x −2 clusters, with M representing a 4d or 5d transition metal, and x denoting the nominal charge. They benchmarked several exchange–correlation functionals by comparing the calculated optical absorption and fluorescence spectra of [Au 13 (dppe) 5 Cl 2 ] 3+ with experimental data. 79 They found that the PBE functional markedly underestimates emission energies, while HSE06 (ref. 99 ) offers the best consistency with experimental observations. Taken together, these studies indicate that no universal exchange–correlation functional has yet been established for ligand-protected gold nanoclusters. The optimal choice remains system-dependent and typically requires benchmarking against experimental spectra or higher-level calculations. Fig. 7. Calculated T 1 → S 0 emission energies of Au@Cu 14 using functionals with different fractions of Hartree–Fock exchange. Open in a new tab In theoretical calculations, electronic excitation properties (such as local and charge-transfer excitations) can be characterized using NTO or hole–electron analysis. Specifically, such analyses are performed to calculate the excited-state properties of clusters via DFT and TDDFT methods combined with optimal exchange–correlation functionals, to further obtain the orbital compositions and transition contributions of the excited-state electronic structure, and ultimately to quantitatively determine the respective contributions of the metal core and ligand shell in the excited-state transitions, thereby elucidating the origin of photoluminescence. For instance, Weerawardene et al. utilized TDDFT to study Au 25 (SR) 18 − clusters and revealed that the emission primarily originates from core-based transitions ( i.e. , excitations out of superatom P orbitals into the lowest two superatom D orbitals). 100 While DFT and TDDFT provide an advantageous compromise between computational efficiency and accuracy, their computational cost escalates rapidly as the size of gold nanoclusters ( e.g. , beyond Au 50 ) increases, limiting the feasibility of simulating larger systems. Concurrently, the electronic and steric influences of protecting ligands markedly affect the photophysical characteristics of gold nanoclusters. 9,34,79 This presents a methodological challenge: simplifying ligand structures to reduce computational demand introduces systematic discrepancies, whereas employing realistic ligand models substantially elevates computational resource requirements. In their study of the dual-emission mechanism of the Au 14 Cd(S-Adm) 12 cluster, 101 Havenridge and co-workers observed that the absorption spectra are highly sensitive to ligands. Simplifying ligands, such as substituting Adm groups with hydrogen or methyl groups, may lead to deviations from the authentic photophysical responses of the actual system. They thus applied a combined approach of TDDFT and tight-binding approximation (TDDFT + TB) 102 to simulate absorption spectra of Au 14 Cd(S-Adm) 12 protected by full ligands. This approach yielded computational efficiency comparable to that of TDDFT while significantly reducing computational effort. Subsequently, they developed an innovative approximation method, termed TDDFT-aas (time-dependent density functional theory approximate auxiliary s-function), which does not depend on tight-binding parameters. 103 Within higher computational efficiency, the calculated absorption spectra of silver and gold nanoparticles using the TDDFT-aas method are in good agreement with the results of TDDFT and TDDFT + TB. Furthermore, the implementation of analytical excited-state gradients within the TDDFT-aas framework facilitates efficient calculation of emission energies for larger nanoclusters. It is worth noting that TDDFT may struggle to accurately describe potential energy surfaces near conical intersections or state crossings, where multi-configurational methods ( e.g. , CASSCF) become essential. In 2016, Wu et al. investigated the photoluminescence mechanisms of ultrasmall ligand-protected gold nanoclusters using the CASPT2//CASSCF approach, which enabled highly accurate predictions of absorption and emission wavelengths in excellent agreement with experimental data. 104 For example, when applied to the cluster [Au 7 (dppp) 4 ] 3+ (dppp = 1,3-bis(diphenylphosphino)propane), the calculated phosphorescence emission energy was 644 nm, matching the experimental value (642 nm) with a deviation of only 2 nm. However, the high computational cost and the need for careful active-space selection limit the applicability of multi-configurational methods to larger or more complex systems. 105 3.2. Simulation of nonradiative processes The PLQY of ligand-protected gold nanoclusters is determined by the competition between radiative and nonradiative decay pathways, serving as a critical parameter for evaluating their photoluminescence efficiency. However, experimental approaches often face challenges in directly resolving the microscopic mechanisms of nonradiative relaxation processes (such as excited-state structural relaxation and electron-nuclear coupling interactions) and in quantitatively assessing their influence on PLQY. Consequently, theoretical calculations become essential for providing deeper mechanistic insight. To theoretically determine the PLQY value, the primary goal is to calculate the radiative rate constant k r and the nonradiative rate constant k nr , 106–108 which in turn allows for a clear determination of the reasons behind low PLQY. Notably, early theoretical studies of ligand-protected gold nanoclusters focused primarily on their absorption and emission spectra, oscillator strengths, emission energies, and radiative rate constants. The radiative rate constant k r is predominantly determined using Einstein's spontaneous emission rate expression ( eqn (1) ), which requires prior calculation of the emission energy and transition dipole moment through methods such as TDDFT. In contrast, modeling nonradiative processes remains a significant challenge in theoretical chemistry and is critical for accurately predicting the PLQY of ligand-protected gold nanoclusters. According to the energy gap law, 109 the IC rate constant k IC increases exponentially as the emission energy E em decreases, which can be expressed as k IC = A e − αE em . However, since vibronic coupling phenomena between electronic states are not explicitly incorporated, this empirical expression provides only a rough estimate and has limited applicability. Additionally, the strong relativistic effect in gold atoms induces significant SOC, which facilitates ISC from singlet to triplet states. Consequently, ligand-protected gold nanoclusters can exhibit phosphorescence, TADF, or dual emission 17–23,50 involving triplet states and the associated ISC and RISC processes. Therefore, a thorough analysis of both ISC and RISC pathways is essential for understanding the luminescence properties and mechanisms of these nanoclusters. Typically, Δ E ST and spin–orbit coupling matrix elements (SOCMEs) are utilized to estimate the competition between ISC and RISC via Fermi's golden rule, thereby helping to identify the emission type of ligand-protected gold nanoclusters. Furthermore, direct evaluation of ISC rate constant k ISC and RISC rate constant k RISC , and their comparison with k r and IC rate constant k IC , enables comprehensive mapping of excited-state decay pathways and facilitates quantitative prediction of excited-state lifetimes and PLQY. 106 However, nonradiative processes, such as IC, ISC, and RISC, involve intricate couplings among electronic and nuclear degrees of freedom, 110,111 making accurate theoretical treatment particularly challenging. According to Fermi's golden rule, the nonradiative transition rate constant k nr governing transitions between two electronic states can be expressed as follows: 2 where 〈 Ψ f | Ĥ ′| Ψ i 〉 represents the nonadiabatic coupling (NAC) Hamiltonian (for IC) or the SOC Hamiltonian (for ISC) between the initial state i and the final state f, ℏ is the reduced Planck constant, and ρ (Δ E ) is the vibrational density of states at the energy gap Δ E . Recent advances by Shuai, 107,112 Marian, 113,114 Santoro, 115 Valiev 116 and others have established general rate-constant formulations that incorporate vibrational overlap and the Duschinsky rotation effect, which have been implemented in computational packages for evaluating nonradiative rate constants. For instance, Xie et al. 51 applied Lin's formalism 117,118 to compute the ISC and RISC rate constants (6.87 × 10 7 s −1 and 2.87 × 10 2 s −1 , respectively). By combining these with the fluorescence and phosphorescence radiative rate constants (1.53 × 10 8 s −1 and 8.52 × 10 3 s −1 , respectively), they provided a theoretical explanation of the dual-emission behavior observed in the Au 42 (PET) 32 cluster. Recently, Li et al. computed nonradiative rate constants for Au 38 S 2 (SR) 20 and Au 30 (SR) 18 with the FCclasses3 program. 49 Based on the computed rate constants and the law of mass action, exciton population evolution diagrams were subsequently constructed to analyze the excited-state decay pathways and emission characteristics of these clusters. In addition to rate-constant calculations, analyses of the Huang–Rhys factors and reorganization energies, which were decomposed into contributions from individual vibrational normal modes, provided further insight into how cluster structure influences nonradiative decay processes. However, the excited-state lifetimes and PLQYs predicted by the current methods tend to be significantly lower than experimental values, primarily due to the neglect of anharmonic vibrational effects. 108,116,119,120 Moreover, these Fermi's golden rule-based approaches require the computation of excited-state vibrational frequencies, which becomes computationally prohibitive for larger clusters. As a practical alternative to Fermi's golden rule-based approaches, Marcus theory has been employed to estimate the ISC and RISC rate constants: 121–124 3 where Δ E ST is the adiabatic energy difference between the S 1 and T n states. k B is the Boltzmann constant and T is the thermodynamic temperature. V SOC is the spin–orbit coupling matrix element (SOCME), and λ is the reorganization energy. Beyond the static excited-state decay models based on Fermi's golden rule discussed above, nonadiabatic molecular dynamics (NA-MD) has emerged as a powerful tool for elucidating the ultrafast excited-state relaxation pathways of ligand-protected gold nanoclusters. 125,126 NA-MD simulations can provide detailed information on the electronic relaxation dynamics, including state-to-state transition times and mechanisms. For instance, NA-MD simulations of the Au 25 (SR) 18 − cluster revealed that the experimentally observed picosecond-scale processes likely originate from core-to-core transitions rather than the previously proposed core-to-half-ring transitions. 102 However, most current NA-MD studies of ligand-protected gold nanoclusters neglect SOC effects, thereby omitting ISC between singlet and triplet states. For ligand-protected gold nanoclusters, reproducing ISC processes that accompany dynamical structural transformations within NA-MD remains highly challenging. This difficulty primarily arises from the need to compute triplet-state energy levels and SOCMEs along nuclear trajectories. Additionally, earlier nonadiabatic simulations of nanoclusters primarily relied on the independent-particle approximation, in which excited states are represented as single-particle excitations constructed from Kohn–Sham orbitals. Within this framework, the state energies and nonadiabatic couplings (NACs) are evaluated using single Slater determinants, thereby neglecting critical carrier–carrier interactions. Recently, Xie and co-workers 127 compared simulation results obtained from both single-particle and many-body (MB) treatments in their NA-MD study of Au 25 (SR) 18 − . Their findings demonstrated that the inclusion of MB effects significantly enhanced the NAC between excited states, resulting in ultrafast IC from higher excited states on a sub-picosecond timescale, which is consistent with previous experimental observations. Future developments of NA-MD simulation methods in ligand-protected gold nanoclusters need to incorporate relativistic effects, ISC processes, and carrier–carrier interactions to more accurately reproduce the excited-state dynamics of ligand-protected gold nanoclusters. 71,72,127–129 4. Mechanistic insights into photoluminescence and excited‐state dynamics Section 3 has already systematically introduced the applicability logic of core methods such as DFT, TDDFT, and NA-MD in the study of ligand-protected gold nanoclusters. Therefore, this section will adopt a “Computational Method → Mechanism Elucidation” framework to systematically summarize the current theoretical research status concerning the photoluminescence origin, excited-state relaxation dynamics, and emission mechanisms, using case studies from typical cluster systems. This summary will explicitly define the core role of different computational methods in elucidating these mechanisms, thereby providing a valuable reference for subsequent theoretical studies and rational structure design. 4.1. Origin of photoluminescence The core–shell structure of ligand-protected gold nanoclusters, consisting of a metal core and a ligand shell (as illustrated in Fig. 8 ), gives rise to diverse PL origins, including core-dominated emission, ligand-involved emission, and charge-transfer processes such as LMCT and MLCT. In some systems, structural relaxation and nonradiative processes (IC or ISC) can reshape the character of the excited state, enabling conversion between initially populated charge-transfer states and a core-dominated excitation state before emission. Furthermore, the intricate interactions between the metal core and ligand shell lead to complex excited-state properties and electronic dynamics, which make experimental identification of the PL origin highly challenging. TDDFT, combined with NTO analysis or electron–hole analysis, has thus emerged as a powerful tool for clarifying excited-state electronic behavior. This method serves as an important tool for elucidating the dominant PL origin and understanding the structural factors that modulate emission by enabling direct comparison among different cluster systems. At present, research on the PL origin in ligand-protected gold nanoclusters can be broadly classified into two major categories. Fig. 8. The photoluminescence origin diagram of ligand-protected gold nanoclusters. Open in a new tab (1) Core-dominated excitation: for the majority of ligand-protected gold nanoclusters, PL originates primarily from core-localized electronic excitation. In this case, the emission is governed by electronic transitions within the metal core, whereas the ligands mainly modulate the orbital energy levels through electron-donating or -withdrawing effects and influence core structural relaxation via steric hindrance. This conclusion has been consistently validated across multiple systems using TDDFT calculations. For example, early transient absorption spectroscopy (TAS) studies on the Au 25 (SR) 18 − cluster suggested that its PL might arise from ligand–metal, charge-transfer, or semi-ring states. 130 To resolve this controversy, Weerawardene et al. employed TDDFT to calculate the excited-state electronic structures of Au 25 (SR) 18 − clusters clearly modified with different R groups (R = H, CH 3 , C 2 H 5 , C 3 H 7 ). 100 Their results clearly demonstrated that the emission originates predominantly from core-based transitions rather than charge-transfer or semi-ring states ( Fig. 9A ), thereby providing a theoretical resolution to the long-standing debate over this system's emission origin. Subsequent TDDFT calculations on systems such as Au 38 (SR) 24 , 131 Au@Cu 14 , 97 Au 22 (C Created by potrace 1.16, written by Peter Selinger 2001-2019 CR) 18 , 132 Au 42 (PET) 32 ( Fig. 9B ), 52 and Au 38 S 2 (S-Adm) 20 (ref. 49 ) further confirmed the dominance of core-localized excitation. Furthermore, Weerawardene et al. 100 illustrated the ligand-dependent modulation of emission energy for Au 25 (SH) 18 − , Au 25 (SCH 3 ) 18 − , and Au 25 (SCH 2 CH 2 CH 3 ) 18 − , and the fluorescence emission energies were calculated to be 0.829 eV, 0.655 eV, and 0.512 eV, respectively. Fig. 9. (A) Kohn–Sham orbitals and orbital energies for the S 1 state of Au 25 (SH) 18 − . Reproduced with permission. 94 Copyright 2016, American Chemical Society. (B) Distributions of hole and electron pairs during the electronic transitions of Au 42 and Au 60 (blue regions represent holes and green regions represent electrons). Reproduced with permission. 53 Copyright 2025, American Chemical Society. (C) Distribution of hole and electron pairs during the T 1 → S 0 transition at the optimized T 1 geometry of Au@Cu 14 and Cl@Cu 14 . Reproduced with permission. 91 Copyright 2025, American Chemical Society. Open in a new tab (2) Charge transfer: in contrast, for certain Au( i ) clusters and specific ligand-protected systems, the ligands participate directly in electronic transitions, giving rise to LMCT. For instance, studies on the PL properties of Au( i ) and Cu( i ) clusters protected by NHC ligands, 133 which exhibit highly efficient TADF, revealed through NTO analysis that the S 1 state in all clusters displayed LMCT characteristics. Similarly, in the Cl-doped Cl@Cu 14 cluster, the T 1 state exhibits charge-transfer features ( Fig. 9C ), 97 whereas the PL origin in the Au@Cu 14 system remains core-localized excitation. Furthermore, He et al. demonstrated that Au doping in AuAg 12 and Au 3 Ag 10 induces significant charge-transfer characteristics. 134 4.2. SOC and emission pathways Ligand-protected gold nanoclusters, which comprise several tens of gold atoms or more, exhibit pronounced SOC owing to the heavy-atom effect of gold. Serving as a pivotal mediato, SOC accelerates ISC between singlet and triplet states and, by regulating the RISC rate, directly dictates the emission pathway, namely fluorescence, phosphorescence, or TADF ( Fig. 10 left). Early studies often overlooked SOC and the associated ISC/RISC processes, which obscured the interpretation of complex photoluminescence behaviors. Here, we summarize advances in elucidating emission types of ligand-protected gold nanoclusters based on quantitative rate-constant calculations and analyses, with particular emphasis on how SOC is incorporated into the models. By computing ISC and RISC rate constants, we clarify the determinative roles of SOC and the Δ E ST in selecting the final emission channel, thereby establishing a “Structural Factors → SOC Effects → Emission Type” regulatory framework ( Fig. 10 right). Fig. 10. Jablonski diagram of fluorescence, TADF, and phosphorescence emission processes (left), and factors influencing ISC and RISC processes (right). Open in a new tab 4.2.1. Influence of SOC and Δ E ST on emission types The SOC effect operates chiefly by opening the phosphorescent radiative pathway and by setting the emission outcome through competition between ISC and RISC. This competition is governed by two key parameters: the Δ E ST and the SOCMEs. A smaller Δ E ST facilitates RISC and favors TADF, whereas a larger Δ E ST suppresses RISC and favors phosphorescence. Larger SOCMEs enhance ISC efficiency and increase the phosphorescent radiative rate. The interplay between Δ E ST and SOCMEs can therefore yield multiple emission types and even dual-emission behavior. This mechanism should be validated with multi-state model simulations and quantitative rate calculations. (1) Differences of emission in isomers: in 2020, the groups of Jin and Wu independently reported two isomeric Au 28 (SR) 20 clusters ( i.e. , Au 28i (SR) 20 and Au 28ii (SR) 20 ). 74,75 Although both possess an identical Au 14 core, differences in their ligand-staple motifs led to markedly different PLQYs. An earlier DFT/TDDFT study by Li et al. 135 proposed an anti-Kasha S 2 fluorescence mechanism but did not incorporate SOC effects or triplet states. More recently, Mazumder et al. combined temperature and oxygen quenching experiments with calculations of ISC and RISC rate constants, 20 revealing the emission type and PL mechanism of two Au 28 (SR) 20 isomers ( Fig. 11A ). Au 28i (SR) 20 exhibited a small Δ E ST (0.126 eV), resulting in comparable T 1 → S 1 RISC and S 1 → T 1 ISC rate constants (2.90 × 10 10 s −1 and 3.53 × 10 10 s −1 , respectively), thus showing both phosphorescence and TADF. In contrast, Au 28ii (SR) 20 has a larger Δ E ST (0.240 eV), which suppresses RISC (9.87 × 10 3 s −1 ) and leads exclusively to phosphorescence via the ISC pathway. Furthermore, rate constant calculations also indicated that the higher-lying T 2 state provides an indirect and efficient nonradiative channel from the excited singlet states to the triplet states, reinforcing the conclusion that SOC directly regulates the emission type. Fig. 11. Simulated mechanism of photoluminescence for (A) Au 28i (SR) 20 and Au 28ii (SR) 20 . Reproduced with permission. 20 Copyright 2024, American Chemical Society. (B) Au 24 , Au 42 , and Au 60 clusters. Reproduced with permission. 52 Copyright 2025, American Chemical Society. (C) Au 5 Ag 11 and Pt 1 Ag 16 . Reproduced with permission. 136 Copyright 2022, Wiley-VCH. (D) Au@Cu 14 and Cl@Cu 14 . Reproduced with permission. 97 Copyright 2025, American Chemical Society. Open in a new tab (2) Three-state modeling of dual emission mechanism: to rationalize the dual emission of Au 42 (PET) 32 , Xie et al. established a three-state kinetic model (S 1 , T 1, and S 0 ) and simulated the ISC and RISC processes. 51 Because S 1 and T 1 share similar localized electronic character, the SOCME between S 1 and T 1 is relatively weak (46.49 cm −1 ), consistent with the El-Sayed rule. In addition, the S 1 state has a large oscillator strength, yielding a high fluorescence radiative rate constant (1.53 × 10 8 s −1 ) that is comparable to the relatively modest S 1 → T 1 ISC rate constant (6.87 × 10 7 s −1 ). Consequently, a fraction of excitons return directly to S 0 via fluorescence, while the remainder populate T 1 through S 1 → T 1 ISC. From T 1 , the subsequent phosphorescence rate (8.52 × 10 3 s −1 ) substantially exceeds the RISC rate (2.87 × 10 2 s −1 ), producing phosphorescence. Overall, the competition among these sequential pathways accounts for the coexistence of fluorescence and phosphorescence in Au 42 (PET) 32 . 4.2.2. Modulation of SOC and Δ E ST Structural factors within a nanocluster, such as central-atom doping, ligand modification, and alloying, tune the spin–orbit landscape by modulating the SOCMEs and the Δ E ST ( Fig. 10 right). They shift the ISC/RISC balance and enable control over the emission channel, including fluorescence, phosphorescence, and TADF. (1) Heavy-atom effect of core atom doping: a comparison of calculations between Au@Cu 14 and Cl@Cu 14 clearly illustrates the mechanism of heavy-atom regulation effect ( Fig. 11D ). 97 Owing to the incorporated Au, Au@Cu 14 shows a substantially larger SOCME value (2.17 cm −1 ) compared to that of Cl@Cu 14 (0.77 cm −1 ). This stronger SOC, together with a smaller Δ E ST (S 1 –T 2 ) of 0.03 eV, results in a faster S 1 → T 2 ISC rate (1.58 × 10 10 s −1 ). As a result, Au@Cu 14 exhibits strong phosphorescence (PLQY = 71.3%) via the S 1 → T 2 → T 1 → S 0 pathway. In contrast, the weaker SOC in Cl@Cu 14 leads to competition between ISC and RISC among S 1 , T 1 , and T 2 , producing the coexistence of TADF and phosphorescence. This comparison clarifies the enhancing role of heavy-atom doping in strengthening SOC. (2) Alloying effects: to investigate how doping concentration influences emission, He et al. performed theoretical calculations on Ag 13 , Au 13 , and the model alloy clusters AuAg 12 and Au 3 Ag 10 , all of which share the identical icosahedral core structure and superatom electronic configuration (1S 2 1P 6 ). 134 For Au 13 , which exhibits strong NIR emission, the large T 2 –T 1 energy gap (0.55 eV) allows T 2 phosphorescence to compete effectively with the T 2 → T 1 IC process, indicating that its emission primarily originates from the T 2 → S 0 transition. In contrast to the localized S 1 and T 1 states found in Au 13 and Ag 13 , the T 1 state in AuAg 12 and Au 3 Ag 10 exhibits pronounced charge-transfer character induced by Au doping ( Fig. 12 ). This modification enhances the SOCMEs (23.97 cm −1 for AuAg 12 and 45.55 cm −1 for Au 3 Ag 10 ) and reduces the Δ E ST relative to Ag 13 (about 0.01 eV for AuAg 12 and 0.03 eV for Au 3 Ag 10 ). These effects accelerate ISC and increase the overall PLQY. Fig. 12. Energy diagrams and images of the hole and electron pairs of (A) Ag 13 , (B) Au 1 Ag 12 , (C) Au 3 Ag 10 and (D) Au 13 . The S r index is defined as the full space integration of a function ( S r ( r )) describing the overlap between electron and hole distributions, and the D index is the distance between a hole and an the electron center of mass. Reproduced with permission. 134 Copyright 2024, Oxford University Press on behalf of the Chinese Academy of Sciences. Open in a new tab Investigation into the analogous dopant effects of various metals on photoluminescence has also been conducted for Au 5 Ag 11 (fluorescence) and Pt 1 Ag 16 (phosphorescence). 136 The computational analysis revealed the photoluminescence mechanisms of Au 5 Ag 11 and PtAg 16 ( Fig. 11C ). For Au 5 Ag 11 , only the T 1 state lies below the S 1 state, with a relatively larger Δ E ST between S 1 and T 1 (0.034 eV) than that of Pt 1 Ag 16 (the Δ E ST between S 1 and T 2 is about 0.002 eV), resulting in a small ISC rate constant that favors radiative emission as fluorescence. In contrast, in Pt 1 Ag 16 , the S 1 and T 2 states are nearly degenerate (Δ E ≈ 0.002 eV), which promotes an efficient ISC from S 1 to T 2 . This facilitates an efficient S 1 → T 2 → T 1 → S 0 decay pathway, ultimately leading to phosphorescence. 136 (3) Size effects: in a theoretical study of one-dimensional ultrathin rod-like ligand-protected gold nanoclusters (Au 24 , Au 42 , and Au 60 ), Luo and co-workers combined rate-constant calculations with multistate model analysis ( Fig. 11B ). 52 They showed that the size indirectly controls the emission type by modulating the SOC effect. The smaller Au 24 cluster exhibits dual fluorescence, attributable to substantial excited-state structural distortion that creates two minima (S 1 -MIN1 and S 1 -MIN2) on the S 1 potential energy surface. These minima yield large Δ E ST values (0.5 eV and 0.81 eV, respectively), which render ISC negligible. By contrast, for the longer rod-like clusters Au 42 and Au 60 , the large transition dipole moments lead to fluorescence radiative rates (1.13 × 10 8 s −1 and 1.71 × 10 8 s −1 , respectively) comparable to ISC rates (1.23 × 10 8 s −1 and 8.68 × 10 8 s −1 , respectively). Together with larger Δ E ST values (0.39 eV and 0.57 eV) that suppress RISC, the result is dual emission comprising fluorescence and phosphorescence. (4) Ligand substituent effects: rate-constant calculations were performed on chiral Au( i )Cu( i ) clusters protected by NHC ligands (R/S-py-X and R/S-ql-X, X = Cl, Br, I) to investigate the effect of the halide substituent ( Fig. 13 ). 133 The results show that as the halogen atomic mass increases, the Δ E ST (S 1 –T 1 ) narrows systematically ( Fig. 13A ): the R-ql-Cl cluster (0.294 eV), R-ql-Br cluster (0.193 eV), and R-ql-I cluster (0.116 eV). This narrowing of the Δ E ST (S 1 –T 1 ) corresponds to an increase in the RISC rate constant: R-ql-I (7.05 × 10 10 s −1 ) > R-ql-Br (4.22 × 10 9 s −1 ) > R-ql-Cl (5.49 ×10 8 s −1 ). Consequently, the R-ql-I cluster exhibits the most efficient TADF owing to its minimal Δ E ST (S 1 –T 1 ). Fig. 13. (A) Energy diagram of R-ql-Cl, R-ql-Br, and R-ql-I indicating TADF and phosphorescence emission processes. (B) The hole and electron pairs for S 1 /S 0 and T 1 /S 0 transitions were obtained by NTO analysis at the optimized S 1 and T 1 geometries of R-ql-Cl, R-ql-Br, and R-ql-I (isovalue of 0.02). Au, yellow; Cu, brown; N, dark blue; C, gray; Cl, blue; Br, green; I, light purple; H, white. Reproduced with permission. 133 Copyright 2023, Springer Nature. Open in a new tab In summary, recent studies establish that SOC, ISC, and RISC are pivotal in determining the photoluminescence of metal nanoclusters. The interaction between singlet and triplet states, which is quantified by SOCMEs and the Δ E ST , governs the competition among fluorescence, phosphorescence, and TADF, thereby shaping the quantum yield and lifetime. Structural factors, including cluster architecture, heteroatom doping, ligand identity and rigidity, and heavy-atom contributions, tune SOCMEs and Δ E ST and also modulate nonradiative channels, producing outcomes ranging from single-channel emission to dual emission. These insights underscore the need for multistate, rate-constant-based models that explicitly incorporate triplet states to explain and predict the diverse luminescent behaviors of metal nanoclusters. 4.3. Excited-state relaxation and nonradiative dynamics Excited-state relaxation dynamics directly determine the cluster's excited-state lifetime and PLQY. These dynamics comprise IC and ISC among excited states, structural relaxation, charge transfer, and electron–hole recombination. Because such processes are not directly observable at the atomic scale, combining NA-MD with DFT/TDDFT enables atomistic elucidation of relaxation pathways, timescales, and nonradiative channels. Current research has successfully revealed how the ligand, metal core structure, and doping effects influence these excited-state relaxation dynamics. By modifying the energy gap and orbital interactions, ligands regulate relaxation rates. TDDFT and NA-MD on Au 25 (SR) 18 − showed that internal relaxation among high-lying singlets (S 6 → S 1 ) occurs on a sub-picosecond timescale, whereas the final S 1 → S 0 decay is about two orders of magnitude slower. 126 Comparing Au 25 (SR) 18 − with different R groups (R = CH 3 , C 2 H 5 , C 3 H 7 , MPA, PET) (MPA = mercaptopropionic acid) revealed S 1 –S 6 relaxation times of 2–18 ps, with MPA producing distinctive dynamics due to a narrower LUMO+1–LUMO+2 gap. 137 For the phosphine-ligated Au 13 (dppe) 5 Cl 2 3+ , which shares the Au 13 core, NA-MD indicated a slower excited-state decay (0.5–20.3 ps) than in thiolate-protected Au 25 (SR) 18 − , consistent with its larger HOMO–LUMO gap. 138 TDDFT and NA-MD on Au 38 (SR) 24 further identified low-frequency vibrations, such as the 24.2 cm −1 bending mode, as key channels for electron–phonon coupling and nonradiative decay, 139 a behavior also observed in Au 18 (SH) 14 . 140 Notably, NA-MD simulations on Au 25 (SR) 18 − by Xie et al. revealed that low-frequency modes, particularly the breathing and stretching modes of Au(core)–Au(core) and Au(core)–Au(ring), are strongly coupled to the electronic transitions. 127 More recent NA-MD work shows that ligands can also modulate carrier lifetimes by altering SOC. 129 For example, thiolate-protected Au 25 (SR) 18 − relaxed in 99 fs, compared with 132 fs for the diphosphine-protected Au 13 3+ cluster. 138 Studies of Au 30 (SR) 18 , Au 36 (SR) 24 , Au 38 S 2 (SR) 20 , and Au 44 (SR) 28 demonstrate that core stacking adjusts excited-state relaxation. 141 Au 30 (SR) 18 with a HCP core exhibits stronger electron–phonon coupling than Au 38 S 2 (SR) 20 with a body-centered cubic (BCC) core, resulting in a shorter lifetime. Additionally, nonadiabatic relaxation also depends on dopants. 142 Yu et al. combined DFT with time-domain NA-MD to examine single-atom doping in MAu 24 (SR) 18 (M = Pd, Pt, Cd, Hg), showing that dopants co-regulate relaxation through changes in the gap, phonon spectrum, and electron–phonon coupling strength. 143 In diphosphine-ligated M@Au 12 clusters, Ir doping enhances core–shell coupling and suppresses recombination, thereby extending the lifetime and confirming the dynamical impact of doping. 144 Overall, real-time NA-MD based on DFT and TDDFT has become indispensable for resolving the excited-state dynamics of ligand-protected gold nanoclusters. These simulations show that larger electronic gaps lengthen lifetimes, stronger electron–phonon coupling accelerates nonradiative decay, and low-frequency vibrations critically mediate energy dissipation. Moreover, subtle structural variations, including isomerism 145 and single-atom doping, can markedly reshape relaxation behavior and photophysical properties. 5. Conclusions and outlook In this review, we summarize recent experimental advances in enhancing the PLQY of ligand-protected gold nanoclusters, with emphasis on strategies that promote radiative decay and suppress nonradiative pathways. Experimentally, increasing the radiative rate constant via HOMO–LUMO gap engineering, achieved through heterometal doping to widen the HOMO–LUMO gap and through enhancement of transition dipole moments, has proven highly effective. Suppressing nonradiative relaxation by rigidifying and constraining ligand motion in core–shell structures has likewise yielded substantial gains in quantum efficiency. Nevertheless, a deeper mechanistic understanding of how transition dipole moments are enhanced is still needed to further improve radiative efficiency. In addition, most ligand-protected gold nanoclusters are soluble only in organic solvents, limiting biomedical applications; consequently, simple and general strategies to convert organo-soluble clusters into water-soluble forms are highly desirable. 64 Theoretically, DFT/TDDFT methods and NA-MD have clarified photoluminescence origins, complex excited-state dynamics, and emission pathways—fluorescence, phosphorescence, TADF, and dual emission—governed by Δ E ST , SOCMEs, and structural relaxation. Even so, improving the computational efficiency of excited-state optimizations and vibrational analyses remains essential, particularly for large clusters. 146,147 Computing radiative and nonradiative rate constants has deepened the understanding of competing decay channels, yet accurately predicting PLQYs is still challenging. Looking ahead, multiscale modeling that incorporates environmental effects such as solvent and matrix interactions will provide more realistic descriptions. 66–70 Moreover, machine-learning approaches trained on high-level quantum-chemical data promise to accelerate the discovery and screening of ligands, dopants, and structural motifs with desirable photophysical properties. 148–154 Author contributions All authors contributed to writing and editing this article. Conflicts of interest There are no conflicts to declare. Acknowledgments The authors acknowledge financial support from NSFC (22373082), the science and technology innovation program of Hunan Province (2023RC1055), the project of innovation team of the ministry of education (IRT_17R90), and the Scientific Research Fund of Hunan Provincial Education Department (22B0154). Biographies Biography Kang Li. Open in a new tab Kang Li is a PhD candidate at Xiangtan University under the supervision of Prof. Yong Pei. His current research interest is focused on the photoluminescence properties of ligand-protected gold nanoclusters. Biography Pu Wang. Open in a new tab Pu Wang received her PhD from Xiangtan University under the supervision of Prof. Yong Pei, and she joined the chemistry faculty of Xiangtan University in 2018. Her current research interests are focused on the structure, properties, and catalytic applications of gold nanoclusters. Biography Yong Pei. Open in a new tab Yong Pei received his bachelor's degree from Xiangtan University in 2001 and his PhD from Nanjing University in 2006. He pursued postdoctoral research at the University of Nebraska-Lincoln from 2006 to 2010. He is now a Professor at Xiangtan University and Hunan Normal University, and his research interests focus on theoretical studies of the structure and optical, catalytic, and electronic properties of coinage metal nanoclusters. Data availability No primary research results, software or code have been included and no new data were generated or analysed as part of this review. References Zheng K. Xie J. Acc. Mater. Res. 2021;2:1104–1116. doi: 10.1021/accountsmr.1c00186. [ DOI ] [ Google Scholar ] Guo Z. Yu G. Zhang Z. Han Y. Guan G. Yang W. Han M. Y. Adv. Mater. 2023:e2206700. doi: 10.1002/adma.202206700. [ DOI ] [ PubMed ] [ Google Scholar ] Huang Y. Chen K. Liu L. Ma H. Zhang X. Tan K. Li Y. Liu Y. Liu C. Wang H. Zhang X.-D. Small. 2023;19:2300145. doi: 10.1002/smll.202300145. [ DOI ] [ PubMed ] [ Google Scholar ] Baghdasaryan A. Dai H. Chem. Rev. 2025;125:5195–5227. doi: 10.1021/acs.chemrev.4c00835. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Yang J. Yang F. Zhang C. He X. Jin R. ACS Mater. Lett. 2022;4:1279–1296. doi: 10.1021/acsmaterialslett.2c00237. [ DOI ] [ Google Scholar ] Yang Z. Yang X. Guo Y. Kawasaki H. ACS Appl. Bio Mater. 2023;6:4504–4517. doi: 10.1021/acsabm.3c00518. [ DOI ] [ PubMed ] [ Google Scholar ] Sharma N. Mohammad W. Le Guével X. Shanavas A. Chem. Biomed. Imaging. 2024;2:462–480. doi: 10.1021/cbmi.4c00021. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Takano S. Hirai H. Nakashima T. Iwasa T. Taketsugu T. Tsukuda T. J. Am. Chem. Soc. 2021;143:10560–10564. doi: 10.1021/jacs.1c05019. [ DOI ] [ PubMed ] [ Google Scholar ] Mutoh K. Yahagi T. Takano S. Kawakita S. Iwasa T. Taketsugu T. Tsukuda T. Nakashima T. Chem. Sci. 2025;16:8240–8246. doi: 10.1039/D5SC00561B. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Kulkarni V. K. Albright E. L. Zeinizade E. Steele E. Chen J. Ding L. Malola S. Takano S. Harrington K. Kwon N. Levchenko T. I. Nambo M. Tsukuda T. Häkkinen H. Stamplecoskie K. Zheng G. Crudden C. M. J. Am. Chem. Soc. 2025;147:4017–4025. doi: 10.1021/jacs.4c12072. [ DOI ] [ PubMed ] [ Google Scholar ] Shi W.-Q. Zeng L. He R.-L. Han X.-S. Guan Z.-J. Zhou M. Wang Q.-M. Science. 2024;383:326–330. doi: 10.1126/science.adk6628. [ DOI ] [ PubMed ] [ Google Scholar ] Luo L. Liu Z. Mazumder A. Jin R. J. Am. Chem. Soc. 2024;146:27993–27997. doi: 10.1021/jacs.4c11703. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Wang X. Zhong Y. Li T. Wang K. Dong W. Lu M. Zhang Y. Wu Z. Tang A. Bai X. Nat. Commun. 2025;16:587. doi: 10.1038/s41467-025-55975-y. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Yang J. S. Zhao Y. J. Li X. M. Dong X. Y. Si Y. B. Xiao L. Y. Hu J. H. Yu Z. Zang S. Q. Angew Chem. Int. Ed. Engl. 2024;63:e202318030. doi: 10.1002/anie.202318030. [ DOI ] [ PubMed ] [ Google Scholar ] Yoshida K. Arima D. Mitsui M. J. Phys. Chem. Lett. 2023;14:10967–10973. doi: 10.1021/acs.jpclett.3c02977. [ DOI ] [ PubMed ] [ Google Scholar ] Zhang W. Luo L. Liu Z. Zhao F. Kong J. Jin R. Luo Y. Zhou M. Sci. Adv. 2025;11:eadx2781. doi: 10.1126/sciadv.adx2781. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Liu Z. Luo L. Kong J. Kahng E. Zhou M. Jin R. Nanoscale. 2024;16:7419–7426. doi: 10.1039/D4NR00677A. [ DOI ] [ PubMed ] [ Google Scholar ] Liu Z. Luo L. Jin R. Adv. Mater. 2024;36:e2309073. doi: 10.1002/adma.202309073. [ DOI ] [ PubMed ] [ Google Scholar ] Luo L. Liu Z. Du X. Jin R. Commun. Chem. 2023;6:22. doi: 10.1038/s42004-023-00819-3. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Mazumder A. Li K. Liu Z. Wang Y. Pei Y. Peteanu L. A. Jin R. ACS Nano. 2024;18:21534–21543. doi: 10.1021/acsnano.4c06702. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Zhu C. Chen Z.-L. Li H. Lu L. Kang X. Xuan J. Zhu M. J. Am. Chem. Soc. 2024;146:23212–23220. doi: 10.1021/jacs.4c05530. [ DOI ] [ PubMed ] [ Google Scholar ] Si W.-D. Zhang C. Zhou M. Tian W.-D. Wang Z. Hu Q. Song K.-P. Feng L. Huang X.-Q. Gao Z.-Y. Tung C.-H. Sun D. Sci. Adv. 2023;9:eadg3587. doi: 10.1126/sciadv.adg3587. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Si W. D. Zhang C. Zhou M. Wang Z. Feng L. Tung C. H. Sun D. Sci. Adv. 2024;10:eadm6928. doi: 10.1126/sciadv.adm6928. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Lin H. Song X. Chai O. J. H. Yao Q. Yang H. Xie J. Adv. Mater. 2024:e2401002. doi: 10.1002/adma.202401002. [ DOI ] [ PubMed ] [ Google Scholar ] Li Q. Zeman C. J. t. Kalkan B. Kirschbaum K. Gianopoulos C. G. Parakh A. Doan D. Lee A. C. Kulikowski J. Schatz G. C. Shen G. Kunz M. Gu X. W. Nano Lett. 2023;23:132–139. doi: 10.1021/acs.nanolett.2c03759. [ DOI ] [ PubMed ] [ Google Scholar ] Zhang L. Fan D. Shi Y. He S. Cui M. Yu H. Zhu M. RSC Adv. 2023;13:16034–16038. doi: 10.1039/D3RA01606D. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Zhou Y. Liao L. Zhuang S. Zhao Y. Gan Z. Gu W. Li J. Deng H. Xia N. Wu Z. Angew. Chem., Int. Ed. 2021;60:8668–8672. doi: 10.1002/anie.202016692. [ DOI ] [ PubMed ] [ Google Scholar ] Wang Y. Gianopoulos C. G. Liu Z. Kirschbaum K. Alfonso D. Kauffman D. R. Jin R. JACS Au. 2024;4:1928–1934. doi: 10.1021/jacsau.4c00152. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Suyama M. Takano S. Tsukuda T. J. Phys. Chem. C. 2020;124:23923–23929. doi: 10.1021/acs.jpcc.0c06765. [ DOI ] [ Google Scholar ] Hirai H. Takano S. Nakashima T. Iwasa T. Taketsugu T. Tsukuda T. Angew Chem. Int. Ed. Engl. 2022;61:e202207290. doi: 10.1002/anie.202207290. [ DOI ] [ PubMed ] [ Google Scholar ] Liu Z. Zhou M. Luo L. Wang Y. Kahng E. Jin R. J. Am. Chem. Soc. 2023;145:19969–19981. doi: 10.1021/jacs.3c06543. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Zhang B. Chen J. Cao Y. Chai O. J. H. Xie J. Small. 2021;17:e2004381. doi: 10.1002/smll.202004381. [ DOI ] [ PubMed ] [ Google Scholar ] Kim A. Zeng C. Zhou M. Jin R. Part. Part. Syst. Charact. 2017;34:1600388. doi: 10.1002/ppsc.201600388. [ DOI ] [ Google Scholar ] Wang Y. Liu Z. Mazumder A. Gianopoulos C. G. Kirschbaum K. Peteanu L. A. Jin R. J. Am. Chem. Soc. 2023;145:26328–26338. doi: 10.1021/jacs.3c09846. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Kim S. Kim H. Lee C. Park I. Kim Y. Moon D. Shim J. H. Ryu S. Park S. S. ACS Nano. 2024;18:29036–29044. doi: 10.1021/acsnano.4c10096. [ DOI ] [ PubMed ] [ Google Scholar ] Chang H. Karan N. S. Shin K. Bootharaju M. S. Nah S. Chae S. I. Baek W. Lee S. Kim J. Son Y. J. Kang T. Ko G. Kwon S.-H. Hyeon T. J. Am. Chem. Soc. 2020;143:326–334. doi: 10.1021/jacs.0c10907. [ DOI ] [ PubMed ] [ Google Scholar ] Luo Z. Yuan X. Yu Y. Zhang Q. Leong D. T. Lee J. Y. Xie J. J. Am. Chem. Soc. 2012;134:16662–16670. doi: 10.1021/ja306199p. [ DOI ] [ PubMed ] [ Google Scholar ] Buschmann D. A. Hirai H. Tsukuda T. Inorg. Chem. Front. 2024;11:6694–6710. doi: 10.1039/D4QI01773K. [ DOI ] [ Google Scholar ] Du X. Higaki T. Ma H. Zhang X. Wang G. Wang H. Jin R. Nanoscale. 2025;17:14408–14414. doi: 10.1039/D5NR01244A. [ DOI ] [ PubMed ] [ Google Scholar ] Luo L. Liu Z. Kong J. Gianopoulos C. G. Coburn I. Kirschbaum K. Zhou M. Jin R. Proc. Natl. Acad. Sci. U. S. A. 2024;121:e2318537121. doi: 10.1073/pnas.2318537121. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Zhou M. Higaki T. Li Y. Zeng C. Li Q. Sfeir M. Y. Jin R. J. Am. Chem. Soc. 2019;141:19754–19764. doi: 10.1021/jacs.9b09066. [ DOI ] [ PubMed ] [ Google Scholar ] Wei J. Halet J. F. Kahlal S. Saillard J. Y. Munoz-Castro A. Inorg. Chem. 2020;59:15240–15249. doi: 10.1021/acs.inorgchem.0c02219. [ DOI ] [ PubMed ] [ Google Scholar ] Fihey A. Maurel F. Perrier A. J. Phys. Chem. C. 2014;118:4444–4453. doi: 10.1021/jp410695v. [ DOI ] [ PubMed ] [ Google Scholar ] Song T. Kong J. Tang S. Cai X. Liu X. Zhou M. Xu W. W. Ding W. Zhu Y. Nano Res. 2023;16:19754–19764. [ Google Scholar ] Hirai H. Takano S. Nakamura T. Tsukuda T. Inorg. Chem. 2020;59:17889–17895. doi: 10.1021/acs.inorgchem.0c00879. [ DOI ] [ PubMed ] [ Google Scholar ] Wang S. Meng X. Das A. Li T. Song Y. Cao T. Zhu X. Zhu M. Jin R. Angew. Chem., Int. Ed. 2014;53:2376–2380. doi: 10.1002/anie.201307480. [ DOI ] [ PubMed ] [ Google Scholar ] Takano S. Tsukuda T. J. Am. Chem. Soc. 2021;143:1683–1698. doi: 10.1021/jacs.0c11465. [ DOI ] [ PubMed ] [ Google Scholar ] Li Q. Zeman C. J. Schatz G. C. Gu X. W. ACS Nano. 2021;15:16095–16105. doi: 10.1021/acsnano.1c04759. [ DOI ] [ PubMed ] [ Google Scholar ] Li K. Wang P. Pei Y. J. Phys. Chem. Lett. 2024;15:9216–9225. doi: 10.1021/acs.jpclett.4c02246. [ DOI ] [ PubMed ] [ Google Scholar ] Luo L. Liu Z. Du X. Jin R. J. Am. Chem. Soc. 2022;144:19243–19247. doi: 10.1021/jacs.2c09107. [ DOI ] [ PubMed ] [ Google Scholar ] Xie X. Y. Cheng K. Q. Chen W. K. Li W. Li Q. Han J. Fang W. H. Cui G. J. Phys. Chem. Lett. 2023;14:10025–10031. doi: 10.1021/acs.jpclett.3c02683. [ DOI ] [ PubMed ] [ Google Scholar ] Luo Y. Li K. Wang P. Pei Y. JACS Au. 2025;5:4593–4603. doi: 10.1021/jacsau.5c00898. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Ma Z. Wang P. Zhou G. Tang J. Li H. Pei Y. J. Phys. Chem. C. 2016;120:13739–13748. doi: 10.1021/acs.jpcc.6b04212. [ DOI ] [ Google Scholar ] Li B. Kang L. Lun Y. Yu J. Song S. Wang Y. Carbon Energy. 2024;6:e547. doi: 10.1002/cey2.547. [ DOI ] [ Google Scholar ] Pei Y. Wang P. Ma Z. Xiong L. Acc. Chem. Res. 2019;52:23–33. doi: 10.1021/acs.accounts.8b00385. [ DOI ] [ PubMed ] [ Google Scholar ] Hossain S. Hirayama D. Ikeda A. Ishimi M. Funaki S. Samanta A. Kawawaki T. Negishi Y. Aggregate. 2022;4:e255. doi: 10.1002/agt2.255. [ DOI ] [ Google Scholar ] Negishi Y. Hashimoto S. Ebina A. Hamada K. Hossain S. Kawawaki T. Nanoscale. 2020;12:8017–8039. doi: 10.1039/D0NR00824A. [ DOI ] [ PubMed ] [ Google Scholar ] Adnan R. H. Madridejos J. M. L. Alotabi A. S. Metha G. F. Andersson G. G. Adv. Sci. 2022;9:e2105692. doi: 10.1002/advs.202105692. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Lei Z. Wan X.-K. Yuan S.-F. Guan Z.-J. Wang Q.-M. Acc. Chem. Res. 2018;51:2465–2474. doi: 10.1021/acs.accounts.8b00359. [ DOI ] [ PubMed ] [ Google Scholar ] Albright E. L. Levchenko T. I. Kulkarni V. K. Sullivan A. I. DeJesus J. F. Malola S. Takano S. Nambo M. Stamplecoskie K. Hakkinen H. Tsukuda T. Crudden C. M. J. Am. Chem. Soc. 2024;146:5759–5780. doi: 10.1021/jacs.3c11031. [ DOI ] [ PubMed ] [ Google Scholar ] Narouz M. R. Osten K. M. Unsworth P. J. Man R. W. Y. Salorinne K. Takano S. Tomihara R. Kaappa S. Malola S. Dinh C. T. Padmos J. D. Ayoo K. Garrett P. J. Nambo M. Horton J. H. Sargent E. H. Hakkinen H. Tsukuda T. Crudden C. M. Nat. Chem. 2019;11:419–425. doi: 10.1038/s41557-019-0246-5. [ DOI ] [ PubMed ] [ Google Scholar ] Luo P. Zhai X. J. Bai S. Si Y. B. Dong X. Y. Han Y. F. Zang S. Q. Angew. Chem., Int. Ed. 2023;62:e202219017. doi: 10.1002/anie.202219017. [ DOI ] [ PubMed ] [ Google Scholar ] Eyyakkandy N. N. Afreen A. Vilangappurath G. Gratious S. Adarsh K. V. Mandal S. J. Phys. Chem. C. 2024;128:18828–18835. doi: 10.1021/acs.jpcc.4c04927. [ DOI ] [ Google Scholar ] Sullivan A. I. Steele E. A. Takano S. Zeinizade E. Chen J. Malola S. Siddhant K. Häkkinen H. Stamplecoskie K. G. Tsukuda T. Zheng G. Crudden C. M. J. Am. Chem. Soc. 2025;147:4230–4238. doi: 10.1021/jacs.4c14240. [ DOI ] [ PubMed ] [ Google Scholar ] Liu Z. Wang Y. Ji W. Ma X. Gianopoulos C. G. Calderon S. Ma T. Luo L. Mazumder A. Kirschbaum K. Dickey E. C. Peteanu L. A. Alfonso D. Jin R. ACS Nano. 2025;19:9121–9131. doi: 10.1021/acsnano.4c18197. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Eskelinen T. Karttunen A. J. Inorg. Chem. 2025;64:9150–9162. doi: 10.1021/acs.inorgchem.5c00761. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Gao Y.-J. Chen W.-K. Wang Z.-R. Fang W.-H. Cui G. Phys. Chem. Chem. Phys. 2018;20:24955–24967. doi: 10.1039/C8CP03657H. [ DOI ] [ PubMed ] [ Google Scholar ] Wang Q. Gao Y.-J. Zhang T.-T. Han J. Cui G. RSC Adv. 2019;9:20786–20795. doi: 10.1039/C9RA02256B. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Liu J. He X. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2023;13:e1650. [ Google Scholar ] Li H. Wang T. Han J. Xu Y. Kang X. Li X. Zhu M. Nat. Commun. 2024;15:5351. doi: 10.1038/s41467-024-49735-7. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Cui G. Thiel W. J. Chem. Phys. 2014;141:124101. doi: 10.1063/1.4894849. [ DOI ] [ PubMed ] [ Google Scholar ] Fang Y.-G. Peng L.-Y. Liu X.-Y. Fang W.-H. Cui G. Comput. Theor. Chem. 2019;1155:90–100. doi: 10.1016/j.comptc.2019.03.025. [ DOI ] [ Google Scholar ] Pyo K. Thanthirige V. D. Kwak K. Pandurangan P. Ramakrishna G. Lee D. J. Am. Chem. Soc. 2015;137:8244–8250. doi: 10.1021/jacs.5b04210. [ DOI ] [ PubMed ] [ Google Scholar ] Chen Y. Zhou M. Li Q. Gronlund H. Jin R. Chem. Sci. 2020;11:8176–8183. doi: 10.1039/D0SC01270J. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Xia N. Yuan J. Liao L. Zhang W. Li J. Deng H. Yang J. Wu Z. J. Am. Chem. Soc. 2020;142:12140–12145. doi: 10.1021/jacs.0c02117. [ DOI ] [ PubMed ] [ Google Scholar ] Dong J. Gan Z. Gu W. You Q. Zhao Y. Zha J. Li J. Deng H. Yan N. Wu Z. Angew Chem. Int. Ed. Engl. 2021;60:17932–17936. doi: 10.1002/anie.202105530. [ DOI ] [ PubMed ] [ Google Scholar ] Deng H. Huang K. Xiu L. Sun W. Yao Q. Fang X. Huang X. Noreldeen H. A. A. Peng H. Xie J. Chen W. Nat. Commun. 2022;13:3381. doi: 10.1038/s41467-022-30760-3. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Liu Z. Lu T. Chen Q. Carbon. 2020;165:461–467. doi: 10.1016/j.carbon.2020.05.023. [ DOI ] [ Google Scholar ] She J. Pei W. Zhou S. Zhao J. J. Phys. Chem. Lett. 2022;13:5873–5880. doi: 10.1021/acs.jpclett.2c01522. [ DOI ] [ PubMed ] [ Google Scholar ] Day P. N. Pachter R. Nguyen K. A. Jin R. J. Phys. Chem. A. 2019;123:6472–6481. doi: 10.1021/acs.jpca.9b02434. [ DOI ] [ PubMed ] [ Google Scholar ] Muniz-Miranda F. Menziani M. C. Pedone A. J. Phys. Chem. A. 2015;119:5088–5098. doi: 10.1021/jp507679f. [ DOI ] [ PubMed ] [ Google Scholar ] Baek H. Moon J. Kim J. J. Phys. Chem. A. 2017;121:2410–2419. doi: 10.1021/acs.jpca.6b11868. [ DOI ] [ PubMed ] [ Google Scholar ] Muniz-Miranda F. Menziani M. C. Pedone A. J. Phys. Chem. C. 2014;118:7532–7544. doi: 10.1021/jp411483x. [ DOI ] [ Google Scholar ] Perdew J. P. Phys. Rev. B:Condens. Matter Mater. Phys. 1986;33:8822–8824. doi: 10.1103/PhysRevB.33.8822. [ DOI ] [ PubMed ] [ Google Scholar ] Perdew J. P. Burke K. Ernzerhof M. Phys. Rev. Lett. 1996;77:3865–3868. doi: 10.1103/PhysRevLett.77.3865. [ DOI ] [ PubMed ] [ Google Scholar ] Tao J. Perdew J. P. Staroverov V. N. Scuseria G. E. Phys. Rev. Lett. 2003;91:146401. doi: 10.1103/PhysRevLett.91.146401. [ DOI ] [ PubMed ] [ Google Scholar ] Perdew J. P. Tao J. Staroverov V. N. Scuseria G. E. J. Chem. Phys. 2004;120:6898–6911. doi: 10.1063/1.1665298. [ DOI ] [ PubMed ] [ Google Scholar ] Adamo C. Barone V. J. Chem. Phys. 1999;110:6158–6170. doi: 10.1063/1.478522. [ DOI ] [ Google Scholar ] Becke A. D. J. Chem. Phys. 1992;96:2155–2160. doi: 10.1063/1.462066. [ DOI ] [ Google Scholar ] Stephens P. J. Devlin F. J. Chabalowski C. F. Frisch M. J. J. Phys. Chem. 1994;98:11623–11627. doi: 10.1021/j100096a001. [ DOI ] [ Google Scholar ] Zhang Y. Yang W. Phys. Rev. Lett. 1998;80:890. doi: 10.1103/PhysRevLett.80.890. [ DOI ] [ Google Scholar ] Staroverov V. N. Scuseria G. E. Tao J. Perdew J. P. J. Chem. Phys. 2003;119:12129–12137. doi: 10.1063/1.1626543. [ DOI ] [ Google Scholar ] Grimme S. J. Phys. Chem. A. 2005;109:3067–3077. doi: 10.1021/jp050036j. [ DOI ] [ PubMed ] [ Google Scholar ] Zhao Y. Truhlar D. G. Theor. Chem. Acc. 2008;120:215–241. [ Google Scholar ] Yanai T. Tew D. P. Handy N. C. Chem. Phys. Lett. 2004;393:51–57. doi: 10.1016/j.cplett.2004.06.011. [ DOI ] [ Google Scholar ] Jacquemin D. Perpète E. A. Ciofini I. Adamo C. Theor. Chem. Acc. 2011;128:127–136. [ Google Scholar ] Luo Y. Wang P. Pei Y. J. Phys. Chem. Lett. 2025;16:3705–3714. doi: 10.1021/acs.jpclett.5c00395. [ DOI ] [ PubMed ] [ Google Scholar ] Qian L. Gao Y. Xu C. Cheng L. J. Phys. Chem. Lett. 2025;16:8683–8689. doi: 10.1021/acs.jpclett.5c02044. [ DOI ] [ PubMed ] [ Google Scholar ] Krukau A. V. Vydrov O. A. Izmaylov A. F. Scuseria G. E. J. Chem. Phys. 2006;125:224106. doi: 10.1063/1.2404663. [ DOI ] [ PubMed ] [ Google Scholar ] Weerawardene K. L. D. M. Aikens C. M. J. Am. Chem. Soc. 2016;138:11202–11210. doi: 10.1021/jacs.6b05293. [ DOI ] [ PubMed ] [ Google Scholar ] Havenridge S. Aikens C. M. J. Chem. Phys. 2021;155:074302. doi: 10.1063/5.0057079. [ DOI ] [ PubMed ] [ Google Scholar ] Rüger R. van Lenthe E. Heine T. Visscher L. J. Chem. Phys. 2016;144:184103. doi: 10.1063/1.4948647. [ DOI ] [ PubMed ] [ Google Scholar ] Wang Y. Havenridge S. Aikens C. M. J. Chem. Phys. 2024;161:024101. doi: 10.1063/5.0213587. [ DOI ] [ PubMed ] [ Google Scholar ] Wu L. Fang W. Chen X. Phys. Chem. Chem. Phys. 2016;18:17320–17325. doi: 10.1039/C6CP02770A. [ DOI ] [ PubMed ] [ Google Scholar ] Veryazov V. Malmqvist P. Å. Roos B. O. Int. J. Quantum Chem. 2011;111:3329–3338. doi: 10.1002/qua.23068. [ DOI ] [ Google Scholar ] Veys K. Bousquet M. H. E. Jacquemin D. Escudero D. J. Chem. Theory Comput. 2023;19:9344–9357. doi: 10.1021/acs.jctc.3c00931. [ DOI ] [ PubMed ] [ Google Scholar ] Shuai Z. Chin. J. Chem. 2020;38:1223–1232. doi: 10.1002/cjoc.202000226. [ DOI ] [ Google Scholar ] Humeniuk A. Buzancic M. Hoche J. Cerezo J. Mitric R. Santoro F. Bonacic-Koutecky V. J. Chem. Phys. 2020;152:054107. doi: 10.1063/1.5143212. [ DOI ] [ PubMed ] [ Google Scholar ] Englman R. Jortner J. Mol. Phys. 1970;18:145–164. doi: 10.1080/00268977000100171. [ DOI ] [ Google Scholar ] Hagai M. Inai N. Yasuda T. Fujimoto K. J. Yanai T. Sci. Adv. 2024;10:eadk3219. doi: 10.1126/sciadv.adk3219. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Penfold T. J. Gindensperger E. Daniel C. Marian C. M. Chem. Rev. 2018;118:6975–7025. doi: 10.1021/acs.chemrev.7b00617. [ DOI ] [ PubMed ] [ Google Scholar ] Peng Q. Niu Y. Shi Q. Gao X. Shuai Z. J. Chem. Theory Comput. 2013;9:1132–1143. doi: 10.1021/ct300798t. [ DOI ] [ PubMed ] [ Google Scholar ] Marian C. M. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2011;2:187–203. [ Google Scholar ] Etinski M. Tatchen J. Marian C. M. J. Chem. Phys. 2011;134:154105. doi: 10.1063/1.3575582. [ DOI ] [ PubMed ] [ Google Scholar ] Cerezo J. Santoro F. J. Comput. Chem. 2023;44:626–643. doi: 10.1002/jcc.27027. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Valiev R. R. Merzlikin B. S. Nasibullin R. T. Kurtzevitch A. Cherepanov V. N. Ramazanov R. R. Sundholm D. Kurten T. Phys. Chem. Chem. Phys. 2023;25:6406–6415. doi: 10.1039/D2CP05275J. [ DOI ] [ PubMed ] [ Google Scholar ] Lin S. H. Xiao W. Z. Dietz W. Phys. Rev. E. 1993;47:3698–3706. doi: 10.1103/PhysRevE.47.3698. [ DOI ] [ PubMed ] [ Google Scholar ] Lin S. H. J. Chem. Phys. 1966;44:3759–3767. doi: 10.1063/1.1726531. [ DOI ] [ Google Scholar ] Mulvihill C. R. Georgievskii Y. Klippenstein S. J. J. Chem. Phys. 2023;159:174104. doi: 10.1063/5.0168612. [ DOI ] [ PubMed ] [ Google Scholar ] Wang Y. Ren J. Shuai Z. J. Chem. Phys. 2021;154:214109. doi: 10.1063/5.0052804. [ DOI ] [ PubMed ] [ Google Scholar ] Ousmanou M. B. Yohana J. Z. Tchangnwa Nya F. Malloum A. Conradie J. Mater. Sci. Semicond. Process. 2025;185:108894. doi: 10.1016/j.mssp.2024.108894. [ DOI ] [ Google Scholar ] Marcus R. A. Rev. Mod. Phys. 1993;65:599–610. [ Google Scholar ] Beljonne D. Shuai Z. Pourtois G. Bredas J. L. J. Phys. Chem. A. 2001;105:3899–3907. doi: 10.1021/jp010187w. [ DOI ] [ Google Scholar ] Köhler A. Bässler H. Mater. Sci. Eng., R. 2009;66:71–109. doi: 10.1016/j.mser.2009.09.001. [ DOI ] [ Google Scholar ] Chen X. Prezhdo O. V. Ma Z. Hou T. Guo Z. Li Y. Phys. Status Solidi B. 2016;253:458–462. doi: 10.1002/pssb.201552405. [ DOI ] [ Google Scholar ] Senanayake R. D. Akimov A. V. Aikens C. M. J. Phys. Chem. C. 2017;121:10653–10662. doi: 10.1021/acs.jpcc.6b09731. [ DOI ] [ Google Scholar ] Xie J. Li W. Pei Y. J. Phys. Chem. Lett. 2025;16:11259–11266. doi: 10.1021/acs.jpclett.5c02926. [ DOI ] [ PubMed ] [ Google Scholar ] Liu X. Y. Chen W. K. Fang W. H. Cui G. J. Chem. Theory Comput. 2023;19:8491–8522. doi: 10.1021/acs.jctc.3c00960. [ DOI ] [ PubMed ] [ Google Scholar ] Zhou Y. Zhou J. Yu X. Pei W. Zhou S. Zhao J. J. Phys. Chem. Lett. 2025;16:8613–8620. doi: 10.1021/acs.jpclett.5c01928. [ DOI ] [ PubMed ] [ Google Scholar ] Miller S. A. Womick J. M. Parker J. F. Murray R. W. Moran A. M. J. Phys. Chem. C. 2009;113:9440–9444. doi: 10.1021/jp9025046. [ DOI ] [ Google Scholar ] Weerawardene K. L. D. M. Guidez E. B. Aikens C. M. J. Phys. Chem. C. 2017;121:15416–15423. doi: 10.1021/acs.jpcc.7b01958. [ DOI ] [ Google Scholar ] Havenridge S. Aikens C. M. J. Phys. Chem. A. 2023;127:9932–9943. doi: 10.1021/acs.jpca.3c04644. [ DOI ] [ PubMed ] [ Google Scholar ] Ma X.-H. Li J. Luo P. Hu J.-H. Han Z. Dong X.-Y. Xie G. Zang S.-Q. Nat. Commun. 2023;14:4121. doi: 10.1038/s41467-023-39802-w. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] He W.-M. Hu J.-H. Cui Y.-J. Li J. Si Y.-B. Wang S.-B. Zhao Y.-J. Zhou Z. Ma L.-F. Zang S.-Q. Natl. Sci. Rev. 2024;11:nwae174. doi: 10.1093/nsr/nwae174. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Li J. Wang P. Pei Y. J. Phys. Chem. Lett. 2022;13:3718–3725. doi: 10.1021/acs.jpclett.2c00539. [ DOI ] [ PubMed ] [ Google Scholar ] Zhu C. Xin J. Li J. Li H. Kang X. Pei Y. Zhu M. Angew. Chem. 2022;134:e202205947. doi: 10.1002/ange.202205947. [ DOI ] [ PubMed ] [ Google Scholar ] Senanayake R. D. Aikens C. M. Phys. Chem. Chem. Phys. 2020;22:5272–5285. doi: 10.1039/c9cp04039k. [ DOI ] [ PubMed ] [ Google Scholar ] Weerawardene K. L. D. M. Pandeya P. Zhou M. Chen Y. Jin R. Aikens C. M. J. Am. Chem. Soc. 2019;141:18715–18726. doi: 10.1021/jacs.9b07626. [ DOI ] [ PubMed ] [ Google Scholar ] Senanayake R. D. Guidez E. B. Neukirch A. J. Prezhdo O. V. Aikens C. M. J. Phys. Chem. C. 2018;122:16380–16388. doi: 10.1021/acs.jpcc.8b03595. [ DOI ] [ Google Scholar ] Senanayake R. D. Aikens C. M. J. Chem. Phys. 2019;151:094702. doi: 10.1063/1.5116902. [ DOI ] [ PubMed ] [ Google Scholar ] Yu X. Pei W. Xu W. W. Zhao Y. Su Y. Zhao J. Inorg. Chem. 2023;62:20450–20457. doi: 10.1021/acs.inorgchem.3c03482. [ DOI ] [ PubMed ] [ Google Scholar ] Pandeya P. Senanayake R. D. Aikens C. M. J. Chem. Phys. 2021;154:184303. doi: 10.1063/5.0045590. [ DOI ] [ PubMed ] [ Google Scholar ] Yu X. Sun Y. Xu W.-w. Fan J. Gao J. Jiang X. Su Y. Zhao J. Nanoscale Horiz. 2022;7:1192–1200. doi: 10.1039/D2NH00281G. [ DOI ] [ PubMed ] [ Google Scholar ] Pei W. Hou L. Yang J. Zhou S. Zhao J. Nanoscale. 2024;16:14081–14088. doi: 10.1039/D4NR02051K. [ DOI ] [ PubMed ] [ Google Scholar ] Sun Y. Yu X. Liu P. Han W. Xu W.-w. Su Y. Zhao J. Nanotechnology. 2022;34:105701. doi: 10.1088/1361-6528/aca80d. [ DOI ] [ PubMed ] [ Google Scholar ] Chen W. K. Fang W. H. Cui G. J. Chem. Phys. 2023;158:044110. doi: 10.1063/5.0129458. [ DOI ] [ PubMed ] [ Google Scholar ] Du J. Liao K. Ma J. Li W. Li S. J. Chem. Theory Comput. 2022;18:7630–7638. doi: 10.1021/acs.jctc.2c00911. [ DOI ] [ PubMed ] [ Google Scholar ] McCandler C. A. Pihlajamäki A. Malola S. Häkkinen H. Persson K. A. ACS Nano. 2024;18:19014–19023. doi: 10.1021/acsnano.4c03094. [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Karmakar A. Dhananjay Karmakar T. Nano Lett. 2025;25:10187–10192. doi: 10.1021/acs.nanolett.5c02218. [ DOI ] [ PubMed ] [ Google Scholar ] Chen T. Li J. Cai P. Yao Q. Ren Z. Zhu Y. Khan S. Xie J. Wang X. Nano Res. 2022;16:4188–4196. doi: 10.1007/s12274-022-5095-7. [ DOI ] [ Google Scholar ] Fang L. Laakso J. Rinke P. Chen X. J. Chem. Phys. 2024;160:094106. doi: 10.1063/5.0180529. [ DOI ] [ PubMed ] [ Google Scholar ] Pihlajamäki A. Hämäläinen J. Linja J. Nieminen P. Malola S. Kärkkäinen T. Häkkinen H. J. Phys. Chem. A. 2020;124:4827–4836. doi: 10.1021/acs.jpca.0c01512. [ DOI ] [ PubMed ] [ Google Scholar ] Tiwari V. Karmakar T. Nano Lett. 2025;25:5940–5946. doi: 10.1021/acs.nanolett.5c01340. [ DOI ] [ PubMed ] [ Google Scholar ] Wang Y. Samarasinghe D. S. N. D. Deng H. Liu B. Aikens C. M. J. Chem. Inf. Model. 2025;65:3892–3902. doi: 10.1021/acs.jcim.4c01495. [ DOI ] [ PubMed ] [ Google Scholar ] Associated Data This section collects any data citations, data availability statements, or supplementary materials included in this article. Data Availability Statement No primary research results, software or code have been included and no new data were generated or analysed as part of this review. Articles from Chemical Science are provided here courtesy of Royal Society of Chemistry ACTIONS View on publisher site PDF (2.9 MB) Cite Collections Permalink PERMALINK Copy RESOURCES Similar articles Cited by other articles Links to NCBI Databases Cite Copy Download .nbib .nbib Format: AMA APA MLA NLM Add to Collections Create a new collection Add to an existing collection Name your collection * Choose a collection Unable to load your collection due to an error Please try again Add Cancel Follow NCBI NCBI on X (formerly known as Twitter) NCBI on Facebook NCBI on LinkedIn NCBI on GitHub NCBI RSS feed Connect with NLM NLM on X (formerly known as Twitter) NLM on Facebook NLM on YouTube National Library of Medicine 8600 Rockville Pike Bethesda, MD 20894 Web Policies FOIA HHS Vulnerability Disclosure Help Accessibility Careers NLM NIH HHS USA.gov Back to Top

Related documents

• State-of-the-art and future developments in dosimetry for SSTR- and PSMA-targeting radioligand therapies.#1016098
• American Association for Thoracic Surgery Quality Gateway: A surgeon case study of its application in adult cardiac surgery for quality assurance.#1016099
• Metabolic Reprogramming of Cancer Stem Cells: Targeting Lipid Flux and Mitochondrial Plasticity to Overcome Therapeutic Resistance.#1016100
• Symptom Bother Is Associated With Medical Care‐Seeking for LUTS in Japanese Women: A Survey of 1800 Women.#1016101