SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Quasi-classical modeling of molecular quantum-dot cellular automata multidriver gates

Ehsan Rahimi* and Shahram Mohammad Nejad

Author affiliations

Nanoptronics Research Center, School of Electrical and Electronic Engineering, Iran University of Science and Technology, Tehran, 16844, Iran

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:274  doi:10.1186/1556-276X-7-274

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/274

Received:29 March 2012
Accepted:25 April 2012
Published:30 May 2012

© 2012 Rahimi and Mohammad Nejad; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Molecular quantum-dot cellular automata (mQCA) has received considerable attention in nanoscience. Unlike the current-based molecular switches, where the digital data is represented by the on/off states of the switches, in mQCA devices, binary information is encoded in charge configuration within molecular redox centers. The mQCA paradigm allows high device density and ultra-low power consumption. Digital mQCA gates are the building blocks of circuits in this paradigm. Design and analysis of these gates require quantum chemical calculations, which are demanding in computer time and memory. Therefore, developing simple models to probe mQCA gates is of paramount importance. We derive a semi-classical model to study the steady-state output polarization of mQCA multidriver gates, directly from the two-state approximation in electron transfer theory. The accuracy and validity of this model are analyzed using full quantum chemistry calculations. A complete set of logic gates, including inverters and minority voters, are implemented to provide an appropriate test bench in the two-dot mQCA regime. We also briefly discuss how the QCADesigner tool could find its application in simulation of mQCA devices.

Electron transfer reactions; Molecular electronics; Molecular gates; Molecular quantum-dots; Quantum cellular automata


Recent advances in molecular electronics on one hand and the limitations of conventional semiconductor devices, on the other hand, have driven a surge of activities towards the realization of molecular devices, circuits, and systems. Achieving the ultimate diminution in size, power consumption, and delay of electronic devices and systems has always been a challenging endeavor of scientists and designers in this field. Due to the prospect of size reduction in electronics offered by molecular-level control of properties, molecular electronics provides means to extend the Moore's law beyond the foreseen limits of small-scale conventional silicon integrated circuits. The small size of molecules allows high device density in the range of 1011 to 1012 devices/cm2[1]. Besides, the chemical self-assembly capacity in manufacturing molecular devices provides many advantages to conventional semiconductor manufacturing technology, including lower manufacturing cost and uniform device reproducibility. Molecular electronics endeavors to use the nonlinear current–voltage characteristics of individual molecules or molecular assemblies as active devices (diodes, transistors, etc.) in electronic circuits. However, the power consumption of molecular current switches at very high frequencies is still a drawback [2]. The π-σ-π mixed-valence type molecules, which provide double-well potentials for electrons, have been proposed and studied by Aviram towards the synthesis of memory, logic, and amplification [3]. Lent proposed using molecules in representing binary information within the molecular quantum-dot cellular automata (mQCA) paradigm [1,4]. Molecular QCA provides an alternative approach to represent and process data, where binary representation lies in the charge configuration within molecules rather than in the on/off states of current switches. A cell in the mQCA model consists of a number of molecular quantum dots (or redox centers) and a few electrons. The electrons tend to occupy antipodal sites as a result of Coulomb repulsion. The Columbic interactions cause electrons to tunnel from one redox center to another in a cell, but not between cells. Thus, it is likely that no current flows, since the neighboring cells are coupled by electrostatic field. Figure 1 depicts two-dot and four-dot mQCA cells and how binary “1” or “0” is represented. The first QCA device was implemented and tested using metal dots at near 0 K [5]. Semiconductor implementation of QCA using GaAs/AlGaAs heterostructure materials has been reported in [6,7] as well. Molecules are good containers for keeping electric charge and mQCA cells have a more promising future to work at room temperature [8]. Nonbonding π or d orbitals of a single molecule (or multiple molecules) can function as quantum dots, where the electric charge is localized in each cell. Synthesis of two-dot and four-dot mixed-valence candidate molecules for mQCA has been reported in [9-14]. Many of these molecules are mixed-valence type and include transition metals to enable fast electron transfer reactions [15]. Molecular QCA gates are the building blocks of circuits in this paradigm. Calculation of the electronic structure of mQCA gates composed of these molecules is challenging, since the number of basis set functions grows exponentially as the number of molecules and atoms are increased. Besides, many of the ab initio methods fail in describing charge distribution in mixed-valence complexes. Therefore, developing semi-classical models to study mQCA gates is of high importance. Currently QCADesigner [16] utilizes nonlinear and two-state approximations to solve metallic-based QCA circuits. Several QCA circuits including combinational as well as sequential circuits have been studied using QCADesigner. Examples are adders, shift registers, RAM, digital data storage, and simple microprocessors [17-26]. In this paper, we center on two-dot mQCA and analyze the validity and accuracy of the two-state model approximation for studying multidriver mQCA gates. This study provides an approach to enhance the QCADesigner tool for simulation of mQCA devices in the future.

