Long-range corrected density functional theory study on static second hyperpolarizabilities of singlet diradical systems

Theoretical Study on the Correlation between Open-Shell Electronic Structures and Third-Order Nonlinear Optical Properties in One-Dimensional Chains of π-Radicals

This paper theoretically investigated the correlation between the open-shell electronic structure and third-order nonlinear optical (NLO) properties of one-dimensional (1D) stacked chains of π-radicals. By employing the finite N-mer models consisting of methyl or phenalenyl radicals with different stacking distances, we evaluated the average and standard deviation of diradical characters yi for N-mer models of π-radicals (yav and ySD). Then, we estimated these diradical characters at the limit of N → ∞. These y-based indices were helpful in discussing the correlation between the open-shell electronic structures and the second hyperpolarizability per dimer at the limit N → ∞, γ∞ for the 1D chains with stacking distance alternation (SDA). The calculated γ∞ values and the polymer/dimer ratio γ∞/γ(N = 2) were enhanced significantly when both the stacking distance and SDA are small. We also found that the spin-unrestricted long-range-corrected (LC-)UBLYP method with the range-separating parameter μ = 0.47 bohr-1 well reproduced the trend of γ∞ of this type of 1D chain estimated at the spin-unrestricted coupled-cluster levels. The present study is expected to contribute to establishing the design guidelines for future high-performance open-shell molecular NLO materials.

Total and orbital density-based analyses of molecules revealing long-range interaction regions

Total and orbital electron densities of molecules are explored for the effect of the long-range correction (LC) for density functional theory (DFT) exchange functionals by comparing to the effect of the ab initio coupled cluster singles and doubles (CCSD) method. Calculating the LC effect on the total electron densities shows that the LC stabilizes the electrons around the long-range interaction regions of kinetic energy density, which are assumed to be electrons other than free electrons and self-interacting electrons, while the CCSD method stabilizes the electrons in the long-range interaction regions in the vertical molecular planes. As a more precise test, the LC effect on orbital densities are compared to the CCSD effect on Dyson orbital densities. Surprisingly, these effects are similar for the unoccupied orbitals, indicating that the LC covers the effects required to reproduce the CCSD Dyson unoccupied orbitals. For exploring the discrepancies between these effects on the occupied orbitals, the photoionization cross sections are calculated as a direct test for the shapes of the HOMOs to investigate the differences between these effects on the occupied orbitals. Consequently, the LC clearly produces the canonical HOMOs close to the CCSD Dyson and experimental ones, except for the HOMO of benzene molecule that mixes with the HOMO - 1 for the CCSD Dyson orbitals. This indicates that the orbital analyses using the photoionization cross sections are available as a direct test for the quality of DFT functionals.

Machine Learning-Based Screening for Potential Singlet Fission Chromophores: The Challenge of Imbalanced Data Sets

Excitation with one photon of a singlet fission (SF) material generates two triplet excitons, thus doubling the solar cell efficiency. Therefore, the SF molecules are regarded as new generation organic photovoltaics, but it is hard to identify them. Recently, it was demonstrated that molecules of low-to-intermediate diradical character (DRC) are potential SF chromophores. This prompts a low-cost strategy for finding new SF candidates by computational high-throughput workflows. We propose a machine learning aided screening for SF entrants based on their DRC. Our data set comprises 469 784 compounds extracted from the PubChem database, structurally rich but inherently imbalanced regarding DRC values. We developed well performing classification models that can retrieve potential SF chromophores. The latter (∼4%) were analyzed by K-means clustering to reveal qualitative structure-property relationships and to extract strategies for molecular design. The developed screening procedure and data set can be easily adapted for applications of diradicaloids in photonics and spintronics.

Theoretical Study on Singlet Fission in Aromatic Diaza s-Indacene Dimers

We theoretically show that diaza (N₂)-substitution to s-indacene with 4n π-electrons, by which the number of π-electrons in N₂-s-indacene amounts to 4n+2, is a new strategy to design efficient singlet fission (SF) molecules. By N₂-substitution, the diradical character and the exchange integral are found to be tuned moderately, leading to satisfying the excitation energy level matching condition for SF with a large triplet excitation energy. On the basis of the effective electronic coupling related to the SF rate, we explore the optimal slip-stack dimer packings for fast SF. Their underlying mechanisms are well understood from the odd-electron density, resonance structure, and frontier orbital distribution, as the functions of the N₂-substituted positions. Furthermore, aromaticities of N₂-s-indacenes are evaluated explicitly on the basis of the magnetically induced current. Although N₂-s-indacenes display strengths of aromaticities similar to that of anthracene, a local decrease in aromaticity is found to correlate to the spatial feature of diradical character, i.e., odd-electron density. The present findings not only newly propose N₂-s-indacenes as feasible SF molecules but also contribute to comprehending the interplay between aromaticity and diradical electronic structures contributing to SF.

