Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We have compared simulations using solutions of Poisson's equation to detailed capacitance measurements on a double quantum dot structure. We tabulate the results and show which cases show good agreement and which do not. The capacitance values are also compared to those calculated by a solution of Laplace's equation. Electron density is plotted and discussed. In order to understand relevant potential barriers we compare simulations at 50 Kelvin to simulations at 15 Kelvin. We show that the charge density does not differ greatly, but that the conduction band potential does. However, a method of estimating the potential at 0 Kelvin based on the charge distribution at 50 Kelvin is shown to be close to the potential at 15 Kelvin. This method was used to estimate potential barriers at 0 Kelvin in two quantum dot structures.

Abstract not provided.

Abstract not provided.

The observation and characterization of a single atom system in silicon is a significant landmark in half a century of device miniaturization, and presents an important new laboratory for fundamental quantum and atomic physics. We compare with multi-million atom tight binding (TB) calculations the measurements of the spectrum of a single two-electron (2e) atom system in silicon - a negatively charged (D-) gated Arsenic donor in a FinFET. The TB method captures accurate single electron eigenstates of the device taking into account device geometry, donor potentials, applied fields, interfaces, and the full host bandstructure. In a previous work, the depths and fields of As donors in six device samples were established through excited state spectroscopy of the D0 electron and comparison with TB calculations. Using self-consistent field (SCF) TB, we computed the charging energies of the D- electron for the same six device samples, and found good agreement with the measurements. Although a bulk donor has only a bound singlet ground state and a charging energy of about 40 meV, calculations show that a gated donor near an interface can have a reduced charging energy and bound excited states in the D- spectrum. Measurements indeed reveal reduced charging energies and bound 2e excited states, at least one of which is a triplet. The calculations also show the influence of the host valley physics in the two-electron spectrum of the donor.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A manager-worker-based parallelization algorithm for Quantum Monte Carlo (QMC-MW) is presented and compared with the pure iterative parallelization algorithm, which is in common use. The new manager-worker algorithm performs automatic load balancing, allowing it to perform near the theoretical maximal speed even on heterogeneous parallel computers. Furthermore, the new algorithm performs as well as the pure iterative algorithm on homogeneous parallel computers. When combined with the dynamic distributable decorrelation algorithm (DDDA) [Feldmann et al., J Comput Chem 28, 2309 (2007)], the new manager-worker algorithm allows QMC calculations to be terminated at a prespecified level of convergence rather than upon a prespecified number of steps (the common practice). This allows a guaranteed level of precision at the least cost. Additionally, we show (by both analytic derivation and experimental verification) that standard QMC implementations are not perfectly parallel as is often claimed.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A finite temperature version of 'exact-exchange' density functional theory (EXX) has been implemented in Sandia's Socorro code. The method uses the optimized effective potential (OEP) formalism and an efficient gradient-based iterative minimization of the energy. The derivation of the gradient is based on the density matrix, simplifying the extension to finite temperatures. A stand-alone all-electron exact-exchange capability has been developed for testing exact exchange and compatible correlation functionals on small systems. Calculations of eigenvalues for the helium atom, beryllium atom, and the hydrogen molecule are reported, showing excellent agreement with highly converged quantumMonte Carlo calculations. Several approaches to the generation of pseudopotentials for use in EXX calculations have been examined and are discussed. The difficult problem of finding a correlation functional compatible with EXX has been studied and some initial findings are reported.

Abstract not provided.

We report the implementation of an iterative scheme for calculating the Optimized Effective Potential (OEP). Given an energy functional that depends explicitly on the Kohn-Sham wave functions, and therefore, implicitly on the local effective potential appearing in the Kohn-Sham equations, a gradient-based minimization is used to find the potential that minimizes the energy. Previous work has shown how to find the gradient of such an energy with respect to the effective potential in the zero-temperature limit. We discuss a density-matrix-based derivation of the gradient that generalizes the previous results to the finite temperature regime, and we describe important optimizations used in our implementation. We have applied our OEP approach to the Hartree-Fock energy expression to perform Exact Exchange (EXX) calculations. We report our EXX results for common semiconductors and ordered phases of hydrogen at zero and finite electronic temperatures. We also discuss issues involved in the implementation of forces within the OEP/EXX approach.