thumbnailFigure 1. Binary representation in the mQCAparadigm. (a) Schematic structure of a two-dot mQCA cell. (b) A two-dot molecule. (c) In two-dot cells, depending on which of the upper (A) or lower (B) quantum-dot is occupied, binary “1” or “0” is represented. (d) In four-dot cells, binary “1” and “0” is represented within the occupation of AB’ or A’B dots correspondingly.


Two-dot molecular QCA test bench

The majority voter (MV) and the inverter (INV) gates [25] are the fundamental building blocks of any circuit in the four-dot QCA architecture. These gates have been schematically shown in Figure 2a, b. Particularly, the MV gate is referred to as the universal QCA gate, since the AND and OR logical operations can be done by this gate, as evident from the truth table shown in Figure 2c.

thumbnailFigure 2. Four-dot QCA gates. (a) The universal MV gate. The majority of the three fixed inputs, which is “0” in this figure, appears at the output as a result of Coulombic interactions and minimum energetics. (b) The inverter gate. (c) Truth table of the MV gate. When d3 = ”0”, the MV gate performs AND logical operations on d1 and d2, and when d3 = ”1” the MV gate functions as a two-input OR gate [25].

Our multidriver minority voter (MinV) gate is composed of m drivers, where m is an odd number, as inputs and one output. Figure 3a schematically illustrates the three-driver MinV model gate in the two-dot mQCA regime. When only one driver (e.g., the d1) is present, the model gate serves as an INV gate (Figure 3b). In the multidriver MinV model gate, all the distances between the centers of the molecules are l, which is equal to the distance between the middle of upper and lower π-bonds as shown in Figure 3c. The inputs of the gates are kept fixed, while the output cells switch to their stable states. To this end, two point charges q and 1-q separated by distance l are used to mimic each input as depicted in Figure 3c.

thumbnailFigure 3. Two-dot QCA gates. (a) Structure of the three-input MinV gate. This gate is composed of three fixed inputs (d1, d2, and d3) and a two-dot molecule (AB) as output. The minority of the three fixed inputs, which is “0” in this figure, appears at the output. (b) When there is only one input, the MinV gate functions as an inverter gate. (c) Two point charges q and 1-q separated by distance l, which is the distance between the middle of upper and lower π-bonds, are used to mimic a fixed input. (d) Truth table of the MinV gate. When d3 = ”0”, the MinV gate performs NAND logical operation on d1 and d2, and when d3 = ”1” the MinV gate functions as a two-input NOR gate.

The MinV gate is an alternative to the MV gate in the four-dot QCA, where the output is inverted. Compared to MV and INV gates in the four-dot architecture, which require 16 and 28 quantum-dots correspondingly, the MinV and INV gates require only 8 and 4 quantum-dots in the two-dot mQCA regime. Consequently, these gates provide a small two-dot mQCA test bench, which make high level quantum chemical calculations feasible. The MinV gate can perform NAND and NOR logical operations, as shown in Figure 3d, and provides a functionally complete logic set to implement any logic function in the two-dot mQCA framework. Additionally, it is possible to implement multi-input (or multidriver) MinV gates, which in turn decrease the total number of gates required to implement a logic circuit. It is important to note that since the MinV gate is not a planar gate, circuits implemented in the two-dot mQCA regime are not planar circuits. We highlight that the practical QCA circuits require clocked-control cells and clocking schemes [21,27-29], which are not addressed in this paper.

Two-state model for molecular QCA gates

The charge configuration in a QCA cell is quantified by the so called ‘polarization’, and is defined as [30]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M1">View MathML</a>


where qAqBqA′ and qB′ are the charges localized at four quantum dots labeled in Figure 1. For two-dot mQCA cells, the polarization is given by the charges qA and qB at the corresponding redox centers in Equation 1. The polarization of a QCA cell varies between −1 and 1, while negative and positive polarizations represent binary “0” and “1”, respectively. In two-dot mQCA cells, the normalized dipole moment of the used two-dot molecule is also identical to the polarization, which is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M2">View MathML</a>


where μα denotes the component of the molecular dipole moment that is parallel to the bridge direction, and the origin is in the middle of the bridge. The dipole moment of an mQCA cell can be obtained through full quantum chemical calculations. An important parameter of a QCA device is the Kink energy (Ek), which is the required energy to excite the system from the ground state to the first excited state. To distinguish a bit value from the thermal environment, Ek must be greater than kBT[31], where T is the operation temperature in Kelvin, and kB is the Boltzmann's constant. The Ek represents the energy cost of cells i and j having opposite polarizations. That is, the electrostatic interaction between all the charges in cells i and j is calculated by [16]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M3">View MathML</a>


where ϵ0 is the permittivity of free space, and ϵr is the relative permittivity of the material system. The Kink energy is then given by [16]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M4">View MathML</a>


where E′i,j and Ei,j denote the electrostatic energy of cells i and j having opposite and same polarizations correspondingly.

