Theoretical Study on the Open-Shell Electronic Structure and Electron Conductivity of [18]Annulene as a Molecular Parallel Circuit Model

Herein, the electron conductivities of [18]annulene and its derivatives are theoretically examined as a molecular parallel circuit model consisting of two linear polyenes. Their electron conductivities are estimated by elastic scattering Green’s function (ESGF) theory and density functional theory (DFT) methods. The calculated conductivity of the [18]annulene does not follow the classical conductivity, i.e., Ohm’s law, suggesting the importance of a quantum interference effect in single molecules. By introducing electron-withdrawing groups into the annulene framework, on the other hand, a spin-polarized electronic structure appears, and the quantum interference effect is significantly suppressed. In addition, the total current is affected by the spin polarization because of the asymmetry in the coupling constant between the molecule and electrodes. From these results, it is suggested that the electron conductivity as well as the quantum interference effect of π-conjugated molecular systems can be designed using their open-shell nature, which is chemically controlled by the substituents.


Introduction
Microfabrication technology for integrated circuits has been developed according to Moore's law proposed in 1965 [1,2].The current gate length of transistors has reached less than 10 nm [3].This miniaturization technology of transistors has increased integration density, which drastically improves the performance of electron devices.On the other hand, it has been suggested that the miniaturization of silicon (Si)-based devices is almost reaching its limit.In addition, energy loss due to leakage currents induced by the miniaturization has been a problem.Recent developments in the fabrication technology of power semiconductors using silicon carbide (SiC) and gallium nitride (GaN) have succeeded in realizing higher efficiency and lower power consumption than those of Si-based semiconductors [4][5][6].As another approach to the problem, molecular electronics, i.e., the use of single molecules for electronic components, has attracted much attention in terms of new nanomaterials as a "bottom-up approach" [7][8][9][10][11].In 1974, a molecular rectifier was first proposed by Aviram and Ratner as molecular electronics [12].
For the use of these molecules as devices, one is required to integrate single-molecule components and to make molecular circuits.Recently, for example, integrated singlemolecule components based on self-assembled monolayers (SAMs) have been proposed, demonstrating a variety of functionalities [44,45].In addition, various molecular devices, such as a biosensor that detects molecular interactions, have been reported [46].Interestingly, it has also been proposed that the non-linear quantum conduction of single molecules can be applied to neuromorphic computers because of its similarity to electrical signals in human brains [47,48].
In recent years, some interesting results on molecular circuits have been reported.For example, Ohm's law, which indicates that the current value of a parallel circuit is equal to the total current through each resistor, is expected to fail in molecular parallel circuits because of quantum interference effects [49].The conductance, G, in nano-scale parallel circuits of two components is given as where G 1 and G 2 are the conductance of respective components and 2 √ G 1 G 2 is the quantum interference term [49].Vazquez et al. indicated that the conductance of a model molecule consisting of two parallel benzenes is larger than twice that of a single benzene, both experimentally and theoretically [50].Recently, some exceptions to Equation (1) have been reported [51,52]; therefore, Equation ( 1) is still under discussion in terms of aromaticity, frontier orbital theory, and orbital interaction [51].Concerning open-shell systems, the relationship between the conductivity of a molecular circuit, quantum interference, and open-shell nature has not yet been discussed.
From those points of view, in this study, we aimed to calculate the electron conductivity of molecular circuit models and to compare this with the conductivity expected from the classical Ohm's law.We focused on [18]annulene as the molecular parallel circuit model, which was assumed to be a parallel circuit that can be divided into two (top and bottom) linear polyenes, as illustrated in Figure 1.The electronic structures and the electron conductivities of the model structures were calculated by density functional theory (DFT) and elastic scattering Green's function (ESGF) methods.

Electron Conductivity Calculation of Single-Molecule
In this study, the electron conductivity of the single molecule, which had discrete orbitals, was calculated using the elastic scattering Green's function proposed by Mujica and Luo [53][54][55].As described below, some model molecules are open-shell systems; therefore, we used the extended Mujica and Luo method for the open-shell systems proposed by Nakanishi et al. [40].In the extended Mujica and Luo method, the Hamiltonian of the system is written as where  ,  , and  are the Hamiltonians of the molecule, the left-sided, and right-sided electrodes, respectively;  is the interaction potential between the molecule and the electrodes; and  , is the interaction between the I-th site of the molecule and the i-th orbital of the electrodes with spin σ (=α or β).Additionally, the transition operator is defined as where  is the Green's function, Figure 1.Molecular parallel circuit models: (a) [18]annulene derivative as the parallel circuit, and (b) two polyenes divided from the [18]annulene derivative.

