A tight binding model for water

We demonstrate for the first time a tight binding model for water incorporating polarizable oxygen atoms. A novel aspect is that we adopt a ``ground up'' approach in that properties of the monomer and dimer only are fitted. Subsequently we make predictions of the structure and properties of hexamer clusters, ice-XI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost two-fold increase in molecular dipole moment as clusters are built up towards the limit of bulk liquid. We concentrate on properties of liquid water, particularly dielectric constant and self diffusion coefficient, which are very well rendered in comparison with experiment. Finally we comment on the question of the contrasting densities of water and ice which is central to an understanding of the subtleties of the hydrogen bond.

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. (C) 2014 AIP Publishing LLC.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Ab initio and tight-binding calculations of noncollinear magnetism in manganese clusters

Combining ab initio and tight-binding calculations, we have studied the noncollinear magnetism in manganese clusters. The oscillations in the per-atom moments observed experimentally are reproduced theoretically. The tendency of antiferromagnetic coupling between near neighbors leads to noncollinear coupling between atoms within the clusters. For clusters containing 12, 13, 15, 19, and 23 atoms, the geometrical structures were optimized from ab initio calculations with collinear coupled spin moments among different atomic sites. For larger clusters such as Mn-36 and Mn-55, the geometries are taken as portions of an fcc structure. Although the local atomic moments have high values close to 4 mu(B), the net moments lie in the range of 0.4-1.2 mu(B)/atom. Taking the noncollinear coupling into account brings the calculated magnetic moments much closer to the experimental results.

Tight-binding model for carbon nanotubes from ab initio calculations