Tougaw and Lent have used a simple Hamiltonian of the extended-Hubbard type to describe the dynamic behavior of four-dot metallic-based QCA nanodevices [32]. Although this Hamiltonian describes the dynamics of the coherent system composed of arrays of four-dot QCA cells elegantly in theory, it is only possible to model the small systems employing this scheme, since the total required number of direct-product basis sets grows exponentially with the number of cells. In other words, an array with N number of four-dot cells and B number of basis sets in each cell requires the total number of direct-product basis sets as [32]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M5">View MathML</a>


By reducing the number of basis sets for each cell and picking up the two orthogonal ones, the Hamiltonian of a four-dot QCA wire can be mapped to Ising model as [32,33]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M6">View MathML</a>


where Ek is the kink energy of four-dot cells, γ is the tunneling energy, and σx and σz are Pauli spin matrices. In semiconductor and metal-dot QCA, the tunneling barriers of the cells are connected to electrodes, and their heights are controlled externally by voltage sources [33]. The steady-state polarization of any cell, j in a block of cells, can be obtained as a solution to the Hartree-Fock intercellular approximation. This approximation decouples the line of N cells into N single cell subsystems and assumes that the cells are only coupled through expectation values of polarizations. The consequent solution is [33],

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M7">View MathML</a>


where Pj is the sum of the polarizations of the neighboring four-dot QCA cells. Equation 7 is currently used in the nonlinear and two-state simulation engine of QCADesigner to solve the metallic-based QCA circuits. It is important to note that mQCA utilizes non-abrupt clocking to reduce the probability of Kink, the property that is not currently present in the QCADesigner as it is based on metallic QCA. In mQCA, the tunneling barriers can be controlled by external electric field [27]. It is demanding to enhance the tool to be able to simulate mQCA circuits. As a primary step towards this end, we present how a similar equation to (7) can be derived directly from the two-state approximation in electron transfer theory [34,35] for two-dot mQCA. We then discuss how these approximations affect the results compared to those obtained from full quantum chemistry calculations.

Equation 8 describes the electron transfer (ET) process in a two-dot mQCA cell, where the two redox centers, A and B, are linked through an intervening bridge, I.

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M8">View MathML</a>


The electronic coupling between the redox centers, which is time independent, is an important factor in the ET process. Within the two-state approximation, the Landau-Zener model [36,37] for avoided crossing of energy surfaces may be applied, where the two diabatic states “1” and “0” denoted by ψa and ψb, and with energies Haa and Hbb interact. The ground state ψ1 and the first excited state of a QCA cell, ψ2, can be related to the diabatic states ψa and ψb by a unitary transformation [34]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M9">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M10">View MathML</a>


In Equations 9 and 10ψ1 and ψ2 are orthonormal, whereas ψa and ψb are not orthogonal in general. The correspondence between diabatic and adiabatic two-state models arises from the secular determinant (Sab= <ψab> is neglected) [38]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M11">View MathML</a>


where E is the adiabatic energy eigenvalue. The η in Equations 9 and 10 satisfies [38]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M12">View MathML</a>


The energy difference between the two diabatic states in the output cell of the MinV gate can be approximated by calculating the difference between the electrostatic energies of the gate for the two output configurations, where the unit charge is localized at sites A and B correspondingly. Using Equation 3, for the “1” and “0” output states (Figure 3a), we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M13">View MathML</a>


Inserting Equation 13 into Equation 12 we have

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M14">View MathML</a>


The Kink energy of two-dot cells can be calculated from Equation 3 and 4 for two neighboring cells as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M15">View MathML</a>


and for each driver, the polarization is defined using Equation 1

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M16">View MathML</a>


Thus, Equation 14 can be rewritten in terms of the Kink energy and the input polarizations

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M17">View MathML</a>


And finally, using Equation 1 and the transformation coefficients in Equations 9 and 10, the output polarization of the MinV gate is obtained

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M18">View MathML</a>


Inserting Equation 17 into Equation 18, we can find the polarization of the output cell of the MinV gate in terms of the polarizations of the inputs straightforwardly as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M19">View MathML</a>


Equation 19 in two-dot mQCA is analogous to Equation 7 in four-dot metal-based QCA, where the tunneling energy γ appears as electronic coupling of redox centers (Hab) in Equation 19. They also imply

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M20">View MathML</a>


The additivity relation in Equation 20 originates from the additivity of electrostatic potential energy in Equation 13 for diabatic states.

Multidriver MinV gates help reduce the number of needed gates for implementation of a logic circuit. Similarly, for an m-input MinV gate, we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M21">View MathML</a>


We refer to Equaiton 21 as the two-state model (TSM) for mQCA gates along this paper. The Ek and the Hab are the only parameters of the TSM. Once the geometrical parameter l is determined, experimentally or from theoretical calculations, the Kink energy can be calculated using Equation 15. The electronic coupling matrix element, Hab, can be calculated using various quantum chemistry techniques [34,38-41] or obtained via spectroscopic experiments, including absorption [42,43], EPR [44], and ultraviolet photoelectron spectroscopy [45]. As we will present, the parameter μ = Ek/2Hab plays an important role in the accuracy of the TSM. It is also the slope of the switching response function at the origin i.e.,

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M22">View MathML</a>


Results and discussion