Diradicalar Character and Ring Stability of Mesoionic Heterocyclic Oxazoles and Thiazoles by Ab Initio Mono and Multi-Reference Methods

Mesoionics are neutral compounds that cannot be represented by a fully covalent or purely ionic structure. Among the possible mesomeric structures of these compounds are the diradical electronic configurations. Theoretical and experimental studies indicate that some mesoionic rings are unstable, which may be related to a significant diradical character, that until then is not quantified. In this work, we investigated the diradical character of four heterocycles: 1,3-oxazol-5-one, 1,3-oxazol-5-thione, 1,3-thiazole-5-one, and 1,3-thiazole-5-thione. The oxazoles are known to be significatively less stable than thiazoles. DFT and ab initio single (B3LYP, MP2, CCSD, and QCISD) and ab initio multi-reference (MR-CISD) methods with three basis sets (6-311+G(d), aug-cc-pVDZ, and aug-cc-pVTZ) were employed to assess the diradical character of the investigated systems, in gas phase and DMSO solvent, from three criteria: (i) HOMO-LUMO energy gap, (ii) determination of energy difference between singlet and triplet wave functions, and (iii) quantification of the most significant diradical character (y₀, determined in the unrestricted formalism). All of the results showed that the diradical character of the investigated systems is very small. However, the calculated electronic structures made it possible to identify the possible origin of the oxazoles instability, which can help the design of mesoionic systems with the desired properties.

Theoretical study on aromatic and open-shell characteristics of carbon nanobelts composed of indeno[1,2-b]fluorene units: dependence on the number of units and charge states

In this study, we theoretically investigate the aromatic and open-shell characteristics of carbon nanobelts (CNBs) composed of five- and six-membered rings. We have designed nanobelts composed of indeno[1,2-b]fluorene ([1,2-b]IF) units, which are referred to as [N]IF-CNB (N: the number of five-membered rings). The number of π-electrons, n π, in neutral [N]IF-CNB is 7N, and thus depending on N and charge states, n π can be 4n + 2 and 4n. Quantum chemical calculations on neutral [6]IF-CNB and [8]IF-CNB and dicationic [8]IF-CNB2+ have revealed that they are expected to exhibit unique aromatic and open-shell characteristics depending on n π, there are several analogies of the electronic structures in [N]IF-CNB to those in [N]annulene. Delocalized and intermediate open-shell electronic structures of [N]IF-CNB are also useful to drastically change the third-order nonlinear optical properties. These results suggest that theoretically designed [N]IF-CNB can be attractive and challenging targets of organic synthesis for realizing novel open-shell functional conjugated macrocycles.

Singlet-Fission-Induced Enhancement of Third-Order Nonlinear Optical Properties of Pentacene Dimers

Using quantum chemical calculations and exciton dynamics simulations, we investigate the static second hyperpolarizability γ [the third-order nonlinear optical (NLO) property at the molecular scale] of slip-stacked pentacene dimer models in the correlated-triplet-pair [1(TT)] state created from the singlet excited state in the singlet fission (SF) process. It is found that the SF induces significant (∼20 times at maximum) enhancement of γ/monomer in the 1(TT) state as compared to that in the singlet ground state. The origin of the remarkable enhancement of γ/monomer is revealed by analyzing the γ density distribution and the intermolecular orbital interaction. Furthermore, we clarify molecular packings suitable for highly efficient SF and largely enhanced γ values of a novel class of SF-induced NLO systems, which have promising potential to surpass the conventional NLO systems.

Structural diradical character

