Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

Molecular orbital theory calculations of the H2O-carbon reaction

Carbon gasification with steam to produce H-2 and CO is an important reaction widely used in industry for hydrogen generation. Although the literature is vast, the. mechanism for the formation of H-2 is still unclear. In particular, little has, been done to investigate the potential of molecular orbital theory to distinguish different mechanism possibilities. In this work, we used molecular orbital theory to demonstrate a favorable energetic pathway where H2O is first physically adsorbed on the virgin graphite surface with negligible change in molecular structure. Chemisorption occurs via O approaching the carbon edge site with one H atom stretching away from the O in the transition state. This is followed by a local minimum. state in which the stretching H is further disconnected from the O atoms and the remaining OH group is still on the carbon edge site. The disconnected H then pivot around the OH group to bond with the H of the OH group and forms H-2. The O atom remaining on the carbon edge site is subsequently desorbed as CO. The reverse occurs when H-2 reacts with the surface oxygen to produce H2O.

A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Ab initio modelling of basal plane oxidation of graphenes and implications for modelling char combustion

Ab initio calculations have been performed to determine the energetics of oxygen atoms adsorbed onto graphene planes and the possible reaction path extracting carbon atorns in the form of carbon monoxide. Front the energetics it is confirmed that this reaction path will not significantly contribute to the gasification of well ordered carbonaceous chars. Modelling results which explore this limit Lire presented. (C) 2002 Elsevier Science Ltd, All rights reserved.

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. All rights reserved.

New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Iminopropadienones R-N=C=C=C=O and carbon suboxide, C3O2. Theoretical and experimental 13C NMR spectra

The C-13 NMR data of five iminopropadienones R-N=C=C=C=O as well as carbon suboxide, C3O2, have been examined theoretically and experimentally. The best theoretical results were obtained using the GIAO/B3LYP/6-31 +G**//MP2/6-31G* level of theory, which reproduces the chemical shifts of the iminopropadienone substituents extremely well while underestimating those of the cumulenic carbons by 5-10 ppm. The computationally faster GIAO/HF/6-31 + G**//B3LYP/6-31 G* level is also adequate. (C) 2004 Elsevier B.V. All rights reserved.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Comparative study of Li, Na, and K adsorptions on graphite by using ab initio method

A comprehensive study has been conducted to compare the adsorptions of alkali metals (including Li, Na, and K) on the basal plane of graphite by using molecular orbital theory calculations. All three metal atoms prefer to be adsorbed on the middle hollow site above a hexagonal aromatic ring. A novel phenomenon was observed, that is, Na, instead of Li or K, is the weakest among the three types of metal atoms in adsorption. The reason is that the SOMO (single occupied molecular orbital) of the Na atom is exactly at the middle point between the HOMO and the LUMO of the graphite layer in energy level. As a result, the SOMO of Na cannot form a stable interaction with either the HOMO or the LUMO of the graphite. On the other hand, the SOMO of Li and K can form a relatively stable interaction with either the HOMO or the LUMO of graphite. Why Li has a relatively stronger adsorption than K on graphite has also been interpreted on the basis of their molecular-orbital energy levels.

The role of Ti as a catalyst for the dissociation of hydrogen on a Mg(0001) surface

In this paper, the dissociative chemisorption of hydrogen on both pure and Ti-incorporated Mg(0001) surfaces are studied by ab initio density functional theory (DFT) calculations. The calculated dissociation barrier of hydrogen molecule on a pure Mg(0001) surface (1.05 eV) is in good agreement with comparable theoretical studies. For the Ti-incorporated Mg(0001) surface, the activated barrier decreases to 0.103 eV due to the strong interaction between the molecular orbital of hydrogen and the d metal state of Ti. This could explain the experimentally observed improvement in absorption kinetics of hydrogen when transition metals have been introduced into the magnesium materials.

The 2.1 angstrom crystal structure of the far-red fluorescent protein HcRed: Inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1 resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu64-Tyr65-Gly66, was observed to adopt both a cis coplanar and a tran. non-coplanar conformation. As a result of these two con formations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the b-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible. 2005 Elsevier Ltd. All rights reserved.

5,6-Dihydroxyindole-2-carboxylic acid: a first principles density functional study

We report first principles density functional calculations for 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and several oxidised forms. DHICA and 5,6-dihydroxyindole (DHI) are believed to be the basic building blocks of the eumelanins. Our results show that carboxylation has a significant effect on the physical properties of the molecules. In particular, the relative stabilities and the highest occupied molecular orbital-lowest unoccupied molecular orbital gaps (calculated with the DeltaSCF method) of the various redox forms are strongly affected. We predict that, in contrast to DHI, the density of unpaired electrons, and hence the ESR signal, in DHICA is negligibly small. (C) 2004 Elsevier B.V. All rights reserved.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.

Chemical and structural disorder in eumelanins: A possible explanation for broadband absorbance

We report the results of an experimental and theoretical study of the electronic and structural properties of a key eumelanin precursor-5,6,-dihydroxyindole-2-carboxylic acid ( DHICA) - and its dimeric forms. We have used optical spectroscopy to follow the oxidative polymerization of DHICA to eumelanin and observe red shifting and broadening of the absorption spectrum as the reaction proceeds. First principles density functional theory calculations indicate that DHICA oligomers ( possible reaction products of oxidative polymerization) have the gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital red-shifted gaps with respect to the monomer. Furthermore, different bonding configurations ( leading to oligomers with different structures) produce a range of gaps. These experimental and theoretical results lend support to the chemical disorder model where the broadband monotonic absorption characteristic of all melanins is a consequence of the superposition of a large number of nonhomogeneously broadened Gaussian transitions associated with each of the components of a melanin ensemble. These results suggest that the traditional model of eumelanin as an amorphous organic semiconductor is not required to explain its optical properties and should be thoroughly reexamined. These results have significant implications for our understanding of the physics, chemistry, and biological function of these important biological macromolecules. Indeed, one may speculate that the robust functionality of melanins in vitro is a direct consequence of its heterogeneity, i.e., chemical disorder is a "low cost" natural resource in these systems