The chemistry of mixed-valence complexes has received considerable attention recently in mQCA device implementation, where the intramolecular electron transfer and charge localization at redox sites are the important key factors. Mixed-valence compounds contain more than one redox site in the same molecule or molecular unit. Simple model molecules for two-dot cells are the π-σ-π mixed-valence types, which were proposed by Aviram and studied later by Hush [3,46]. In the Aviram's model molecule (1, 4-diallyl butane cation), the two π-allyl groups form two redox centers and are connected by a σ-butyl bridge. One of the allyl groups is a neutral radical, while the other one is anionic (or cationic). The possibility of charge localization in some mixed-valence mQCA candidate molecules has been examined theoretically as well as experimentally [15,47,48]. Advances in quantum chemistry in the past half century provide reliable methods to explore the electronic structure of molecules; however, many of the ab initio techniques fail in describing charge distribution in mixed-valence complexes. The unrestricted Hartree-Fock method overestimates the charge localization due to the neglect of electron correlation effects [49]. In the density functional theory (DFT) method, the exchange potential defined in hybrid functional leads to underestimation of charge localization [47,49]. The complete active space self-consistent field (CASSCF) method [50,51] is believed to be the most reliable for describing charge distribution in mixed-valence complexes [50]. However, the multi-determinant CASSCF calculations scale with the system size, which makes this method highly demanding in computer time and memory. The number of Slater determinants has factorial dependence on both the number of active electrons and particularly on the number of active orbitals generating many-electron configurations (full configuration interaction (CI) within the active space). This is much more significant than any dependence on the number of one-electron basis functions. The number of Slater determinants in a full CI calculation is given by:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M23">View MathML</a>


where M is the number of active orbitals, Nα and Nβ are the numbers of active electrons with α- and β-spins, respectively, and the quantities in parentheses are binomial coefficients:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M24">View MathML</a>


We present full quantum chemistry calculations of the steady-state output polarization of the universal MinV model gate serving as INV and three-input MinV gates. The results based on full quantum chemical calculations are compared to the results obtained from the TSM. The π-σ-π mixed-valence type molecules, descended from Aviram's original idea are analyzed. These molecules include 1, 6-heptadiene, 1, 8-nonadiene, and 1, 4-diallyl butane radical cations, which will be referred to as molecule 1, molecule 2, and molecule 3 in this paper, respectively (Figure 4). We optimized the geometry of these monocations using the DFT/B3LYP method. The dot-dot distance, l, in these molecules is between 0.5 to 0.8 Å. Bistability and electron localizability of these molecules have been studied in [3,46,49].

thumbnailFigure 4. Geometry of the molecules we used in our calculations. 1, 6-heptadiene, 1, 8-nonadiene, and 1, 4-diallyl butane radical monocations from left to right. In the first two molecules, coordinates of the three highlighted atoms have been set to xy plane, while the coordinates of the central carbon have been set to origin. For 1, 4-diallyl butane, the origin has been set in the middle of the bond between the two central carbon atoms.

Koopmans' theorem [52] has found extensive application in calculation of the ET matrix element, Hab for symmetric molecules. Under the two-state approximation, Hab is related to adiabatic energies of the ground and first excited state (E1 and E2) as [34]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M25">View MathML</a>


When no driver is present or the sum of the input drivers is zero, Po = 0; and from Equation 18 it is clear that Cos2η = 0, thus

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M26">View MathML</a>


According to the one electron Koopmans' theorem, the ionization potential of the highest occupied molecular orbital (HOMO) and HOMO-1 can be expressed in terms of the molecular orbital (MO) energies, i.e., [38-41]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M27">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M28">View MathML</a>


Inserting Equations 27 and 28 into Equation 26, the electronic coupling element is obtained in terms of the MO energies as [38-41]

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M29">View MathML</a>


The state-averaged CASSCF (SA/CASSCF) method [53] can be used to calculate the electronic coupling element of asymmetric molecules. One can obtain Hab by calculating the energies of the ground and first excited states and use them in Equation 26 within the two-state approximation. We have calculated the electronic coupling elements using both methods. The calculations for molecule 1 and 2 were based on SA/CASSCF(34), and the calculations for molecule 3 were based on SA/CASSCF(56). In 1, 4-diallyl butane cation, the allyl π-bonds are aromatic, and the active space is extended to five electrons in six orbitals. The calculated ET matrix elements have been compiled in Table 1. The Kink energies of the molecules have also been calculated using Equation 15. Table 1 includes all required parameters of the TSM. All calculations reported here were performed in Gaussian 09 [54], using 6-31 G(d) basis set. Various basis sets have been extensively tested to examine the basis set dependency of the results. Application of larger basis sets did not significantly influence the energy difference between the ground state and the first excited state. The results from the two methods are in good agreement. The ET matrix element, Hab decays exponentially with dot-dot distance, l[55]. The dot-dot distance for molecule 3 is less than that of molecule 2; however, the Hab has been remarkably decreased. This is due to the symmetry of this molecule and the aromatic bonds of the radical allyls. The geometrical parameter, l, and the type of head groups, play an important role in determining the ET matrix element. To obtain a more accurate electronic coupling element, the overlap integral, Sab, should be taken into account as described in the work of Farazdel et al [56]. Aviram [3] has obtained a negligible overlap integral for molecule 3 with dot-dot distance of 7 Å. In mQCA, the electron transfer drama should have a little effect on the geometric parameters [4]. Consequently, candidate molecules should possess fast electron transfer reactions, and the relaxation of nuclear degrees of freedom can be ignored. Table 1 also lists the changes in the head groups' π-bonds (Δζ) as a consequence of ET reactions. It is seen that ET reactions in molecule 3 should be faster compared to the other two molecules.

