Time-dependent Density Functional Theory Research Articles

Theoretical Insights into the Effect of Different Numbers of Thiophene Groups on Hydrogen Bond Interaction and Excited-State Intramolecular Proton-Transfer Process for Flavonoid Derivatives.

In this study, we systematically explored the impact of varying the number of thiophene groups on the hydrogen bond interaction and excited-state intramolecular proton-transfer (ESIPT) processes in flavonoid derivatives (STF, DTF, and TTF) using the density functional theory and time-dependent density functional theory methods. Initially, a thorough analysis of the optimized geometric structures revealed that the intramolecular hydrogen bond in the S1 state is enhanced and gradually weakened as the number of thiophene groups increases. To gain a deeper understanding of the hydrogen bond interaction, topological analysis, interaction region indicator scatter plots, and isosurface plots were employed. These images provide further insights that align with the structural analysis. Additionally, we observed a red-shift in the electronic spectra (absorption and fluorescence spectra), which is primarily attributed to the narrowing of the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, as elucidated by the frontier molecular orbitals. Furthermore, a combined analysis between the hole-electron distribution and the transition density matrix heat map shows that electron excitation involves the unidirectional charge-transfer mechanism. In the end, by conducting relaxed potential energy curve scans, we found that an increase in the number of thiophene groups elevates the energy barrier for ESIPT, making it more challenging for the reaction. In summary, our study underscores the vital effect of thiophene-substituted numbers in modulating the ESIPT process, which is able to provide valuable insights for the design and synthesis of desired organic fluorescent probes in the future.

Intramolecular Charge Transfer and Stimuli-Responsive Emission in Cholesterol-Appended Phenothiazine-Cyanostyryl-Based Donor-Acceptor Systems.

Organic fluorescent molecules have received considerable attention owing to their various optoelectronic applications. Herein, we report the design and synthesis of two cholesterol-functionalized cyanostyrene-phenothiazine-based D-π-A systems that are emissive in both the solution and solid states. The newly synthesized cholesterol-appended phenothiazine-cyanostyrene diads PTCS-1 and PTCS-2 vary in the N-alkylation of phenothiazine, respectively, with─octyl and─hexyl chains. Both molecules are highly fluorescent and show reasonably good quantum yields in nonpolar solvents because of twisted intramolecular charge transfer (TICT). The molecules exhibit aggregation-induced emission in the solid state. Due to the presence of flexible alkyl chains in the phenothiazine and cholesterol moieties, PTCS-1 and PTCS-2 show mechanochromic luminescence switching in response to external shear stress and emission recovery under methanol vapor. Powder X-ray diffraction studies prove that the emission switching on the applied stimuli in both PTCS-1 and PTCS-2 is attributed to the reversible transformation between the crystalline and amorphous states. Time-dependent density functional theory (TD-DFT) studies are carried out to gain insight into the ICT interactions. TD-DFT analysis at the TD-M06-2X/def2-TZVP level further revealed that in both molecules, the lowest unoccupied molecular orbital (LUMO) + 2, LUMO, highest occupied molecular orbital (HOMO), and HOMO - 1 orbitals are responsible for the charge transfer interactions. These ICT interactions are identified as π-π* type interactions.

Absorption and fluorescence spectroscopy of cold proflavine ions isolated in the gas phase.

