Analytical potential energy function and spectroscopic parameters for the ground and excited states of NaH

set of energies at different internuclear distances for the ground electronic state and two excited electronic states of NaH molecule have been calculated using valence internally contracted multireference configuration interaction(MRCI) including Davidson correction and three basis sets. Then, a potential energy curve (PEC) for each state was determined by extrapolating MRCI energies to the complete basis sets limit. Based on the PECs, accurate vibrational energy levels and rotational constants were determined. The computational PECs are were fitted to analytical potential energy functions using the Murrell-Sorbie potential function. Then, accurate spectroscopic parameters were calculated. Compared with experimental results, values obtained with the basis set extrapolation yield a potential energy curve that gives accurate vibrational energy levels, rotational constants and spectroscopic parameters for the NaH molecule. (C) 2004 Elsevier B.V. All rights reserved.

Golden-rule treatment on the ClO/ClO+ electron-transfer system

The structures, properties and electron transfer reactivity of the ClO/ClO+ coupling system are studied in this paper at ab initio (HF and MP2) levels and the density functional theory (DFT: B3LYP, B3P86, B3PW91) levels employing 6311 + G(3df) basis set and on the basis of the golden-rule of the time-dependent perturbation theory. Investigations indicate that the results got from the B3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated by using the B3LYP/6-311 + G(3df) method, and then the electron transfer rates are determined at this level. The electronic coupling matrix element of EC.6 is very small, only 0.03 kcal/mol, while that of EC.7 is the biggest, being 12.41 kcal/mol, the corresponding electron transfer rate is also the fastest among these seven encounter complexes. The averaged electron transfer rate is about 1.672 X 10(11) M-1 s(-1). It is indicated that the structures optimized by B3LYP method are more reliable than the results got from the other four methods. It also testified that the electronic coupling matrix element is the vital factor that significantly affects the electron transfer rate. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical study on the ClO/ClO- system electron-transfer reactivity by the golden-rule

The structures, properties and electron transfer reactivity of the ClO/ClO- coupling system are studied in this paper at ab initio (UHF and UMP2) levels and the Density Functional Theory (DFT: UB3LYP, UB3P86, UB3PW91) levels employing 6311 + G(3df) basis set and on the basis of the Golden-rule of the time-dependent perturbation theory. Investigations indicate that the results obtained using the UB3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. For this coupling system, six stable coupling modes have been found which correspond to six different encounter complexes and denote six different electron transfer mechanism: four O-O directly linked structures (one collinear: D-h, one anti-parallel: C-s, two twist: C-2) and two Cl-O linked structures (cis- and anti- C-s structures). The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated for the electron transfer reactions via these six different mechanism at the UB3LYP/6-311 + G(3df) level, and then the electron transfer rates are determined at the same level. The most favorable coupling mode to the electron transfer is the anti-parallel mechanism. The averaged electron transfer rate is about 5.58 X 10(11) M-1 s(-1). It is also implied that the B3LYP method can give more reasonable results for the electron transfer reactivity of this system. (C) 2003 Elsevier B.V. All rights reserved.

Bonding and correlation analysis of various SiCO isomers