Table 1. Two-state model parameters for the used molecules

INV gates

The INV gate in two-dot mQCA is the nucleus of all other gates. Once its operation and switching properties are clearly understood, the properties of more intricate structures such as multidriver MinV gates can be derived from extrapolating the results obtained from the inverters, based on the additivity relation (Equation 20). The analysis of inverters can be extended to explain the behavior of more complex gates, which in turn form the building blocks for modules such as adders, multipliers, and processors. Table 2 compares the output polarizations of the INV gates, obtained from full quantum chemistry calculations and the TSM. The normalized dipole moments (Equation 2) of the monocations adjacent to fixed inputs (point charges as fixed drivers) have been calculated based on SA/CASSCF(3,4) for molecule 1 and 2, and SA/CASSCF(5,6) for molecule 3. The root mean square errors (RMSE) of the results obtained from the two methods have been calculated. The RMSE decreases with the increase of the μ parameter (or decrease of the ET matrix element, Hab), determining the degree of agreement between the results. The saturation polarization of the output is also dependant on the μ parameter. For INV gates, it is obtained by setting the sum of the input drivers to one in Equation 21 as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/274/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/274/mathml/M30">View MathML</a>


Table 2. INV gates

Equation 30 shows that the saturation polarization of the output increases with the increase of μ. This is also evident from the results in Table 2.

Two-driver devices

In the model MinV gate, the number of input drivers (m) should be odd. No logic operation is performed when m is even. However, neglecting the logic, the two-driver device is an appropriate small model system for studying the additivity relation, and how the accuracy of the TSM is influenced by the number of drivers. Here, the MinV gate is probed when only the two input drivers d1 and d3 are present (Figure 3a). We have calculated the normalized dipole moments of the gates' outputs based on the SA/CASSCF calculations. The output polarizations have also been calculated by the TSM. The results obtained from the two methods are compiled in Table 3, which are in good agreement. The conclusions from analysis of the INV gates can be extrapolated to two-driver devices as well. Compared to INV gates, the increase in the RMSE of the two-driver devices, composed of molecule 1 or 2, is mainly attributed to the asymmetric head groups. In other words, the effect of d1 on the head groups is different from that of d3, where Pd1= Pd3. In molecule 3, the allyl head groups are symmetric, and the TSM error mainly arises from the classical approximation of the intercellular interactions. It is important to note that the output polarization of the two-driver devices can be calculated by employing the additivity relation on the output polarizations of the INV gates. The additivity relation has been validated for the SA/CASSCF method as well. Through full quantum chemistry calculations, the output polarizations of the two-driver devices were obtained. We also used the results of the INV gates, P(Pd1+ Pd2) from Table 2, to examine the additivity relation for SA/CASSCF method. As expected, RMSE is highly dependent on the symmetry of the head groups. Unlike molecule 1 and 2, for the case of molecule 3, each driver has an exactly same effect on the head allyl groups, which leads to smaller RMSE. We also highlight that employing additivity relation decreases the computational cost of SA/CASSCF calculations. Table 2 and Table 3 also show how the accuracy of the TSM is affected by the number of input drivers. It is seen that RMSE is approximately doubled when the number of input drivers is scaled up by a factor of two.

Table 3. Two-driver devices

Three-input MinV gates

The multidriver MinV gate is a universal gate in the two-dot mQCA. The output polarizations of these gates with three fixed input drivers are shown in Table 4. This table also shows that the output polarizations of the MinV gates can be obtained from extrapolating the INV output polarizations using the additivity relation. For the model molecules, considering the spatial location of the d2, the effect of d2 on the head groups is different from the same effect of d1 and d3, while in the TSM, drivers with equal polarizations have same effects on the head groups and are treated the same. Quantum chemical calculations show that despite the equal sum of the input polarizations, the output polarizations are not equal particularly when the sum of the drivers is zero. Table 4 also displays that the SA/CASSCF method returns different output polarizations, while the sum of the input polarizations is zero. This is the main reason of the decrease in the accuracy of the results obtained from the TSM for MinV gates. Ignoring these points by avoiding the null state logic occurrence, the two-state approximation results are fairly in good agreement with the quantum chemical calculations. Table 2 and Table 4 show how the accuracy of the two-state model is decreased with the number of drivers. It is seen that RMSE is tripled when the number of input drivers is scaled up by a factor of three.

Table 4. Three-driverMinVgates