Proflavine, a fluorescent cationic dye with strong absorption in the visible, has been proposed as a potential contributor to diffuse interstellar bands (DIBs). To investigate this hypothesis, it is essential to examine the spectra of cold and isolated ions for comparison. Here, we report absorption spectra of proflavine ions, trapped in a liquid-nitrogen-cooled ion trap filled with helium-buffer gas, as well as fluorescence spectra to provide further information on the intrinsic photophysics. We find absorption- and fluorescence-band maxima at 434.2 ± 0.1 and 434.7 ± 0.3nm, corresponding to a Stokes shift of maximum 48cm-1, which indicates minor differences between ground-state and excited-state geometries. Based on time-dependent density functional theory, we assign the emitting state to S2 as its geometry closely resembles that of S0, whereas the S1 geometry differs from that of S0. As a result, simulated spectra involving S1 exhibit long Franck-Condon progressions, absent in the experimental spectra. The latter displays well-resolved vibrational features, assigned to transitions involving in-plane vibrational modes where the vibrational quantum number changes by one. Dominant transitions are associated with vibrations localized on the NH2 moieties. Experiments repeated at room temperature yield broader spectra with maxima at 435.5 ± 1nm (absorption) and 438.0 ± 1nm (fluorescence). We again conclude that prevalent fluorescence arises from S2, i.e., anti-Kasha behavior, in agreement with previous work. Excited-state lifetimes are 5 ± 1ns, independent of temperature. Importantly, we exclude the possibility that a narrow DIB at 436.4nm originates from cold proflavine cations as the band is redshifted compared to our absorption spectra.

Non-adiabatic excited-state time-dependent GW molecular dynamics (TDGW) satisfying extended Koopmans' theorem: An accurate description of methane photolysis.

There is a longstanding difficulty that time-dependent density functional theory relying on adiabatic local density approximation is not applicable to the electron dynamics, for example, for an initially excited state, such as in photochemical reactions. To overcome this, we develop non-adiabatic excited-state time-dependent GW molecular dynamics (TDGW) on the basis of the extended quasiparticle theory. Replacing Kohn-Sham orbitals/energies with correlated, interacting quasiparticle orbitals/energies allows the full correspondence to the excited-state surfaces and corresponding total energies, with satisfying extended Koopmans' theorem. We demonstrate the power of TDGW using methane photolysis, CH4→CH3•+H, an important initiation reaction for combustion/pyrolysis and hydrogen production of methane. We successfully explore several possible pathways and show how this reaction dynamics is captured accurately through simultaneously time-tracing all quasiparticle levels. TDGW scales as O(NB3-4), where NB is the number of basis functions, which is distinctly advantageous to performing dynamics using configuration interaction and coupled cluster methods.

Molecular engineering of benzothiadiazole core based non-fullerene acceptors to tune the optoelectronic properties of perovskite solar cells

Non-fullerene acceptors exhibit durable absorption capacity, and exceptional power conversion efficiency that make them appealing photovoltaic material for use as hole transport materials in perovskite solar cells (PSCs). To interpret the benzothiadiazole as the central-core unit, five novel acceptors (U1-U5) have been successfully fabricated by end-capped engineering. Their electrical and optical performance have been investigated theoretically. Utilizing density functional theory (DFT) and time-dependent DFT computations, we have examined the photo-physical and optoelectronic characteristics. U4 exhibits 3.80 eV optical band gap with an absorption (λmax) of 670.58 nm. Lowest excitation energy (1.85 eV), greater LHE (0.98) and improved PCE (18.31 eV) validate the U4 designed molecule as a brilliant photovoltaic material for PSCs. Outstanding charge transfer at the donor–acceptor junction has also been revealed by thorough analysis of U4/PTB7-Th. Therefore, substitution of alternate functional groups in terminal moieties is efficient for effective PSCs that can be used as photon harvesting materials.

Electronic structure and optical spectral analysis of the MnO42- anion with consideration of site and Jahn-Teller distortion.

The ground state of 3d1 MnO42- was studied by density functional theory (DFT) and complete active space self-consistent field (CASSCF) methods in terms of a variety of molecular point group structures to ascertain the site and Jahn-Teller (JT) distortion effect. Modeling results from UB3LYP/6-31+G(d) calculations with natural bond orbital analysis show the four Mn-O bonds are coordinate covalent. The one-electron matrix elements from CASSCF(AILFT) (abinitio ligand field theory) with a second order perturbation treatment were used to calculate the parameters of the angular overlap model. These allowed the calculation of JT stabilization energies and first and second order JT coupling constants for MnO42- with C2v and D2d symmetry. Absorption spectra and excited state transition energies were calculated assuming the lattice distorted geometry was C2v with time-dependent DFT (TDDFT) using the unrestricted coulomb-attenuating UCAM-B3LYP density functional and with the spin-adapted spin-flip DFT using the UBH&HLYP density functional, both with the AUG-CC-PVTZ basis set. The mean absolute deviation of the calculations from experimental excited states for the ligand field and ligand-to-metal charge transfer (LMCT) bands was better than 0.1eV. Several new assignments for LMCT excited states were made on the basis of the TDDFT excitation results.

Exploring electron donor and acceptor effects: DFT analysis of ESIPT/GSIPT in 2-(oxazolinyl)-phenols for photophysical and luminophore enhancement.

Understanding excited-state intramolecular proton transfer (ESIPT) is essential for designing organic molecules to enhance photophysical and luminophore properties in the development of optoelectronic devices. In this context, an attempt has been made to understand the impact of substituents on the ESIPT process of 2-(oxazolinyl)-phenol. Electron donating (EDG: -NH2, -OCH3, and -CH3) and electron withdrawing (EWG: -Cl, -Br, -COOH, -CF3, -CN, and -NO2) substitutions have been computationally designed and screened through density functional theory (DFT) and time-dependent density-functional theory (TDDFT) calculations. Furthermore, the ground state intramolecular proton transfer and ESIPT mechanisms of these designed luminophores are explored using the transition state theory. The results reveal that molecules with EDG show higher absorption and emission peaks than molecules with EWG and also indicate that the mobility of charge carriers in 2-(oxazolinyl)-phenol derivatives is significantly influenced by substituents. We found that the EWGs decrease the reorganization energy and increase the vertical ionization potential and electron affinity values, as well as the highest occupied molecular orbital-lowest unoccupied molecular orbital gap, compared to the EDG substituted molecules. Significantly, the excited state (S1) of the keto emission (K) form shows notably larger values for the EDG substitutions. The intersystem crossing pathway efficiency weakens with reduced spin-orbit coupling matrix element in the enol form with electron-donating substituents and vice versa in the keto form during S1-T3 transitions. Our research links intramolecular proton transfers and triplet generation, making these substituted molecules appealing for optoelectronic devices. Introducing EDGs, such as -NH2, boosts the ESIPT reaction in 2-(oxazolinyl)-phenol. This study guides designing ESIPT emitters with unique photophysical properties.

Fission fragment distributions within time-dependent density functional theory

Fission fragment distributions within time-dependent density functional theory

Extending the limit of LR-TDDFT on two different approaches: Numerical algorithms and new Sunway heterogeneous supercomputer

First-principles time-dependent density functional theory (TDDFT) is a powerful tool to accurately describe the excited-state properties of molecules and solids in condensed matter physics, computational chemistry, and materials science. However, a perceived drawback in TDDFT calculations is its ultrahigh computational cost O(N5∼N6) and large memory usage O(N4) especially for plane-wave basis set, confining its applications to large systems containing thousands of atoms. Here, we present a massively parallel implementation of linear-response TDDFT (LR-TDDFT) and accelerate LR-TDDFT in two different aspects: (1) numerical algorithms on the X86 supercomputer and (2) optimizations on the heterogeneous architecture of the new Sunway supercomputer. Furthermore, we carefully design the parallel data and task distribution schemes to accommodate the physical nature of different computation steps. By utilizing these two different methods, our implementation can gain an overall speedup of 10x and 80x and efficiently scales to large systems up to 4096 and 2744 atoms within dozens of seconds.

Effects of Nanowire Doping on Plasmon-Enhanced N2 Dissociation.

Doping a transition metal element into plasmonic systems can tune the optical properties of the system, which will potentially facilitate the plasmon-enhanced catalytic process. In this study, we applied the linear-response time-dependent density functional theory (LR-TDDFT) method with real-time electron dynamics and mean-field Ehrenfest dynamics methods to computationally investigate the effects of doping silver nanowires on plasmon-enhanced N2 dissociation. We calculated the absorption spectra for different doped systems, applied an external electric field to the system, and performed mean-field Ehrenfest dynamics to examine how plasmonic excitation will affect the N2 activation or dissociation. In addition, we also studied how the transition metal dopant affects the system's electronic structure and potential energy surface.