The structural properties for various SiCO isomers in the singlet and triplet states have been investigated using CASSCF methods with a 6-311 +G* basis set and also using three DFT and MP2 with same basis set for those systems except for the linear singlet state. The detailed bonding character is discussed, and the state-state correlations and the isomerization mechanism are also determined. Results indicate that there are four different isomers for each spin state, and for all isomers, the triplet state is more stable than the corresponding singlet state. The most stable is the linear SiCO ((3)Sigma(-)) species and may be refer-red to the ground state. At the CASSCF-MP2(full)/6-311+G* level, the state-state energy separations of the other triplet states relative to the ground state are 43.2 (cyclic), 45.2 (linear SiOC), and 75.6 kcal/mol (linear CSiO), respectively, whereas the triplet-singlet state excitation energies for each configuration are 17.3 (linear SiCO), 2.2 (cyclic SiCO), 10.2 (linear SiOC), and 18.5 kcal/mol (linear CSiO), respectively. SiCo ((3)Sigma(-)) may be classified as silene (carbonylsilene), and its COdelta- moiety possesses CO- property. The dissociation energy of the ground state is 42.5 kcal/mol at the CASSCF-MP2(full)/6-311+G* level and falls within a range of 36.5-41.5 kcal/mol at DFT level, and of 23.7-28.9 kcal/mol at the wave function-correlated level, whereas the vertical IP is 188.8 kcal/mol at the CASSCF-MP2(full)/6-311+G* level and is very close to the first IP of Si atom. Three linear isomers (SiCO, SiOC, and CSiO) have similar structural bonding character. SiOC may be referred to the iso-carbonyl Si instead of the aether compound, whereas the CSiO isomer may be considered as the combination of C (the analogue of Si) with SiO (the analogue of CO). The bonding is weak for all linear species, and the corresponding potential energy surfaces are flat, and thus these linear molecules are facile. Another important isomer is of cyclic structure, it may be considered as the combination of CO with Si by the side pi bond. This structure has the smallest triplet state-singlet state excitation energy (similar to2.2 kcal/mol); the C-O bonds are longer, and the corresponding vibrational frequencies are significantly smaller than those of the other linear species. This cyclic species is not classified as an epoxy compound. State-state correlation analysis and the isomerization pathway searches have indicated that there are no direct correlations among three linear structures for each spin state, but they may interchange by experiencing two transition states and one cyclic intermediate. The easiest pathway is to break the Si-O bond to go to the linear SiCO, but its inverse process is very difficult. The most difficult process is to break the C-O bond and to go to the linear CSiO.

Synthesis, photophysical properties, and theoretical studies on pyrrole-containing bromo Re(I) complex

Two bromo rhenium(I) carbonyl complexes with the formula of [Re(CO)(3)(L)Br], where L = 1,10-phenanthroline (Phen-Re) and 5-(1H-pyrrol-1-yl)-1,10-phenanthroline (Pyph-Re), were successfully synthesized with the aim to analyze the effect of the pyrrole (Py) moiety on the photophysical properties of Pyph-Re. It was found that the triplet metal-to-ligand charge-transfer d pi (Re) --> pi*(N-N) emission of Phen-Re and Pyph-Re centered at ca. 527 nm with the luminescence quantum yield (LQY) of 0.015 and ca. 578 nm with the LQY of 0.011, respectively. At the same time, the geometrical structures of the ground state and the absorption spectral properties of Phen-Re and Pyph-Re were also calculated with the 6-31G* basis set employed on C, H, N, O, and Br atoms, and LANL2DZ adopted on Re atom.

First Principles Investigation on the Ultra-Incompressible and Hard TaN

The structural, electronic, and mechanical properties of TaN were investigated by use of the density functional theory (DFT). Eight structures were considered, i.e.. hexagonal WC TaN, NiAs, wurtzite, and CoSn structures. cubic NaCl. zinc-blende and CsCl structures. The results indicate that TaN in TaN-type structure is the most stable at ambient conditions among the considered structures. Above 5 GPa, TaN in WC-type structure becomes energetically the most stable phase. They are also stable both thermodynamically and mechanically. TaN in WC-type has the largest shear Modulus 243 GPa and large bulk modulus 337 GPa among the considered structures. The Volume compressibility is slightly larger than diamond, but smaller than c-BN at pressures from 0 to 100 GPa. The compressibility along the c axis is smaller than the linear compressibility of both diamond and c-BN.

Ground state of lutetium dimer by density functional methods

Bond distances, vibrational frequencies and dissociation energies for the ground state of Lu-2 were studied by density functional methods B3LYP, B3PW91, BLYP, BHLYP, BP86, B3P86, MPW1PW91, PBE1PBE and SVWN with CEP-121G and SDD basis sets. Singlet state is predicted to be the most stable. CEP-121G has a better overall performance than SDD. At CEP-121G basis set, all density functional methods used in this study perform well in reproducing the spectroscopic constants.

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYP

Theoretical investigation of LaC2 and LaC2+ clusters

LaC2 (with doublet and quartet states) and LaC2+ (with singlet and triplet states) cluster have been studied by using the B3LYP (Becke three-parameter/Lee-Yang-Parr) density functional method and the HF (Hartree-Fock) method with LANLIDZ basis set. For each cluster, four possible isomers in C-2v, C-s, C-proportional to v and D-proportional to h symmetries have been investigated. The results indicate that structures in C-s symmetry are local minima in all cases and, in most cases (particularly for high spin states), our initial guess in C-s symmetry converges to structures in C-2v symmetry. For the isomers in C-2v, C-proportional to v and D-proportional to h symmetries, local minima were found to be dependent on the method and spin state. The two clusters may also exist as linear chains. The ordering of the binding energies for the isomers in all spin states is C-s similar to C-2v < C-proportional to v < D-proportional to h. The ionization potential of LaC2 is reported as well. (C) 1998 Elsevier Science B.V.

Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC3n+ (n = 0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter-Lee-Yang-Parr) density functional method. The basis set is Dunning/ Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C-2 upsilon symmetry, the other two are Linear chains with C-infinity upsilon symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC3+, doublet and quartet for LaC3 and LaC32+, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C-2 upsilon symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. (C) 1998 John Wiley & Sons, Inc.

Theoretical study on adsorption of Na+ and Na+(H2O)(n) (n=1-6) on a clean Si(111) surface

To model the adsorption of Na+ in aqueous solution on the semiconductor surface, the interactions of Na+ and Na+(H2O)(n) (n = 1-6) with a clean Si(111) surface were investigated by using hybrid density functional theory (B3LYP) and Moller-Plesset second-order perturbation (MP2) methods. The Si(111) surface was described with Si8H12, Si16H20, and Si22H21 Cluster models. The effect of the basis set superposition error (BSSE) was taken into account by applying the counterpoise (CP) correction. The calculated results indicated that the interactions between the Na+ cation and the dangling bonds of the Si(111) surface are primarily electrostatic with partial orbital interactions. The magnitude of the binding energies depends weakly on the adsorption sites and the size of the clusters. When water molecules are present, the interaction between the Nal and Si(I 11) surfaces weakens and the binding energy has the tendency to saturate. On a Si22H21 cluster described surface, the optimized Na+-surface distance for Na+(H2O)(5) adsorbed at on-top site is 4.16 angstrom and the CP-corrected binding energy (MP2) is -35.4 kJ/mol, implying a weakly adsorption of hydrated Na+ cation on clean Si(111) surface.

A Connection Between GRBF and MLP

Both multilayer perceptrons (MLP) and Generalized Radial Basis Functions (GRBF) have good approximation properties, theoretically and experimentally. Are they related? The main point of this paper is to show that for normalized inputs, multilayer perceptron networks are radial function networks (albeit with a non-standard radial function). This provides an interpretation of the weights w as centers t of the radial function network, and therefore as equivalent to templates. This insight may be useful for practical applications, including better initialization procedures for MLP. In the remainder of the paper, we discuss the relation between the radial functions that correspond to the sigmoid for normalized inputs and well-behaved radial basis functions, such as the Gaussian. In particular, we observe that the radial function associated with the sigmoid is an activation function that is good approximation to Gaussian basis functions for a range of values of the bias parameter. The implication is that a MLP network can always simulate a Gaussian GRBF network (with the same number of units but less parameters); the converse is true only for certain values of the bias parameter. Numerical experiments indicate that this constraint is not always satisfied in practice by MLP networks trained with backpropagation. Multiscale GRBF networks, on the other hand, can approximate MLP networks with a similar number of parameters.

Linear Object Classes and Image Synthesis from a Single Example Image