A reliable first-principles description of singlet diradical character is essential for predicting nonlinear optical and magnetic properties of molecules. As diradical and closed-shell electronic structures differ in their distribution of single, double, triple, and aromatic bonds, modeling electronic diradical character requires accurate bond-length patterns, in addition to accurate absolute bond lengths. We therefore introduce structural diradical character, which we suggest as an additional measure for comparing first-principles calculations with experimental data. We employ this measure to identify suitable exchange-correlation functionals for predicting the bond length patterns and electronic diradical character of a biscobaltocene with the potential for photoswitchable nonlinear optical activity. Out of four popular approximate exchange-correlation functionals with different exact-exchange admixtures (BP86, TPSS, B3LYP, TPSSh), the two hybrid functionals TPSSh and B3LYP perform best for diradical bond length patterns, with TPSSh being best for the organometallic validation systems and B3LYP for the organic ones (for which the D3 dispersion correction was included). Still, none of the functionals is suitable for correctly describing relative bond lengths across the range of molecules studied, so that none can be recommended for predictive studies of (potential) diradicals without reservation. © 2018 Wiley Periodicals, Inc.

Theoretical Study on Open-Shell Singlet Character and Second Hyperpolarizabilities in Cofacial π-Stacked Dimers Composed of Weak Open-Shell Antiaromatic Porphyrins

From the analysis based on the broken-symmetry density functional theory (DFT) calculations, we in this study propose a strategy to enhance the open-shell characters and third-order nonlinear optical (NLO) properties of π-stacked dimers composed of antiaromatic molecules with weak open-shell characters. For this purpose, we here constructed cofacial π-stacked dimer models composed of aromatic and antiaromatic Ni^II porphyrins in order to examine the π-π stacking distance (R) dependence of the diradical characters (y) and static second hyperpolarizabilities (γ). The antiaromatic porphyrin dimers are found to have intermediate y around R∼3.3 Å, the result of which originates in the unique intermolecular interactions between the antiaromatic monomers. Static γ along the stacking direction of such antiaromatic porphyrin dimers with intermediate diradical characters are shown to be enhanced significantly as compared to those of the isolated monomers and the aromatic porphyrin dimers.

Exploring the novel donor-nanotube archetype as an efficient third-order nonlinear optical material: asymmetric open-shell carbon nanotubes

Contrary to the enormous number of previous studies on carbon nanotubes (CNTs), herein, we realized the origin of the intrinsic open-shell diradical character and second hyperpolarizability γ using a broken symmetry approach. This study was inspired by our recent findings (S. Muhammad, et al., Nanoscale, 2016, 8, 17998 and Nakano, et al., J. Phys. Chem. C, 2016, 120, 1193). We performed structural modifications through a unique asymmetric donor-nanotube framework, which led to a novel paradigm of modified CNTs with tunable open-shell diradical character and remarkably superior NLO response properties. Interestingly, asymmetry and diradical character were found to be the crucial factors to modulate the second hyperpolarizability γ. We initially performed a comparative analysis of the diradical characters and γ amplitudes of boron nitride nanotubes (BNNTs) and CNTs possessing significant ionic characters and covalent characters, respectively. The basic findings for these simple configurations were further extended to the donor-acceptor CNT paradigm, which finally led to excellent asymmetric donor-CNT configurations with remarkably larger γ amplitudes. Furthermore, among the CNTs, finite length zigzag CNT(6,0)3 were modified with different donor-acceptor configurations. Interestingly, for the first time, unique donor-nanotube configurations [1,4-(NH2)2CNT-(6,0)3 and 1,4-(NH2)2CNT-(6,0)5] were found; they showed significantly robust γ amplitudes as large as 2519 × 103 and 4090 × 103 a.u. at the LC-UBLYP(μ = 0.33)/6-31G* level of theory. Additionally, several molecular level insights have been obtained for these novel donor-nanotube configurations using their odd electron densities, molecular electrostatic maps, densities of states and γ density analyses to highlight the realization of these novel materials for highly efficient optical and NLO applications.

Tunability of Open-Shell Character, Charge Asymmetry, and Third-Order Nonlinear Optical Properties of Covalently Linked (Hetero)Phenalenyl Dimers