Here we present a parametrized tight-binding (TB) model to calculate the band structure of single-wall carbon nanotubes (SWNTs). On the basis of ab initio calculations we fit the band structure of nanotubes of different radii with results obtained with an orthogonal TB model to third neighbors, which includes the effects of orbital hybridization by means of a reduced set of parameters. The functional form for the dependence of these parameters on the radius of the tubes can be used to interpolate appropriate TB parameters for different SWNTs and to study the effects of curvature on their electronic properties. Additionally, we have shown that the model gives an appropriate description of the optical spectra of SWNTs, which can be useful for a proper assignation of SWNTs` specific chirality from optical absorption experiments.

Análise da potencialidade a condução do BDT através da densidade eletrônica de estados via tight binding e desenvolvimento de software (B3J)

Neste trabalho nos propomos a fazer um estudo acerca da potencialidade de condução eletrônica no polímero BDT (1,3-benzoditiol 4H-ciclopenta[2,1-b:3,4b’]). O estudo usual de polímeros conjugados é feito de modo a obter sua densidade de estados com diversos tipos e níveis de dopagem. O método de Huckel é o mais utilizado e se baseia na separabilidade das ligações sigma e pi que é possível quando a molécula estudada é plana. Os polímeros conjugados são em sua maioria planos e estão inseridos nesta aproximação. O monômero do BDT apresenta sua geometria fora do plano por apresentar ligações com orbitais sp3. Para contornar esse problema foi desenvolvido o programa B3J, que considera todos os orbitais de valencia (s, px, py e pz). O programa B3J calcula a densidade de estados de sistemas poliméricos. O estudo das bandas do BDT foi feito com este software. Calculamos a densidade de estados do sistema neutro e com diversos níveis de dopagem, com distribuição aleatória e ordenada dos defeitos, dopagem do tipo n e do tipo p. O comportamento do quadrado do coeficiente da expansão da função de onda foi obtido para polímeros de até 20 monômeros. Estes cálculos foram obtidos com geometrias dos métodos AM1 e PM3. Obtivemos os espectros de absorção de oligômeros a fim de inferir seu comportamento para um polímero. Foram utilizados cálculos de otimização de geometria através dos métodos semi-empíricos AM1 e PM3 e ZINDO/S e o método DFT. Em outro objetivo desta monografia há o estudo do aproveitamento de tetrâmeros de BDT como dispositivos eletrônicos. Tais oligômeros foram otimizados em diversos valores de potencial elétrico, com a inserção em suas cadeias de moléculas doadoras e aceitadoras para induzir um aumento no momento de dipolo da mesma.

Peptide nucleic acids: Mechanism of tight binding and application in artificial nuclease

Peptide nucleic acids (PNA) are mimics of nucleic acids with a peptidic backbone. Duplexes and triplexes formed between PNA and DNA or RNA possess remarkable thermal stability, they are resistant to nuclease cleavage and can better discriminate mismatches. Understanding the mechanism for the tight binding between PNA and oligonucleotides is important for the design and development of better PNA-based drugs.^ We have performed molecular dynamics (MD) simulations of 8-mer PNA/DNA duplex and two analogous duplexes with chiral modification of PNA strand (D- or L-Alanine modification). MD simulations were performed with explicit water and Na$\sp{+}$ counter ions. The 1.5-ns simulations were carried out with AMBER using periodic boundary and particle mesh Ewald summation. The point charges for PNA monomers were derived from fitting electrostatic potentials, obtained from ab initio calculation, to atomic centers using RESP. Derived charges reveal significantly altered charge distribution on the PNA bases and predict the Watson-Crick H-bonds involving PNA to be stronger. Results from NMR studies investigating H-bond interactions between DNA-DNA and DNA-PNA base pairs in non-polar environment are consistent with this prediction. MD simulations demonstrated that the PNA strand is more flexible than the DNA strand in the same duplex. That this flexibility might be important for the duplex stability is tested by introducing modification into the PNA backbones. Results from MD simulation revealed dramatically altered structures for the modified PNA-DNA duplexes. Consistent with previous NMR results, we also found no intrachain hydrogen bonds between O7$\sp\prime$ and N1$\sp\prime$ of the neighboring residues in our MD study. Our study reveals that in addition to the lack of charge repulsion, stronger Watson-Crick hydrogen bonds together with flexible backbone are important factors for the enhanced stability of the PNA-DNA duplex.^ In a related study, we have developed an application of Gly-Gly-His-(Gly)$\sb3$-PNA conjugate as an artificial nuclease. We were able to demonstrate cleavage of single stranded DNA at a single site upon Ni(II) binding to Gly-Gly-His tripeptide and activation of nuclease with monoperoxyphthalic acid. ^

Cyclopropane-1,1-dicarboxylate is a slow-, tight-binding inhibitor of rice ketol-acid reductoisomerase

Ketol-acid reductoisomerase (EC 1.1.1.86) catalyses the second reaction in the biosynthesis of the branched-chain amino acids. The reaction catalyzed consists of two stages, the first of which is an alkyl migration from one carbon atom to its neighbour. The likely transition state is therefore a cyclopropane derivative, and cyclopropane-1,1-dicarboxylate(CPD) has been reported to inhibit the Escherichia coli enzyme. In addition, this compound causes the accumulation of the substrate of ketol-acid reductoisomerase in plants. Here, we investigate the inhibition of the purified rice enzyme. The cDNA was cloned, and the recombinant protein was expressed in E. coli, purified and characterized kinetically. The purified enzyme is strongly inhibited by cyclopropane-1,1-dicarboxylate, with an inhibition constant of 90 nM. The inhibition is time-dependent and this is due to the low rate constants for formation (2.63 X 10(5) M-1 min(-1)) and dissociation (2.37 x 10(-2) min(-1)) of the enzyme-inhibitor complex. Other cyclopropane derivatives are much weaker inhibitors while dimethylmalonate is moderately effective. (c) 2004 Elsevier Ireland Ltd. All rights reserved.

A tight-binding analysis of the electronic properties of III-nitride semiconductors

This thesis divides into two distinct parts, both of which are underpinned by the tight-binding model. The first part covers our implementation of the tight-binding model in conjunction with the Berry phase theory of electronic polarisation to probe the atomistic origins of spontaneous polarisation and piezoelectricity as well as attempting to accurately calculate the values and coefficients associated with these phenomena. We first develop an analytic model for the polarisation of a one-dimensional linear chain of atoms. We compare the zincblende and ideal wurtzite structures in terms of effective charges, spontaneous polarisation and piezoelectric coefficients, within a first nearest neighbour tight-binding model. We further compare these to real wurtzite structures and conclude that accurate quantitative results are beyond the scope of this model but qualitative trends can still be described. The second part of this thesis deals with implementing the tight-binding model to investigate the effect of local alloy fluctuations in bulk AlGaN alloys and InGaN quantum wells. We calculate the band gap evolution of Al1_xGaxN across the full composition range and compare it to experiment as well as fitting bowing parameters to the band gap as well as to the conduction band and valence band edges. We also investigate the wavefunction character of the valence band edge to determine the composition at which the optical polarisation switches in Al1_xGaxN alloys. Finally, we examine electron and hole localisation in InGaN quantum wells. We show how the built-in field localises the carriers along the c-axis and how local alloy fluctuations strongly localise the highest hole states in the c-plane, while the electrons remain delocalised in the c-plane. We show how this localisation affects the charge density overlap and also investigate the effect of well width fluctuations on the localisation of the electrons.

Benchmarking Density Functional Based Tight-Binding for Silver and Gold Materials: From Small Clusters to Bulk

International audience

Aharonov-Bohm oscillation and chirality effect in optical activity of single-wall carbon nanotubes

We study the Aharonov-Bohm effect in the optical phenomena of single-wall carbon nanotubes (SWCN) and also their chirality dependence. Especially, we consider the natural optical activity as a proper observable and derive its general expression based on a comprehensive symmetry analysis, which reveals the interplay between the enclosed magnetic flux and the tubule chirality for arbitrary chiral SWCN. A quantitative result for this optical property is given by a gauge invariant tight-binding approximation calculation to stimulate experimental measurements.

Plasmons in vertically coupled InAs/GaAs quantum dots

We investigate plasmon excitations in a quantum wire that consists of an infinite one-dimensional array of vertically coupled InAs/GaAs strained quantum dots (QDs). The research is carried out in the framework of random-phase approximation using effective-mass theory. Our formalism is capable of studying plasmons with strong tunneling among QDs, which frustrate the conventionally adopted tight-binding approximation. Based on this formalism, a systematic study on the intraminiband or intrasubband plasmon in vertically coupled InAs/GaAs strained QDs is presented. It is found that an increase of the dot spacing will inevitably reduce the plasmon energy. In contrast, the role of dot height is relatively complex and depends on the dot spacing. The results demonstrate the possibility to engineer collective excitations in low dimensional systems by simply changing their geometric configuration.

THE CHARGE-DISTRIBUTION AND THE ELECTRIC-FIELD GRADIENT IN THE HIGH-TEMPERATURE SUPERCONDUCTOR TLBA2CAN-1CUNO2N+3 (N = 1 TO 3)

The net charges at atoms in the high-temperature superconductor TlBa2Can-1CunO2n+3 (n = 1 to 3) are calculated by means of the tight-binding approximation based on the EHMO method. The results indicate that the charge distribution in this kind of compounds possesses a specially layered arrangement. An insulating Ba-Ba layer is inserted between the Cu-O layer and the Tl-O layer. There may exist a weak coupling between the Cu-O layer and the Tl-O layer through the interaction of the same O(2) atom with both the Cu atom and the Tl atom. The existence of the Ca in the compounds can cause the valence fluctuation at the Cu atom. The calculated electric field gradients at atoms implies that the conducting electron or hole may move in the Cu-O layer, which is closest to the Tl-O layer, along the a-b plane.