Axial coordinated manganeseporphyrin/tetraazaporphyrin-4-(10-phenylanthracen-9-yl)pyridine dyads: self-assembly, structure and spectral properties in ground and excited states.

Self-assembly of new donor-acceptor systems based on (5,10,15,20-tetraphenylporphinato)manganese(III)/(5,10,15,20-tetra-4-tert-butylphenylporphinato)manganese(III)/(octakis(4-tert-butylphenyl)tetraazaporphinato)manganese(III) acetate ((AcO)MnTPP/(AcO)MnTBPP/(AcO)MnTAP) and 4-(10-phenylanthracen-9-yl)pyridine (PyAn) was studied using fluorescence spectroscopy and mass spectrometry. It was found that the coordination complexes of 1:1 composition (dyads) are formed in toluene. The spectral properties, the chemical structures and redox behavior of the dyads were described using 1H NMR, IR, EPR spectroscopy and cyclic voltammetry, respectively. The dynamic processes and the characteristics in the excited state of the dyads were obtained using the femtosecond transient absorption spectroscopy method. Density functional theory (DFT), time-dependent DFT methods were used to elucidate the dyad electronic structures and to establish the differences in their frontier molecular orbitals. The analysis of the lambda parameter and the distance of hole-pair interaction was indicated more favorable charge transfer between the macrocycle and the axial PyAn fragment in (AcO)(PyAn)MnTAP. The calculated values of the zero-field splitting parameters D and E/D, together with the g tensors of the lowest spin-orbit state for (AcO)MnTPP and (AcO)(PyAn)MnTPP were obtained using the combination of DFT and Multireference Perturbation Theory (CASSCF/NEVPT2) simulations. The data obtained develop the fundamental basis in the field of photovoltaics and show the prospects for the study of molecular systems of this class.

Photophysical Studies of a Zr(IV) Complex with Two Pyrrolide-Based Tetradentate Schiff Base Ligands.

The reaction of two equivalents of N,N'-bis(2-pyrrolylmethylidene)-1,2-phenylenediamine (H2bppda) with tetrabenzylzirconium provided the air- and moisture-stable eight-coordinate complex Zr(bppda)2. Temperature-dependent steady-state and time-resolved emission spectroscopy established weak photoluminescence (ΦPL = 0.4% at 293 K) by a combination of prompt fluorescence and thermally activated delayed fluorescence (TADF) upon visible light excitation at and around room temperature. TADF emission is strongly quenched by 3O2 and shows highly temperature-sensitive emission lifetimes of hundreds of microseconds. The lifetime of the lowest energy singlet excited state, S1, was established by transient absorption spectroscopy and shows rapid deactivation (τ = 142 ps) by prompt fluorescence and intersystem crossing to the triplet state, T1. Time-dependent density functional theory (TD-DFT) calculations predict moderate ligand-to-metal charge transfer (LMCT) contributions of 25-30% for the S1 and T1 states. A comparison of Zr(bppda)2 to related zirconium pyridine dipyrrolide complexes, Zr(PDP)2, revealed important electronic structure changes due to the eight-coordinate ligand environment in Zr(bppda)2, which were correlated to differences in the photophysical properties between the two compound classes.

Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States.

The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.

Inverse Hypercorroles.

Ground-state and time-dependent density functional theory (TDDFT) calculations with the long-range-corrected, Coulomb-attenuating CAMY-B3LYP exchange-correlation functional and large, all-electron STO-TZ2P basis sets have been used to examine the potential "inverse hypercorrole" character of meso-p-nitrophenyl-appended dicyanidocobalt(III) corrole dianions. The effect is most dramatic for 5,15-bis(p-nitrophenyl) derivatives, where it manifests itself in intense NIR absorptions. The 10-aryl groups in these complexes play a modulatory role, as evinced by experimental UV-visible spectroscopic and electrochemical data for a series of 5,15-bis(p-nitrophenyl) dicyanidocobalt(III) corroles. TDDFT (CAMY-B3LYP) calculations ascribe these features clearly to a transition from the corrole's a2u-like HOMO (retaining the D4h irrep used for metalloporphyrins) to a nitrophenyl-based LUMO. The outward nature of this transition contrasts with the usual phenyl-to-macrocycle direction of charge transfer transitions in many hyperporphyrins and hypercorroles; thus, the complexes studied are aptly described as inverse hypercorroles.