Tunability of the open-shell character, charge asymmetry, and third-order nonlinear optical (NLO) properties of covalently linked (hetero)phenalenyl dimers are investigated by using the density functional theory method. By changing the molecular species X and substitution position (i, j) for the linker part, a variety of intermonomer distances R and relative alignments between the phenalenyl dimers can be realized from the geometry optimizations, resulting in a wide-range tuning of diradical character y and charge asymmetry. It is found that the static second hyperpolarizabilities along the stacking direction, γ_yyyy , are one-order enhanced for phenalenyl dimer systems exhibiting intermediate y, a feature that is in good agreement with the "y-γ correlation". By replacing the central carbon atoms of the phenalenyl rings with a boron or a nitrogen, we have also designed covalently linked heterophenalenyl dimers. The introduction of such a charge asymmetry to the open-shell systems, which leads to closed-shell ionic ground states, is found to further enhance the γ_yyyy values of the systems having longer intermonomer distance R with intermediate ionic character, that is, charge asymmetry. The present results demonstrate a promising potential of covalently linked NLO dimers with intermediate open-shell/ionic characters as a new building block of highly efficient NLO systems.

Unravelling the dissociation pathways of acetic acid upon electron transfer in potassium collisions: experimental and theoretical studies

Electron transfer in alkali-molecule collisions with gas phase acetic acid and its deuterated analogues resulting in OH^- formation requires considerable internal rearrangement in the temporary negative ion. At a collision energy well above the threshold of negative ion formation, electron transfer from potassium to CH₃COOH/CH₃COOD and CD₃COOH results not only in H transfer from CH₃ to COOH/COOD, but also in H release from COOH and subsequent rearrangement to eliminate OH^-. These processes are also investigated by theoretical post-Hartree-Fock and DFT calculations. The combination of both studies reveals that the most favourable intermediate mechanism occurs via diol formation. Such intramolecular H transfer is reported here for the first time in the context of electron transfer induced dissociation experiments in alkali-molecule collisions. A comprehensive fragmentation study is presented and dissociation mechanisms are suggested.

Third-Order Nonlinear Optical Properties of One-Dimensional Quinoidal Oligothiophene Derivatives Involving Phenoxyl Groups

The diradical characters (y) and third‐order nonlinear optical (NLO) properties of open‐shell quinoidal oligothiophene derivatives with phenoxyl groups, and the corresponding reduced (hydrogenated)‐state oligomers, are investigated by using the broken‐symmetry density functional theory method. The oxidized (dehydrogenated) states are predicted to have an open‐shell singlet ground state and their y values increase with the number of units. Static second hyperpolarizabilities (γ) of the open‐shell oligomers with intermediate y are shown to be enhanced significantly compared with those of the closed‐shell analogues. Furthermore, owing to the effective diradical distances, the γ values of open‐shell oligomers are found to exceed that of s‐indaceno[1,2,3‐cd;5,6,7‐c′d′]diphenalene, which is known as an organic molecule with the largest two‐photon absorption cross‐section in this size of the pure hydrocarbons. This feature extends the range of efficient open‐shell third‐order NLO materials to a novel class of one‐dimensional conjugated oligomers with redox‐based high tunability of third‐order NLO properties.

Theoretical Study on the Open-Shell Singlet Nature and the Second Hyperpolarizabilities of Corannulene Derivatives with Two Phenoxyl Radicals

Using the spin-unrestricted density functional theory method, we investigate the interplay between the diradical character y and second hyperpolarizabilities γ (the third-order nonlinear optical (NLO) properties at the molecular scale) of corannulene derivatives with two phenoxyl radicals. This molecule in the singlet state exhibits intermediate y and thus displays a significantly larger γ value than the triplet state and the closed-shell bis-phenol analogue. We also examine the planar molecules involving a coronene moiety in place of the curved corannulene. The intermediate y and large γ values of the corannulene systems are found to originate not from their curved skeleton but from the equilibrium between benzenoid/quinoid resonance forms due to delocalization of the radical electrons of the terminal phenoxyl rings. The longitudinal γ value of the singlet state is found to be comparable to that of s-indaceno[1,2,3-cd;5,6,7-c'd']diphenalene, which is known to be one of the organic molecules with the largest two-photon absorption cross section in this size of pure hydrocarbons. The present system is thus expected to be a promising candidate for highly efficient open-shell NLO molecules.

A theoretical study on quasi-one-dimensional open-shell singlet ladder oligomers: multi-radical nature, aromaticity and second hyperpolarizability

Quasi-one-dimensional open-shell singlet ladder oligomers composed of indenofluorene-like units show strong correlation between multiradical nature, aromaticity and second hyperpolarizability.

Theoretical Study on the Second Hyperpolarizailities of Oligomeric Systems Composed of Carbon and Silicon π-Structures