Molecular QCA gates are the building blocks of more complex modules. Probing molecular devices requires quantum chemical calculations, which are challenging as the molecular system grows in size. A semi-classical model was derived directly from the two-state approximation in the ET theory, serving as a device for studying mQCA gates. This model is very similar to the two-state model which is currently the core of the QCADesigner simulation engine for solving circuits based on metallic QCA. The range of applications and limitations of this model for mQCA gates was investigated carefully. The parametric TSM can be used to study more complex mQCA gates composed of practical candidate mixed-valence molecules, where exploiting the SA/CASSCF method is of high computational cost. A complete set of logic gates were implemented within the two-dot mQCA framework. These gates include INV and MinV gates, which provide a small molecular test bench, making further analysis by quantum chemistry methods, particularly SA/CASSCF, practical. The INV gate was studied as a nucleus of all other gates. It was also presented that output polarizations of all other gates can be derived from extrapolating the results obtained from inverters based on the additivity relation. We compared the results obtained from the TSM to those obtained from SA/CASSCF calculations for INV and MinV gates. The degree of agreement between the TSM and quantum chemical calculations is highly dependent on the μ parameter and the symmetry of the head groups. Additionally, application of the additivity relation for CASSCF method can in turn reduce the computational cost. It is important to note that we did not address questions of surface attachment, input/output, clocked control, layout, and patterning, which are the requirements of a practical QCA system. Moreover, we did not consider the relaxation of nuclear degrees of freedom associated with electron transfer. It is presented that for mQCA, the electron localization and Coulombic interactions play the key roles, and nuclear positions can be considered frozen (nuclear relaxation even assists charge localization) [4]. Although we limited our focus on the two-dot mQCA, it merits highlighting that the model can also be used for four-dot cells, since they can be considered as double two-dot cells. Our focus was on the mQCA gates as building blocks of circuits. The two-state model may be applied to simulate mQCA circuits as well, as it is currently used iteratively for simulation of metallic QCA circuits in the QCADesigner. However, to determine the additive error resulting from exploiting the two-state model for solving mQCA circuits, further quantum chemical calculations on the mQCA clocked circuits composed of several molecules are required, which are extremely challenging at the time, and have not been addressed in this paper.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ER developed the theories and carried out the quantum chemical calculations. SMN supervised the project. All authors read and approved the final manuscript.