Longitudinal Extension of Double π-Helix Enables Near-Infrared Amplified Dissymmetry and Chiroptical Response.

Near-infrared (NIR) circularly polarized light absorbing or emitting holds great promise for highly sensitive and precise bioimaging, biosensing, and photodetectors. Aiming at designing NIR chiral molecular systems with amplified dissymmetry and robust chiroptical response, herein, we present a series of double π-helical dimers with longitudinally extended π-entwined substructures via Ullmann or Yamamoto homocoupling reactions. Circular dichroism (CD) spectra revealed an approximate linear bathochromic shift with the rising number of naphthalene subunits, indicating a red to NIR chiroptical response. Particularly, the terrylene diimide-entwined dimers exhibited the strongest CD intensities, with the maximal |Δε| reaching up to 393 M-1 cm-1 at 666 nm for th-TDI[2]; and a record-high chiroptical response (|ΔΔε|) between the neutral and dianionic species of 520 M-1 cm-1 at 833 nm for th-TDI[2]Cl was achieved upon further reduction to its dianionic state. Time-dependent density functional theory (TDDFT) calculations suggested that the pronounced intensification of the CD spectra originated from a simultaneous enhancement of both electric (μ) and magnetic (m) transition dipole moments, ultimately leading to an overall increase in the rotatory strength (R). Notably, the circularly polarized luminescence (CPL) brightness (BCPL) reached 77 M-1 cm-1 for th-TDI[2]Cl, among the highest values reported for NIR-CPL emitters. Furthermore, all chiral dianions exhibited excellent air stability under ambient conditions with half-life times of up to 10 days in N-methylpyrrolidone (NMP), which is significant for future biological applications and chiroptic switches.

The density–density response function in time-dependent density functional theory: Mathematical foundations and pole shifting

The density–density response function in time-dependent density functional theory: Mathematical foundations and pole shifting

Versatile Aggregation Induced Emissive Molecular Probes Incorporating different Donor/Acceptor groups for Tunable Multiple Color Fluorescence

Excited state intramolecular proton transfer (ESIPT) and aggregation induced emission active multicolour fluorescent probes for bioimaging applications has grown popularly in recent years. Herein, we report a series of different donor/acceptor units (Cl, NO2 and O-CH3) bearing anthracene-salicylaldehyde based fluorescent probes and their antibacterial activity and fluorescence imaging of E. Coli bacteria and zebrafish embryos. All the compounds (ASA1, ASA2, ASA3 and ASA4) synthesized via facile one pot condensation reaction of anthracene hydrazide with Salicylaldehyde, 5-Chlorosalicylaldehyde, 5-Nitrosalicylaldehyde, 5-methoxysalicylaldehyde and characterized using all the possible spectroscopic techniques. The noteworthy tuneable emission wavelength shifts from 471 nm to 625 nm and ESIPT characteristics of the compounds were demonstrated through photophysical experiments and single crystal X-ray diffraction analysis. All the experimental results were completely analysed theoretically on the ground as well as the excited state using density functional theory (DFT) and time dependant density functional theory (TD-DFT) calculations. Compounds ASA1-ASA4 showed fascinating different colour fluorescent (Green, Yellow and Orange) aggregation induced emission features in the THF-H2O mixture. AIE characteristics of ASA1, ASA2, ASA3 and ASA4 were recognized using absorption, emission, fluorescence quantum yield (φ), DLS and SEM analysis. ASA1, ASA2, ASA3 and ASA4 showed excellent antibacterial activity and ROS generation against E. Coli bacteria. Further, owing to the multicolour fluorescent property of the compounds successfully applied for fluorescent imaging in E. Coli bacteria and live zebrafish embryos.