To explore the prospect of molecules involving silicon-silicon multiple bonds as nonlinear optical molecular systems, the relationship between the structure and the second hyperpolarizabilities γ of the oligomeric systems composed of carbon and silicon π-structures is investigated using the density functional theory method. It is found that these compounds indicate intramolecular charge transfer (ICT) from the silicon units to the carbon units together with nonzero diradical characters. The γ values of these compounds are shown to be 2–13 times as large as those of the carbon analogs. Although asymmetric carbon and silicon π-systems exhibit comparable enhancement to the corresponding symmetric systems, donor-π-donor structures exhibit remarkable enhancement of γ despite of their both-end short silicon π-chain moieties (donor units). Further analysis using the odd electron and γ densities clarifies that the intermediate diradical character also contributes to the enhancement of γ. These results predict that even short π-conjugated silicone moieties can cause remarkable enhancement of γ by introducing them into π-conjugated hydrocarbon structures.

Intermolecular modulation of IR intensities in the solid state. The role of weak interactions in polyethylene crystal: A computational DFT study

Density functional theory calculations with periodic boundary conditions are exploited to study the infrared spectrum of crystalline polyethylene. Spectral changes lead by the intermolecular packing in the orthorhombic three-dimensional crystal are discussed by means of a careful comparison with calculations carried out for an isolated polymer chain in the all-trans conformation, described as an ideal one-dimensional crystal. The results are analyzed in the framework of the "oligomer approach" through the modelling of the IR spectrum of n-alkanes of different lengths. The study demonstrates that a relevant absorption intensity modulation of CH₂ deformation transitions takes place in the solid state. This finding suggests a new interpretation for the experimental evidences collected in the past by means of IR intensity measurement during thermal treatment. Moreover, the comparison between calculations for 3-D crystal and for the isolated polyethylene chain (1-D crystal) allows to put in evidence the effect of the local electric field on the computed infrared intensities. This observation provides guidelines for the comparison between infrared absorption intensities predicted for an isolated unit and for a molecule belonging to a crystal, through the introduction of suitable correction factors based on the refraction index of the material and depending on the dimensionality of such units (0D-molecule; 1D-polymer; 2D-slab).

Challenging compounds for calculating molecular second hyperpolarizabilities: the triplet state of the trimethylenemethane diradical and two derivatives

The second hyperpolarizability γ of trimethylenemethane (TMM) and two 1,3-dipole derivatives (NXA and OXA) in their triplet ground state has been evaluated at the UCCSD(T) level with the d-aug-cc-pVDZ extended basis set, highlighting that γ decreases from TMM to NXA and OXA, following the opposite order of their permanent dipole moments. These results are then used to benchmark a broad range of levels of approximation. So, the UMP2, UMP4, and UCCSD methods can be used to characterize γ of TMM and NXA but not of OXA. In that case, the large field-induced charge transfer contribution is difficult to handle using the MPn methods and only the UCCSD method provides values close to the UCCSD(T) reference. Turning to the performance of DFT with typical exchange-correlation functionals, the UM06-2X functional, which contains 54% of HF exchange, performs very well with a maximum of 4.5% of difference with respect to the reference values. On the other hand, employing less HF exchange leads to an overestimation of the responses whereas range-separated hybrids generally underestimate the second hyperpolarizabilities. Finally, the use of spin-projected methods for these 1,3-dipole triplet molecules has a little impact since the spin contamination is almost negligible.

Nonlinear optical properties in open-shell molecular systems

For more than 30 years, nonlinear optical (NLO) properties of molecular systems have been actively studied both theoretically and experimentally due to their potential applications in photonics and optoelectronics. Most of the NLO molecular systems are closed-shell species, while recently open-shell molecular species have been theoretically proposed as a new class of NLO systems, which exhibit larger NLO properties than the traditional closed-shell NLO systems. In particular, the third-order NLO property, the second hyperpolarizability γ, was found to be strongly correlated to the diradical character y, which is a quantum-chemically defined index of effective bond weakness or of electron correlation: the γ values are enhanced in the intermediate y region as compared to the closed-shell (y = 0) and pure open-shell (y = 1) domains. This principle has been exemplified by accurate quantum-chemical calculations for polycyclic hydrocarbons including graphene nanoflakes, multinuclear transition-metal complexes, main group compounds, and so on. Subsequently, some of these predictions have been substantiated by experiments, including two-photon absorption. The fundamental mechanism of the y-γ correlation has been explained by using a simple two-site model and the valence configuration interaction method. On the basis of this y-γ principle, several molecular design guidelines for controlling γ have been proposed. They consist in tuning the diradical characters through chemical modifications of realistic open-shell singlet molecules. These results open a new path toward understanding the structure-NLO property relationships and toward realizing a new class of highly efficient NLO materials. WIREs Comput Mol Sci 2016, 6:198-210. doi: 10.1002/wcms.1242.