ER was affiliated with Norwegian University of Science and Technology. Calculations presented in this work have been carried out on Stallo. ER also thanks Professor Sven Larsson for many enlightening discussions on electron transfer theory at Chalmers University of Technology.


  1. Lent CS: Bypassing the transistor paradigm.

    Science 2000, 288:1597-1599. Publisher Full Text OpenURL

  2. Niemier M, Kogge P, Murphy R, Rodrigues A, Dysart T, Frost S: Data flow in molecular QCA: logic can “sprint”, but the memory wall can still be a "hurdle". University of Notre Dame: Technical report; 2005. OpenURL

  3. Aviram A: Molecules for memory, logic, and amplification.

    J Am Chem Soc 1988, 110:5687-5692. Publisher Full Text OpenURL

  4. Lent CS, Isaksen B, Lieberman M: Molecular quantum-dot cellular automata.

    J Am Chem Soc 2003, 125:1056-1063. PubMed Abstract | Publisher Full Text OpenURL

  5. Orlov AO, Amlani I, Kummamuru RK, Ramasubramaniam R, Toth G, Lent CS, Bernstein GH, Snider GL: Experimental demonstration of clocked single-electron switching in quantum-dot cellular automata.

    Appl Phys Lett 2000, 77:295-297. Publisher Full Text OpenURL

  6. Smith C: Realization of quantum-dot cellular automata using semiconductor quantum dots.

    Superlatt Microst 2003, 34:195-203. Publisher Full Text OpenURL

  7. Gardelis S, Smith C, Cooper J, Ritchie D, Linfield E, Jin Y: Evidence for transfer of polarization in a quantum-dot cellular automata cell consisting of semiconductor quantum dots.

    Phys Rev B 2003, 67:033302. OpenURL

  8. Lent CS, Isaksen B: Clocked molecular quantum-dot cellular automata.

    IEEE T Electron Dev 2003, 50:1890-1896. Publisher Full Text OpenURL

  9. Qi H, Sharma S, Li Z, Snider GL, Orlov AO, Lent CS, Fehlner TP: Molecular quantum cellular automata cells. electric field driven switching of a silicon surface bound array of vertically oriented two-dot molecular quantum cellular automata.

    J Am Chem Soc 2003, 125:15250-15259. PubMed Abstract | Publisher Full Text OpenURL

  10. Jiao J, Long GJ, Grandjean F, Beatty AM, Fehlner TP: Building blocks for the molecular expression of quantum cellular automata. Isolation and characterization of a covalently bonded square array of two ferrocenium and two ferrocene complexes.

    J Am Chem Soc 2003, 125:7522-7523. PubMed Abstract | Publisher Full Text OpenURL

  11. Jiao J, Long GJ, Rebbouh L, Grandjean F, Beatty AM, Fehlner TP: Properties of a mixed-valence (feII)2(feIII)2 square cell for utilization in the quantum cellular automata paradigm for molecular electronics.

    J Am Chem Soc 2005, 127:17819-17831. PubMed Abstract | Publisher Full Text OpenURL

  12. Qi H, Gupta A, Noll BC, Snider GL, Lu Y, Lent C, Fehlner TP: Dependence of field switched ordered arrays of dinuclear mixed-valence complexes on the distance between the redox centers and the size of the counterions.

    J Am Chem Soc 2005, 127:15218-15227. PubMed Abstract | Publisher Full Text OpenURL

  13. Li Z, Fehlner TP: Molecular QCA cells. 2. Characterization of an unsymmetrical dinuclear mixed-valence complex bound to a Au surface by an organic linker.

    Inorg Chem 2003, 42:5715-5721. PubMed Abstract | Publisher Full Text OpenURL

  14. Lu Y, Lent CS: Theoretical study of molecular quantum-dot cellular automata.

    J Comput Elec 2005, 4:115-118. Publisher Full Text OpenURL

  15. Braun-Sand SB, Wiest O: Theoretical studies of mixed-valence transition metal complexes for molecular computing.

    J Phys Chem A 2002, 107:285-291. OpenURL

  16. Walus K, Dysart TJ, Jullien GA, Budiman RA: QCADesigner: a rapid design and simulation tool for quantum-dot cellular automata.

    IEEE T Nanotechnol 2004, 3:26-31. Publisher Full Text OpenURL

  17. Vetteth A, Walus K, Dimitrov SV, Jullien GA: Quantum-dot cellular automata carry-look-ahead adder and barrel shifter. Richardson, TX: IEEE Conf on Emerging Telecom Tech; 2002. OpenURL

  18. Frost SE, Rodrigues AF, Janiszewski AW, Rausch RT, Kogge PM: Memory in Motion: a Study of Storage Structures in QCA.

    1st Workshop on Non-Silicon Computation 2002. OpenURL

  19. Niemier MT, Kogge PM: Logic in wire: using quantum dots to implement a microprocessor.

    ICECS '99 6th IEEE Int Conf Circuits and Systems 1999, 3:1211-1215. PubMed Abstract | Publisher Full Text OpenURL

  20. Walus K, Vetteth A, Jullien GA, Dimitrov VS: RAM design using quantum-dot cellular automata.

    Nanotechnology Conf and Trade Show 2003, 2:160-163. OpenURL

  21. Lent CS, Tougaw PD: A device architecture for computing with quantum dots.

    Proc of the IEEE 1997, 85:541-557. Publisher Full Text OpenURL

  22. Rahimi E, Nejad SM: Secure clocked QCA logic for implementation of quantum cryptographic processors.

    Int Conf on Applied Electronics 2009, 217-220. OpenURL

  23. Rahimi E, Nejad SM: Quantum-Dot Cellular ROM: A nano-scale level approach to digital data storage.

    6th IEEE Int Conf on Communication Systems, Networks and Digital Signal Processing 2008, 618-621. OpenURL

  24. Rahimi E, Nejad SM: A novel architecture for quantum-dot cellular ROM.

    7th IEEE Int Conf on Communication Systems, Networks and Digital Signal Processing 2010, 347-350. OpenURL

  25. Tougaw PD, Lent CS: Logical devices implemented using quantum cellular automata.

    J Appl Phys 1994, 75:1818-1825. Publisher Full Text OpenURL

  26. Mohammad Nejad S, Rahimi E: QCA: The prospective technology for digital telecommunication systems. In Nanotechnology for Telecommunications. CRC Press; 2010:275-307. OpenURL

  27. Hennessy K, Lent CS: Clocking of molecular quantum-dot cellular automata.

    J Vac Sci Technol B 2001, 19:1752-1755. OpenURL

  28. Orlov AO, Kummamuru R, Ramasubramaniam R, Lent CS, Bernstein GH, Snider GL: Clocked quantum-dot cellular automata shift register.

    Surf Sci 2003, 532–535:1193-1198. OpenURL

  29. Karim F, Walus K, Ivanov A: Analysis of field-driven clocking for molecular quantum-dot cellular automata based circuits.

    J Comput Elec 2010, 9:16-30. Publisher Full Text OpenURL

  30. Lent CS, Tougaw PD, Porod W, Bernstein GH: Quantum cellular automata.

    Nanotechnology 1993, 4:49. Publisher Full Text OpenURL

  31. Timler J, Lent CS: Power gain and dissipation in quantum-dot cellular automata.

    J Appl Phys 2002, 91:823-831. Publisher Full Text OpenURL

  32. Tougaw PD, Lent CS: Dynamic behavior of quantum cellular automata.

    J Appl Phys 1996, 80:4722. Publisher Full Text OpenURL

  33. Géza T: Correlation and coherence in quantum-dot cellular automata. PhD Thesis. University of Notre Dame; 2000. OpenURL

  34. Newton MD: Quantum chemical probes of electron-transfer kinetics: the nature of donor-acceptor interactions.

    Chem Rev 1991, 91:767-792. Publisher Full Text OpenURL

  35. Newton MD: Control of electron transfer kinetics: models for medium reorganization and donor–acceptor coupling.

    Adv Chem Phys 2007, 106:303-375. OpenURL

  36. Zener C: Non-adiabatic crossing of energy levels.

    Proc R Soc A 1932, 137:696-702. Publisher Full Text OpenURL

  37. Zener C: Discussion of excited diatomic molecules by external perturbations.

    Proc R Soc A 1933, 140:660-668. Publisher Full Text OpenURL

  38. Larsson S, Braga M: Transfer of mobile electrons in organic molecules.

    Chem Phys 1993, 176:367-375. Publisher Full Text OpenURL

  39. Li X-Y, Tang X-S, Xiao S-Q, He F-C: Application of Koopmans’ theorem in evaluating electron transfer matrix element of long-range electron transfer.

    J Mol Struct (THEOCHEM) 2000, 530:49-58. Publisher Full Text OpenURL

  40. Li X-Y, Tang X-S, He F-C: Electron transfer in poly (p-phenylene) oligomers: effect of external electric field and application of Koopmans theorem.

    Chem Phys 1999, 248:137-146. Publisher Full Text OpenURL

  41. Rodriguez-Monge L, Larsson S: Conductivity in polyacetylene.

    IV. Ab initio calculations for a two-site model for electron transfer between allyl anion and allyl. Int J Quantum Chem 1997, 61:847-857. OpenURL

  42. Hush NS: Intervalence-transfer absorption. Part 2. Theoretical considerations and spectroscopic data.

    Prog Inorg Chem Vol 2007, 8:391-444. OpenURL

  43. Hush NS: Homogeneous and heterogeneous optical and thermal electron transfer.

    Electrochim Acta 1968, 13:1005-1023. Publisher Full Text OpenURL

  44. Bailey SE, Zink JI, Nelsen SF: Contributions of symmetric and asymmetric normal coordinates to the intervalence electronic absorption and resonance raman spectra of a strongly coupled p-phenylenediamine radical cation.

    J Am Chem Soc 2003, 125:5939-5947. PubMed Abstract | Publisher Full Text OpenURL

  45. Coropceanu V, Gruhn NE, Barlow S, Lambert C, Durivage JC, Bill TG, Nöll G, Marder SR, Brédas J-L: Electronic couplings in organic mixed-valence compounds: the contribution of photoelectron spectroscopy.

    J Am Chem Soc 2004, 126:2727-2731. PubMed Abstract | Publisher Full Text OpenURL

  46. Hush NS, Wong AT, Bacskay GB, Reimers JR: Electron and energy transfer through bridged systems. 6. Molecular switches: the critical field in electric field activated bistable molecules.

    J Am Chem Soc 1990, 112:4192-4197. Publisher Full Text OpenURL

  47. Cabrero J, Calzado CJ, Maynau D, Caballol R, Malrieu JP: Metal−ligand delocalization in magnetic orbitals of binuclear complexes.

    J Phys Chem A 2002, 106:8146-8155. Publisher Full Text OpenURL

  48. Lu Y, Quardokus R, Lent CS, Justaud F, Lapinte C, Kandel SA: Charge localization in isolated mixed-valence complexes: an STM and theoretical study.

    J Am Chem Soc 2010, 132:13519-13524. PubMed Abstract | Publisher Full Text OpenURL

  49. Lu Y, Lent CS: A metric for characterizing the bistability of molecular quantum-dot cellular automata.

    Nanotechnology 2008, 19:155703. PubMed Abstract | Publisher Full Text OpenURL

  50. Malmqvist P-Å, Roos BO: The CASSCF state interaction method.

    Chem Phys Lett 1989, 155:189-194. Publisher Full Text OpenURL

  51. Yamamoto N, Vreven T, Robb MA, Frisch MJ, Bernhard Schlegel H: A direct derivative MC-SCF procedure.

    Chem Phys Lett 1996, 250:373-378. Publisher Full Text OpenURL

  52. Koopmans T: Über die zuordnung von wellenfunktionen und eigenwerten zu den einzelnen elektronen eines atoms.

    Physica 1934, 1:104-113. Publisher Full Text OpenURL

  53. Ruedenberg K, Cheung LM, Elbert ST: MCSCF optimization through combined use of natural orbitals and the brillouin–levy–berthier theorem.

    Int J Quantum Chem 1979, 16:1069-1101. Publisher Full Text OpenURL

  54. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, et al.:

    Gaussian 09. Revision B.01. 2009. OpenURL

  55. Newton MD, Cave RJ: Molecular Control of Electron and Hole Transfer Processes: Theory and Applications.

    Mol Electron 1997., 73 OpenURL

  56. Farazdel A, Dupuis M, Clementi E, Aviram A: Electric-field induced intramolecular electron transfer in spiro.pi.-electron systems and their suitability as molecular electronic devices. A theoretical study.

    J Am Chem Soc 1990, 112:4206-4214. Publisher Full Text OpenURL