Controlling Chirogenic Effects in Porphyrin Based Supramolecular Systems: Theoretical Analysis Versus Experimental Observations.

Electronic circular dichroism (ECD) spectroscopy is a widely employed method for studying chiral analysis, requiring the presence of a chromophore close to a chiral centre. Porphyrinoids are found to be one of the best chromophoric systems serving for this purpose and enabling the application of ECD spectroscopy for chirality determination across diverse classes of organic compounds. Consequently, it is crucial to understand the induction mechanisms of ECD in the porphyrin-based complexes. The present study explores systematically the influence of secondary chromophores, bonded to an achiral zinc porphyrin or to chiral guest molecules, on the B-region of ECD spectra using the time-dependent density functional theory (TD-DFT) calculations. The study analyses the impact of change in both the conformation of achiral porphyrin (host) and change in position and conformation of chiral organic molecule (guest) on the B-band of ECD spectra (energy, intensity, sign of Cotton effect). Finally, conclusions made on model complexes are applied to published experimental data, contributing to a deeper understanding of various factors influencing ECD spectra in chiral systems. In addition, a computer program aimed to help rationalise ECD spectra by visualizing corresponding orbital energies, rotatory strengths, electric and magnetic transition moments, and angles between them, is presented.

Theoretical insights into long-range coupling of electron-hole pairs in TCTA–PO-T2T exciplex

This work reports theoretical investigations concerning long-range coupling of electron-hole pairs in spatially separated exciplex (SSE) systems via Density Functional Theory and Time-Dependent Density Functional Theory. Based on TCTA−PO-T2T parent exciplex, where 4,4′,4″-Tris(carbazol-9-yl)triphenylamine (TCTA) serves as donor (D) and 2,4,6-Tris[3-(diphenylphosphinyl)phenyl]-1,3,5-triazine (PO-T2T) as acceptor (A), SSEs are constructed by intentionally tuning the D−A distance. The calculation results demonstrate that all SSEs, even if the D−A distance is over 14 Å, exhibit the clear and strong charge transfer character in terms of DCT, qCT, and t indexes. As the D−A distance increases, the energy gap of SSEs is increased, resulting in the blueshift of emission spectra. Calculation and experiment results show a good consistency, indicating that our model could well describe SSEs. Meanwhile, the reduced energy gap between CT and 3LE and more degenerate states could boost reverse intersystem crossing via vibronic coupling and hyperfine coupling, eventually improving the electroluminescent performance of SSEs. Our study suggests that manipulation of relative energy alignment via controlling D−A distance not only promotes the properties of exciplexes, but also offers guidance for designing thermally activated delayed fluorescence emitters through-space charge transfer.

Ultrafast hydrogen production in boron/oxygen-codoped graphitic carbon nitride revealed by nonadiabatic dynamics simulations.

Graphitic carbon nitride (g-C3N4 or GCN) shows promise in photocatalytic water splitting, despite facing the challenge of rapid electron-hole recombination. In this study, we investigated the influence of boron/oxygen codoping on the photocatalytic performance of GCN systems for hydrogen generation. First-principles calculations and nonadiabatic molecular dynamics (NAMD) simulations were employed to reveal that the recombination time of photogenerated carriers could be increased by 16% to 64% in the codoped systems compared to the pristine GCN. The time-dependent density functional theory (TDDFT) scheme was utilized to select energy windows and initiate dynamics in cluster models of B/O co-doped heptazine with water molecules. Notably, we observed efficient direct photodissociation of hydrogen atoms from water molecules within 60 fs and proton hops within the hydrogen-bonded network within 80 fs in the co-doped system, diverging from the previously proposed mechanism for pristine heptazine in NAMD simulations. This discovery underscores the significant role of faster proton-coupled electron transfer (PCET) reactions and rapid radiationless relaxation in achieving high photocatalytic efficiency in water splitting. Our work enhances the understanding of the internal mechanism of highly efficient photocatalysts for water splitting and provides a new design strategy for doped GCN.