Open-Shell Singlet Nature and σ-/π-Conjugation Effects on the Third-Order Nonlinear Optical Properties of Si Chains: Polysilane and Poly(disilene-1,2-diyl)

We theoretically investigate the open-shell singlet nature and σ-/π-conjugation effects on the longitudinal second hyperpolarizabilities (γ) of one-dimensional chains involving silicon-silicon bonds, that is, polysilane and poly(disilene-1,2-diyl), by comparison with their carbon analogues, polyethylene and polyacetylene, respectively. It is found that poly(disilene-1,2-diyl) has less bond length alternation than polyacetylene and that σ-conjugation of polysilane is less effective on the enhancement of γ than π-conjugation of polyacetylene, whereas π-conjugation of poly(disilene-1,2-diyl) indicates a more than 20 times greater enhancement of the γ than that of polyacetylene, which is known to be a typical nonlinear optical molecule with large γ, for one-dimensional chains involving 3-5 double bonds. Further theoretical analyses of poly(disilene-1,2-diyl) reveal the σ- and π-electrons contribute negatively and positively to the γ, respectively. The latter contribution is significantly larger than the former and thus causes the remarkable enhancement of γ amplitudes due to the emerging open-shell singlet nature in the long π-conjugation length.

PPV Polymerization through the Gilch Route: Diradical Character of Monomers

Despite various studies on the polymerization of poly(p-phenylene vinylene) (PPV) through different precursor routes, detailed mechanistic knowledge on the individual reaction steps and intermediates is still incomplete. The present study aims to gain more insight into the radical polymerization of PPV through the Gilch route. The initial steps of the polymerization involve the formation of a p-quinodimethane intermediate, which spontaneously self-initiates through a dimerization process leading to the formation of diradical species; chain propagation ensues on both sides of the diradical or chain termination occurs by the formation of side products, such as [2.2]paracyclophanes. Furthermore, different p-quinodimethane systems were assessed with respect to the size of their aromatic core as well as the presence of heteroatoms in/on the conjugated system. The nature of the aromatic core and the specific substituents alter the electronic structure of the p-quinodimethane monomers, affecting the mechanism of polymerization. The diradical character of the monomers has been investigated with several advanced methodologies, such as spin-projected UHF, CASSCF, CASPT2, and DMRG calculations. It was shown that larger aromatic cores led to a higher diradical character in the monomers, which in turn is proposed to cause rapid initiation.

Diradical character and nonlinear optical properties of buckyferrocenes: focusing on the use of suitably modified fullerene fragments

Buckyferrocenes, made up of suitably modified fullerenes fragments having an uninterrupted π-electron network, show intermediate diradical character and robustly large third-order nonlinear optical responses.

Computational investigation on redox-switchable nonlinear optical properties of a series of polycyclic p-quinodimethane molecules

The polycyclic p-quinodimethanes are proposed to be the novel candidates of the high-performance nonlinear optical (NLO) materials because of their large third order polarizabilities (γ). We investigate the switchable NLO responses of a series of polycyclic p-quinodimethanes with redox properties by employing the density functional theory (DFT). The polycyclic p-quinodimethanes are forecasted to exhibit obvious pure diradical characters because of their large y 0 index (the y 0 index is a value between 0 [closed-shell state] and 1 [pure biradical state]). The γ values of these polycyclic p-quinodimethanes and their corresponding one-electron and two-electron reduced/oxidized species are calculated by the (U)BHandHLYP method. The γ values of polycyclic p-quinodimethanes and their corresponding one-electron reduced species are all positive and significantly different. The large differences of the γ values are due to a change in the transition energy and are related to the different delocalization of the spin density, which demonstrates that the NLO switching is more effective on one-electron reduction reactions. Therefore, the study on these polycyclic p-quinodimethanes provides a guideline for a molecular design of highly efficient NLO switching.