Electron Conductivity Calculation of Single-Molecule
In this study, the electron conductivity of the single molecule, which had discrete orbitals, was calculated using the elastic scattering Green's function proposed by Mujica and Luo [53][54][55].As described below, some model molecules are open-shell systems; therefore, we used the extended Mujica and Luo method for the open-shell systems proposed by Nakanishi et al. [40].In the extended Mujica and Luo method, the Hamiltonian of the system is written as where H M , H L , and H R are the Hamiltonians of the molecule, the left-sided, and rightsided electrodes, respectively; U is the interaction potential between the molecule and the electrodes; and γ iI,σ is the interaction between the I-th site of the molecule and the i-th orbital of the electrodes with spin σ (=α or β).Additionally, the transition operator is defined as where G is the Green's function, Here, z is a complex variable.Assuming that the electrodes only interact directly with the end-sites, i.e., site 1 and N of the molecule (see Figure 2) [53][54][55], the transition matrix can be written as , (10) where φ η σ is the eigenstate of the total Hamiltonian H (H ).Since the end-sites are defined as the anchor atoms connected to the electrodes, the eigenstate, η, that overlaps with the end-sites only contributes to matrix elements in the Green's function.
Here, ϕ η σ is approximated by the orbitals obtained from the Kohn-Sham equation for the finite systems consisting of the molecule sandwiched between the electrodes (H [56].The conduction current density, i LR σ , of the system under the applied voltage V D by the electrodes is given as where the constants e, m e , k B , and ℏ are the elementary charge, the electron mass, the Boltzmann constant, and the reduced Plank constant, respectively.T is the temperature of the system and f σ (E) is the Fermi distribution.E F,σ is the Fermi energy, defined as the intermediate value between the orbital energies of the HOMO and the LUMO of the extended molecule consisting of the molecule and the electrodes.Assuming that the interactions between the different scattering channels are negligible because the spacing between the molecular orbitals is large enough, then the transition probability is written as where Γ 2 η,σ is the spin-dependent escape rate determined by Fermi's golden rule, Here, 1 ϕ  [37,56].Thus, the γ between the left electrode and site 1 is written as where V L,σ (LUMO) is the interaction between the LUMO of the molecule and the HOMO of the electrode at spin σ, and ∆E σ,LUMO is the energy difference between them.∆E σ,HOMO-LUMO is the HOMO-LUMO gap of the extended molecule and d 2 1,σ (LUMO) is the ratio of the electron density at end-site 1 to that of the entire molecule.In other words, V L,σ (LUMO) expresses the strength of the junction between the molecule and the electrodes, d * 1,σ (LUMO) represents the degree of transmission from the electrode to the molecule, and their product, γ L1,σ (LUMO), determines the electron conductivity derived from the electron transfer between the electrodes and the molecule.Equations ( 15)-( 17) can also be applied to the interaction between the end-site N and the right electrode with subscripts N and R.
In this way, the current can be calculated by considering transition probabilities based on tunneling through the discrete molecular orbitals and their energies without an external electric field.The total conduction is I LR = Ai LR , where A is the effective injection area where the electrons transfer from the electrodes.Here, it is approximated that A ≈ πr 2 s , where r s is the radius of a sphere with the same volume as an electron, r s = [3/(4πn)] 1/3 , and n is the density of the electron, n = (2m e E F,σ ) 3/2 / 3π 2 ℏ 3 [37].
Notably, the relationship between the electron conductivities and the quantum interference of the derivatives is discussed through the following formula, instead of Equation ( 1) under the same voltage in this paper.

Computational Details
In this study, [18]annulene derivatives were used as models of the molecular parallel circuit, which was divided into two linear polyenes (top: Polyene A; bottom: Polyene B), as illustrated in Figure 1.The [18]annulene derivatives and the two polyenes were connected to the gold electrodes via thioketone groups, as illustrated in Figure 2. The anchor sulfur atoms were assumed to be connected to the bridge-site of the Au(111) surface, and the electrodes were approximated as gold dimers [40].Several [18]annulene derivatives that involve electron-donating/withdrawing groups (X) were also examined together with the non-substituted (NS) derivative.Since some computational models show open-shell nature, as explained below, density functional theory (DFT) calculations with the spinunrestricted approximation that enabled the spin polarization to be performed involved the static electron correlation effect.In addition, BHandHLYP [57], which is reported to be effective for spin-polarized open-shell systems, was utilized as the functional for all calculations [58,59].
These [18]annulene derivatives were first geometrically optimized by the BHandHLYP/ 6-31+G*(S), 6-31G*(other atoms) [57,[60][61][62] level of theory, and then they were divided into the two polyenes (i.e., Polyene A and Polyene B).Only the hydrogen atoms added at the boundaries of the polyene models were optimized after the partition.The Cartesian coordinates of these model molecules are summarized in Table S1 in the Supplementary Materials.For these models, the electronic structures were obtained by BHandHLYP/LANL2DZ [63] (Au), 6-31+G* (S) and 6-31G* (other atoms) levels of theory.For all models, the charge and the spin states were assumed to be neutral and singlet, respectively.
The electron conductivity was calculated through ESGF methods using the DFT results, assuming that the temperature was 300 K and the intramolecular electron conduction only occurred through the molecular orbitals from LUMO+9 to HOMO−9 [40].All DFT calculations were performed in the gas phase using Gaussian09 Rev. D01 [64]; the electron conductivity was simulated using our self-developed program.
The electron conductivity was calculated through ESGF methods using the DFT results, assuming that the temperature was 300 K and the intramolecular electron conduc tion only occurred through the molecular orbitals from LUMO+9 to HOMO−9 [40].Al DFT calculations were performed in the gas phase using Gaussian09 Rev. D01 [64]; the electron conductivity was simulated using our self-developed program.

Conductivity of [18]Annulene: The Molecular Parallel Circuit Model
Figure 3 shows the current-voltage (I-V) characteristics of [18]annulene ( annulene ) The current increased non-linearly with the increase in voltage, showing typical quantum electron conductivity characteristics.For the comparison, a total current of two independ ent polyenes (Polyene A and Polyene B) ( Polyene A +  Polyene B =  A+B ) and an estimated cur rent for quantum systems, i.e.,  A+B plus the quantum interference (QI) term ( Polyene A +  Polyene B + 2  Polyene A  Polyene B =  A+B+Q ) were also plotted.As shown in the figure, the di rectly calculated I-V characteristics of the [18]annulene were larger than the simple sum ( A+B ), but smaller than the value including the QI term ( A+B+Q ).To consider the difference in the conductivity between the [18]annulene and the polyenes ( annulene and  A+B ,  A+B+Q ) several parameters, i.e., total spin angular momentum (〈 〉) and the coupling constants such as  ,  , are summarized in Table 1.In addition, Site-overlap values of each mo lecular orbital are also shown in Table 2.The Site-overlap values indicate that the mos effective orbital to conductivity is HOMO−5 in annulene, whose electron distribution is illustrated in Figure 4a.

Conductivity of [18]Annulene: The Molecular Parallel Circuit Model
Figure 3 shows the current-voltage (I-V) characteristics of [18]annulene (I annulene ).The current increased non-linearly with the increase in voltage, showing typical quantum electron conductivity characteristics.For the comparison, a total current of two independent polyenes (Polyene A and Polyene B) (I Polyene A + I Polyene B = I A+B ) and an estimated current for quantum systems, i.e., I A+B plus the quantum interference (QI) term (I Polyene A + I Polyene B + 2 I Polyene A I Polyene B = I A+B+Q ) were also plotted.As shown in the figure, the directly calculated I-V characteristics of the [18]annulene were larger than the simple sum (I A+B ), but smaller than the value including the QI term (I A+B+Q ).To consider the difference in the conductivity between the [18]annulene and the polyenes (I annulene and I A+B , I A+B+Q ), several parameters, i.e., total spin angular momentum (⟨ Ŝ2 ⟩) and the coupling constants such as γ L1 , γ NR , are summarized in Table 1.In addition, Site-overlap values of each molecular orbital are also shown in Table 2.The Site-overlap values indicate that the most effective orbital to conductivity is HOMO−5 in annulene, whose electron distribution is illustrated in Figure 4a.
sults, assuming that the temperature was 300 K and the intramolecular electron conduction only occurred through the molecular orbitals from LUMO+9 to HOMO−9 [40].All DFT calculations were performed in the gas phase using Gaussian09 Rev. D01 [64]; the electron conductivity was simulated using our self-developed program.

Conductivity of [18]Annulene: The Molecular Parallel Circuit Model
Figure 3 shows the current-voltage (I-V) characteristics of [18]annulene ( annulene ).The current increased non-linearly with the increase in voltage, showing typical quantum electron conductivity characteristics.For the comparison, a total current of two independent polyenes (Polyene A and Polyene B) ( Polyene A +  Polyene B =  A+B ) and an estimated current for quantum systems, i.e.,  A+B plus the quantum interference (QI) term ( Polyene A +  Polyene B + 2  Polyene A  Polyene B =  A+B+Q ) were also plotted.As shown in the figure, the directly calculated I-V characteristics of the [18]annulene were larger than the simple sum ( A+B ), but smaller than the value including the QI term ( A+B+Q ).To consider the difference in the conductivity between the [18]annulene and the polyenes ( annulene and  A+B ,  A+B+Q ), several parameters, i.e., total spin angular momentum (〈 〉) and the coupling constants such as  ,  , are summarized in Table 1.In addition, Site-overlap values of each molecular orbital are also shown in Table 2.The Site-overlap values indicate that the most effective orbital to conductivity is HOMO−5 in annulene, whose electron distribution is illustrated in Figure 4a.[18]annulene ( annulene ), the sum of independent polyenes ( A+B ), and the values including the quantum interference (QI) term ( A+B+Q ).[18]annulene (I annulene ), the sum of independent polyenes (I A+B ), and the values including the quantum interference (QI) term (I A+B+Q ).

Table 1. Summary of ⟨
Ŝ2 ⟩ and coupling constants of the [18]annulene and the divided independent polyenes.In this table, α and β represent α and β orbitals, respectively.Other parameters related to the transition probability, such as orbital energy, Fermi energy, and so on, are shown in Tables S2 and S3, and Figure S1 in the Supplementary Materials.
[ On the other hand, one can find significant difference in ⟨ Ŝ2 ⟩ values between [18]annulene and polyenes.The ⟨ Ŝ2 ⟩ values represent whether the molecule is the closed-shell [66].From the calculated ⟨ Ŝ2 ⟩ values in Table 1, the [18]annulene is found to be a closed-shell system, whereas the divided polyenes are weak open-shell systems.Moreover, according to Tables 1 and 2, the coupling constants and Site-overlap values of α and β orbitals are the same in the [18]annulene, while they are different in the divided polyenes.In the closed-shell systems, the α and β orbital distributions are spatially similar to each other.On the other hand, in the open-shell systems, the α and β orbitals show spatially different distributions, indicating spin-polarized electronic states.In other words, α and β orbitals tend to localize at the left and the right sides of the polyenes, respectively.As a result, the parameters of the α and β orbitals are the same in the [18]annulene that has the closed-shell electronic structure, in contrast, the divided polyenes show different parameter values between the α and β orbitals due to their open-shell electronic structures.Therefore, the results indicate that the electron is not easily transmitted through the molecules in the polyenes because of the spin polarization.Indeed, the parameters related to the conductivity are different between the left and right sides (γ L1,σ and γ NR,σ ); thus, the conductivity of the polyenes is decreased in comparison with that of the annulene.On the other hand, I A+B+Q that fully involves the QI term becomes larger than I annulene , suggesting that the quantum interference must be examined partially in the systems with open-shell nature.

Effects of Substituents
As explained in the previous section, the electron conductivity and the QI term are affected by the open-shell nature of the molecules.The result suggests a possibility that the electron conductivity of molecules can be controlled by changing the degree of spin polarization with the introduction of the substituent groups (X).From this point of view, the relationship between the open-shell nature and the conductivity is examined by introducing electron-donating (X = -OCH3, -OH) and -withdrawing (X = -CN, -NO2) substituents to the [18]annulene and the divided polyenes.Here, a pair of substituents were introduced symmetrically to carbon atoms with a larger distribution in HOMO−5, which dominantly contributes to the Site-overlap value in the non-substituted (NS) annulene, as shown in Figure 4a. Figure 4b shows the current values of the substituted annulenes and polyenes at 1.0 V.The introduction of the substituents significantly decreases the conductivity for both annulenes and polyenes, regardless of the electron-donating/withdrawing groups.However, the relationship in conductivity between the annulenes and the sum of two polyenes is different between the electron-donating and -withdrawing groups.In the case of the electron-donating groups, the current values of the annulenes are close to the total current values of the polyenes, including the QI term ( annulene ≈  A+B+Q ), whereas they become comparable to the simple sum of the current values of the two polyenes in the case of the withdrawing groups ( annulene ≈  A+B ).This suggests that the effect of the QI term depends on the π electronic states in the annulene ring.
To further understand the effect, we examined the 〈 〉 values.In the case of the NS model and its derivatives with electron-donating groups (-OCH3 and -OH), the calculated 〈 〉 values are zero, although their polyenes show non-zero values, as summarized in Table 3.Therefore, the NS annulene and its derivatives with electron-donating groups are closed-shell systems, while the corresponding polyenes are open-shell systems.In contrast, in the case of the derivatives with electron-withdrawing groups (-CN and -NO2), both the derivatives and their polyenes show non-zero 〈 〉 values, indicating the openshell electronic structures.In order to confirm the effect of the open-shell nature in detail,

Effects of Substituents
As explained in the previous section, the electron conductivity and the QI term are affected by the open-shell nature of the molecules.The result suggests a possibility that the electron conductivity of molecules can be controlled by changing the degree of spin polarization with the introduction of the substituent groups (X).From this point of view, the relationship between the open-shell nature and the conductivity is examined by introducing electron-donating (X = -OCH 3 , -OH) and -withdrawing (X = -CN, -NO 2 ) substituents to the [18]annulene and the divided polyenes.Here, a pair of substituents were introduced symmetrically to carbon atoms with a larger distribution in HOMO−5, which dominantly contributes to the Site-overlap value in the non-substituted (NS) annulene, as shown in Figure 4a. Figure 4b shows the current values of the substituted annulenes and polyenes at 1.0 V.The introduction of the substituents significantly decreases the conductivity for both annulenes and polyenes, regardless of the electron-donating/withdrawing groups.However, the relationship in conductivity between the annulenes and the sum of two polyenes is different between the electron-donating and -withdrawing groups.In the case of the electron-donating groups, the current values of the annulenes are close to the total current values of the polyenes, including the QI term (I annulene ≈ I A+B+Q ), whereas they become comparable to the simple sum of the current values of the two polyenes in the case of the withdrawing groups (I annulene ≈ I A+B ).This suggests that the effect of the QI term depends on the π electronic states in the annulene ring.
To further understand the effect, we examined the ⟨ Ŝ2 ⟩ values.In the case of the NS model and its derivatives with electron-donating groups (-OCH 3 and -OH), the calculated ⟨ Ŝ2 ⟩ values are zero, although their polyenes show non-zero values, as summarized in Table 3.Therefore, the NS annulene and its derivatives with electron-donating groups are closed-shell systems, while the corresponding polyenes are open-shell systems.In contrast, in the case of the derivatives with electron-withdrawing groups (-CN and -NO 2 ), both the derivatives and their polyenes show non-zero ⟨ Ŝ2 ⟩ values, indicating the open-shell electronic structures.In order to confirm the effect of the open-shell nature in detail, the γ L1 , γ NR and the Site-overlap totals of each derivative are summarized in Table 4.In addition, the parameters related to the transition probability, such as the Siteoverlap of each molecular orbital, the HOMO-LUMO gap, coupling constants, etc., are summarized in Table S2; the orbital energy is shown in Table S3 and Figure S1 in the By summarizing these results, we can suggest that: (i) In the case where annulene derivatives have closed-shell electronic structures, the QI term becomes significant (I annulene ≈ I A+B+Q ).(ii) In the case where the annulene derivatives have open-shell electronic structures, the QI term is suppressed; therefore, the current of the annulenes is close to the sum of the current values of two polyenes (I annulene ≈ I A+B ).This is similar to Ohmic behavior.In other words, the spin polarization, i.e., a separation between α and β orbitals, avoids the quantum interference in the molecular parallel circuit, and it can be realized by the introduction of an electron-withdrawing group.Moreover, in such spin-polarized systems, the total current also decreased due to the asymmetry in the coupling constants and Site-overlap values.
Finally, a reason why the introduction of donating groups as well as both the electronwithdrawing groups decreases the conductivity in comparison with the NS model is explained here.As mentioned above, the total current of the annulene derivatives with the electron-withdrawing groups decreases due to the spin polarization.In contrast, a difference in the molecular orbital distribution, especially HOMO−5, becomes important for models with electron-donating groups.The distribution in HOMO−5 of the NS model and substituted annulenes of the electron-donating groups are shown in Figure 5.The introduction of the electron-donating groups suppresses a degree of delocalization of HOMO−5 that contributes to the conductivity.This decrease is explained as follows: the π orbital energies of the annulene ring are shifted by the substitution effect, so that the energy difference between the π orbitals of the annulene ring and junction moieties (-S-Au 2 ) becomes larger.As a result, the Site-overlap values become small due to a decrease in conjugation between the ring and junction moieties.
(i) Sum of HOMO−9 and LUMO+9.(ii) The small differences in Site-overlap between Polyenes A and B are considered to originate from the convergence accuracy of the molecular orbitals (here, 10 −7 a.u. in electron density).
By summarizing these results, we can suggest that: (i) In the case where annulene derivatives have closed-shell electronic structures, the QI term becomes significant ( annulene ≈  A+B+Q ).(ii) In the case where the annulene derivatives have open-shell electronic structures, the QI term is suppressed; therefore, the current of the annulenes is close to the sum of the current values of two polyenes ( annulene ≈  A+B ).This is similar to Ohmic behavior.In other words, the spin polarization, i.e., a separation between α and β orbitals, avoids the quantum interference in the molecular parallel circuit, and it can be realized by the introduction of an electron-withdrawing group.Moreover, in such spin-polarized systems, the total current also decreased due to the asymmetry in the coupling constants and Site-overlap values.
Finally, a reason why the introduction of donating groups as well as both the electron-withdrawing groups decreases the conductivity in comparison with the NS model is explained here.As mentioned above, the total current of the annulene derivatives with the electron-withdrawing groups decreases due to the spin polarization.In contrast, a difference in the molecular orbital distribution, especially HOMO−5, becomes important for models with electron-donating groups.The distribution in HOMO−5 of the NS model and substituted annulenes of the electron-donating groups are shown in Figure 5.The introduction of the electron-donating groups suppresses a degree of delocalization of HOMO−5 that contributes to the conductivity.This decrease is explained as follows: the π orbital energies of the annulene ring are shifted by the substitution effect, so that the energy difference between the π orbitals of the annulene ring and junction moieties (-S-Au2) becomes larger.As a result, the Site-overlap values become small due to a decrease in conjugation between the ring and junction moieties.

Conclusions
In this study, we examined the electron conductivity of [18]annulene and its derivatives with electron-donating/-withdrawing groups as models for molecular parallel circuits.In the case of the NS model and its derivatives with electron-donating groups, the electronic structures are closed-shell.Furthermore, their current values are larger than those of the sum of the two divided polyenes, indicating the importance of quantum interference effect in these systems.In the derivatives with the electron-withdrawing groups, however, the quantum interference effect is suppressed by the spin polarization of the open-shell nature.
Thus, it is suggested that the quantum interference effect of the π-conjugated molecular parallel circuit can be designed with an open-shell nature, which is chemically controlled by the substituent groups.The open-shell systems with spin polarization due to the strong electron correlation have been a focus because of their applicability to functional materials [67,68].However, a relationship between their open-shell nature and conductivity has not been clarified.Concerning the electron conductivity, our group has examined the relationship between electronic/spin structures and the electron conductivity of onedimensional polynuclear nickel complexes through theoretical calculations [41][42][43].Such extended metal atom chains (EMACs) are considered to be promising candidates for the molecular wire [69,70]; however, we found that they become candidates for single molecular switches controlled by an external magnetic field just by introducing the openshell nature [41][42][43].Thus, open-shell systems have great potential for the development of highly functional single molecular devices, which cannot be realized by closed-shell systems.Recently, it has been suggested that brain-like computers can be realized by combining molecular circuits and machine leaning [48].Furthermore, quantum interference has been proposed to be important for making devices [71].Therefore, for the rational design of high-performance molecular circuits, we need to determine the relationship between the open-shell nature, electron conductivity, and quantum interference.
In this study, we found a possibility that the quantum interference effect in the molecular circuit can also be controlled by chemical substituents using the open-shell nature.To the best of our knowledge, this is the first report that discusses a relationship between quantum interference and the open-shell nature.

η σ and ϕ η σ N
denote the site-orbital overlap matrix elements between the end-sites and the extended molecule.Their product 1 , which represents the delocalization of the molecular orbitals of the extended molecule; γ is the coupling constant, as explained below.Luo et al. proposed that the occupied molecular orbitals of electrodes interact with the LUMO of the molecule based on the frontier orbital theory

Figure 2 .
Figure 2. Calculated model systems.Electrodes are approximated by the gold dimers.The labels 1 and N represent the sulfur atoms at end-sites; L and R are the left and right electrodes, respectively (see the equations presented in Section 2.1).The distances between the middle point of the Au dimer and the sulfur atom were fixed at 2.30 Å, and the Au-Au distance was assumed to be 2.88 Å[37,65]

Figure 3 .
Figure 3. I-V characteristics of the[18]annulene ( annulene ), the sum of independent polyenes ( A+B ) and the values including the quantum interference (QI) term ( A+B+Q ).

Figure 2 .
Figure 2. Calculated model systems.Electrodes are approximated by the gold dimers.The labels 1 and N represent the sulfur atoms at end-sites; L and R are the left and right electrodes, respectively (see the equations presented in Section 2.1).The distances between the middle point of the Au dimer and the sulfur atom were fixed at 2.30 Å, and the Au-Au distance was assumed to be 2.88 Å[37,65].

Figure 2 .
Figure 2. Calculated model systems.Electrodes are approximated by the gold dimers.The labels 1 and N represent the sulfur atoms at end-sites; L and R are the left and right electrodes, respectively (see the equations presented in Section 2.1).The distances between the middle point of the Au dimer and the sulfur atom were fixed at 2.30 Å, and the Au-Au distance was assumed to be 2.88 Å[37,65].

Figure 3 .
Figure3.I-V characteristics of the[18]annulene (I annulene ), the sum of independent polyenes (I A+B ), and the values including the quantum interference (QI) term (I A+B+Q ).

Figure 4 .
Figure 4. (a) Molecular orbital distribution of HOMO−5 of the [18]annulene extended molecule.The isovalue of the distribution is 0.02.Substituents are introduced at the carbons indicated by arrows.(b) Current of the [18]annulene derivatives (I annulene ) and the divided polyenes (I A+B = I Polyene A + I Polyene B , I A+B+Q = I Polyene A + I Polyene B + 2 I Polyene A I Polyene B ) at 1.0 V.

Figure 5 .
Figure 5. Electron distribution in HOMO−5 of (a) NS annulenes, (b) OCH3-substituted annulenes, and (c) OH-substituted annulenes, where the isovalues are 0.02.Other molecular orbitals of all models are summarized in Figure S2 in the Supplementary Materials.

Figure 5 .
Figure 5. Electron distribution in HOMO−5 of (a) NS annulenes, (b) OCH 3 -substituted annulenes, and (c) OH-substituted annulenes, where the isovalues are 0.02.Other molecular orbitals of all models are summarized in Figure S2 in the Supplementary Materials.

Table 2 .
[18]-overlap values of each orbital for the[18]annulene and their divided independent polyenes.Total values of each orbital are also indicated in the same table.

Table 4 .
Summary of the coupling constants and total of Site-overlap (i),(ii) of each of the extended molecules.(a) OCH 3 -substituted, (b) OH-substituted, (c) CN-substituted, and (d) NO 2 -substituted molecules.α and β in the table represent α and β orbitals, respectively.Sum of HOMO−9 and LUMO+9.(ii) The small differences in Site-overlap between Polyenes A and B are considered to originate from the convergence accuracy of the molecular orbitals (here, 10 −7 a.u. in electron density).