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.

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.

12C2H2-Ar van der Waals complex

New theoretical and experimental results on the acetylene-Ar van der Waals complex are presented and the literature is reviewed. New ab initio calculations at the MP2 level were performed using large basis sets with diffuse functions and taking into account the basis set superposition error. It was found that the structure of acetylene is not significantly altered by the complexation and that its vibrational frequencies are only slightly lowered. Finally, it was observed that the calculated properties of the complex (structure, vibrational spectrum, bond dissociation energy) are not sensitive to the structure imposed on acetylene. Experimentally, acetylene-Ar was produced in a supersonic expansion under experimental conditions corresponding to 9 K rotational temperature. Thanks to the performances of CW-CRDS detection, the Ka = 0 ← 1, 1 ← 0, and 2 ← 1 sub- bands of the v1 + v3 band could be recorded and resolved and most of their lines assigned. Upper-state rotational constants were fitted, however not including the upper Ka = 2 state, which shows K-doubling the opposite of the expected. The Lorentzian width of most line profiles sets the mean lifetime to some 7.5 ns. Local perturbations affecting line positions and/or line widths are demonstrated. Additional series of lines tentatively attributed to acetylene-Ar are discussed.© 2009 American Chemical Society.

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C 2H2-Kr and C2H2-Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n.1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C2H2 moiety corresponds to distances of 4.064 and 4.229Å, and angles of 65.22° and 68.67° for C 2H2-Kr and C2H2-Xe, respectively. The interaction energy of both complexes is estimated to be -151.88 (1.817 kJ mol-1) and -182.76 cm-1 (2.186 kJ mol-1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed. © 2012 Taylor and Francis.

Equilibrium structure of the hydrogen bonded dimer H2O ⋯ HF

The equilibrium structure of the hydrogen bonded complex H2O HF has been calculated ab initio using the CCSD(T) method with basis sets up to sextuple- quality with diffuse functions and taking into account the basis set superposition error correction. The calculations carried out confirm the importance of diffuse functions and of counterpoise correction to obtain an accurate geometry. The most important point is that the basis set convergence is extremely slow and, for this reason an accurate ab initio structure requires a very large basis set. Nevertheless, the ab initio structure is significantly different from the experimental r0 and rm structures. Analysis of the basis set convergence and of the approximations used for the determination of the experimental structures indicates that the ab initio structure is expected to be more reliable.

Comparison of the experimental, semi-experimental and ab initio equilibrium structures of acetylene: Influence of relativisitic effects and of the diagonal Born-Oppenheimer corrections

The equilibrium structure of acetylene (also named ethyne) has been reinvestigated to resolve the small discrepancies noted between different determinations. The size of the system as well as the large amount of available experimental data provides the quite unique opportunity to check the magnitude and relevance of various contributions to equilibrium structure as well as to verify the accuracy of experimental results. With respect to pure theoretical investigation, quantum-chemical calculations at the coupled-cluster level have been employed together with extrapolation to the basis set limit, consideration of higher excitations in the cluster operator, inclusion of core correlation effects as well as relativistic and diagonal Born-Oppenheimer corrections. In particular, it is found that the extrapolation to the complete basis set limit, the inclusion of higher excitations in the electronic-correlation treatment and the relativistic corrections are of the same order of magnitude. It also appears that a basis set as large as a core-valence quintuple-zeta set is required for accurately accounting for the inner-shell correlation contribution. From a pure experimental point of view, the equilibrium structure has been determined using very accurate rotational constants recently obtained by a global analysis (that is to say that all non-negligible interactions are explicitely included in the Hamiltonian matrix) of rovibrational spectra. Finally, a semi-experimental equilibrium structure (where the equilibrium rotational constants are obtained from the experimental ground state rotational constants and computed rovibrational corrections) has been obtained from the available experimental ground-state rotational constants for ten isotopic species corrected for computed vibrational corrections. Such a determination led to the revision of the ground-state rotational constants of two isotopologues, thus showing that structural determination is a good method to identify errors in experimental rotational constants. The three structures are found in a very good agreement, and our recommended values are rCC 120.2958(7) pm and rCH 106.164(1) pm. © 2011 American Institute of Physics.

Energy levels and transition probabilities for nitrogen-like Fe xx

Energies of the 700 lowest levels in Fe XX have been obtained using the multiconfiguration Dirac-Fock method. Configuration interaction method on the basis set of transformed radial orbitals with variable parameters taking into account relativistic corrections in the Breit-Pauli approximation was used to crosscheck our presented results. Transition probabilities, oscillator and line strengths are presented for electric dipole (E1), electric quadrupole (E2) and magnetic dipole (M1) transitions among these levels. The total radiative transition probabilities from each level are also provided. Results are compared with data compiled by NIST and with other theoretical work.

Multi-state CDW approximation to electron capture during slow p-H collisions

The first complete multi-state CDW close coupling calculations which use a fully normalized basis set are performed. The results obtained at impact energies in the region of 10 keV for total and n = 2 capture cross sections are in reasonably good accord with experiment despite the fact that only the ground states of both species and the n = 2 states of the projectile are incorporated into the model. The theory has significant advantages over other atomic and molecular expansions which may require extensive bases to obtain similar accuracy.

Prediction of the formation and stabilities of energetic salts and ionic liquids based on ab initio electronic structure calculations

A computational approach to predict the thermodynamics for forming a variety of imidazolium-based salts and ionic liquids from typical starting materials is described. The gas-phase proton and methyl cation acidities of several protonating and methylating agents, as well as the proton and methyl cation affinities of many important methyl-, nitro-, and cyano- substituted imidazoles, have been calculated reliably by using the computationally feasible DFT (B3LYP) and MP2 (extrapolated to the complete basis set limit) methods. These accurately calculated proton and methyl cation affinities of neutrals and anions are used in conjunction with an empirical approach based on molecular volumes to estimate the lattice enthalpies and entropies of ionic liquids, organic solids, and organic liquids. These quantities were used to construct a thermodynamic cycle for salt formation to reliably predict the ability to synthesize a variety of salts including ones with potentially high energetic densities. An adjustment of the gas phase thermodynamic cycle to account for solid- and liquid-phase chemistries provides the best overall assessment of salt formation and stability. This has been applied to imidazoles (the cation to be formed) with alkyl, nitro, and cyano substituents. The proton and methyl cation donors studied were as follows: HCl, HBr, HI, (HO)(2)SO2, HSO3CF3 (TfOH), and HSO3(C6H4)CH3 (TsOH); CH3Cl, CH3Br, CH3I, (CH3O)(2)SO2, CH3SO3CF3 (TfOCH3) and CH3SO3(C6H4)CH3 (TsOCH3). As substitution of the cation with electron-withdrawing groups increases, the triflate reagents appear to be the best overall choice as protonating and methylating agents. Even stronger alkylating agents should be considered to enhance the chances of synthetic success. When using the enthalpies of reaction for the gas-phase reactants (eq 6) to form a salt, a cutoff value of - 13 kcal mol(-1) or lower (more negative) should be used as the minimum value for predicting whether a salt can be synthesized.

Lowest energy potentials of Li2 + F accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Electronic states and spin-forbidden cooling transitions of AlH and AlF

The feasibility of laser cooling AlH and AlF is investigated using ab initio quantum chemistry. All the electronic states corresponding to the ground and lowest two excited states of the Al atom are calculated using multi-reference configuration interaction (MRCI) and the large AV6Z basis set for AlH. The smaller AVQZ basis set is used to calculate the valence electronic states of AlF. Theoretical Franck-Condon factors are determined for the A(1)Pi -> X(1)Sigma(+) transitions in both radicals and found to agree with the highly diagonal factors found experimentally, suggesting computational chemistry is an effective method for screening suitable laser cooling candidates. AlH does not appear to have a transition quite as diagonal as that in SrF (which has been laser cooled) but the A(1)Pi -> X(1)Sigma(+) transition transition of AlF is a strong candidate for cooling with just a single laser, though the cooling frequency is deep in the UV. Furthermore, the a (3)Pi -> X(1)Sigma(+) transitions are also strongly diagonal and in AlF is a practical method for obtaining very low final temperatures around 3 mu K.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

Ab initio potentials of F + Li-2 accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li-2 + F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li-2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.