The optimized effective potential (OEP) method allows orbital-dependent functionals to be used in density functional theory (DFT), which, in particular, allows exact exchange formulations of the exchange energy to be used in DFT calculations. Because the exact exchange is inherently self-interaction correcting, the resulting OEP calculations have been found to yield superior band-gaps for condensed-phase systems. Here we apply these methods to the isolated atoms He and Be, and compare to high quality experiments and calculations to demonstrate that the orbital energies accurately reproduce the excited state spectrum for these species. These results suggest that coupling the exchange-only OEP calculations with proper (orbital-dependent or other) correlation functions might allow quantitative accuracy from DFT calculations.

Abstract not provided.

We derive here the form for the exact exchange energy density for a density that decays with Gaussian-type behavior at long range. This functional is intermediate between the B88 and the PW91 exchange functionals. Using this modified functional to match the form expected for Gaussian densities, we propose the X3LYP extended functional. We find that X3LYP significantly outperforms Becke three parameter Lee-Yang-Parr (B3LYP) for describing van der Waals and hydrogen bond interactions, while performing slightly better than B3LYP for predicting heats of formation, ionization potentials, electron affinities, proton affinities, and total atomic energies as validated with the extended G2 set of atoms and molecules. Thus X3LYP greatly enlarges the field of applications for density functional theory. In particular the success of X3LYP in describing the water dimer (with Re and De within the error bars of the most accurate determinations) makes it an excellent candidate for predicting accurate ligand-protein and ligand-DNA interactions.

Molecular compounds-comprised of mechanically interlocked components-such as rotaxanes and catenanes can be designed to display readily controllable internal movements of one component with respect to the other. Since theweak noncovalent bonding interactions that contribute to the template-directed synthesis of such compounds live on between the components thereafter, they can be activated such that the components move in either a linear fashion (rotaxanes) or a rotary manner (catenanes). These molecules can be activated by switching the recognition elements off and on between components chemically, electrically, or optically, such that they perform motions reminiscent of the moving parts in macroscopic machines. This review will highlight how the emergence ofthe mechanical bond in chemistry during the last two decades has brought with it a real prospect of integrating a bottom-up approach, based on molecular design and micro- and nanofabrication, to construct molecular electronic devices that store information at very high densities using minimal power. Although most of the research reported in this review on switchable catenanes and rotaxanes has been carried out in the context of solution-phase mechanical processes, recent results demonstrate that relative mechanical movements between the components in interlocked molecules can be stimulated (a) chemically in Langmuir and Langmuir-Blodgett films, (b) electrochemically as self-assembled monolayers on gold, and (c) electronically within the settings of solid-state devices. Not only has reversible, electronically driven switching been observed in devices incorporating a bistable [2]catenane, but a crosspoint random access memory circuit has been fabricated using an amphiphilic, bistable [2]rotaxane. The experiments provide strong evidence that switchable catenanes and rotaxanes operate mechanically in a soft-matter environment and can withstand simple device-processing steps. Studies on single-walled carbon nanotubes used as one of the electrodes in molecular switch tunnel junctions have revealed that interfacial chemical interactions involving electrodes containing carbon, silicon, and oxygen are good choices when carrying out molecular electronics on the class of rotaxane- and catenane-based molecules reported in this review. This conclusion is supported by differential conductance measurements (at 4K) made with single-molecule transistors using the break-junction method. It transpires that the electronic transport properties in such devices are more sensitive to the chemical nature of the molecule-electrode contacts than the details of the molecules' electronic structure away from the contacts. This result has profound implications for molecular electronics and highlights the importance of also considering the molecules and the electrodes as an integrated system. It all adds up to an integrated systems-oriented approach to nanotechnology that finds its inspiration in the transfer of concepts like molecular recognition from the life sciences into materials science and provides a model for how, in principle, to transfer elements of traditional chemistry to technology platforms that are being developed on the nanoscale. Before there can be any serious prospect of a technology, there has to be some good, sound science in the making. Molecular electronics is very much in its infancy and, as such, it can be expected to give rise to a great deal of intellectually stimulating science before it stands half a chance of becoming a viable companion to silicon-based technology.