The need to generate new views of a 3D object from a single real image arises in several fields, including graphics and object recognition. While the traditional approach relies on the use of 3D models, we have recently introduced techniques that are applicable under restricted conditions but simpler. The approach exploits image transformations that are specific to the relevant object class and learnable from example views of other "prototypical" objects of the same class. In this paper, we introduce such a new technique by extending the notion of linear class first proposed by Poggio and Vetter. For linear object classes it is shown that linear transformations can be learned exactly from a basis set of 2D prototypical views. We demonstrate the approach on artificial objects and then show preliminary evidence that the technique can effectively "rotate" high- resolution face images from a single 2D view.

RABBIT: A Compiler for SCHEME

We have developed a compiler for the lexically-scoped dialect of LISP known as SCHEME. The compiler knows relatively little about specific data manipulation primitives such as arithmetic operators, but concentrates on general issues of environment and control. Rather than having specialized knowledge about a large variety of control and environment constructs, the compiler handles only a small basis set which reflects the semantics of lambda-calculus. All of the traditional imperative constructs, such as sequencing, assignment, looping, GOTO, as well as many standard LISP constructs such as AND, OR, and COND, are expressed in macros in terms of the applicative basis set. A small number of optimization techniques, coupled with the treatment of function calls as GOTO statements, serve to produce code as good as that produced by more traditional compilers. The macro approach enables speedy implementation of new constructs as desired without sacrificing efficiency in the generated code. A fair amount of analysis is devoted to determining whether environments may be stack-allocated or must be heap-allocated. Heap-allocated environments are necessary in general because SCHEME (unlike Algol 60 and Algol 68, for example) allows procedures with free lexically scoped variables to be returned as the values of other procedures; the Algol stack-allocation environment strategy does not suffice. The methods used here indicate that a heap-allocating generalization of the "display" technique leads to an efficient implementation of such "upward funargs". Moreover, compile-time optimization and analysis can eliminate many "funargs" entirely, and so far fewer environment structures need be allocated at run time than might be expected. A subset of SCHEME (rather than triples, for example) serves as the representation intermediate between the optimized SCHEME code and the final output code; code is expressed in this subset in the so-called continuation-passing style. As a subset of SCHEME, it enjoys the same theoretical properties; one could even apply the same optimizer used on the input code to the intermediate code. However, the subset is so chosen that all temporary quantities are made manifest as variables, and no control stack is needed to evaluate it. As a result, this apparently applicative representation admits an imperative interpretation which permits easy transcription to final imperative machine code. These qualities suggest that an applicative language like SCHEME is a better candidate for an UNCOL than the more imperative candidates proposed to date.

Investigation of numerical atomic orbitals for first-principles calculations of the electronic and transport properties of silicon nanowire structures

This thesis is focused on the application of numerical atomic basis sets in studies of the structural, electronic and transport properties of silicon nanowire structures from first-principles within the framework of Density Functional Theory. First we critically examine the applied methodology and then offer predictions regarding the transport properties and realisation of silicon nanowire devices. The performance of numerical atomic orbitals is benchmarked against calculations performed with plane waves basis sets. After establishing the convergence of total energy and electronic structure calculations with increasing basis size we have shown that their quality greatly improves with the optimisation of the contraction for a fixed basis size. The double zeta polarised basis offers a reasonable approximation to study structural and electronic properties and transferability exists between various nanowire structures. This is most important to reduce the computational cost. The impact of basis sets on transport properties in silicon nanowires with oxygen and dopant impurities have also been studied. It is found that whilst transmission features quantitatively converge with increasing contraction there is a weaker dependence on basis set for the mean free path; the double zeta polarised basis offers a good compromise whereas the single zeta basis set yields qualitatively reasonable results. Studying the transport properties of nanowire-based transistor setups with p+-n-p+ and p+-i-p+ doping profiles it is shown that charge self-consistency affects the I-V characteristics more significantly than the basis set choice. It is predicted that such ultrascaled (3 nm length) transistors would show degraded performance due to relatively high source-drain tunnelling currents. Finally, it is shown the hole mobility of Si nanowires nominally doped with boron decreases monotonically with decreasing width at fixed doping density and increasing dopant concentration. Significant mobility variations are identified which can explain experimental observations.