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.

Time delay between photoemission from the 2p and 2s subshells of neon

The R-matrix incorporating time (RMT) method is a method developed recently for solving the time-dependent Schrödinger equation for multielectron atomic systems exposed to intense short-pulse laser light. We have employed the RMT method to investigate the time delay in the photoemission of an electron liberated from a 2p orbital in a neon atom with respect to one released from a 2s orbital following absorption of an attosecond xuv pulse. Time delays due to xuv pulses in the range 76-105 eV are presented. For an xuv pulse at the experimentally relevant energy of 105.2 eV, we calculate the time delay to be 10.2±1.3 attoseconds (as), somewhat larger than estimated by other theoretical calculations, but still a factor of 2 smaller than experiment. We repeated the calculation for a photon energy of 89.8 eV with a larger basis set capable of modeling correlated-electron dynamics within the neon atom and the residual Ne ion. A time delay of 14.5±1.5 as was observed, compared to a 16.7±1.5 as result using a single-configuration representation of the residual Ne⁺ ion. 

Recoil-induced electronic excitation and ionization in one-and two-electron ions

The electronic redistribution of an ion or atom induced by a sudden recoil of the nucleus occurring during the emission or capture of a neutral particle is theoretically investigated. For one-electron systems, analytical expressions are derived for the electronic transition probabilities to bound and continuum states. The quality of a B-spline basis set approach is evaluated from a detailed comparison with the analytical results. This numerical approach is then used Io study the dynamics of two-electron systems (neutral He and Ne ) using correlated wavefunctions for both the target and daughter ions. The total transition probabilities to discrete states, autoionizing states and direct single- and double-ionization probabilities are calculated from the pseudospectra. Sum rules for transition probabilities involving an initial bound state and a complete final series are discussed.

Modeling study of woody biomass:Interactions of cellulose, hemicellulose, and lignin

Lignocellulosic biomass pretreatment and the subsequent thermal conversion processes to produce solid, liquid, and gas biofuels are attractive solutions for today's energy challenges. The structural study of the main components in biomass and their macromolecular complexes is an active and ongoing research topic worldwide. The interactions among the three main components, cellulose, hemicellulose, and lignin, are studied in this paper using electronic structure methods, and the study includes examining the hydrogen bond network of cellulose-hemicellulose systems and the covalent bond linkages of hemicellulose-lignin systems. Several methods (semiempirical, Hartree-Fock, and density functional theory) using different basis sets were evaluated. It was shown that theoretical calculations can be used to simulate small model structures representing wood components. By comparing calculation results with experimental data, it was concluded that B3LYP/6-31G is the most suitable basis set to describe the hydrogen bond system and B3LYP/6-31G(d,p) is the most suitable basis set to describe the covalent system of woody biomass. The choice of unit model has a much larger effect on hydrogen bonding within cellulose-hemicellulose system, whereas the model choice has a minimal effect on the covalent linkage in the hemicellulose-lignin system. © 2011 American Chemical Society.

A Fast Algorithm for Sparse Support Vector Machines for Mobile Computing Applications

This research presents a fast algorithm for projected support vector machines (PSVM) by selecting a basis vector set (BVS) for the kernel-induced feature space, the training points are projected onto the subspace spanned by the selected BVS. A standard linear support vector machine (SVM) is then produced in the subspace with the projected training points. As the dimension of the subspace is determined by the size of the selected basis vector set, the size of the produced SVM expansion can be specified. A two-stage algorithm is derived which selects and refines the basis vector set achieving a locally optimal model. The model expansion coefficients and bias are updated recursively for increase and decrease in the basis set and support vector set. The condition for a point to be classed as outside the current basis vector and selected as a new basis vector is derived and embedded in the recursive procedure. This guarantees the linear independence of the produced basis set. The proposed algorithm is tested and compared with an existing sparse primal SVM (SpSVM) and a standard SVM (LibSVM) on seven public benchmark classification problems. Our new algorithm is designed for use in the application area of human activity recognition using smart devices and embedded sensors where their sometimes limited memory and processing resources must be exploited to the full and the more robust and accurate the classification the more satisfied the user. Experimental results demonstrate the effectiveness and efficiency of the proposed algorithm. This work builds upon a previously published algorithm specifically created for activity recognition within mobile applications for the EU Haptimap project [1]. The algorithms detailed in this paper are more memory and resource efficient making them suitable for use with bigger data sets and more easily trained SVMs.