We use the density functional theory and x-ray and neutron diffraction to investigate the crystal structures and reaction mechanisms of intermediate phases likely to be involved in decomposition of the potential hydrogen storage material LiAlH{sub 4}. First, we explore the decomposition mechanism of monoclinic LiAlH4 into monoclinic Li{sub 3}AlH{sub 6} plus face-centered cubic (fcc) Al and hydrogen. We find that this reaction proceeds through a five-step mechanism with an overall activation barrier of 36.9 kcal/mol. The simulated x ray and neutron diffraction patterns from LiAlH{sub 4} and Li{sub 3}AlH{sub 6} agree well with experimental data. On the other hand, the alternative decomposition of LiAlH{sub 4} into LiAlH2 plus H2 is predicted to be unstable with respect to that through Li{sub 3}AlH{sub 6}. Next, we investigate thermal decomposition of Li{sub 3}AlH{sub 6} into fcc LiH plus Al and hydrogen, occurring through a four-step mechanism with an activation barrier of 17.4 kcal/mol for the rate-limiting step. In the first and second steps, two Li atoms accept two H atoms from AlH{sub 6} to form the stable Li-H-Li-H complex. Then, two sequential H2 desorption steps are followed, which eventually result in fcc LiH plus fcc Al and hydrogen: Li{sub 3}AlH{sub 6}(monoclinic) {yields} 3 LiH(fcc) + Al(fcc) + 3/2 H{sub 2} is endothermic by 15.8 kcal/mol. The dissociation energy of 15.8 kcal/mol per formula unit compares to experimental enthalpies in the range of 9.8-23.9 kcal/mol. Finally, we explore thermal decomposition of LiH, LiH(s) + Al(s) {yields} LiAl(s) + 1/2 H{sub 2}(g) is endothermic by 4.6 kcal/mol. The B32 phase, which we predict as the lowest energy structure for LiAl, shows covalent bond characters in the Al-Al direction. Additionally, we determine that transformation of LiH plus Al into LiAlH is unstable with respect to transformation of LiH through LiAl.

We use quantum mechanics to characterize the structure and current-voltage performance of the Stoddart-Heath rotaxane-based programmable electronic switch. We find that the current when the ring is on the DNP is 37?58 times the current when the ring is on the TTF, in agreement with experiment (ratio of 10?100). This establishes the basis for iterative experimental?theoretical efforts to optimize systems for molecule-based electronics which we illustrate by predicting the effect of adding a group such as CN to the rotaxane.

The mechanism of hydroarylation of olefins by a homogeneous Ph-Ir(acac){sub 2}(L) catalyst is elucidated by first principles quantum mechanical methods (DFT), with particular emphasis on activation of the catalyst, catalytic cycle, and interpretation of experimental observations. On the basis of this mechanism, we suggest new catalysts expected to have improved activity. Initiation of the catalyst from the inert trans-form into the active cis-form occurs through a dissociative pathway with a calculated {Delta}H(0 K){sub {+-}} = 35.1 kcal/mol and {Delta}G(298 K){sub {+-}} = 26.1 kcal/mol. The catalytic cycle features two key steps, 1,2-olefin insertion and C?H activation via a novel mechanism, oxidative hydrogen migration. The olefin insertion is found to be rate determining, with a calculated {Delta}H(0 K){sub {+-}} = 27.0 kcal/mol and {Delta}G(298 K){sub {+-}} = 29.3 kcal/mol. The activation energy increases with increased electron density on the coordinating olefin, as well as increased electron-donating character in the ligand system. The regioselectivity is shown to depend on the electronic and steric characteristics of the olefin, with steric bulk and electron withdrawing character favoring linear product formation. Activation of the C?H bond occurs in a concerted fashion through a novel transition structure best described as an oxidative hydrogen migration. The character of the transition structure is seven coordinate Ir{sup V}, with a full bond formed between the migrating hydrogen and iridium. Several experimental observations are investigated and explained: (a) The nature of L influences the rate of the reaction through a ground-state effect. (b) The lack of {beta}-hydride products is due to kinetic factors, although {beta}-hydride elimination is calculated to be facile, all further reactions are kinetically inaccessible. (c) Inhibition by excess olefin is caused by competitive binding of olefin and aryl starting materials during the catalytic cycle in a statistical fashion. On the basis of this insertion-oxidative hydrogen transfer mechanism we suggest that electron-withdrawing substituents on the acac ligands, such as trifluoromethyl groups, are good modifications for catalysts with higher activity.

Abstract not provided.

Abstract not provided.

{sm_bullet}HF/DFT are one-particle approximation to the Schrodinger equation {sm_bullet} The one-particle, mean field approaches are what lead to the nonlinear eigenvalue problem {sm_bullet} DFT includes a parameterized XC functional that reproduces many-electron effects -Very accurate ground state structures and energies - Problematic for excited states, band gaps