Elliptical orbital studies of H21 and H2 molecules /|nS. K. Gupta. -- 260 St. Catharines, Ont. : [s. n.],

Calculations are performed on the \S <:Jd ground states of d ' + the H and HC) molecules using a basis set of non-integral ~ ~ I elliptical orbitals. Different variational wavefunctions constructed i- for H~ involved one parameter to three par~~eter variation. In order to l"'educe the ntunber of parameters in most commonly 0- used basis orbitals set, the importance of the term (,+~) Y\ over the term ;u 'Where n is a variational pararneter and the value of cr may be given by boundary condition or cusp condition is outlined in Chapters II and III. It is found that the two parameter -+ wavefunction for H~ including the ternl (~+~) , a- given by the bound~ condition, gives lower variational energies than any wavefunction published to date for small and moderate internuclear separations. c;. In order to find out the importance of the term (I +~ ) Y\ over ;U for the two electron problem, the variational energy is computed for the H~ molecule from unrestricted two parameter closed shell wavefunctions including the term U+ft)

An analysis of the methyl rotation dynamics in the So (x' A) and T (a A) states of thioacetaldehyde from Pyrolysis jet spectra /

Jet-cooled, laser-induced phosphorescence excitation spectra (LIP) of thioacetaldehyde CH3CHS, CH3CDS, CD3CHS and CD3CDS have been observed over the region 15800 - 17300 cm"^ in a continuous pyrolysis jet. The vibronic band structure of the singlet-triplet n -* n* transition were attributed to the strong coupling of the methyl torsion and aldehydic hydrogen wagging modes . The vibronic peaks have been assigned in terms of two upper electronic state (T^) vibrations; the methyl torsion mode v^g, and the aldehydic hydrogen wagging mode v^^. The electronic origin O^a^ is unequivocally assigned as follows: CH3CHS (16294.9 cm"'' ), CH3CDS (16360.9 cm"'' ), CD3CHS (16299.7 cm"^ ), and CD3CDS (16367.2 cm"'' ). To obtain structural and dynamical information about the two electronic states, potential surfaces V(e,a) for the 6 (methyl torsion) and a (hydrogen wagging) motions were generated by ab initio quantum mechanical calculations with a 6-3 IG* basis in which the structural parameters were fully relaxed. The kinetic energy coefficients BQ(a,e) , B^(a,G) , and the cross coupling term B^(a,e) , were accurately represented as functions of the two active coordinates, a and 9. The calculations reveal that the molecule adopts an eclipsed conformation for the lower Sq electronic state (a=0°,e=0"') with a barrier height to internal rotation of 541.5 cm"^ which is to be compared to 549.8 cm"^ obtained from the microwave experiment. The conformation of the upper T^ electronic state was found to be staggered (a=24 . 68° ,e=-45. 66° ) . The saddle point in the path traced out by the aldehyde wagging motion was calculated to be 175 cm"^ above the equilibrium configuration. The corresponding maxima in the path taken by methyl torsion was found to be 322 cm'\ The small amplitude normal vibrational modes were also calculated to aid in the assignment of the spectra. Torsional-wagging energy manifolds for the two states were derived from the Hamiltonian H(a,e) which was solved variationally using an extended two dimensional Fourier expansion as a basis set. A torsionalinversion band spectrum was derived from the calculated energy levels and Franck-Condon factors, and was compared with the experimental supersonic-jet spectra. Most of the anomalies which were associated with the interpretation of the observed spectrum could be accounted for by the band profiles derived from ab initio SCF calculations. A model describing the jet spectra was derived by scaling the ab initio potential functions. The global least squares fitting generates a triplet state potential which has a minimum at (a=22.38° ,e=-41.08°) . The flatter potential in the scaled model yielded excellent agreement between the observed and calculated frequency intervals.

Application of Reptation Quantum Monte Carlo and Related Methods to LiH

This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.

Experimental and Computational Studies on Copper and Hydrobenzoin Derivatives in Catalytic Asymmetric Reactions

The use of theory to understand and facilitate catalytic enantioselective organic transformations involving copper and hydrobenzoin derivatives is reported. Section A details the use of theory to predict, facilitate, and understand a copper promoted amino oxygenation reaction reported by Chemler et al. Using Density Functional Theory (DFT), employing the hybrid B3LYP functional and a LanL2DZ/6-31G(d) basis set, the mechanistic details were studied on a N-tosyl-o-allylaniline and a [alpha]-methyl-[gamma]-alkenyl sulfonamide substrate. The results suggest the N-C bond formation proceeds via a cisaminocupration, and not through a radical-type mechanism. Additionally, the origin of diastereoselection observed with [alpha]-methyl-[gamma]-alkenyl sulfonamide arises from avoidance of unfavourable steric interactions between the methyl substituent and the N -protecting group. Section B details the computationally guided, experimental investigation of two hydrobenzoin derivatives as ligands/ catalysts, as well as the attempted synthesis of a third hydrobenzoin derivative. The bis-boronic acid derived from hydrobenzoin was successful as a Lewis acid catalyst in the Bignielli reaction and the Conia ene reaction, but provided only racemic products. The chiral diol derived from hydrobenzoin successfully increased the rate of the addition of diethyl zinc to benzaldehyde in the presence of titanium tetraisopropoxide, however poor enantioinduction was obseverved. Notably, the observed reactivity was successfully predicted by theoretical calculations.

Méthode de calcul à N-corps basée sur la G0W0 : étude du couplage électron-phonon dans le C60 et développement d’une approche accélérée pour matériaux organiques

La présente thèse porte sur les limites de la théorie de la fonctionnelle de la densité et les moyens de surmonter celles-ci. Ces limites sont explorées dans le contexte d'une implémentation traditionnelle utilisant une base d'ondes planes. Dans un premier temps, les limites dans la taille des systèmes pouvant être simulés sont observées. Des méthodes de pointe pour surmonter ces dernières sont ensuite utilisées pour simuler des systèmes de taille nanométrique. En particulier, le greffage de molécules de bromophényle sur les nanotubes de carbone est étudié avec ces méthodes, étant donné l'impact substantiel que pourrait avoir une meilleure compréhension de ce procédé sur l'industrie de l'électronique. Dans un deuxième temps, les limites de précision de la théorie de la fonctionnelle de la densité sont explorées. Tout d'abord, une étude quantitative de l'incertitude de cette méthode pour le couplage électron-phonon est effectuée et révèle que celle-ci est substantiellement plus élevée que celle présumée dans la littérature. L'incertitude sur le couplage électron-phonon est ensuite explorée dans le cadre de la méthode G0W0 et cette dernière se révèle être une alternative substantiellement plus précise. Cette méthode présentant toutefois de sévères limitations dans la taille des systèmes traitables, différents moyens théoriques pour surmonter ces dernières sont développés et présentés dans cette thèse. La performance et la précision accrues de l'implémentation résultante laissent présager de nouvelles possibilités dans l'étude et la conception de certaines catégories de matériaux, dont les supraconducteurs, les polymères utiles en photovoltaïque organique, les semi-conducteurs, etc.

Enhancement of Microwave Properties of Planar Filters and Antennas using Photonic Bandgap (PBG) Structures

In this thesis, we explore the design, computation, and experimental analysis of photonic crystals, with a special emphasis on structures and devices that make a connection with practically realizable systems. First, we analyze the propenies of photonic-crystal: periodic dielectric structures that have a band gap for propagation. The band gap of periodically loaded air column on a dielectric substrate is computed using Eigen solvers in a plane wave basis. Then this idea is extended to planar filters and antennas at microwave regime. The main objectives covered in this thesis are:• Computation of Band Gap origin in Photonic crystal with the abet of Maxwell's equation and Bloch-Floquet's theorem • Extension of Band Gap to Planar structures at microwave regime • Predict the dielectric constant - synthesized dieletric cmstant of the substrates when loaded with Photonic Band Gap (PBG) structures in a microstrip transmission line • Identify the resonant characteristic of the PBG cell and extract the equivalent circuit based on PBG cell and substrate parameters for microstrip transmission line • Miniaturize PBG as Defected Ground Structures (DGS) and use the property to be implemented in planar filters with microstrip transmission line • Extended the band stop effect of PBG / DGS to coplanar waveguide and asymmetric coplanar waveguide. • Formulate design equations for the PBG / DGS filters • Use these PBG / DGS ground plane as ground plane of microstrip antennas • Analysis of filters and antennas using FDID method

Theoretical study of the fluorinating power of SF6 to vanadocene

Màster Oficial en Química Teòrica i Computacional Curs: 2008-2009, Director: Juan J. Novoa Vide

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested