Ab inito investigation of hydrodesulferization
Tilson, J.L.; Marshall, C.L.; Brenner, J.R.
1997-12-31
We utilize MPP and vector computers to model the interaction of large sulfur-containing species bonded with hydrodesulfurization (HDS) catalysts. This work is made possible by the availability large aggregate memory, parallel computers. The ability of modern non-traditional computers to solve large-scale scientific problems has been demonstrated. This success is accomplished, in part, by access to portable low- and user-level software tools which exhibit good control over NUMA. ab initio SCF methods are used to characterize the energies of adsorption of poly-aromatic, sulfur-containing hydrocarbons onto a series of molecular-based models of MoS2. These hydrocarbon include thiophene (TP), benzothiophene (BT) and dibenzothiophene (DBT) which are representative of heavy crude feedstocks. Our work attempts to ascertain if a consistent model of adsorption onto a MoS2 surface exists and to determine how the electronic and steric configuration of each species affect their energies of interaction with the metal surface.
Ab-inito calculation of energy level alignment and vacuum level shift at CuPc/C60 interfaces
NASA Astrophysics Data System (ADS)
Sai, Na; Zhu, Xiaoyang; Chelikowsky, James; Leung, Kevin
2012-02-01
The alignment of the donor and acceptor enegy levels is of crucial importance for organic photovotaic performance. We investigate the interfaical electronic structure and energy level alignment of copper phthalocyanine (CuPc)/fullerene (C60) using ab-inito density functional theory calculations including van der Waals interactions and hybrid density functionals. We show that energy level alignment critically depends on the standing-up and lying-down orientation of the CuPc molecules relative to C60 at the interface. We calculate the magnitude of the interface dipole at different molecular orientations and compare them to the vacuum level shift observed in photoemission spectroscopy. The validity of existing theoretical models which invoke charge transfer on this organic interface will be discussed in light of our predictions. Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Deparment of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Ab inito molecular-dynamics study of EC decomposition process on Li2O2 surfaces
NASA Astrophysics Data System (ADS)
Ando, Yasunobu; Ikeshoji, Tamio; Otani, Minoru
2015-03-01
We have simulated electrochemical reactions of the EC molecule decomposition on Li2O2 substrate by ab initio molecular dynamics combined with the effective screening medium method. EC molecules adsorb onto the peroxide spontaneously. We find through the analysis of density of states that the adsorption state is stabilized by hybridization of the sp2 orbital and the surface states of the Li2O2. After adsorption, EC ring opens, which leads to the decomposition of the peroxide and the formation of a carboxy group. This kind of alkyl carbonates formed on the Li2O2 substrate was found in experiments actually Nanosystem Research Institute, AIST; ESICB, Kyoto University
Ab initio quantum chemistry: Methodology and applications
Friesner, Richard A.
2005-01-01
This Perspective provides an overview of state-of-the-art ab initio quantum chemical methodology and applications. The methods that are discussed include coupled cluster theory, localized second-order Moller–Plesset perturbation theory, multireference perturbation approaches, and density functional theory. The accuracy of each approach for key chemical properties is summarized, and the computational performance is analyzed, emphasizing significant advances in algorithms and implementation over the past decade. Incorporation of a condensed-phase environment by means of mixed quantum mechanical/molecular mechanics or self-consistent reaction field techniques, is presented. A wide range of illustrative applications, focusing on materials science and biology, are discussed briefly. PMID:15870212
Quantum phases of AB 2 fermionic chains
NASA Astrophysics Data System (ADS)
Murcia-Correa, L. S.; Franco, R.; Silva-Valencia, J.
2016-02-01
A fermionic chain is a one-dimensional system with fermions that interact locally and can jump between sites in the lattice, in particular an AB n chain type, where A and B are sites that exhibit a difference in energy level of Δ and site B is repeated n-times, such that the unit cell has n +1 sites. A limit case of this model, called the ionic Hubbard model (n = 1), has been widely studied due to its interesting physics and applications. In this paper, we study the ground state of an AB 2 chain, which describes the material R 4[Pt 2(P 2O5H2)4X] · nH 2 O. Specifically, we consider a filling with two electrons per unit cell, and using the density matrix renormalization group method we found that the system exhibits the band insulator and Mott correlated insulator phases, as well as an intermediate phase between them. For couplings of Δ = 2,10 and 20, we estimate the critical points that separate these phases through the structure factor and the energy gap in the sector of charge and spin, finding that the position of the critical point rises as a function of Δ.
Ab Initio Computation of the Energies of Circular Quantum Dots
Lohne, M. Pedersen; Hagen, Gaute; Hjorth-Jensen, M.; Kvaal, S.; Pederiva, F.
2011-01-01
We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.
The implementation of ab initio quantum chemistry calculations on transporters.
Cooper, M D; Hillier, I H
1991-06-01
The RHF and geometry optimization sections of the ab initio quantum chemistry code, GAMESS, have been optimized for a network of parallel microprocessors, Inmos T800-20 transputers, using both indirect and direct SCF techniques. The results indicate great scope for implementation of such codes on small parallel computer systems, very high efficiencies having been achieved, particularly in the cases of direct SCF and geometry optimization with large basis sets. The work, although performed upon one particular parallel system, the Meiko Computing Surface, is applicable to a wide range of parallel systems with both shared and distributed memory. PMID:1919615
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
Makhov, Dmitry V.; Shalashilin, Dmitrii V.; Glover, William J.; Martinez, Todd J.
2014-08-07
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
NASA Astrophysics Data System (ADS)
Makhov, Dmitry V.; Glover, William J.; Martinez, Todd J.; Shalashilin, Dmitrii V.
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics.
Makhov, Dmitry V; Glover, William J; Martinez, Todd J; Shalashilin, Dmitrii V
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. PMID:25106573
Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.
Ab initio quantum dynamics using coupled-cluster.
Kvaal, Simen
2012-05-21
The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given. PMID:22612082
Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2004-01-01
We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
Ab initio quantum transport calculations using plane waves
NASA Astrophysics Data System (ADS)
Garcia-Lekue, A.; Vergniory, M. G.; Jiang, X. W.; Wang, L. W.
2015-08-01
We present an ab initio method to calculate elastic quantum transport at the nanoscale. The method is based on a combination of density functional theory using plane wave nonlocal pseudopotentials and the use of auxiliary periodic boundary conditions to obtain the scattering states. The method can be applied to any applied bias voltage and the charge density and potential profile can either be calculated self-consistently, or using an approximated self-consistent field (SCF) approach. Based on the scattering states one can straightforwardly calculate the transmission coefficients and the corresponding electronic current. The overall scheme allows us to obtain accurate and numerically stable solutions for the elastic transport, with a computational time similar to that of a ground state calculation. This method is particularly suitable for calculations of tunneling currents through vacuum, that some of the nonequilibrium Greens function (NEGF) approaches based on atomic basis sets might have difficulty to deal with. Several examples are provided using this method from electron tunneling, to molecular electronics, to electronic devices: (i) On a Au nanojunction, the tunneling current dependence on the electrode-electrode distance is investigated. (ii) The tunneling through field emission resonances (FERs) is studied via an accurate description of the surface vacuum states. (iii) Based on quantum transport calculations, we have designed a molecular conformational switch, which can turn on and off a molecular junction by applying a perpendicular electric field. (iv) Finally, we have used the method to simulate tunnel field-effect transistors (TFETs) based on two-dimensional transition-metal dichalcogenides (TMDCs), where we have studied the performance and scaling limits of such nanodevices and proposed atomic doping to enhance the transistor performance.
Microsolvation of methyl hydrogen peroxide: Ab initio quantum chemical approach
NASA Astrophysics Data System (ADS)
Kulkarni, Anant D.; Rai, Dhurba; Bartolotti, Libero J.; Pathak, Rajeev K.
2009-08-01
Methyl hydrogen peroxide (MHP), one of the simplest organic hydroperoxides, is a strong oxidant, with enhanced activity in aqueous ambience. The present study investigates, at the molecular level, the role of hydrogen bonding that is conducive to cluster formation of MHP with water molecules from its peroxide end, with the methyl group remaining hydrophobic for up to five water molecules. Ab initio quantum chemical computations on MHP⋯(H2O)n, [n =1-5] are performed at second order Møller-Plesset (MP2) perturbation theory employing the basis sets 6-31G(d,p) and 6-311++G(2d,2p) to study the cluster formation of MHP with water molecules from its peroxide end and hydrophobic hydration due to the methyl group. Successive addition of water molecules alters the hydrogen bonding pattern, which leads to changes in overall cluster geometry and in turn to IR vibrational frequency shifts. Molecular co-operativity in these clusters is gauged directly through a detailed many-body interaction energy analysis. Molecular electrostatic potential maps are shown to have a bearing on predicting further growth of these clusters, which is duly corroborated through sample calculations for MHP⋯(H2O)8. Further, a continuum solvation model calculation for energetically stable clusters suggests that this study should serve as a precursor for pathways to aqueous solvation of MHP.
Heats of Segregation of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
Ab inito study of Ag-related defects in ZnO
NASA Astrophysics Data System (ADS)
Wan, Qixin; Xiong, Zhihua; Li, Dongmei; Liu, Guodong
2008-12-01
Using first-principles calculations, we investigated the structure and electronic properties of Ag-related defects in ZnO. The calculation results indicate that AgZn behaves as acceptor. Simultaneously, by comparing the formation energy and electronic structure of Ag-related defects in ZnO, Oi-AgZn behaves as acceptor in Ag-doped ZnO and it is better to gain p-type ZnO. However, Hi-AgZn complex has the lowest formation energy. Thus, the formation of the other point defects is greatly suppressed by the formation of Hi in Ag-doped ZnO. Moreover, the H atoms can be easily dissociated from hydrogen-passivated complexes by post-annealing at moderate temperatures, thus, codoping Ag with H may be a good method to achieve p-type in Ag-doped ZnO.
Quantum fluctuations and isotope effects in ab initio descriptions of water
Wang, Lu; Markland, Thomas E.; Ceriotti, Michele
2014-09-14
Isotope substitution is extensively used to investigate the microscopic behavior of hydrogen bonded systems such as liquid water. The changes in structure and stability of these systems upon isotope substitution arise entirely from the quantum mechanical nature of the nuclei. Here, we provide a fully ab initio determination of the isotope exchange free energy and fractionation ratio of hydrogen and deuterium in water treating exactly nuclear quantum effects and explicitly modeling the quantum nature of the electrons. This allows us to assess how quantum effects in water manifest as isotope effects, and unravel how the interplay between electronic exchange and correlation and nuclear quantum fluctuations determine the structure of the hydrogen bond in water.
Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo
Zen, Andrea; Luo, Ye Mazzola, Guglielmo Sorella, Sandro; Guidoni, Leonardo
2015-04-14
Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J. P.; Maris, P.; Honkanen, H.; Li, J.; Shirokov, A. M.; Brodsky, S. J.; Harindranath, A.
2009-12-17
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually, we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J.P.; Maris, P.; Shirokov, A.M.; Honkanen, H.; li, J.; Brodsky, S.J.; Harindranath, A.; Teramond, G.F.de; /Costa Rica U.
2009-08-03
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually,we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
Ab Initio study of multiple exciton generation in layered structure quantum dots
NASA Astrophysics Data System (ADS)
Zhang, Zhiyong; Zimmerman, Paul; Cui, Yi; Musgrave, Charles
2011-03-01
Multiple Exciton Generation (MEG) can potentially increase the photovoltaic conversion efficiency significantly and has been reported in a large number of systems and has been extensively studies theoretically and experimentally. Here we report our study of the MEG process in inorganic layered structure quantum dots using high level Ab Initio methods that are capable of electronic states of multi-exciton in character. Our results show that multiple states that are of multi-exciton character exist in quantum dots and different mechanisms govern the MEG process in quantum dots: (1) MEG through an internal crossing mechanism from a optically active state to an optically dark multi-exciton state, as in the singlet fission process of pentacene; and (2) direct multi-exciton generation through an optically active excited state. We also discuss detailed structure evolution of quantum dots, from stable molecular like structures of various shapes and sizes, to larger quantum dots of bulk like bonding motifs with distinctive surface structures and illustrate the correlation between structure and the multi-exciton states.
NASA Astrophysics Data System (ADS)
John, Christopher; Spura, Thomas; Habershon, Scott; Kühne, Thomas D.
2016-04-01
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Amorphous Ge quantum dots embedded in crystalline Si: ab initio results.
Laubscher, M; Küfner, S; Kroll, P; Bechstedt, F
2015-10-14
We study amorphous Ge quantum dots embedded in a crystalline Si matrix through structure modeling and simulation using ab initio density functional theory including spin-orbit interaction and quasiparticle effects. Three models are generated by replacing a spherical region within diamond Si by Ge atoms and creating a disordered bond network with appropriate density inside the Ge quantum dot. After total-energy optimisations of the atomic geometry we compute the electronic and optical properties. We find three major effects: (i) the resulting nanostructures adopt a type-I heterostructure character; (ii) the lowest optical transitions occur only within the Ge quantum dots, and do not involve or cross the Ge-Si interface. (iii) for larger amorphous Ge quantum dots, with diameters of about 2.0 and 2.7 nm, absorption peaks appear in the mid-infrared spectral region. These are promising candidates for intense luminescence at photon energies below the gap energy of bulk Ge. PMID:26402441
Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides
NASA Astrophysics Data System (ADS)
Sharma, Vinit; Krogel, Jaron T.; Kent, P. R. C.; Reboredo, Fernando A.
One of the critical scientific challenges of contemporary research is to obtain an accurate theoretical description of the electronic properties of strongly correlated systems such as transition metal oxides and rare-earth compounds, since state-of-art ab-initio methods based on approximate density functionals are not always sufficiently accurate. Quantum Monte Carlo (QMC) methods, which use statistical sampling to evaluate many-body wave functions, have the potential to answer this challenge. Owing to the few fundamental approximations made and the direct treatment of electron correlation, QMC methods are among the most accurate electronic structure methods available to date. We assess the accuracy of the diffusion Monte Carlo method in the case of rocksalt manganese oxide (MnO). We study the electronic properties of this strongly-correlated oxide, which has been identified as a suitable candidate for many applications ranging from catalysts to electronic devices. ``This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.'' Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides.
Renison, C Alicia; Fernandes, Kyle D; Naidoo, Kevin J
2015-07-01
This article describes an extension of the quantum supercharger library (QSL) to perform quantum mechanical (QM) gradient and optimization calculations as well as hybrid QM and molecular mechanical (QM/MM) molecular dynamics simulations. The integral derivatives are, after the two-electron integrals, the most computationally expensive part of the aforementioned calculations/simulations. Algorithms are presented for accelerating the one- and two-electron integral derivatives on a graphical processing unit (GPU). It is shown that a Hartree-Fock ab initio gradient calculation is up to 9.3X faster on a single GPU compared with a single central processing unit running an optimized serial version of GAMESS-UK, which uses the efficient Schlegel method for s- and l-orbitals. Benchmark QM and QM/MM molecular dynamics simulations are performed on cellobiose in vacuo and in a 39 Å water sphere (45 QM atoms and 24843 point charges, respectively) using the 6-31G basis set. The QSL can perform 9.7 ps/day of ab initio QM dynamics and 6.4 ps/day of QM/MM dynamics on a single GPU in full double precision. © 2015 Wiley Periodicals, Inc. PMID:25975864
NASA Astrophysics Data System (ADS)
Tokizaki, Chihiro; Yoshida, Takahiko; Takayanagi, Toshiyuki
2016-05-01
The cyclooctatetraene (COT) anion has a stable D4h structure that is similar to the transition state configurations of the neutral C-C bond-alternation (D4h ↔ D8h ↔ D4h) and ring-inversion (D2d ↔ D4h ↔ D2d) unimolecular reactions. The previously measured photodetachment spectrum of COT- revealed the reaction dynamics in the vicinity of the two transition states on the neutral potential energy surface. In this work, the photodetachment spectrum is calculated quantum mechanically on ab initio-level potential energy surfaces within a three degree-of-freedom reduced-dimensionality model. Very good agreement has been obtained between theory and experiment, providing reliable interpretations for the experimental spectrum. A detailed picture of the reactive molecular dynamics of the COT unimolecular reaction in the transition state region is also discussed.
NASA Astrophysics Data System (ADS)
Matsuda, Yuki
This dissertation focuses on ab-initio quantum mechanical calculations of nanoelectronics in three research topics: contact resistance properties of carbon nanotubes and graphenes (Chapters 1 through 3), electrical properties of carbon nanotubes (Chapter 4) and silicon nanowires (Chapter 5). Through all the chapters, the aim of the research is to provide useful guidelines for experimentalists. Chapter 1 presents the contact resistance of metal electrode-carbon nanotube and metal electrode-graphene interfaces for various deposited metals, based on first-principles quantum mechanical density functional and matrix Green's function methods. Chapters 2 and 3 describe inventive ways to enhance contact resistance properties as well as mechanical stabilities using "molecular anchors" (Chapter 2) or using "end-contacted" (or end-on) electrodes (Chapter 3). Chapters 1 through 3 also provide useful guidelines for nanotube assembly process which is one of the main obstacles in nanoelectronics. Chapter 4 shows accurate and detailed band structure properties of single-walled carbon nanotubes using B3LYP hybrid functional, which are critical parameters in determining the electronic properties such as small band gaps (˜0.1 eV) and effective masses. Chapter 5 details both structural and electronic properties of silicon nanowires. These results lead to the findings controlling the diameter and surface coverage by adsorbates (e.g., hydrogen) of silicon nanowires can be effectively used to optimize their properties for various applications. All the theoretical results are compared with other theoretical studies and experimental data. Notably, electronic studies using B3LYP show excellent agreement with experimental studies quantitatively, which previous quantum mechanical calculations had failed. These studies show how quantum mechanical predictions of complex phenomena can be effectively investigated computationally in nanomaterials and nanodevices. Given the difficulty, expense
NASA Astrophysics Data System (ADS)
Cuong, Nguyen Tien; Mizuta, Hiroshi; Cong, Bach Thanh; Otsuka, Nobuo; Chi, Dam Hieu
2012-09-01
Graphene is a promising candidate as a material used in nano-scale devices because of recent developments in advanced experimental techniques. Motivated by recent successful fabrications of U-shaped graphene channel transistors by using the gallium focused ion beam technology, we have performed ab-initio calculations to investigate the electronic properties and quantum transport in U-shaped graphene nanoribbons. The electronic properties are calculated using a numerical atomic orbital basis set in the framework of the density functional theory. The transport properties are investigated using the non-equilibrium Green's function method. The transmission spectra of U-shaped graphenes are analyzed in order to reveal the quantum transport of the systems. We found that the graphene nanoribbons tend to open a band gap when U-shaped structures are formed in both armchair and zigzag cases. The geometrical structures of U-shaped GNRs had enormous influences on the electron transport around the Fermi energy due to the formation of quasi-bound states at zigzag edges. The obtained results have provided valuable information for designing potential nano-scale devices based on graphenes.
NASA Astrophysics Data System (ADS)
Vitali, Ettore; Shi, Hao; Qin, Mingpu; Zhang, Shiwei
The possibility of calculating dynamical correlation functions from first principles provides a unique opportunity to explore the manifold of the excited states of a quantum many-body system. Such calculations allow us to predict interesting physical properties like spectral functions, excitation spectra and charge and spin gaps, which are more difficult to access from usual equilibrium calculations. We address the ab-initio calculation of dynamical Green functions and two-body correlation functions in the Auxiliary-field Quantum Monte Carlo method, using the two-dimensional Hubbard model as an example. When the sign problem is not present, an unbiased estimation of imaginary time correlation functions is obtained. We discuss in detail the complexity and the stability of the calculations. Moreover, we propose a new approach which is expected to be very useful when dealing with dilute systems, e.g. for cold gases, allowing calculations with a very favorable complexity in the system size. Supported by NSF, DOE SciDAC, and Simons Foundation.
NASA Astrophysics Data System (ADS)
Zeng, Xiancheng; Hu, Hao; Hu, Xiangqian; Cohen, Aron J.; Yang, Weitao
2008-03-01
Electron transfer (ET) reactions are one of the most important processes in chemistry and biology. Because of the quantum nature of the processes and the complicated roles of the solvent, theoretical study of ET processes is challenging. To simulate ET processes at the electronic level, we have developed an efficient density functional theory (DFT) quantum mechanical (QM)/molecular mechanical (MM) approach that uses the fractional number of electrons as the order parameter to calculate the redox free energy of ET reactions in solution. We applied this method to study the ET reactions of the aqueous metal complexes Fe(H2O)62+/3+ and Ru(H2O)62+/3+. The calculated oxidation potentials, 5.82 eV for Fe(II/III) and 5.14 eV for Ru(II/III), agree well with the experimental data, 5.50 and 4.96 eV, for iron and ruthenium, respectively. Furthermore, we have constructed the diabatic free energy surfaces from histogram analysis based on the molecular dynamics trajectories. The resulting reorganization energy and the diabatic activation energy also show good agreement with experimental data. Our calculations show that using the fractional number of electrons (FNE) as the order parameter in the thermodynamic integration process leads to efficient sampling and validate the ab initio QM/MM approach in the calculation of redox free energies.
Ab initio path-integral molecular dynamics and the quantum nature of hydrogen bonds
NASA Astrophysics Data System (ADS)
Yexin, Feng; Ji, Chen; Xin-Zheng, Li; Enge, Wang
2016-01-01
The hydrogen bond (HB) is an important type of intermolecular interaction, which is generally weak, ubiquitous, and essential to life on earth. The small mass of hydrogen means that many properties of HBs are quantum mechanical in nature. In recent years, because of the development of computer simulation methods and computational power, the influence of nuclear quantum effects (NQEs) on the structural and energetic properties of some hydrogen bonded systems has been intensively studied. Here, we present a review of these studies by focussing on the explanation of the principles underlying the simulation methods, i.e., the ab initio path-integral molecular dynamics. Its extension in combination with the thermodynamic integration method for the calculation of free energies will also be introduced. We use two examples to show how this influence of NQEs in realistic systems is simulated in practice. Project supported by the National Natural Science Foundation of China (Grant Nos. 11275008, 91021007, and 10974012) and the China Postdoctoral Science Foundation (Grant No. 2014M550005).
Marsalek, Ondrej; Markland, Thomas E
2016-02-01
Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost. PMID:26851913
Conte, Riccardo; Aspuru-Guzik, Alán; Ceotto, Michele
2013-10-17
A time-dependent semiclassical approach for vibrational spectra calculations is shown to describe deep tunneling splittings, resonances, and quantum frequencies in multidimensional multiwell systems, by propagating a very limited number of classical trajectories. The approach is tested on ammonia by evolving eight trajectories on a full-dimensional PES. Quantum effects are reproduced, and results are in good agreement with time-independent quantum calculations. All the features are maintained when ab initio "on-the-fly" dynamics is adopted, thus demonstrating that precomputation of the PES can be avoided. The approach overcomes the typical scaling issues of quantum mechanical techniques without introducing any simplifications nor reductions of dimensionality of the problem. The proposed methodology is promising for further applications to systems of major complexity. PMID:26705583
NASA Astrophysics Data System (ADS)
Philipp, Dean Michael
Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree-Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by ``breaking'' bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. The mixed modeling presented here can be used for single point energy calculations and geometry optimizations. Results from tests of the method to find relative conformational energies and geometries of alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations.
Communication: Towards ab initio self-energy embedding theory in quantum chemistry
NASA Astrophysics Data System (ADS)
Lan, Tran Nguyen; Kananenka, Alexei A.; Zgid, Dominika
2015-12-01
The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green's function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.
Baima, Jacopo; Zelferino, Alessandro; Olivero, Paolo; Erba, Alessandro; Dovesi, Roberto
2016-01-21
Quantum-mechanical ab initio calculations are performed to elucidate the vibrational spectroscopic features of a common irradiation-induced defect in diamond, i.e. the neutral vacancy. Raman spectra are computed analytically through a Coupled-Perturbed-Hartree-Fock/Kohn-Sham approach as a function of both different defect spin states and defect concentration. The experimental Raman features of defective diamond located in the 400-1300 cm(-1) spectral range, i.e. below the first-order line of pristine diamond at 1332 cm(-1), are well reproduced, thus corroborating the picture according to which, at low damage densities, this spectral region is mostly affected by non-graphitic sp(3) defects. No peaks above 1332 cm(-1) are found, thus ruling out previous tentative assignments of different spectral features (at 1450 and 1490 cm(-1)) to the neutral vacancy. The perturbation introduced by the vacancy to the thermal nuclear motion of carbon atoms in the defective lattice is discussed in terms of atomic anisotropic displacement parameters (ADPs), computed from converged lattice dynamics calculations. PMID:26686374
Communication: Towards ab initio self-energy embedding theory in quantum chemistry.
Lan, Tran Nguyen; Kananenka, Alexei A; Zgid, Dominika
2015-12-28
The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green's function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces. PMID:26723581
NASA Astrophysics Data System (ADS)
Erba, A.; Ferrabone, M.; Baima, J.; Orlando, R.; Rérat, M.; Dovesi, R.
2013-02-01
The vibration spectrum of single-walled zigzag boron nitride (BN) nanotubes is simulated with an ab initio periodic quantum chemical method. The trend towards the hexagonal monolayer (h-BN) in the limit of large tube radius R is explored for a variety of properties related to the vibrational spectrum: vibration frequencies, infrared intensities, oscillator strengths, and vibration contributions to the polarizability tensor. The (n,0) family is investigated in the range from n = 6 (24 atoms in the unit cell and tube radius R = 2.5 Å) to n = 60 (240 atoms in the cell and R = 24.0 Å). Simulations are performed using the CRYSTAL program which fully exploits the rich symmetry of this class of one-dimensional periodic systems: 4n symmetry operators for the general (n,0) tube. Three sets of infrared active phonon bands are found in the spectrum. The first one lies in the 0-600 cm-1 range and goes regularly to zero when R increases; the connection between these normal modes and the elastic and piezoelectric constants of h-BN is discussed. The second (600-800 cm-1) and third (1300-1600 cm-1) sets tend regularly, but with quite different speed, to the optical modes of the h-BN layer. The vibrational contribution of these modes to the two components (parallel and perpendicular) of the polarizability tensor is also discussed.
Communication: Towards ab initio self-energy embedding theory in quantum chemistry
Lan, Tran Nguyen; Kananenka, Alexei A.; Zgid, Dominika
2015-12-28
The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green’s function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.
Ab Initio Quantum Mechanical/Molecular Mechanical Studies of Histone Modifying Enzymes
NASA Astrophysics Data System (ADS)
Zhang, Yingkai
Histone proteins that form the nucleosome core are subject to a variety of post-translational transformations. These histone modifications make up the histone code which extends the information in the genetic code and is emerging as an essential mechanism to regulate gene expression. In spite of a current flurry of significant advances in experimental studies, there has been little theoretical understanding regarding how enzymes generate or remove these modifications. Very recently, we have made excellent progresses in investigating two such important histone-modifying enzyme families: zinc-dependent histone deacetylases (HDACs) and histone lysine methyltransferases (HKMTs). Our studies on a histonedeacetylase- like protein HDLP suggested a novel catalytic mechanism. The simulations on HKMT SET7/9 have characterized the histone lysine methylation reaction and elucidated the origin of enzyme catalysis. Our computational approaches centered on the pseudobond ab initio quantum mechanical/molecular mechanical (QM/MM) method, which allows for accurate modeling of the chemistry at the reaction active site while properly including the effects of the protein environment
NASA Astrophysics Data System (ADS)
Smith, Grant D.; Jaffe, Richard L.; Yoon, Do. Y.
1998-06-01
High-level ab initio quantum chemistry calculations are shown to predict conformer populations of 1,2-dimethoxypropane and 5-methoxy-1,3-dioxane that are consistent with gas-phase NMR vicinal coupling constant measurements. The conformational energies of the cyclic ether 5-methoxy-1,3-dioxane are found to be consistent with those predicted by a rotational isomeric state (RIS) model based upon the acyclic analog 1,2-dimethoxypropane. The quantum chemistry and RIS calculations indicate the presence of strong attractive 1,5 C(H 3)⋯O electrostatic interactions in these molecules, similar to those found in 1,2-dimethoxyethane.
Wiebeler, Christian; Schumacher, Stefan
2014-09-11
Photochromism is a light-induced molecular process that is likely to find its way into future optoelectronic devices. In further optimization of photochromic materials, light-induced conversion efficiencies as well as reaction times can usually only be determined once a new molecule was synthesized. Here we use nonadiabatic ab initio molecular dynamics to study the electrocyclic reaction of diarylethenes, comparing normal- and inverse-type systems. Our study highlights that reaction quantum yields can be successfully predicted in accord with experimental findings. In particular, we find that inverse-type diarylethenes show a significantly higher reaction quantum yield and cycloreversion on times typically as short as 100 fs. PMID:25140609
Wong, Kim F.; Sonnenberg, Jason L.; Paesani, Francesco; Yamamoto, Takeshi; Vaníček, Jiří; Zhang, Wei; Schlegel, H. Bernhard; Case, David A.; Cheatham, Thomas E.; Miller, William H.; Voth, Gregory A.
2010-01-01
The rates of intramolecular proton transfer are calculated on a full-dimensional reactive electronic potential energy surface that incorporates high level ab initio calculations along the reaction path and by using classical Transition State theory, Path-Integral Quantum Transition State Theory, and the Quantum Instanton approach. The specific example problem studied is malonaldehyde. Estimates of the kinetic isotope effect using the latter two methods are found to be in reasonable agreement with each other. Improvements and extensions of this practical, yet chemically accurate framework for the calculations of quantized, reactive dynamics are also discussed. PMID:21116485
The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations
Sellier, J.M. Dimov, I.
2014-09-15
The aim of ab-initio approaches is the simulation of many-body quantum systems from the first principles of quantum mechanics. These methods are traditionally based on the many-body Schrödinger equation which represents an incredible mathematical challenge. In this paper, we introduce the many-body Wigner Monte Carlo method in the context of distinguishable particles and in the absence of spin-dependent effects. Despite these restrictions, the method has several advantages. First of all, the Wigner formalism is intuitive, as it is based on the concept of a quasi-distribution function. Secondly, the Monte Carlo numerical approach allows scalability on parallel machines that is practically unachievable by means of other techniques based on finite difference or finite element methods. Finally, this method allows time-dependent ab-initio simulations of strongly correlated quantum systems. In order to validate our many-body Wigner Monte Carlo method, as a case study we simulate a relatively simple system consisting of two particles in several different situations. We first start from two non-interacting free Gaussian wave packets. We, then, proceed with the inclusion of an external potential barrier, and we conclude by simulating two entangled (i.e. correlated) particles. The results show how, in the case of negligible spin-dependent effects, the many-body Wigner Monte Carlo method provides an efficient and reliable tool to study the time-dependent evolution of quantum systems composed of distinguishable particles.
NASA Astrophysics Data System (ADS)
Ran, Shi-Ju
2016-05-01
In this work, a simple and fundamental numeric scheme dubbed as ab initio optimization principle (AOP) is proposed for the ground states of translational invariant strongly correlated quantum lattice models. The idea is to transform a nondeterministic-polynomial-hard ground-state simulation with infinite degrees of freedom into a single optimization problem of a local function with finite number of physical and ancillary degrees of freedom. This work contributes mainly in the following aspects: (1) AOP provides a simple and efficient scheme to simulate the ground state by solving a local optimization problem. Its solution contains two kinds of boundary states, one of which play the role of the entanglement bath that mimics the interactions between a supercell and the infinite environment, and the other gives the ground state in a tensor network (TN) form. (2) In the sense of TN, a novel decomposition named as tensor ring decomposition (TRD) is proposed to implement AOP. Instead of following the contraction-truncation scheme used by many existing TN-based algorithms, TRD solves the contraction of a uniform TN in an opposite way by encoding the contraction in a set of self-consistent equations that automatically reconstruct the whole TN, making the simulation simple and unified; (3) AOP inherits and develops the ideas of different well-established methods, including the density matrix renormalization group (DMRG), infinite time-evolving block decimation (iTEBD), network contractor dynamics, density matrix embedding theory, etc., providing a unified perspective that is previously missing in this fields. (4) AOP as well as TRD give novel implications to existing TN-based algorithms: A modified iTEBD is suggested and the two-dimensional (2D) AOP is argued to be an intrinsic 2D extension of DMRG that is based on infinite projected entangled pair state. This paper is focused on one-dimensional quantum models to present AOP. The benchmark is given on a transverse Ising
Revealing halogen bonding interactions with anomeric systems: an ab initio quantum chemical studies.
Lo, Rabindranath; Ganguly, Bishwajit
2015-02-01
A computational study has been performed using MP2 and CCSD(T) methods on a series of O⋯X (X=Br, Cl and I) halogen bonds to evaluate the strength and characteristic of such highly directional noncovalent interactions. The study has been carried out on a series of dimeric complexes formed between interhalogen compounds (such as BrF, BrCl and BrI) and oxygen containing electron donor molecule. The existence and consequences of the anomeric effect of the electron donor molecule has also been investigated through an exploration of halogen bonding interactions in this halogen bonded complexes. The ab initio quantum chemical calculations have been employed to study both the nature and directionality of the halogen molecules toward the sp(3) oxygen atom in anomeric systems. The presence of anomeric nO→σ*CN interaction involves a dominant role for the availability of the axial and equatorial lone pairs of donor O atom to participate with interhalogen compounds in the halogen-bonded complexes. The energy difference between the axial and equatorial conformers with interhalogen compounds reaches up to 4.60 kJ/mol, which however depends upon the interacting halogen atoms and its attaching atoms. The energy decomposition analysis further suggests that the total halogen bond interaction energies are mainly contributed by the attractive electrostatic and dispersion components. The role of substituents attached with the halogen atoms has also been evaluated in this study. With the increase of halogen atom size and the positive nature of σ-hole, the halogen atom interacted more with the electron donor atom and the electrostatic contribution to the total interaction energy enhances appreciably. Further, noncovalent interaction (NCI) studies have been carried out to locate the noncovalent halogen bonding interactions in real space. PMID:25522359
Kanno, Manabu; Ito, Yuta; Shimakura, Noriyuki; Koseki, Shiro; Kono, Hirohiko; Fujimura, Yuichi
2015-01-21
We theoretically investigated the mechanism of ultrafast nonradiative transition through conical intersections in photoexcited pyrazine by ab initio quantum dynamical calculations. This work was motivated by the recent theoretical and experimental studies that presented conflicting results: the former is the on-the-fly semiclassical surface hopping calculation combined with the time-dependent density functional theory, which showed that nonadiabatic transitions from the optically bright S2 ((1)B(2u), ππ*) state to the optically dark S3 ((1)A(u), nπ*) and S4 ((1)B(2g), nπ*) states take place predominantly at the initial stage of electronic relaxation [U. Werner et al., Chem. Phys., 2008, 349, 319]; the latter is the pump-probe photoelectron spectroscopic measurement, which reported the S2 lifetime (22 ± 3 fs) of nonradiative decay to the almost dark S1 ((1)B(3u), nπ*) state [Y.-I. Suzuki et al., J. Chem. Phys., 2010, 132, 174302]. We constructed adiabatic and diabatic potential energy surfaces of these ππ* and nπ* states using the multireference configuration interaction method and calculated their diabatic couplings within two-dimensional subspaces spanned by selected ground-state normal coordinates. Contrary to the surface hopping study, our nuclear wave packet simulations demonstrated that nonadiabatic transitions to the S3 and S4 states are so small that the conventional two-state (S1 and S2) picture is valid. Ultrafast internal conversion of pyrazine, which is deemed to proceed with a 22 fs lifetime, in fact consists of three consecutive steps: (i) the wave packet excited to the S2 state travels toward the S2-S1 conical intersection in 10 fs, (ii) the nonadiabatic transition to the S1 state progresses at a rapid rate corresponding to a transient lifetime of 7 fs, and (iii) intramolecular vibrational energy redistribution occurs in the S1 state in about 80 fs after optical excitation. To verify this prediction, time-resolved experiments with a
Moin, Syed Tarique; Hofer, Thomas S.; Weiss, Alexander K. H.; Rode, Bernd M.
2013-07-07
Ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) were successfully applied to Cu(II) embedded in water to elucidate structure and to understand dynamics of ligand exchange mechanism. From the simulation studies, it was found that using an extended large quantum mechanical region including two shells of hydration is required for a better description of the dynamics of exchanging water molecules. The structural features characterized by radial distribution function, angular distribution function and other analytical parameters were consistent with experimental data. The major outcome of this study was the dynamics of exchange mechanism and reactions in the first hydration shell that could not be studied so far. The dynamical data such as mean residence time of the first shell water molecules and other relevant data from the simulations are close to the results determined experimentally. Another major characteristic of hydrated Cu(II) is the Jahn-Teller distortion which was also successfully reproduced, leading to the final conclusion that the dominating aqua complex is a 6-coordinated species. The ab initio QMCF-MD formalism proved again its capabilities of unraveling even ambiguous properties of hydrated species that are far difficult to explore by any conventional quantum mechanics/molecular mechanics (QM/MM) approach or experiment.
Xie, Binbin; Liu, Lihong; Cui, Ganglong; Fang, Wei-Hai; Cao, Jun; Feng, Wei; Li, Xin-qi
2015-11-21
In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N2CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N2CO photodissociation at λ > 335 nm is an ultrafast process and the two C-N bonds are broken in a stepwise way, giving birth to CO and N2 as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C-N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes. PMID:26590527
Influence of pressure on the properties of GaN/AlN multi-quantum wells - Ab initio study
NASA Astrophysics Data System (ADS)
Strak, Pawel; Sakowski, Konrad; Kaminska, Agata; Krukowski, Stanislaw
2016-06-01
Pressure dependence of physical properties of GaN/AlN multi-quantum wells (MQWs) was investigated using ab intio calculations. The influence of pressure was divided into two main contributions: pressure affecting the properties of GaN and AlN bulk semiconductors and an influence on systems of polar quantum wells deposited on various substrates. An influence of hydrostatic, uniaxial, and tetragonal strain on the crystallographic structure, polarization (piezoelectricity), and the bandgap of the bulk systems is assessed using ab initio calculations. It was shown that when a partial relaxation of the structure is assumed, the tetragonal strain may explain an experimentally observed reduction of pressure coefficients for polar GaN/AlN MQWs. The MQWs were also simulated directly using density functional theory (DFT) calculations. A comparison of these two approaches confirmed that nonlinear effects induced by the tetragonal strain related to lattice mismatch between the substrates and the polar MQWs systems are responsible for a drastic decrease of the pressure coefficients of photoluminescence (PL) energy experimentally observed in polar GaN/AlGaN MQWs.
ab initio based tight-binding investigation of quantum spin Hall effect in InAs/GaSb quantum wells
NASA Astrophysics Data System (ADS)
Wu, Quansheng; Soluyanov, Alexey; Troyer, Matthias
Quantum spin Hall state is a topologically non-trivial quantum state, which can be used for designing various quantum devices including those potentially useful for quantum computing. Type-II InAs/GaSb semiconductor quantum well was predicted to realize this state of matter. In this work, we systematically investigate topological properties of this system using symmetrized Wannier-based tight-binding models. The model parameters are derived from first-principles hybrid functional calculations, which capture the right band gap and effective masses of both InAs and GaSb. By varying the thickness of InAs and GaSb layers, three possible phases are obtained: normal insulator, quantum spin Hall insulator, and semimetal, allowing us to optimize the growth conditions for the quantum spin Hall phase realization. Most importantly, we identify optimal growth directions, showing that a significant increase of the topological gap can be obtained by growing the quantum well in the [111]-direction. Phase diagrams are obtained for different layer thicknesses and growth directions. Effects of strain and applied electric fields are also discussed.
NASA Astrophysics Data System (ADS)
Dieterich, Johannes M.; Werner, Hans-Joachim; Mata, Ricardo A.; Metz, Sebastian; Thiel, Walter
2010-01-01
Energy and free energy barriers for acetaldehyde conversion in aldehyde oxidoreductase are determined for three reaction pathways using quantum mechanical/molecular mechanical (QM/MM) calculations on the solvated enzyme. Ab initio single-point QM/MM energies are obtained at the stationary points optimized at the DFT(B3LYP)/MM level. These ab initio calculations employ local correlation treatments [LMP2 and LCCSD(T0)] in combination with augmented triple- and quadruple-zeta basis sets, and the final coupled cluster results include MP2-based corrections for basis set incompleteness and for the domain approximation. Free energy perturbation (FEP) theory is used to generate free energy profiles at the DFT(B3LYP)/MM level for the most important reaction steps by sampling along the corresponding reaction paths using molecular dynamics. The ab initio and FEP QM/MM results are combined to derive improved estimates of the free energy barriers, which differ from the corresponding DFT(B3LYP)/MM energy barriers by about 3 kcal mol-1. The present results confirm the qualitative mechanistic conclusions from a previous DFT(B3LYP)/MM study. Most favorable is a three-step Lewis base catalyzed mechanism with an initial proton transfer from the cofactor to the Glu869 residue, a subsequent nucleophilic attack that yields a tetrahedral intermediate (IM2), and a final rate-limiting hydride transfer. The competing metal center activated pathway has the same final step but needs to overcome a higher barrier in the initial step on the route to IM2. The concerted mechanism has the highest free energy barrier and can be ruled out. While confirming the qualitative mechanistic scenario proposed previously on the basis of DFT(B3LYP)/MM energy profiles, the present ab initio and FEP QM/MM calculations provide corrections to the barriers that are important when aiming at high accuracy.
Walker, Brent; Michaelides, Angelos
2010-11-01
The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation. PMID:21054031
DFT and ab initio quantum chemical studies on p-cyanobenzoic acid
NASA Astrophysics Data System (ADS)
Arjunan, V.; Rani, T.; Varalakshmy, L.; Mohan, S.; Tedlamelekot, F.
2011-05-01
The Fourier transform infrared (FTIR) and FT-Raman spectra of p-cyanobenzoic acid (CBA) have been recorded in the range 4000-400 and 4000-100 cm -1, respectively. The complete vibrational assignment and analysis of the fundamental modes of the compound were carried out using the observed FTIR and FT-Raman data. The vibrational frequencies determined experimentally were compared with theoretical wavenumbers obtained from ab initio HF and DFT-B3LYP gradient calculations employing 6-31G**, 6-311++G** and cc-pVTZ basis sets for the optimised geometry of the compound. The geometry and normal modes of vibration obtained from the HF and DFT methods are in good agreement with the experimental data. The normal coordinate analysis was also carried out with ab initio force fields utilising Wilson's FG matrix method. The interactions of cyano and carboxylic acid groups with the skeletal vibrational modes were investigated.
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-21
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H{sub 2}O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0–4000 cm{sup −1} is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
NASA Astrophysics Data System (ADS)
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-01
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H2O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0-4000 cm-1 is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
Ab Initio Studies of Shock-Induced Chemical Reactions of Inter-Metallics
NASA Astrophysics Data System (ADS)
Zaharieva, Roussislava; Hanagud, Sathya
2009-06-01
Shock-induced and shock assisted chemical reactions of intermetallic mixtures are studied by many researchers, using both experimental and theoretical techniques. The theoretical studies are primarily at continuum scales. The model frameworks include mixture theories and meso-scale models of grains of porous mixtures. The reaction models vary from equilibrium thermodynamic model to several non-equilibrium thermodynamic models. The shock-effects are primarily studied using appropriate conservation equations and numerical techniques to integrate the equations. All these models require material constants from experiments and estimates of transition states. Thus, the objective of this paper is to present studies based on ab initio techniques. The ab inito studies, to date, use ab inito molecular dynamics. This paper presents a study that uses shock pressures, and associated temperatures as starting variables. Then intermetallic mixtures are modeled as slabs. The required shock stresses are created by straining the lattice. Then, ab initio binding energy calculations are used to examine the stability of the reactions. Binding energies are obtained for different strain components super imposed on uniform compression and finite temperatures. Then, vibrational frequencies and nudge elastic band techniques are used to study reactivity and transition states. Examples include Ni and Al.
Fano effect in an AB interferometer with a quantum dot side-coupled to a single Majorana bound state
NASA Astrophysics Data System (ADS)
Zeng, Qi-Bo; Chen, Shu; Lü, Rong
2016-02-01
We study the conductance and interference effects through an AB interferometer with an embedded quantum dot (QD) side-coupled to a single Majorana bound state (MBS) by using non-equilibrium Green's function method. The energy levels appearing in the QD are calculated by diagonalizing the Hamiltonian of the embedded QD-MBS system. When the single QD energy level ɛ0 is set to 0, there are three discrete energy levels in the QD appearing at around ω = 0, ±√{ ɛM2 + 2λ2 } due to the coupling with MBS where ɛM is the coupling strength between the two MBSs at the two ends of the nanowire and λ is the coupling strength between the MBS and the QD. Asymmetric Fano lineshapes are found around these levels in the conductance due to the interference between electrons traversing through different paths. The phase shift of electrons through the QD changes from π / 2 to - π / 2 at each of these three energy values. However, the phase does not vary smoothly between these three energy levels but shows severe changes from - π / 2 to π / 2 at ω = ±√{ ɛM2 +λ2 }. As a comparison, we also study the similar AB interferometer in which the QD-MBS system is replaced by a normal QD-QD system or a simple single QD system, which shows only two or one Fano peak and the phase shifts from π / 2 to - π / 2 only at the Fano peaks. These differences reflect the distinct influences of Majorana bound state on the transport properties of AB interferometer.
Giese, Timothy J; York, Darrin M
2016-06-14
A new approach for performing Particle Mesh Ewald in ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with extended atomic orbital basis sets is presented. The new approach, the Ambient-Potential Composite Ewald (CEw) method, does not perform the QM/MM interaction with Mulliken charges nor electrostatically fit charges. Instead the nuclei and electron density interact directly with the MM environment, but in a manner that avoids the use of dense Fourier transform grids. By performing the electrostatics with the underlying QM density, the CEw method avoids self-consistent field instabilities that have been encountered with simple charge mapping procedures. Potential of mean force (PMF) profiles of the p-nitrophenyl phosphate dissociation reaction in explicit solvent are computed from PBE0/6-31G* QM/MM molecular dynamics simulations with various electrostatic protocols. The CEw profiles are shown to be stable with respect to real-space Ewald cutoff, whereas the PMFs computed from truncated and switched electrostatics produce artifacts. PBE0/6-311G**, AM1/d-PhoT, and DFTB2 QM/MM simulations are performed to generate two-dimensional PMF profiles of the phosphoryl transesterification reactions with ethoxide and phenoxide leaving groups. The semiempirical models incorrectly produce a concerted ethoxide mechanism, whereas PBE0 correctly produces a stepwise mechanism. The ab initio reaction barriers agree more closely to experiment than the semiempirical models. The failure of Mulliken-charge QM/MM-Ewald is analyzed. PMID:27171914
Ab initio quantum chemistry in parallel-portable tools and applications
Harrison, R.J.; Shepard, R. ); Kendall, R.A. )
1991-01-01
In common with many of the computational sciences, ab initio chemistry faces computational constraints to which a partial solution is offered by the prospect of highly parallel computers. Ab initio codes are large and complex (O(10{sup 5}) lines of FORTRAN), representing a significant investment of communal effort. The often conflicting requirements of portability and efficiency have been successfully resolved on vector computers by reliance on matrix oriented kernels. This proves inadequate even upon closely-coupled shared-memory parallel machines. We examine the algorithms employed during a typical sequence of calculations. Then we investigate how efficient portable parallel implementations may be derived, including the complex multi-reference singles and doubles configuration interaction algorithm. A portable toolkit, modeled after the Intel iPSC and the ANL-ACRF PARMACS, is developed, using shared memory and TCP/IP sockets. The toolkit is used as an initial platform for programs portable between LANS, Crays and true distributed-memory MIMD machines. Timings are presented. 53 refs., 4 tabs.
Luo, Ye Sorella, Sandro; Zen, Andrea
2014-11-21
We present a systematic study of a recently developed ab initio simulation scheme based on molecular dynamics and quantum Monte Carlo. In this approach, a damped Langevin molecular dynamics is employed by using a statistical evaluation of the forces acting on each atom by means of quantum Monte Carlo. This allows the use of an highly correlated wave function parametrized by several variational parameters and describing quite accurately the Born-Oppenheimer energy surface, as long as these parameters are determined at the minimum energy condition. However, in a statistical method both the minimization method and the evaluation of the atomic forces are affected by the statistical noise. In this work, we study systematically the accuracy and reliability of this scheme by targeting the vibrational frequencies of simple molecules such as the water monomer, hydrogen sulfide, sulfur dioxide, ammonia, and phosphine. We show that all sources of systematic errors can be controlled and reliable frequencies can be obtained with a reasonable computational effort. This work provides convincing evidence that this molecular dynamics scheme can be safely applied also to realistic systems containing several atoms.
Cvitas, Marko T; Althorpe, Stuart C
2009-04-23
We describe a quantum wave packet method for computing the state-to-state quantum dynamics of 4-atom AB + CD --> ABC + D reactions. The approach is an extension to 4-atom reactions of a version of the reactant-product decoupling (RPD) approach, applied previously to 3-atom reactions ( J. Chem. Phys. 2001, 114 , 1601 ). The approach partitions the coordinate space of the reaction into separate reagent, strong-interaction, and product regions, using a system of artificial absorbing and reflecting potentials. It employs a partitioned version of the split-operator propagator, which is more efficient than partitioning the (exact) time-dependent Schrodinger equation. The wave packet bounces off a reflecting potential in the entrance channel, which generates a source term; this is transformed efficiently from reagent to product Jacobi coordinates by exploiting some simple angular momentum properties. The efficiency and accuracy of the method is demonstrated by numerical tests on the benchmark OH + H(2) --> H(2)O + H reaction. PMID:19298045
Iodine-polyphenylacetylene charge-transfer complex: an ab initio quantum-chemical assessment
NASA Astrophysics Data System (ADS)
Andreocci, M. V.; Bossa, M.; Furlani, A.; Polzonetti, G.; Russo, M. V.
1991-07-01
The ab initio MO-LCAO-HF method has been used to calculate the electronic structure of the iodine-polyphenylacetylene charge-transfer complex (PPAI 2). Two models have been considered for the PPA molecule: a simple one containing two phenyl groups and a more realistic one containing six phenyl groups. The calculations give automatically the charge separation between I 5 and the polymer, and show that the total charge separation can be less than 1 e at short distances. The computed charges at the energy minimum have been succesfully introduced into the curve fitting of the I-3d 5/2 core level spectrum of PPAI 2 films, giving good agreement between experimental and theoretical results.
Xin, Xukai; Li, Bo; Jung, Jaehan; Yoon, Young Jun; Biswas, Rana; Lin, Zhiqun
2014-07-24
Quantum dot-sensitized solar cells (QDSSCs) have emerged as a promising solar architecture for next-generation solar cells. The QDSSCs exhibit a remarkably fast electron transfer from the quantum dot (QD) donor to the TiO2 acceptor with size quantization properties of QDs that allows for the modulation of band energies to control photoresponse and photoconversion efficiency of solar cells. In order to understand the mechanisms that underpin this rapid charge transfer, the electronic properties of CdSe and PbSe QDs with different sizes on the TiO2 substrate are simulated using a rigorous ab initio density functional method. Our method capitalizes on localized orbitalmore » basis set, which is computationally less intensive. Quite intriguingly, a remarkable set of electron bridging states between QDs and TiO2 occurring via the strong bonding between the conduction bands of QDs and TiO2 is revealed. Such bridging states account for the fast adiabatic charge transfer from the QD donor to the TiO2 acceptor, and may be a general feature for strongly coupled donor/acceptor systems. All the QDs/TiO2 systems exhibit type II band alignments, with conduction band offsets that increase with the decrease in QD size. This facilitates the charge transfer from QDs donors to TiO2 acceptors and explains the dependence of the increased charge transfer rate with the decreased QD size.« less
NASA Astrophysics Data System (ADS)
Hagiwara, Yohsuke; Tateno, Masaru
2010-10-01
We review the recent research on the functional mechanisms of biological macromolecules using theoretical methodologies coupled to ab initio quantum mechanical (QM) treatments of reaction centers in proteins and nucleic acids. Since in most cases such biological molecules are large, the computational costs of performing ab initio calculations for the entire structures are prohibitive. Instead, simulations that are jointed with molecular mechanics (MM) calculations are crucial to evaluate the long-range electrostatic interactions, which significantly affect the electronic structures of biological macromolecules. Thus, we focus our attention on the methodologies/schemes and applications of jointed QM/MM calculations, and discuss the critical issues to be elucidated in biological macromolecular systems.
Grid-based methods for biochemical ab initio quantum chemical applications
Colvin, M.E.; Nelson, J.S.; Mori, E.
1997-01-01
A initio quantum chemical methods are seeing increased application in a large variety of real-world problems including biomedical applications ranging from drug design to the understanding of environmental mutagens. The vast majority of these quantum chemical methods are {open_quotes}spectral{close_quotes}, that is they describe the charge distribution around the nuclear framework in terms of a fixed analytic basis set. Despite the additional complexity they bring, methods involving grid representations of the electron or solvent charge can provide more efficient schemes for evaluating spectral operators, inexpensive methods for calculating electron correlation, and methods for treating the electrostatic energy of salvation in polar solvents. The advantage of mixed or {open_quotes}pseudospectral{close_quotes} methods is that they allow individual non-linear operators in the partial differential equations, such as coulomb operators, to be calculated in the most appropriate regime. Moreover, these molecular grids can be used to integrate empirical functionals of the electron density. These so-called density functional methods (DFT) are an extremely promising alternative to conventional post-Hartree Fock quantum chemical methods. The introduction of a grid at the molecular solvent-accessible surface allows a very sophisticated treatment of a polarizable continuum solvent model (PCM). Where most PCM approaches use a truncated expansion of the solute`s electric multipole expansion, e.g. net charge (Born model) or dipole moment (Onsager model), such a grid-based boundary-element method (BEM) yields a nearly exact treatment of the solute`s electric field. This report describes the use of both DFT and BEM methods in several biomedical chemical applications.
Gurka, D F; Titus, R; Robins, K; Wong, A; Wurrey, C J; Durig, J R; Shen, Z; Burkhard, L P
1996-12-01
Computational chemistry, in conjunction with gas chromatography/mass spectrometry/Fourier transform infrared spectrometry (GC/MS/FT-IR), was used to tentatively identify seven tetrachlorobutadiene (TCBD) isomers detected in an environmental sample. Computation of the TCBD infrared spectra was performed with the Gaussian quantum chemistry software. The Hartree-Fock/6-31G* level of theory was employed, with IR frequencies scaled by a standard factor of 0.89. This approach shows great promise as a means of characterizing or confirming environmental analyte identifications when standard spectra, or pure standards required to measure standard spectra, are unavailable. PMID:21619333
Ivanov, Sergei D; Grant, Ian M; Marx, Dominik
2015-09-28
With the goal of computing quantum free energy landscapes of reactive (bio)chemical systems in multi-dimensional space, we combine the metadynamics technique for sampling potential energy surfaces with the ab initio path integral approach to treating nuclear quantum motion. This unified method is applied to the double proton transfer process in the formic acid dimer (FAD), in order to study the nuclear quantum effects at finite temperatures without imposing a one-dimensional reaction coordinate or reducing the dimensionality. Importantly, the ab initio path integral metadynamics technique allows one to treat the hydrogen bonds and concomitant proton transfers in FAD strictly independently and thus provides direct access to the much discussed issue of whether the double proton transfer proceeds via a stepwise or concerted mechanism. The quantum free energy landscape we compute for this H-bonded molecular complex reveals that the two protons move in a concerted fashion from initial to product state, yet world-line analysis of the quantum correlations demonstrates that the protons are as quantum-uncorrelated at the transition state as they are when close to the equilibrium structure. PMID:26429008
Ivanov, Sergei D. Grant, Ian M.; Marx, Dominik
2015-09-28
With the goal of computing quantum free energy landscapes of reactive (bio)chemical systems in multi-dimensional space, we combine the metadynamics technique for sampling potential energy surfaces with the ab initio path integral approach to treating nuclear quantum motion. This unified method is applied to the double proton transfer process in the formic acid dimer (FAD), in order to study the nuclear quantum effects at finite temperatures without imposing a one-dimensional reaction coordinate or reducing the dimensionality. Importantly, the ab initio path integral metadynamics technique allows one to treat the hydrogen bonds and concomitant proton transfers in FAD strictly independently and thus provides direct access to the much discussed issue of whether the double proton transfer proceeds via a stepwise or concerted mechanism. The quantum free energy landscape we compute for this H-bonded molecular complex reveals that the two protons move in a concerted fashion from initial to product state, yet world-line analysis of the quantum correlations demonstrates that the protons are as quantum-uncorrelated at the transition state as they are when close to the equilibrium structure.
NASA Astrophysics Data System (ADS)
Ivanov, Sergei D.; Grant, Ian M.; Marx, Dominik
2015-09-01
With the goal of computing quantum free energy landscapes of reactive (bio)chemical systems in multi-dimensional space, we combine the metadynamics technique for sampling potential energy surfaces with the ab initio path integral approach to treating nuclear quantum motion. This unified method is applied to the double proton transfer process in the formic acid dimer (FAD), in order to study the nuclear quantum effects at finite temperatures without imposing a one-dimensional reaction coordinate or reducing the dimensionality. Importantly, the ab initio path integral metadynamics technique allows one to treat the hydrogen bonds and concomitant proton transfers in FAD strictly independently and thus provides direct access to the much discussed issue of whether the double proton transfer proceeds via a stepwise or concerted mechanism. The quantum free energy landscape we compute for this H-bonded molecular complex reveals that the two protons move in a concerted fashion from initial to product state, yet world-line analysis of the quantum correlations demonstrates that the protons are as quantum-uncorrelated at the transition state as they are when close to the equilibrium structure.
NASA Astrophysics Data System (ADS)
Hara Chakravarty, Krishna; Thomsan, Kaj
2013-04-01
Through both field and laboratory studies it has been well established that wett-ability and flow rate in carbonate and sandstone reservoirs during water flooding is dictated by the sub-surface chemistry connected to stability of the water film amidst the oil phase and the rock surface, but there is no mathematical model to explain the flow rate and wett-ability mechanism on insertion of a low saline water in the rock system and its interaction with the hydrocarbons. Here, a model for calculating the feasible surface reactions has been developed based on free energy profiles of the different combinations of water, metal (Calcium, strontium and Magnesium) ions and hydrocarbons over a temperature regime of 273 K to 373 K .The approach uses Density Functional Theory based Local Density Approximation on molecular dynamic simulation using standard plane wave pseudopotential scheme. Herein we observe a two shell model during the interaction of water molecule with free energy of -2.8 and -2.5 kJ/mol at 2.4 and 4.2 Å heights from the calcite surface. Mg2+-water complex also retains a two layer model, while Ca2+ and Sr2+ metal ion-water complex show one major layer model of the free energy profile. The metal ion-water cluster adsorbs as an inner sphere complex directly above a surface carbonate group at a height of 3.9, 4.1, and 4.2 Å for calcium, strontium and magnesium respectively at a temperature of 273 K ,which is almost at the same height of the outer shell(4.2 Å) of bulk water profile and thus provides a local energy well for interaction among the complexes. Moreover at this height Mg2+ shows significantly more stability than water or other metal ion complex and it retains a similar pattern of energy profile at temperatures upto 373 K while outer shells loose their structure at higher temperatures in other metal ion complexes. A series of calculation were conducted to achieve the most stable polar hydrocarbon-metal ion-water complex for each metal individually .A free energy profile was also developed for polar hydrocarbon-Metal ion-water complex (CH3COO(-)-M2+-[H2O]n).Herein the two layer model was observed with Ca2+ and Sr2+ while the inner shell was still inconsistent with the 4.2 ( i.e. 3.9 to 4.4 ) Å local energy well model, the outer shell was observed near 6 Å. While with Mg2+ one major shell free energy model was obtained with several minor variations with changing temperature. Here the Ca2+ complex was significantly more stable being the only complex with a negative energy profile around 4.2 Å. Further these observations were together investigated to propose a wett-ability and flow rate model at on the surface described as: Mg2+ + Calcite +(5-6)H2O - -> Mg2+[H2O]5-6- Calcite + eV Calcite- [H2O]4 -Ca2+-CH3COO(-) + eV + SO42-+ (Ca2+,Sr2+) -> Calcite- [H2O] + Calcite- Ca2+- SO42-+ (Ca2+,Sr2+)-(-)OOCCH3 A large part of the work also includes the categorical rejection of many other possible structures, complexes and clusters formations. The model abides to all the experimentally observed points both in field and lab in wett ability alteration, and enhanced oil recovery, and through these energy variation the model categorically explains the effect of temperature on the flow behavior. This model provides a new view to the activity of Mg2+ in surface chemistry for carbonate rock structures in general and can have application in study of erosion, corrosion and other fields as well. This is the first comprehensive mathematical model to extensively explain the low salinity behaviors observed in enhanced oil recovery from a microscopic level to the best of our knowledge.
Nekrashevich, S. S. Gritsenko, V. A.; Klauser, R.; Gwo, S.
2010-10-15
Charge transfer {Delta}Q = 0.35e at the Si-N bond in silicon nitride is determined experimentally using photoelectron spectroscopy, and the ionic formula of silicon nitride Si{sub 3}{sup +1.4}N{sub 4}{sup -1.05} is derived. The electronic structure of {alpha}-Si{sub 3}N{sub 4} is studied ab initio using the density functional method. The results of calculations (partial density of states) are compared with experimental data on X-ray emission spectroscopy of amorphous Si{sub 3}N{sub 4}. The electronic structure of the valence band of amorphous Si{sub 3}N{sub 4} is studied using synchrotron radiation at different excitation energies. The electron and hole effective masses m{sub e}{sup *} {approx} m{sub h}{sup *} {approx} 0.5m{sub e} are estimated theoretically. The calculated values correspond to experimental results on injection of electrons and holes into silicon nitride.
Miller, J.; Miaskiewicz, K.; Osman, R.
1993-12-01
Studies of ring-saturated pyrimidine base lesions are used to illustrate an integrated modeling approach that combines quantum-chemical calculations with molecular dynamics simulation. Electronic-structure calculations on the lesions in Isolation reveal strong conformational preferences due to interactions between equatorial substituents to the pyrimidine ring. Large distortions of DNA should result when these interactions force the methyl group of thymine to assume an axial orientation, as is the case for thymine glycol but not for dihydrothymine. Molecular dynamics simulations of the dodecamer d(CGCGAATTCGCG){sub 2} with and without a ring-saturated thymine lesion at position T7 support this conclusion. Implications of these studies for recognition of thymine lesions by endonuclease III are also discussed.
Ab-Initio Hamiltonian Approach to Light Nuclei And to Quantum Field Theory
Vary, J.P.; Honkanen, H.; Li, Jun; Maris, P.; Shirokov, A.M.; Brodsky, S.J.; Harindranath, A.; de Teramond, G.F.; Ng, E.G.; Yang, C.; Sosonkina, M.; /Ames Lab
2012-06-22
Nuclear structure physics is on the threshold of confronting several long-standing problems such as the origin of shell structure from basic nucleon-nucleon and three-nucleon interactions. At the same time those interactions are being developed with increasing contact to QCD, the underlying theory of the strong interactions, using effective field theory. The motivation is clear - QCD offers the promise of great predictive power spanning phenomena on multiple scales from quarks and gluons to nuclear structure. However, new tools that involve non-perturbative methods are required to build bridges from one scale to the next. We present an overview of recent theoretical and computational progress with a Hamiltonian approach to build these bridges and provide illustrative results for the nuclear structure of light nuclei and quantum field theory.
NASA Astrophysics Data System (ADS)
Liu, Haiyan; Lu, Zhenyu; Cisneros, G. Andrés; Yang, Weitao
2004-07-01
The determination of reaction paths for enzyme systems remains a great challenge for current computational methods. In this paper we present an efficient method for the determination of minimum energy reaction paths with the ab initio quantum mechanical/molecular mechanical approach. Our method is based on an adaptation of the path optimization procedure by Ayala and Schlegel for small molecules in gas phase, the iterative quantum mechanical/molecular mechanical (QM/MM) optimization method developed earlier in our laboratory and the introduction of a new metric defining the distance between different structures in the configuration space. In this method we represent the reaction path by a discrete set of structures. For each structure we partition the atoms into a core set that usually includes the QM subsystem and an environment set that usually includes the MM subsystem. These two sets are optimized iteratively: the core set is optimized to approximate the reaction path while the environment set is optimized to the corresponding energy minimum. In the optimization of the core set of atoms for the reaction path, we introduce a new metric to define the distances between the points on the reaction path, which excludes the soft degrees of freedom from the environment set and includes extra weights on coordinates describing chemical changes. Because the reaction path is represented by discrete structures and the optimization for each can be performed individually with very limited coupling, our method can be executed in a natural and efficient parallelization, with each processor handling one of the structures. We demonstrate the applicability and efficiency of our method by testing it on two systems previously studied by our group, triosephosphate isomerase and 4-oxalocrotonate tautomerase. In both cases the minimum energy paths for both enzymes agree with the previously reported paths.
Ab initio quantum Monte Carlo simulations of the uniform electron gas without fixed nodes
NASA Astrophysics Data System (ADS)
Groth, S.; Schoof, T.; Dornheim, T.; Bonitz, M.
2016-02-01
The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in the warm dense matter regime, e.g., dense plasmas and laser excited solids. Also, the quality of density functional theory calculations crucially relies on the availability of accurate data for the exchange-correlation energy. Recently, results for N =33 spin-polarized electrons at high density, rs=r ¯/aB≲4 , and low temperature have been obtained with the configuration path integral Monte Carlo (CPIMC) method [T. Schoof et al., Phys. Rev. Lett. 115, 130402 (2015), 10.1103/PhysRevLett.115.130402]. To achieve these results, the original CPIMC algorithm [T. Schoof et al., Contrib. Plasma Phys. 51, 687 (2011), 10.1002/ctpp.201100012] had to be further optimized to cope with the fermion sign problem (FSP). It is the purpose of this paper to give detailed information on the manifestation of the FSP in CPIMC simulations of the UEG and to demonstrate how it can be turned into a controllable convergence problem. In addition, we present new thermodynamic results for higher temperatures. Finally, to overcome the limitations of CPIMC towards strong coupling, we invoke an independent method—the recently developed permutation blocking path integral Monte Carlo approach [T. Dornheim et al., J. Chem. Phys. 143, 204101 (2015), 10.1063/1.4936145]. The combination of both approaches is able to yield ab initio data for the UEG over the entire density range, above a temperature of about one half of the Fermi temperature. Comparison with restricted path integral Monte Carlo data [E. W. Brown et al., Phys. Rev. Lett. 110, 146405 (2013), 10.1103/PhysRevLett.110.146405] allows us to quantify the systematic error arising from the free particle nodes.
Řezáč, Jan; Hobza, Pavel
2014-08-12
Hydrogen fluoride dimer is a perfect model system for studying hydrogen bonding. Its size makes it possible to apply the most advanced theoretical methods available, yet it is a full-featured complex of molecules with nontrivial electronic structure and dynamic properties. Moreover, the dissociation energy of the HF dimer has been measured experimentally with an unparalleled accuracy of ±1 cm(-1)(Bohac et al. J. Chem. Phys. 1992, 9, 6681). In this work, we attempt to reproduce it by purely ab initio means, using advanced quantum-mechanical computational methods free of any empiricism. The purpose of this study is to demonstrate the capabilities of today's computational chemistry and to point out its limitations by identifying the contributions that introduce the largest uncertainty into the result. The dissociation energy is calculated using a composite scheme including large basis set CCSD(T) calculations, contributions of higher excitations up to CCSDTQ, relativistic and diagonal Born-Oppenheimer corrections and anharmonic vibrational calculations. The error of the calculated dissociation energy is 0.07 kcal/mol (25 cm(-1), 2.5%) when compared to the experiment. The major part of this error can be attributed to the inaccuracy of the calculations of the zero-point vibrational energy. PMID:26588277
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined. PMID:25106565
NASA Astrophysics Data System (ADS)
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-07
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
NASA Astrophysics Data System (ADS)
Brinck, Nils Tore
The objective of this study has been to investigate the use of computed molecular properties in predicting and interpreting intermolecular interactions. The molecular properties have been calculated rigorously from ab initio wave functions. We have found the electrostatic potential to be a good tool for the analysis of nonbonding intermolecular interactions. It is demonstrated that the calculated electrostatic potentials around carbon-halogen bonds can be used to explain the directional preferences of halogen interactions in crystals. We also show that the orientation of the molecules in weak gas phase complexes between dihalogens and Lewis bases can be rationalized from their electrostatic potentials. However an analysis of the bonding in boron trifluoride and boron trichloride and their complexes with ammonia indicates that the relative stabilities of these complexes are dictated by charge transfer rather than electrostatics. The higher binding affinity for boron trichloride compared to boron trifluoride is explained by the higher charge capacity of the former. This is contrary to the commonly accepted explanation, which is based on backbonding. A local charge separation index has been defined from surface electrostatic potentials. This provides a measure of local polarity, even for molecules with zero dipole moments. Average local ionization energies computed on molecular surfaces permit predictions of relative reactivities of various sites toward electrophilic attack. There is a very good relationship between {rm p}K_{a} and minima in |{I}({bf r}), designated as |{I}_{S,min }, for a series of azines and azoles. Excellent correlations have also been found between { rm p}K_{a} values for a variety of carbon, oxygen and nitrogen acids and the |{I}_{S,min} of their conjugate bases. A study of the electrostatic potentials and average local ionization energies of the V-VII hydrides of the first three rows of the periodic table and their anions demonstrates the
Heats of Segregation of BCC Metals Using Ab Initio and Quantum Approximate Methods
NASA Technical Reports Server (NTRS)
Good, Brian; Chaka, Anne; Bozzolo, Guillermo
2003-01-01
Many multicomponent alloys exhibit surface segregation, in which the composition at or near a surface may be substantially different from that of the bulk. A number of phenomenological explanations for this tendency have been suggested, involving, among other things, differences among the components' surface energies, molar volumes, and heats of solution. From a theoretical standpoint, the complexity of the problem has precluded a simple, unified explanation, thus preventing the development of computational tools that would enable the identification of the driving mechanisms for segregation. In that context, we investigate the problem of surface segregation in a variety of bcc metal alloys by computing dilute-limit heats of segregation using both the quantum-approximate energy method of Bozzolo, Ferrante and Smith (BFS), and all-electron density functional theory. In addition, the composition dependence of the heats of segregation is investigated using a BFS-based Monte Carlo procedure, and, for selected cases of interest, density functional calculations. Results are discussed in the context of a simple picture that describes segregation behavior as the result of a competition between size mismatch and alloying effects
Xin, Xukai; Li, Bo; Jung, Jaehan; Yoon, Young Jun; Biswas, Rana; Lin, Zhiqun
2014-07-24
Quantum dot-sensitized solar cells (QDSSCs) have emerged as a promising solar architecture for next-generation solar cells. The QDSSCs exhibit a remarkably fast electron transfer from the quantum dot (QD) donor to the TiO_{2} acceptor with size quantization properties of QDs that allows for the modulation of band energies to control photoresponse and photoconversion efficiency of solar cells. In order to understand the mechanisms that underpin this rapid charge transfer, the electronic properties of CdSe and PbSe QDs with different sizes on the TiO_{2} substrate are simulated using a rigorous ab initio density functional method. Our method capitalizes on localized orbital basis set, which is computationally less intensive. Quite intriguingly, a remarkable set of electron bridging states between QDs and TiO_{2} occurring via the strong bonding between the conduction bands of QDs and TiO_{2} is revealed. Such bridging states account for the fast adiabatic charge transfer from the QD donor to the TiO_{2} acceptor, and may be a general feature for strongly coupled donor/acceptor systems. All the QDs/TiO_{2} systems exhibit type II band alignments, with conduction band offsets that increase with the decrease in QD size. This facilitates the charge transfer from QDs donors to TiO_{2} acceptors and explains the dependence of the increased charge transfer rate with the decreased QD size.
Mills, Jeffrey D; Ben-Nun, Michal; Rollin, Kyle; Bromley, Michael W J; Li, Jiabo; Hinde, Robert J; Winstead, Carl L; Sheehy, Jeffrey A; Boatz, Jerry A; Langhoff, Peter W
2016-08-25
Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions. The general theory in its various forms incorporates the early semiempirical atoms- and diatomics-in-molecules approaches of Moffitt, Ellison, Tully, Kuntz, and others in a comprehensive mathematical setting, and generalizes the developments of Eisenschitz, London, Claverie, and others addressing electron permutation symmetry adaptation issues, completing these early attempts to treat van der Waals and chemical forces on a common basis. Exact expressions are obtained for molecular Hamiltonian matrices and for associated energy eigenvalues as sums of separate atomic and interaction-energy terms, similar in this respect to the forms of classical force fields. The latter representation is seen to also provide a long-missing general definition of the energies of individual atoms and of their interactions within molecules and matter free from subjective additional constraints. A computer code suite is described for calculations of the many-electron atomic eigenspectra and
Ab initio quantum chemical investigation of arsenic sulfide molecular diversity from As4S6 and As4
NASA Astrophysics Data System (ADS)
Kyono, Atsushi
2013-10-01
The structural diversity of arsenic sulfide molecules in compositions between As4S6 and As4 was investigated using ab initio quantum chemical calculations. The As4S6 molecule consists of four trigonal pyramid coordinations of As atoms bonding to three S atoms. In the As4S5 composition, only one type of molecular configuration corresponds to an uzonite-type molecule. In the As4S4 composition, two molecular configurations exist with realgar-type and pararealgar-type molecules. Three molecular configurations are in the As4S3 composition. The first configuration comprises trigonal pyramidal As atom coordinations of two types: bonding to two S atoms and one As atom, and bonding to one S atom and two As atoms. The second is the molecular configuration of dimorphite. The third comprises trigonal pyramidal As atom coordinations of two types: bonding to three As atoms, and bonding to one As atom and two S atoms. The As4S2 composition allows molecular configurations of two types. One is comprised of trigonal pyramidal As atom configurations of one type bonding to two As atoms and one S atom. The other comprises trigonal pyramidal As atom coordinations of three types: bonding to two S atoms and one As atoms, bonding to one S atom and two As atoms, and bonding to three As atoms. The As4S molecule has trigonal pyramidal As atom coordinations of two types: bonding to one S atom and two As atoms, and bonding to three As atoms. The As4S composition permits only one molecular configuration, which suggests that the mineral duranusite comprises the As4S molecular geometry. In all, ten molecular configurations are predicted in the molecular hierarchy of the arsenic sulfide binary system. The simulated Raman spectral profiles are helpful in searching for undiscovered arsenic sulfide minerals.
NASA Astrophysics Data System (ADS)
Akin-Ojo, Omololu; Song, Yang; Wang, Feng
2008-08-01
A new method called adaptive force matching (AFM) has been developed that is capable of producing high quality force fields for condensed phase simulations. This procedure involves the parametrization of force fields to reproduce ab initio forces obtained from condensed phase quantum-mechanics/molecular-mechanics (QM/MM) calculations. During the procedure, the MM part of the QM/MM is iteratively improved so as to approach ab initio quality. In this work, the AFM method has been tested to parametrize force fields for liquid water so that the resulting force fields reproduce forces calculated using the ab initio MP2 and the Kohn-Sham density functional theory with the Becke-Lee-Yang-Parr (BLYP) and Becke three-parameter LYP (B3LYP) exchange correlation functionals. The AFM force fields generated in this work are very simple to evaluate and are supported by most molecular dynamics (MD) codes. At the same time, the quality of the forces predicted by the AFM force fields rivals that of very expensive ab initio calculations and are found to successfully reproduce many experimental properties. The site-site radial distribution functions (RDFs) obtained from MD simulations using the force field generated from the BLYP functional through AFM compare favorably with the previously published RDFs from Car-Parrinello MD simulations with the same functional. Technical aspects of AFM such as the optimal QM cluster size, optimal basis set, and optimal QM method to be used with the AFM procedure are discussed in this paper.
NASA Technical Reports Server (NTRS)
Wang, Dunyou
2003-01-01
A time-dependent wave-packet approach is presented for the quantum dynamics study of the AB+CDE reaction system for zero total angular momentum. A seven-degree-of-freedom calculation is employed to study the chemical reaction of H2+C2H yields H + C2H2 by treating C2H as a linear molecule. Initial state selected reaction probabilities are presented for various initial ro-vibrational states. This study shows that vibrational excitation of H2 enhances the reaction probability, whereas the excitation of C2H has only a small effect on the reactivity. An integral cross section is also reported for the initial ground states of H2 and C2H. The theoretical and experimental results agree with each other very well when the calculated seven dimensional results are adjusted to account for the lower transition state barrier heights found in recent ab initio calculations.
Li, Junjie; Li, Xiaohu; Iyengar, Srinivasan S
2014-06-10
We discuss a multiconfigurational treatment of the "on-the-fly" electronic structure within the quantum wavepacket ab initio molecular dynamics (QWAIMD) method for coupled treatment of quantum nuclear effects with electronic structural effects. Here, multiple single-particle electronic density matrices are simultaneously propagated with a quantum nuclear wavepacket and other classical nuclear degrees of freedom. The multiple density matrices are coupled through a nonorthogonal configuration interaction (NOCI) procedure to construct the instantaneous potential surface. An adaptive-mesh-guided set of basis functions composed of Gaussian primitives are used to simplify the electronic structure calculations. Specifically, with the replacement of the atom-centered basis functions positioned on the centers of the quantum-mechanically treated nuclei by a mesh-guided band of basis functions, the two-electron integrals used to compute the electronic structure potential surface become independent of the quantum nuclear variable and hence reusable along the entire Cartesian grid representing the quantum nuclear coordinates. This reduces the computational complexity involved in obtaining a potential surface and facilitates the interpretation of the individual density matrices as representative diabatic states. The parametric nuclear position dependence of the diabatic states is evaluated at the initial time-step using a Shannon-entropy-based sampling function that depends on an approximation to the quantum nuclear wavepacket and the potential surface. This development is meant as a precursor to an on-the-fly fully multireference electronic structure procedure embedded, on-the-fly, within a quantum nuclear dynamics formalism. We benchmark the current development by computing structural, dynamic, and spectroscopic features for a series of bihalide hydrogen-bonded systems: FHF(-), ClHCl(-), BrHBr(-), and BrHCl(-). We find that the donor-acceptor structural features are in good
Long, Run; Prezhdo, Oleg V
2015-07-01
Hybrid organic/inorganic polymer/quantum dot (QD) solar cells are an attractive alternative to the traditional cells. The original, simple models postulate that one-dimensional polymers have continuous energy levels, while zero-dimensional QDs exhibit atom-like electronic structure. A realistic, atomistic viewpoint provides an alternative description. Electronic states in polymers are molecule-like: finite in size and discrete in energy. QDs are composed of many atoms and have high, bulk-like densities of states. We employ ab initio time-domain simulation to model the experimentally observed ultrafast photoinduced dynamics in a QD/polymer hybrid and show that an atomistic description is essential for understanding the time-resolved experimental data. Both electron and hole transfers across the interface exhibit subpicosecond time scales. The interfacial processes are fast due to strong electronic donor-acceptor, as evidenced by the densities of the photoexcited states which are delocalized between the donor and the acceptor. The nonadiabatic charge-phonon coupling is also strong, especially in the polymer, resulting in rapid energy losses. The electron transfer from the polymer is notably faster than the hole transfer from the QD, due to a significantly higher density of acceptor states. The stronger molecule-like electronic and charge-phonon coupling in the polymer rationalizes why the electron-hole recombination inside the polymer is several orders of magnitude faster than in the QD. As a result, experiments exhibit multiple transfer times for the long-lived hole inside the QD, ranging from subpicoseconds to nanoseconds. In contrast, transfer of the short-lived electron inside the polymer does not occur beyond the first picosecond. The energy lost by the hole on its transit into the polymer is accommodated by polymer's high-frequency vibrations. The energy lost by the electron injected into the QD is accommodated primarily by much lower-frequency collective and
NASA Astrophysics Data System (ADS)
Yung, Y. Y.; Choi, B. H.; Tang, K. T.
1980-01-01
Three dimensional quantum mechanical calculations are carried out for the reactive scattering of D+H2→DH+H on the ab initio potential energy surface calculated by Liu and Siegbahn and fitted by Truhlar and Horowitz. The differential and total cross sections as well as the S matrix elements are obtained from the adiabatic distorted wave method. Threshold energy, cross sections and product distributions over final states are all in good agreement with experimental measurements. Results are also compared with the corresponding ones obtained on the Porter-Karplus and the Yates-Lester semi-empirical surfaces.
Hua -Gen Yu; Han, Huixian; Guo, Hua
2016-03-29
Vibrational energy levels of the ammonium cation (NH4+) and its deuterated isotopomers are calculated using a numerically exact kinetic energy operator on a recently developed nine-dimensional permutation invariant semiglobal potential energy surface fitted to a large number of high-level ab initio points. Like CH4, the vibrational levels of NH4+ and ND4+ exhibit a polyad structure, characterized by a collective quantum number P = 2(v1 + v3) + v2 + v4. As a result, the low-lying vibrational levels of all isotopomers are assigned and the agreement with available experimental data is better than 1 cm–1.
Yamaji, Youhei
2015-12-31
Recently, condensed-matter ab initio approaches to strongly correlated electrons confined in crystalline solids have been developed and applied to transition-metal oxides and molecular conductors. In this paper, an ab initio scheme based on constrained random phase approximations and localized Wannier orbitals is applied to a spin liquid candidate Na{sub 2}IrO{sub 3} and is shown to reproduce experimentally observed specific heat.
Xie, Binbin; Liu, Lihong; Cui, Ganglong; Fang, Wei-Hai; Cao, Jun; Feng, Wei; Li, Xin-qi
2015-11-21
In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N{sub 2}CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N{sub 2}CO photodissociation at λ > 335 nm is an ultrafast process and the two C—N bonds are broken in a stepwise way, giving birth to CO and N{sub 2} as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C—N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes.
NASA Astrophysics Data System (ADS)
Farahani, Pooria; Lundberg, Marcus; Karlsson, Hans O.
2013-11-01
The SN2 substitution reactions at phosphorus play a key role in organic and biological processes. Quantum molecular dynamics simulations have been performed to study the prototype reaction Cl-+PH2Cl→ClPH2+Cl-, using one and two-dimensional models. A potential energy surface, showing an energy well for a transition complex, was generated using ab initio electronic structure calculations. The one-dimensional model is essentially reflection free, whereas the more realistic two-dimensional model displays involved resonance structures in the reaction probability. The reaction rate is almost two orders of magnitude smaller for the two-dimensional compared to the one-dimensional model. Energetic errors in the potential energy surface is estimated to affect the rate by only a factor of two. This shows that for these types of reactions it is more important to increase the dimensionality of the modeling than to increase the accuracy of the electronic structure calculation.
Kishi, Ryohei; Fujii, Hiroaki; Minami, Takuya; Shigeta, Yasuteru; Nakano, Masayoshi
2015-01-22
In this study, we apply the ab initio molecular orbital - configuration interaction based quantum master equation (MOQME) approach to the calculation and analysis of the dynamic first hyperpolarizabilities (β) of asymmetric π-conjugated molecules. In this approach, we construct the excited state models by the ab initio configuration interaction singles method. Then, time evolutions of system reduced density matrix ρ(t) and system polarization p(t) are calculated by the QME approach. Dynamic β in the second harmonic generation is calculated based on the nonperturbative definition of nonlinear optical susceptibility, using the frequency domain system polarization p(ω). Spatial contributions of electrons to β are analyzed based on the dynamic hyperpolarizability density map, which visualizes the second-order response of charge density oscillating with a frequency of 2ω. We apply the present method to the calculation of the dynamic β of a series of donor/acceptor substituted polyene oligomers, and then discuss the applicability of the MOQME method to the calculation and analysis of dynamic NLO properties of molecular systems.
NASA Astrophysics Data System (ADS)
Passler, Peter P.; Rode, Bernd M.
2015-10-01
Hydration properties of trivalent thulium and ytterbium ions in aqueous solution are investigated via quantum mechanical charge field molecular dynamics (QMCF-MD) simulations. The QMCF-MD formalism is a special type of QM/MM simulation, where the chemically most relevant part of the system - in this case the ion with its first and second hydration shells - is treated by quantum mechanics. The mean ionsbnd O distances and the average coordination numbers of the first hydration shells are compared with experimental EXAFS data. Mean ligand residence times, vibrational frequencies and force constants were evaluated to characterise the dynamics of the systems.
NASA Astrophysics Data System (ADS)
Reyes-Lillo, Sebastian E.; Rangel, Tonatiuh; Bruneval, Fabien; Neaton, Jeffrey B.
2016-07-01
The Ruddlesden-Popper (RP) homologous series Srn +1TinO3 n +1 provides a useful template for the study and control of the effects of dimensionality and quantum confinement on the excited state properties of the complex oxide SrTiO3. We use ab initio many-body perturbation theory within the G W approximation and the Bethe-Salpeter equation approach to calculate quasiparticle energies and absorption spectra of Srn +1TinO3 n +1 for n =1 -5 and ∞ . Our computed direct and indirect optical gaps are in excellent agreement with spectroscopic measurements. The calculated optical spectra reproduce the main experimental features and reveal excitonic structure near the gap edge. We find that electron-hole interactions are important across the series, leading to significant exciton binding energies that increase for small n and reach a value of 330 meV for n =1 , a trend attributed to increased quantum confinement. We find that the lowest-energy singlet exciton of Sr2TiO4 (n =1 ) localizes in the two-dimensional plane defined by the TiO2 layer, and we explain the origin of its localization.
NASA Astrophysics Data System (ADS)
Stepanian, S. G.; Reva, I. D.; Radchenko, E. D.; Latajka, Z.; Wierzejewska, M.; Ratajczak, H.
1999-06-01
The molecular structure of 1-methylaminophosphinic acid (AMPA) was investigated with the matrix isolation IR spectroscopy and ab initio calculations performed with RHF, MP2, MP3, MP4(DQ), MP4(SDQ) and MP4(SDTQ) methods. Three pseudopotential basis sets designed as CEP-31G were used in the calculations: Basis Set I-CEP-31G with the d-functions on phosphorus; Basis Set II-CEP-31G with the d-functions on all heavy atoms; Basis Set III-CEP-31G with the d-functions on all heavy atoms and p-functions on hydrogens. Four stable molecular and four stable zwitterion conformers of aminophosphinic acid were found via ab initio calculations. According to the calculations, molecular conformers are always more stable than the zwitterion conformers, irrespective of the basis set size and level of theory. This result is in good agreement with matrix IR spectrum of the AMPA. The presence of the bands of OH stretching and NH 2 bending vibrations and the absence of the bands of POO - and NH 3+ vibrations are the evidence of molecular structure of AMPA in the isolated state. An increased number of vibrational bands is found in the IR spectrum. It is explained by the high conformation lability of AMPA molecules which is related to very low barrier of rotation about C-P bond. The IR spectrum is actually determined by multiple sites of AMPA molecule packed in the Ar crystal, which considerably increases the number of bands in the IR spectrum.
Topological states in α -Sn and HgTe quantum wells: A comparison of ab initio results
NASA Astrophysics Data System (ADS)
Küfner, Sebastian; Bechstedt, Friedhelm
2015-01-01
Both α -Sn and HgTe are expected to have similar topological properties because of their inverted band structure and zero-gap character. We investigate how the different crystal symmetries and the bonding to barrier materials act to the quantum phase transition versus the thickness of the corresponding quantum well (QW) structures. They are simulated by (SnSn)N(CdTe)M and (HgTe)N(CdTe)M(110 ) superlattices. Their electronic structures and eigenstates are studied by means of first-principles calculations using the modified Becke-Johnson exchange-correlation functional and spin-orbit interaction. Significant differences are observed for the two QW materials. A topological transition between trivial insulator and quantum spin Hall phase together with the formation of topologically protected edge states are observed in the case of HgTe QWs, while these features are considerably modified for α -Sn . The different behaviors are discussed in the light of the different symmetry, spin-orbit interaction, and interface bonding. For a better understanding of the influence of the interface electrostatics also results for (HgTe)N(InSb)M(110 ) systems are discussed.
Aziz, Saadullah G.; Alyoubi, Abdulrahman O.; Elroby, Shaaban A.; Osman, Osman I.; Hilal, Rifaat H.
2015-01-01
The present study aims at a fundamental understanding of bonding characteristics of the C–Br and O–Br bonds. The target molecular systems are the isomeric CH3OBr/BrCH2OH system and their decomposition products. Calculations of geometries and frequencies at different density functional theory (DFT) and Hartree–Fock/Møller–Plesset (HF/MP2) levels have been performed. Results have been assessed and evaluated against those obtained at the coupled cluster single-double (Triplet) (CCSD(T)) level of theory. The characteristics of the C–Br and O–Br bonds have been identified via analysis of the electrostatic potential, natural bond orbital (NBO), and quantum theory of atoms in molecules (QTAIM). Analysis of the electrostatic potential (ESP) maps enabled the quantitative characterization of the Br σ-holes. Its magnitude seems very sensitive to the environment and the charge accumulated in the adjacent centers. Some quantum topological parameters, namely ∇2ρ, ellipticity at bond critical points and the Laplacian bond order, were computed and discussed. The potential energy function for internal rotation has been computed and Fourier transformed to characterize the conformational preferences and origin of the barriers. NBO energetic components for rotation about the C–Br and O–Br bonds as a function of torsion angle have been computed and displayed. PMID:25815595
NASA Astrophysics Data System (ADS)
Yamashita, Takefumi; Takatsuka, Kazuo
2007-02-01
The infrared spectrum of phenol-water cationic cluster, [PhOH•H2O]+, taken by Sawamura et al. [J. Phys. Chem. 100, 8131 (1996)] is puzzling in that the peak due to the stretching mode of the phenolic OH (3657cm-1 for a neutral monomer and 3524cm-1 for PhOH•H2O) seemingly disappears and instead an extremely broad tail extending down to 2900cm-1 is observed. The present authors theoretically ascribe this anomalous spectrum to an inhomogeneous broadening of the OH stretching peak caused by the hydrogen bond, the strength of which has been greatly enhanced by ionization of the phenyl ring. Indeed they estimate that the peak position is at 2300cm-1 and the spectral width can become as wide as 1000cm-1 at the cluster energy of 32kcal/mol. This surprisingly wide broadening can be generic in hydrogen-bond systems, which in turn is useful to study the nature of the hydrogen-bond assisted dynamics in various systems such as those in DNA and proteins. To study the present system quantitatively, the authors have developed an ab initio mixed quantum-classical method, in which the nuclear motions on an adiabatic ab initio potential surface are treated such that only the OH stretching motion is described quantum mechanically, while all the other remaining modes are treated classically with on-the-fly scheme. This method includes the implementation of many numerical methodologies, which enables it to deal with a relatively large molecular system. With this theoretical method, the authors analyze the present anomalous broadening in a great detail. In particular, they suggest that one can extract direct information about the hydrogen-bond dynamics with respect to the clear correlation between the vibrational excitation energy of the OH stretching and intermolecular distance by means of a time-resolved infrared spectroscopy: Reflecting the slow and wide-range variation of the intermolecular distance of the relevant hydrogen bond, the time-resolved spectrum is predicted to vary
NASA Astrophysics Data System (ADS)
Ramya, T.; Gunasekaran, S.; Ramkumaar, G. R.
2015-10-01
The experimental and theoretical spectra of (S)-2-Oxopyrrolidin-1-yl Butanamide (S2OPB) were studied. FT-IR and FT-Raman spectra of S2OPB in the solid phase were recorded and analyzed in the range 4000-450 and 5000-50 cm-1 respectively. The structural and spectroscopic analyses of S2OPB were calculated using ab initio Hartree Fock (HF) and density functional theory calculations (B3PW91, B3LYP) with 6-31G(d,p) basis set. A complete vibrational interpretation has been made on the basis of the calculated Potential Energy Distribution (PED). The HF, B3LYP and B3PW91 methods based NMR calculation has been used to assign the 1H NMR and 13C NMR chemical shift of S2OPB. Comparative study on UV-Vis spectral analysis between the experimental and theoretical (B3PW91, B3LYP) methods and the global chemical parameters and local descriptor of reactivity through the Fukui function were performed. Finally the thermodynamic properties of S2OPB were calculated at different temperatures and the corresponding relations between the properties and temperature were also studied.
Ramya, T; Gunasekaran, S; Ramkumaar, G R
2015-10-01
The experimental and theoretical spectra of (S)-2-Oxopyrrolidin-1-yl Butanamide (S2OPB) were studied. FT-IR and FT-Raman spectra of S2OPB in the solid phase were recorded and analyzed in the range 4000-450 and 5000-50 cm(-1) respectively. The structural and spectroscopic analyses of S2OPB were calculated using ab initio Hartree Fock (HF) and density functional theory calculations (B3PW91, B3LYP) with 6-31G(d,p) basis set. A complete vibrational interpretation has been made on the basis of the calculated Potential Energy Distribution (PED). The HF, B3LYP and B3PW91 methods based NMR calculation has been used to assign the (1)H NMR and (13)C NMR chemical shift of S2OPB. Comparative study on UV-Vis spectral analysis between the experimental and theoretical (B3PW91, B3LYP) methods and the global chemical parameters and local descriptor of reactivity through the Fukui function were performed. Finally the thermodynamic properties of S2OPB were calculated at different temperatures and the corresponding relations between the properties and temperature were also studied. PMID:25956325
Mewes, Stefanie A; Mewes, Jan-Michael; Dreuw, Andreas; Plasser, Felix
2016-01-28
Excitonic effects play a fundamental role in the photophysics of organic semiconductors such as poly(para phenylene vinylene) (PPV). The emergence of these effects is examined for PPV oligomers based on high level ab initio excited-state calculations. The computed many-body wavefunctions are subjected to our recently developed exciton analysis protocols to provide a qualitative and quantitative characterization of excitonic effects. The discussion is started by providing high-level benchmark calculations using the algebraic-diagrammatic construction for the polarization propagator in third order of perturbation theory (ADC(3)). These calculations support the general adequacy of the computationally more efficient ADC(2) method in the case of singly excited states but also reveal the existence of low-energy doubly excited states. In a next step, a series of oligomers with chains of two to eight phenyl rings is studied at the ADC(2) level showing that the confinement effects are dominant for small oligomers, while delocalized exciton bands emerge for larger systems. In the case of the largest oligomer, the first twenty singlet and triplet excited states are computed and a detailed analysis in terms of the Wannier and Frenkel models is presented. The presence of different Wannier bands becomes apparent, showing a general trend that exciton sizes are lowered with increasing quasi-momentum within the bands. PMID:26700493
NASA Astrophysics Data System (ADS)
Prasetyo, Niko; Armunanto, Ria
2016-05-01
Structures and dynamics of Ag+ in 18.6% aqueous ammonia have been studied using Quantum Mechanical Charge Field Molecular Dynamics (QMCF-MD) simulation at the Hartree-Fock (HF) level theory employing LANL2DZ ECP basis set for Ag+ and Dunning DZP for solvent molecules. Structural properties are in excellent agreement with previous QM/MM and experiments studies. [Ag(NH3)2(H2O)3]+ was found as dominant species during simulation time. For 20 ps of simulation time, a labile first solvation shell was observed with both fast ammonia and water ligands exchanges. QMCF-MD framework describes first solvation shell more labile than conventional QM/MM MD simulation.
The excited J = 01 Σu+ levels of D2: Measurements and ab initio quantum defect study
NASA Astrophysics Data System (ADS)
Glass-Maujean, M.; Jungen, Ch.; Dickenson, G. D.; de Oliveira, N.; Ubachs, W.
2016-02-01
The DESIRS beamline of the SOLEIL synchrotron facility, equipped with a vacuum ultraviolet Fourier-transform spectrometer has been used to measure P (N″ = 1) (N -N″ = - 1) absorption transitions of the D2 molecule. Some 44 P-lines were assigned and their transition frequencies determined up to excitation energies of 134,000 cm-1 above the ground state, thereby extending the earlier work by various authors, and considerably improving the spectral accuracy (<0.1 cm-1). The assignments have been aided by first principles multichannel quantum defect theory (MQDT) calculations. These calculations also provide predictions of the autoionization widths of the upper levels which agree well with the observed resonance widths.
NASA Astrophysics Data System (ADS)
Wouters, Sebastian; Poelmans, Ward; Ayers, Paul W.; Van Neck, Dimitri
2014-06-01
The density matrix renormalization group (DMRG) has become an indispensable numerical tool to find exact eigenstates of finite-size quantum systems with strong correlation. In the fields of condensed matter, nuclear structure and molecular electronic structure, it has significantly extended the system sizes that can be handled compared to full configuration interaction, without losing numerical accuracy. For quantum chemistry (QC), the most efficient implementations of DMRG require the incorporation of particle number, spin and point group symmetries in the underlying matrix product state (MPS) ansatz, as well as the use of so-called complementary operators. The symmetries introduce a sparse block structure in the MPS ansatz and in the intermediary contracted tensors. If a symmetry is non-abelian, the Wigner-Eckart theorem allows to factorize a tensor into a Clebsch-Gordan coefficient and a reduced tensor. In addition, the fermion signs have to be carefully tracked. Because of these challenges, implementing DMRG efficiently for QC is not straightforward. Efficient and freely available implementations are therefore highly desired. In this work we present CheMPS2, our free open-source spin-adapted implementation of DMRG for ab initio QC. Around CheMPS2, we have implemented the augmented Hessian Newton-Raphson complete active space self-consistent field method, with exact Hessian. The bond dissociation curves of the 12 lowest states of the carbon dimer were obtained at the DMRG(28 orbitals, 12 electrons, DSU(2) = 2500)/cc-pVDZ level of theory. The contribution of 1 s core correlation to the X1Σg+ bond dissociation curve of the carbon dimer was estimated by comparing energies at the DMRG(36o, 12e, DSU(2) = 2500)/cc-pCVDZ and DMRG-SCF(34o, 8e, DSU(2) = 2500)/cc-pCVDZ levels of theory.
NASA Astrophysics Data System (ADS)
Churakov, S. V.; Parrinello, M.
2003-04-01
The carbonation of brucite (Mg(OH)2) has been considered as a potential technology for cleaning industrial carbon dioxide waste. The kinetics of the reaction Mg(OH)2 + CO2 -> MgCO3 + H2O have been studied experimentally at 573°C by Bearat at al. [1]. Their experiments suggest that the carbonation of magnesium hydroxide proceeds by the reaction Mg(OH)2 -> MgO + H2O followed by the adsorption of CO2 molecules on the dehydrated brucite surface. Due to the large difference in volumes between Mg(OH)2 and MgO, dehydration causes the formation of dislocations and cracks, allowing water molecules to leave the brucite surface and facilitating the advance of the carbonation front in the bulk solid. The detailed mechanism of this process is however unknown. We used the Car-Parrinello ab initio molecular dynamics method to study the structure and dynamics of the (0001), (1-100) and (11-20) surfaces of brucit and calculated the enthalpy and activation barrier of H2O nucleation and dehydration on different surfaces. The results obtained are in agreement with previous studies of brucite dehydration by Masini and Bernasconi [2]. The reactive Car-Parrinello molecular dynamics method [3] has been applied to investigate the detailed mechanism of the dehydration-carbonation reaction at the (1-100) interface of brucite with the gas phase. Based on the results of our MD simulations and the calculated enthalpy of CO2 adsorption on the dehydrated brucite surfaces we propose a mechanism for the dehydration/carbonation reaction. [1] Bearat H, McKelvy MJ, Chizmeshya AVG, Sharma R, Carpenter RW (2002) J. Amer. Ceram. Soc. 85(4):742 [2] Masini P and Bernasconi M (2001) J. Phys. Cond. Mat. 13: 1-12 [3] Iannuzzi M, Laio A and Parrinello M (2003) Phys. Rev. Lett. (submitted)
NASA Astrophysics Data System (ADS)
de la Mora, Pablo; Cosio-Castañeda, Carlos; Martinez-Anaya, Oliver; Morales, Francisco; Tavizon, Gustavo
2016-09-01
In this work, a theoretical study of the electrical properties of the Bi2-ySryIr2O7 (Bi2-ySryIr2O16O2) α-pyrochlore-type solid solution is presented. Quantum ab initio DFT(WIEN2k) calculations were performed in order to understand the electrical resistivity changes associated to the Bi substitution by Sr in this system. The main crystallographic modification associated to this substitution is the x position of the 48f oxygen (x, 1/8, 1/8) (O1); this substitution substantially modifies the Bi/Sr-O1 and Ir-O1 atomic distances, increasing the former and diminishing the latter. Experimentally, the Bi2-ySryIr2O7 samples are metallic and the electrical resistivity increases with the Sr content. Electronic structure calculations for Bi2Ir2O7 and BiSrIr2O7 show that, regardless of structural changes, there is only a small change of electrical conductivity with the Sr substitution, and the experimentally observed increase of the resistivity can be explained in terms of a larger impact on the electronic structure of both; the Sr 'impurities' as well as of the thermal Sr oscillations.
NASA Astrophysics Data System (ADS)
Bandekar, Jagdeesh; Klima, Suzanne
1992-10-01
Ab initio quantum chemical computations were carried out on (a) dimethyl ether, (b) N-methyl formamide, (c) dimethyl ether- N-methyl formamide complex, and (d) N-methyl formamide dimer to compute the strengths of hydrogen bonds (H-bonds) between the NH groups and CO and ether COC groups. The basis set used was the 3-21G set of the GAUSSIAN 80 program obtained from QCPE, Bloomington, IN. Variations in the strengths of these two H-bonds with the N . . O distance (where O is either carbonyl or ether group oxygen) were studied and found to be similar in behavior. The strength of the "ether" hydrogen bond is computed to be 10.32 kcal mol -1, which is quite significant compared to the value of 10.11 kcal mol -1 for the more accepted "carbonyl" hydrogen bond. The "ether" hydrogen bond is found to be directional, specific and non-negligible. Work with two more basis sets has indicated that the results so obtained are not dependent on their choice. Possible importance of such a hydrogen bond in polyurethanes, inhalation anesthetics, and depsi-peptides is indicated.
Wood, Geoffrey P F; Sreedhara, Alavattam; Moore, Jamie M; Wang, John; Trout, Bernhardt L
2016-05-12
An assessment of the mechanisms of (•)OH and (•)OOH radical-mediated oxidation of tryptophan was performed using density functional theory calculations and ab initio plane-wave Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics simulations. For the (•)OH reactions, addition to the pyrrole ring at position 2 is the most favored site with a barrierless reaction in the gas phase. The subsequent degradation of this adduct through a H atom transfer to water was intermittently observed in aqueous-phase molecular dynamics simulations. For the (•)OOH reactions, addition to the pyrrole ring at position 2 is the most favored pathway, in contrast to the situation in the model system ethylene, where concerted addition to the double bond is preferred. From the (•)OOH position 2 adduct QM/MM simulations show that formation of oxy-3-indolanaline occurs readily in an aqueous environment. The observed transformation starts from an initial rupture of the O-O bond followed by a H atom transfer with the accompanying loss of an (•)OH radical to solution. Finally, classical molecular dynamics simulations were performed to equate observed differential oxidation rates of various tryptophan residues in monoclonal antibody fragments. It was found that simple parameters derived from simulation correlate well with the experimental data. PMID:27082439
Chen, Liuyang; Shao, Kejie; Chen, Jun; Yang, Minghui; Zhang, Dong H
2016-05-21
This work performs a time-dependent wavepacket study of the H2 + C2H → H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H↔H + C2H2, H + C2H2 → HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H → H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values. PMID:27208951
Song, Hongwei; Lu, Yunpeng; Li, Jun; Yang, Minghui; Guo, Hua
2016-04-28
An initial state selected time-dependent wave packet method is applied to study the dynamics of the OH + CHD3 reaction with a six-dimensional model on a newly developed full-dimensional ab initio potential energy surface (PES). This quantum dynamical (QD) study is complemented by full-dimensional quasi-classical trajectory (QCT) calculations on the same PES. The QD results indicate that both translational energy and the excitation of the CH stretching mode significantly promote the reaction while the excitation of the umbrella mode has a negligible effect on the reactivity. For this early barrier reaction, interestingly, the CH stretching mode is more effective than translational energy in promoting the reaction except at very low collision energies. These QD observations are supported by QCT results. The higher efficacy of the CH stretching model in promoting this early barrier reaction is inconsistent with the prediction of the naively extended Polanyi's rules, but can be rationalized by the recently proposed sudden vector projection model. PMID:27131546
NASA Astrophysics Data System (ADS)
Song, Hongwei; Lu, Yunpeng; Li, Jun; Yang, Minghui; Guo, Hua
2016-04-01
An initial state selected time-dependent wave packet method is applied to study the dynamics of the OH + CHD3 reaction with a six-dimensional model on a newly developed full-dimensional ab initio potential energy surface (PES). This quantum dynamical (QD) study is complemented by full-dimensional quasi-classical trajectory (QCT) calculations on the same PES. The QD results indicate that both translational energy and the excitation of the CH stretching mode significantly promote the reaction while the excitation of the umbrella mode has a negligible effect on the reactivity. For this early barrier reaction, interestingly, the CH stretching mode is more effective than translational energy in promoting the reaction except at very low collision energies. These QD observations are supported by QCT results. The higher efficacy of the CH stretching model in promoting this early barrier reaction is inconsistent with the prediction of the naively extended Polanyi's rules, but can be rationalized by the recently proposed sudden vector projection model.
NASA Astrophysics Data System (ADS)
de la Mora, Pablo; Cosio-Castañeda, Carlos; Martinez-Anaya, Oliver; Morales, Francisco; Tavizon, Gustavo
2016-09-01
In this work, a theoretical study of the electrical properties of the Bi2-ySryIr2O7 (Bi2-ySryIr2O16O2) α-pyrochlore-type solid solution is presented. Quantum ab initio DFT(WIEN2k) calculations were performed in order to understand the electrical resistivity changes associated to the Bi substitution by Sr in this system. The main crystallographic modification associated to this substitution is the x position of the 48f oxygen (x, 1/8 , 1/8 ) (O1); this substitution substantially modifies the Bi/Sr-O1 and Ir-O1 atomic distances, increasing the former and diminishing the latter. Experimentally, the Bi2-ySryIr2O7 samples are metallic and the electrical resistivity increases with the Sr content. Electronic structure calculations for Bi2Ir2O7 and BiSrIr2O7 show that, regardless of structural changes, there is only a small change of electrical conductivity with the Sr substitution, and the experimentally observed increase of the resistivity can be explained in terms of a larger impact on the electronic structure of both; the Sr 'impurities' as well as of the thermal Sr oscillations.
Liu, Jianbo; Miller, William H.; Fanourgakis, G. S.; Xantheas, Sotiris S.; Imoto, Sho; Saito, Shinji
2011-12-28
The dynamical properties of liquid water play an important role in many processes in Nature. In this paper we focus on the infrared (IR) absorption spectrum of liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with the local Gaussian approximation (LGA) [Liu and Miller, J. Chem. Phys. 131, 074113 (2009)] and an ab initio based, flexible, polarizable Thole-type model (TTM3-F) [Fanourgakis and Xantheas, J. Chem. Phys. 128, 074506 (2008)]. Although the LSC-IVR (LGA) gives the exact result for the isolated 3-dimensional shifted harmonic stretching model, it yields a blue-shifted peak position for the more realistic anharmonic stretching potential. By using the short time information of the LSCIVR correlation function, however, it is shown how one can obtain more accurate results for the position of the stretching peak. Due to the physical decay in the condensed phase system, the LSC-IVR (LGA) is a good and practical approximate quantum approach for the IR spectrum of liquid water. The present results offer valuable insight into future attempts to improve the accuracy of the TTM3-F potential in reproducing the IR spectrum of liquid water.
NASA Astrophysics Data System (ADS)
Chen, Liuyang; Shao, Kejie; Chen, Jun; Yang, Minghui; Zhang, Dong H.
2016-05-01
This work performs a time-dependent wavepacket study of the H2 + C2H → H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H↔H + C2H2, H + C2H2 → HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H → H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values.
NASA Astrophysics Data System (ADS)
Dornheim, T.; Groth, S.; Schoof, T.; Hann, C.; Bonitz, M.
2016-05-01
In a recent publication [S. Groth et al., Phys. Rev. B 93, 085102 (2016), 10.1103/PhysRevB.93.085102], we have shown that the combination of two complementary quantum Monte Carlo approaches, namely configuration path integral Monte Carlo [T. Schoof et al., Phys. Rev. Lett. 115, 130402 (2015), 10.1103/PhysRevLett.115.130402] and permutation blocking path integral Monte Carlo [T. Dornheim et al., New J. Phys. 17, 073017 (2015), 10.1088/1367-2630/17/7/073017], allows for the accurate computation of thermodynamic properties of the spin-polarized uniform electron gas over a wide range of temperatures and densities without the fixed-node approximation. In the present work, we extend this concept to the unpolarized case, which requires nontrivial enhancements that we describe in detail. We compare our simulation results with recent restricted path integral Monte Carlo data [E. W. Brown et al., Phys. Rev. Lett. 110, 146405 (2013), 10.1103/PhysRevLett.110.146405] for different energy contributions and pair distribution functions and find, for the exchange correlation energy, overall better agreement than for the spin-polarized case, while the separate kinetic and potential contributions substantially deviate.
Many-body ab initio diffusion quantum Monte Carlo applied to the strongly correlated oxide NiO
Mitra, Chandrima; Krogel, Jaron T.; Santana, Juan A.; Reboredo, Fernando A.
2015-10-28
We present a many-body diffusion quantum Monte Carlo (DMC) study of the bulk and defect properties of NiO. We find excellent agreement with experimental values, within 0.3%, 0.6%, and 3.5% for the lattice constant, cohesive energy, and bulk modulus, respectively. The quasiparticle bandgap was also computed, and the DMC result of 4.72 (0.17) eV compares well with the experimental value of 4.3 eV. Furthermore, DMC calculations of excited states at the L, Z, and the gamma point of the Brillouin zone reveal a flat upper valence band for NiO, in good agreement with Angle Resolved Photoemission Spectroscopy results. To study defect properties, we evaluated the formation energies of the neutral and charged vacancies of oxygen and nickel in NiO. A formation energy of 7.2 (0.15) eV was found for the oxygen vacancy under oxygen rich conditions. For the Ni vacancy, we obtained a formation energy of 3.2 (0.15) eV under Ni rich conditions. These results confirm that NiO occurs as a p-type material with the dominant intrinsic vacancy defect being Ni vacancy.
NASA Astrophysics Data System (ADS)
Zeng, Xiancheng; Hu, Hao; Hu, Xiangqian; Yang, Weitao
2009-04-01
A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
Zeng Xiancheng; Hu Hao; Hu Xiangqian; Yang Weitao
2009-04-28
A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids 'on-the-fly' QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
Bankura, Arindam; Chandra, Amalendu
2015-01-28
The dynamics of proton transfer (PT) through hydrogen bonds in a two-dimensional water layer confined between two graphene sheets at room temperature are investigated through ab initio and quantum-classical simulations. The excess proton is found to be mostly solvated as an Eigen cation where the hydronium ion donates three hydrogen bonds to the neighboring water molecules. In the solvation shell of the hydronium ion, the three coordinated water molecules with two donor hydrogen bonds are found to be properly presolvated to accept a proton. Although no hydrogen bond needs to be broken for transfer of a proton to such presolvated water molecules from the hydronium ion, the PT rate is still found to be not as fast as it is for one-dimensional chains. Here, the PT is slowed down as the probability of finding a water with two donor hydrogen bonds in the solvation shell of the hydronium ion is found to be only 25%-30%. The hydroxide ion is found to be solvated mainly as a complex anion where it accepts four H-bonds through its oxygen atom and the hydrogen atom of the hydroxide ion remains free all the time. Here, the presolvation of the hydroxide ion to accept a proton requires that one of its hydrogen bonds is broken and the proton comes from a neighboring water molecule with two acceptor and one donor hydrogen bonds. The coordination number reduction by breaking of a hydrogen bond is a slow process, and also the population of water molecules with two acceptor and one donor hydrogen bonds is only 20%-25% of the total number of water molecules. All these factors together tend to slow down the hydroxide ion migration rate in two-dimensional water layers compared to that in three-dimensional bulk water.
Long, Run; Prezhdo, Oleg V
2011-11-30
Following recent experiments [Science 2010, 328, 1543; PNAS 2011, 108, 965], we report an ab initio nonadiabatic molecular dynamics (NAMD) simulation of the ultrafast photoinduced electron transfer (ET) from a PbSe quantum dot (QD) into the rutile TiO(2) (110) surface. The system forms the basis for QD-sensitized semiconductor solar cells and demonstrates that ultrafast interfacial ET is instrumental for achieving high efficiencies in solar-to-electrical energy conversion. The simulation supports the observation that the ET successfully competes with energy losses due to electron-phonon relaxation. The ET proceeds by the adiabatic mechanism because of strong donor-acceptor coupling. High frequency polar vibrations of both QD and TiO(2) promote the ET, since these modes can rapidly influence the donor-acceptor state energies and coupling. Low frequency vibrations generate a distribution of initial conditions for ET, which shows a broad variety of scenarios at the single-molecule level. Compared to the molecule-TiO(2) interfaces, the QD-TiO(2) system exhibits pronounced differences that arise due to the larger size and higher rigidity of QDs relative to molecules. Both donor and acceptor states are more delocalized in the QD system, and the ET is promoted by optical phonons, which have relatively low frequencies in the QD materials composed of heavy elements. In contrast, in molecular systems, optical phonons are not thermally accessible under ambient conditions. Meanwhile, TiO(2) acceptor states resemble surface impurities due to the local influence of molecular chromophores. At the same time, the photoinduced ET at both QD-TiO(2) and molecule-TiO(2) interfaces is ultrafast and occurs by the adiabatic mechanism, as a result of strong donor-acceptor coupling. The reported state-of-the-art simulation generates a detailed time-domain atomistic description of the interfacial ET process that is fundamental to a wide variety of applications. PMID:22007727
NASA Astrophysics Data System (ADS)
Belyayev, S. T.
2013-06-01
In 1947 I became a second-year student at Moscow State University's Physics and Engineering Department, where a part of the week's classes were taught at base organizations. Our group's base was the future Kurchatov Institute, at that time known as the mysterious "Laboratory N^circ 2," and later as LIPAN.
NASA Astrophysics Data System (ADS)
De Fazio, Dario; de Castro-Vitores, Miguel; Aguado, Alfredo; Aquilanti, Vincenzo; Cavalli, Simonetta
2012-12-01
In this work we critically revise several aspects of previous ab initio quantum chemistry studies [P. Palmieri et al., Mol. Phys. 98, 1835 (2000);, 10.1080/00268970009483387 C. N. Ramachandran et al., Chem. Phys. Lett. 469, 26 (2009)], 10.1016/j.cplett.2008.12.035 of the HeH_2^+ system. New diatomic curves for the H_2^+ and HeH+ molecular ions, which provide vibrational frequencies at a near spectroscopic level of accuracy, have been generated to test the quality of the diatomic terms employed in the previous analytical fittings. The reliability of the global potential energy surfaces has also been tested performing benchmark quantum scattering calculations within the time-independent approach in an extended interval of energies. In particular, the total integral cross sections have been calculated in the total collision energy range 0.955-2.400 eV for the scattering of the He atom by the ortho- and para-hydrogen molecular ion. The energy profiles of the total integral cross sections for selected vibro-rotational states of H_2^+ (v = 0, …,5 and j = 1, …,7) show a strong rotational enhancement for the lower vibrational states which becomes weaker as the vibrational quantum number increases. Comparison with several available experimental data is presented and discussed.
Szalay, Péter G; Holka, Filip; Fremont, Julien; Rey, Michael; Peterson, Kirk A; Tyuterev, Vladimir G
2011-03-01
The aim of the study was to explore the limits of ab initio methods towards the description of excited vibrational levels up to the dissociation limit for molecules having more than two electrons. To this end a high level ab initio potential energy function was constructed for the four-electron LiH molecule in order to accurately predict a complete set of bound vibrational levels corresponding to the electronic ground state. It was composed from: (a) an ab initio non-relativistic potential obtained at the MR-CISD level including size-extensivity corrections and quintuple-sextuple ζ extrapolation of the basis, (b) MVD relativistic corrections obtained at icMR-CISD/cc-pwCV5Z level, and (c) DBOC obtained at the MR-CISD/cc-pwCVTZ level. Finally, the importance of non-adiabatic effects was also tested by using atomic masses in the vibrational kinetic energy operator. The calculated vibrational levels were compared with those obtained from experimental data [J. A. Coxon and C. S. Dickinson, J. Chem. Phys., 2004, 121, 9378]. Our best estimate of the potential curve results in vibrational energies with a RMS deviation of only ∼1 cm(-1) for the entire set of all empirically determined vibrational levels known so far. These results represent a drastic improvement over previous theoretical predictions of vibrational levels of (7)LiH up to dissociation, D(0), which was predicted to be 19,594 cm(-1). PMID:21180724
NASA Astrophysics Data System (ADS)
Szalay, Péter G.; Holka, Filip; Fremont, Julien; Rey, Michael; Tyuterev, Vladimir G.
2011-06-01
The aim of the study was to explore the limits of initio methods towards the description of excited vibrational levels up to the dissociation limit for molecules having more than two electrons. To this end a high level ab initio potential energy function was constructed for the four-electron LiH molecule in order to accurately predict a complete set of bound vibrational levels corresponding to the electronic ground state. It was composed from: a) an ab initio non-relativistic potential obtained at the MR-CISD level including size-extensivity corrections and quintuple-sextuple ζ extrapolation of the basis, b) MVD (Mass-velocity-Darwin) relativistic corrections obtained at icMR-CISD/cc-pwCV5Z level, and c) DBOC (Diagonal Born-Oppenheimer correction) obtained at the MR-CISD/cc-pwCVTZ level. Finally, the importance of non-adiabatic effects was also tested by using atomic masses in the vibrational kinetic energy operator and by calculation of non-adiabatic coupling by ab initio methods. The calculated vibrational levels were compared with those obtained from experimental data [J.A. Coxon and C.S. Dickinson, J. Chem. Phys., 2004, 121, 9378]. Our best estimate of the potential curve results in vibrational energies with a RMS deviation of only ˜1 wn\\ for the entire set of all empirically determined vibrational levels known so far. These results represent a drastic improvement over previous theoretical predictions of vibrational levels of ^7LiH up to dissociation, D_0, which was predicted to be 19594 Cm-1. In addition, rotational levels have also been calculated. The RMS deviation between our ab initio calculations and empirical results by Coxon and Dickinson for rotational spacings Δ E = E(v, J = 1)-E(v, J = 0) over all available vibrational states of ^7LiH from v = 0 to v= 20 is 0.010 wn (with nuclear masses) and 0.006 wn (with atomic masses). Note that for high vibrational states with v > 6 this falls within the uncertainty of the measurements.
Hu Mei; Liu Xinguo; Tan Ruishan; Li Hongzheng; Xu Wenwu
2013-05-07
A new global potential energy surface for the ground electronic state (1{sup 2}A Prime ) of the Ar+H{sub 2}{sup +}{yields}ArH{sup +}+H reaction has been constructed by multi-reference configuration interaction method with Davidson correction and a basis set of aug-cc-pVQZ. Using 6080 ab initio single-point energies of all the regions for the dynamics, a many-body expansion function form has been used to fit these points. The quantum reactive scattering dynamics calculations taking into account the Coriolis coupling (CC) were carried out on the new potential energy surface over a range of collision energies (0.03-1.0 eV). The reaction probabilities and integral cross sections for the title reaction were calculated. The significance of including the CC quantum scattering calculation has been revealed by the comparison between the CC and the centrifugal sudden approximation calculation. The calculated cross section is in agreement with the experimental result at collision energy 1.0 eV.
van Genderen, E.; Clabbers, M. T. B.; Das, P. P.; Stewart, A.; Nederlof, I.; Barentsen, K. C.; Portillo, Q.; Pannu, N. S.; Nicolopoulos, S.; Gruene, T.; Abrahams, J. P.
2016-01-01
Until recently, structure determination by transmission electron microscopy of beam-sensitive three-dimensional nanocrystals required electron diffraction tomography data collection at liquid-nitrogen temperature, in order to reduce radiation damage. Here it is shown that the novel Timepix detector combines a high dynamic range with a very high signal-to-noise ratio and single-electron sensitivity, enabling ab initio phasing of beam-sensitive organic compounds. Low-dose electron diffraction data (∼0.013 e− Å−2 s−1) were collected at room temperature with the rotation method. It was ascertained that the data were of sufficient quality for structure solution using direct methods using software developed for X-ray crystallography (XDS, SHELX) and for electron crystallography (ADT3D/PETS, SIR2014). PMID:26919375
NASA Astrophysics Data System (ADS)
Pandey, Prasenjit; Chakraborty, Tanmoy; Mukherjee, Asok K.
2013-10-01
Ab initio theory at the HF/6-311G(d,p) level has been used to compute the hydrogen bonding thermodynamics in bulk liquid ethanol. Inter-cluster hydrogen bonding is assumed to mimic the H-bonding in bulk ethanol. Rotation of the clusters has been neglected, but translational and vibrational motions are taken into account for calculating bulk thermodynamic parameters. Results are well in agreement with an earlier report [J. Chem. Phys. 116, 4212 (2002)]. For a more accurate dipole moment of monomer, MP2/6-311++G(d,p) calculation was done. Use of the computed thermodynamic data in a statistical model yields the Kirkwood-Frohlich correlation factor and the dielectric constant of ethanol (21.0) close to the experimental value, 24.3 at 298 K.
Erba, Alessandro; Maul, Jefferson; Demichelis, Raffaella; Dovesi, Roberto
2015-05-01
The thermochemical behavior of α-Al2O3 corundum in the whole temperature range 0-2317 K (melting point) and under pressures up to 12 GPa is predicted by applying ab initio methods based on the density functional theory (DFT), the use of a local basis set and periodic-boundary conditions. Thermodynamic properties are treated both within and beyond the harmonic approximation to the lattice potential. In particular, a recent implementation of the quasi-harmonic approximation, in the Crystal program, is here shown to provide a reliable description of the thermal expansion coefficient, entropy, constant-volume and constant-pressure specific heats, and temperature dependence of the bulk modulus, nearly up to the corundum melting temperature. This is a remarkable outcome suggesting α-Al2O3 to be an almost perfect quasi-harmonic crystal. The effect of using different computational parameters and DFT functionals belonging to different levels of approximations on the accuracy of the thermal properties is tested, providing a reference for further studies involving alumina polymorphs and, more generally, quasi-ionic minerals. PMID:25864537
NASA Astrophysics Data System (ADS)
Tiwald, P.; Gräfe, S.; Burgdörfer, J.; Wirtz, L.
2013-12-01
We study the non-adiabatic charge transfer dynamics during the collision of a slow proton with a lithium fluoride surface employing a quantum-chemistry based dynamics approach. The surface is modeled by an Li5F1 + H+ cluster embedded in a large matrix of point charges. Going beyond the adiabatic (or Born-Oppenheimer) approximation, we apply multi-reference configuration-interaction methods that allow for the calculation of ground and excited states of the embedded cluster as well as of the non-adiabatic couplings between them. This information serves as input for the determination of the neutralization probability of a proton scattered off a LiF surface using Tully's semi-classical surface hopping algorithm.
Ab initio infrared and Raman spectra
NASA Astrophysics Data System (ADS)
Fredkin, Donald R.; Komornicki, Andrew; White, Steven R.; Wilson, Kent R.
1983-06-01
We discuss several ways in which molecular absorption and scattering spectra can be computed ab initio, from the fundamental constants of nature. These methods can be divided into two general categories. In the first, or sequential, type of approach, one first solves the electronic part of the Schrödinger equation in the Born-Oppenheimer approximation, mapping out the potential energy, dipole moment vector (for infrared absorption) and polarizability tensor (for Raman scattering) as functions of nuclear coordinates. Having completed the electronic part of the calculation, one then solves the nuclear part of the problem either classically or quantum mechanically. As an example of the sequential ab initio approach, the infrared and Raman rotational and vibrational-rotational spectral band contours for the water molecule are computed in the simplest rigid rotor, normal mode approximation. Quantum techniques are used to calculate the necessary potential energy, dipole moment, and polarizability information at the equilibrium geometry. A new quick, accurate, and easy to program classical technique involving no reference to Euler angles or special functions is developed to compute the infrared and Raman band contours for any rigid rotor, including asymmetric tops. A second, or simultaneous, type of ab initio approach is suggested for large systems, particularly those for which normal mode analysis is inappropriate, such as liquids, clusters, or floppy molecules. Then the curse of dimensionality prevents mapping out in advance the complete potential, dipole moment, and polarizability functions over the whole space of nuclear positions of all atoms, and a solution in which the electronic and nuclear parts of the Born-Oppenheimer approximation are simultaneously solved is needed. A quantum force classical trajectory (QFCT) molecular dynamic method, based on linear response theory, is described, in which the forces, dipole moment, and polarizability are computed quantum
Skerencak-Frech, Andrej; Maiwald, Martin; Trumm, Michael; Froehlich, Daniel R; Panak, Petra J
2015-02-16
The complexation of Cm(III) with oxalate is studied in aqueous solution as a function of the ligand concentration, the ionic strength (NaCl), and the temperature (T = 20–90 °C) by time-resolved laser fluorescence spectroscopy (TRLFS) and quantum chemical calculations. Four complex species ([Cm(Ox)n](3–2n), n = 1, 2, 3, 4) are identified, and their molar fractions are determined by peak deconvolution of the emission spectra. The conditional log K′n(T) values of the first three complexes are calculated and extrapolated to zero ionic strength with the specific ion interaction theory approach. The [Cm(Ox)4](5–) complex forms only at high temperatures. Thus, the log K4(0)(T) value was determined at T > 60 °C. The log K1(0)(25 °C) = 6.86 ± 0.02 decreases by 0.1 logarithmic units in the studied temperature range. The log K2(0)(25 °C) = 4.68 ± 0.09 increases by 0.35, and log K3(0)(25 °C) = 2.11 ± 0.05 increases by 0.37 orders of magnitude. The log Kn(0)(T) (n = 1, 2, 3) values are linearly correlated with the reciprocal temperature. Thus, their temperature dependencies are fitted with the linear Van’t Hoff equation yielding the standard reaction enthalpy (ΔrHm(0)) and standard reaction entropy (ΔrSm(0)) of the stepwise formation of the [Cm(Ox)n](3–2n) species (n = 1, 2, 3). Furthermore, the binary ion–ion interaction coefficients of the four Cm(III) oxalate species with Cl(–)/Na(+) are determined. The binding energies, bond lengths, and bond angles of the different Cm(III) oxalate complexes are calculated in the gas phase as well as in a box containing 1000 H2O molecules by ab inito calculations and molecular dynamics simulations, respectively. PMID:25646935
Kubas, Adam; Blumberger, Jochen; Hoffmann, Felix; Heck, Alexander; Elstner, Marcus; Oberhofer, Harald
2014-03-14
We introduce a database (HAB11) of electronic coupling matrix elements (H{sub ab}) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H{sub ab} values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.
NASA Astrophysics Data System (ADS)
Vaks, V. G.
2013-06-01
I had the good fortune to be a student of A. B. Migdal - AB, as we called him in person or in his absence - and to work in the sector he headed at the Kurchatov Institute, along with his other students and my friends, including Vitya Galitsky, Spartak Belyayev and Tolya Larkin. I was especially close with AB in the second half of the 1950s, the years most important for my formation, and AB's contribution to this formation was very great. To this day, I've often quoted AB on various occasions, as it's hard to put things better or more precisely than he did; I tell friends stories heard from AB, because these stories enhance life as AB himself enhanced it; my daughter is named Tanya after AB's wife Tatyana Lvovna, and so on. In what follows, I'll recount a few episodes in my life in which AB played an important or decisive role, and then will share some other memories of AB...
NASA Astrophysics Data System (ADS)
Thiessen, P. A.; Treder, H.-J.
Jedes initium wird durch experimenta crucis zum eventus. Jedes theoretisch interpretierbare ex-eventu-Resultat führt auf ein neues Initium. Gerade dies ist die gemeinsame Aussage von Atomistik, Quantenmechanik und Relativitätstheorie.Translated AbstractAb initio vel ex eventu. IIEvery initium becomes an eventus by experimenta crucis. Every theoretically interpretable ex-eventu result leads to a new initium. Right this is the joint assertion of atomism, quantum mechanics, and relativity.
NASA Astrophysics Data System (ADS)
Wu, Xiangyang
1999-07-01
The heterocyclic amine 2-amino-3-methylimidazo (4, 5-f) quinoline (IQ) is one of a number of carcinogens found in barbecued meat and fish. It induces tumors in mammals and is probably involved in human carcinogenesis, because of great exposure to such food carcinogens. IQ is biochemically activated to a derivative which reacts with DNA to form a covalent adduct. This adduct may deform the DNA and consequently cause a mutation. which may initiate carcinogenesis. To understand this cancer initiating event, it is necessary to obtain atomic resolution structures of the damaged DNA. No such structures are available experimentally due to synthesis difficulties. Therefore, we employ extensive molecular mechanics and dynamics calculations for this purpose. The major IQ-DNA adduct in the specific DNA sequence d(5'G1G2C G3CCA3') - d(5'TGGCGCC3') with IQ modified at G3 is studied. The d(5'G1G2C G3CC3') sequence has recently been shown to be a hot-spot for mutations when IQ modification is at G3. Although this sequence is prone to -2 deletions via a ``slippage mechanism'' even when unmodified, a key question is why IQ increases the mutation frequency of the unmodified DNA by about 104 fold. Is there a structural feature imposed by IQ that is responsible? The molecular mechanics and dynamics program AMBER for nucleic acids with the latest force field was chosen for this work. This force field has been demonstrated to reproduce well the B-DNA structure. However, some parameters, the partial charges, bond lengths and angles, dihedral parameters of the modified residue, are not available in the AMBER database. We parameterized the force field using high level ab initio quantum calculations. We created 800 starting conformations which uniformly sampled in combination at 18° intervals three torsion angles that govern the IQ-DNA orientations, and energy minimized them. The most important structures are abnormal; the IQ damaged guanine is rotated out of its standard B
Energy Science and Technology Software Center (ESTSC)
2007-05-01
A hierarchical, modular modeling environment for hybrid simulations of sequential-modular, systems dynamics, discrete-event, and agent-based paradigms Omega-AB models contain a hierarchically-defined module tree that specifies the execution logic for the simulation, and a multi-network graph that defines the environment within which the simulation occurs. Modules are the fundamental buildinig blocks of an Omega-AB model and can define anything from a basic mathematical operation to a complex behavioral response model. Modules rely on the "plug-in" conceptmore » which allows developers to build independent module libraries that are gathered, linked, and instantiated by the Omega-AB engine at run time. Inter-module communication occurs through two complimentary systems: pull-based "ports" for general computation patterns and push-based "plugs" for event processing. The simulation environment is an abstract graph of nodes and links. Agents (module sub-trees headed up by an Agent module) reside at nodes and relate to their neighbors through typed links. To facilitate the construction and visualization of complex, interacting networks with dramatically different structure, Omega-AB provides a system for organizing the nodes into hierarchica trees that describe "slices" of the overall network.« less
Ab initio infrared and Raman spectra
NASA Technical Reports Server (NTRS)
Fredkin, D. R.; White, S. R.; Wilson, K. R.; Komornicki, A.
1983-01-01
It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.
Ab Initio Studies of Stratospheric Ozone Depletion Chemistry
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Head-Gordon, Martin; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
An overview of the current understanding of ozone depletion chemistry, particularly with regards the formation of the so-called Antarctic ozone hole, will be presented together with an outline as to how ab initio quantum chemistry can be used to further our understanding of stratospheric chemistry. The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results will be shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.
Ables, E.; Bionta, R.; Olson, H.; Ott, L.; Parker, E.; Wright, D.; Wuest, C
1996-02-01
After investigating a number of materials, we discovered that an ABS plastic doped with a conducting polymer performs well as the resistive electrode in a narrow gap RPC (resistive plate chamber). Operating in the streamer mode, we find efficiencies of 90-96% with low noise and low strip multiplicities. We have also studied a variety of operating gases and found that a mixture containing SF{sub 6}, a non-ozone depleting gas, argon and isobutane gives good streamer mode performance, even with isobutane concentrations of 20% or less.
NASA Astrophysics Data System (ADS)
Abrams, Daniel S.
This thesis describes several new quantum algorithms. These include a polynomial time algorithm that uses a quantum fast Fourier transform to find eigenvalues and eigenvectors of a Hamiltonian operator, and that can be applied in cases (commonly found in ab initio physics and chemistry problems) for which all known classical algorithms require exponential time. Fast algorithms for simulating many body Fermi systems are also provided in both first and second quantized descriptions. An efficient quantum algorithm for anti-symmetrization is given as well as a detailed discussion of a simulation of the Hubbard model. In addition, quantum algorithms that calculate numerical integrals and various characteristics of stochastic processes are described. Two techniques are given, both of which obtain an exponential speed increase in comparison to the fastest known classical deterministic algorithms and a quadratic speed increase in comparison to classical Monte Carlo (probabilistic) methods. I derive a simpler and slightly faster version of Grover's mean algorithm, show how to apply quantum counting to the problem, develop some variations of these algorithms, and show how both (apparently distinct) approaches can be understood from the same unified framework. Finally, the relationship between physics and computation is explored in some more depth, and it is shown that computational complexity theory depends very sensitively on physical laws. In particular, it is shown that nonlinear quantum mechanics allows for the polynomial time solution of NP-complete and #P oracle problems. Using the Weinberg model as a simple example, the explicit construction of the necessary gates is derived from the underlying physics. Nonlinear quantum algorithms are also presented using Polchinski type nonlinearities which do not allow for superluminal communication. (Copies available exclusively from MIT Libraries, Rm. 14- 0551, Cambridge, MA 02139-4307. Ph. 617-253-5668; Fax 617-253-1690.)
Nuclear binding near a quantum phase transition
NASA Astrophysics Data System (ADS)
Lee, Dean
2016-03-01
I review recent ab initio results by the Nuclear Lattice Effective Field Theory Collaboration showing that nature lies close to a quantum phase transition between an alpha-particle gas and nuclear liquid. I discuss the control parameter of this transition and the implications for clustering in nuclei and improving ab initio nuclear structure calculations.
Ab initio determination of light hadron masses.
Dürr, S; Fodor, Z; Frison, J; Hoelbling, C; Hoffmann, R; Katz, S D; Krieg, S; Kurth, T; Lellouch, L; Lippert, T; Szabo, K K; Vulvert, G
2008-11-21
More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab initio calculation of the masses of protons, neutrons, and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega-electron volts are used to extrapolate to the physical point, with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties. PMID:19023076
Ab Initio Modeling of Molecular Radiation
NASA Technical Reports Server (NTRS)
Jaffe, Richard; Schwenke, David
2014-01-01
Radiative emission from excited states of atoms and molecules can comprise a significant fraction of the total heat flux experienced by spacecraft during atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, some of this radiative flux can be absorbed by molecular constituents in the boundary layer that are formed by the ablation process. Ab initio quantum mechanical calculations are carried out to predict the strengths of these emission and absorption processes. This talk will describe the methods used in these calculations using, as examples, the 4th positive emission bands of CO and the 1g+ 1u+ absorption in C3. The results of these calculations are being used as input to NASA radiation modeling codes like NeqAir, HARA and HyperRad.
Ab initio melting curve of osmium
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Burakovsky, N.; Preston, D. L.
2015-11-01
The melting curve of osmium up to a pressure P of 500 GPa is obtained from an extensive suite of ab initio quantum molecular dynamics (QMD) simulations using the Z method. The ab initio P =0 melting point of Os is 3370 ±75 K; this range encompasses all of the available data in the literature and corroborates the conclusion of J. W. Arblaster [Platinum Metals Rev. 49, 166 (2005)], 10.1595/147106705X70264 that the melting temperature of pure Os is 3400 ±50 K and that the 3300 K typically quoted in the literature is the melting point of impure Os. The T =0 equation of state (EOS) of Os and the P dependence of the optimized c /a ratio for the hexagonal unit cell, both to pressures ˜900 GPa, are obtained in the ab initio approach as validation of its use. Although excellent agreement with the available experimental data (P ≲80 GPa) is found, it is the third-order Birch-Murnaghan EOS with B0'=5 rather than the more widely accepted B0'=4 that describes the QMD data to higher pressures, in agreement with the more recent experimental EOS by Godwal et al. The theoretical melting curve of Os obtained earlier by Joshi et al. is shown to be inconsistent with our QMD results, and the possible reason for this discrepancy is suggested. Regularities in the melting curves of Os and five other third-row transition metals (Ta, W, Re, Pt, Au) could be used to estimate the currently unknown melting curves of Hf and Ir.
NASA Astrophysics Data System (ADS)
Tapia, O.; Andres, J.; Aullo, J. M.; Bränden, C.-I.
1985-11-01
The electronic mechanisms of a model hydride transfer reaction are theoretically studied with ab inito RHF and UHF SCF MO procedures at the 4-31G basis set level and analytical gradient methods. The model system describes the reduction of cyclopropenyl cation to cyclopropene by the oxidation of lithium hydride to lithium cation. The molecular fragments corresponding to the asymptotic reactive channels characterizing the stepwise mechanisms currently discussed in the literature have been characterized. The binding energy between the fragments is estimated within a simple electrostatic approximate scheme. The results show that a hydride-ion mechanism is a likely pathway for this particular system. The system is thereafter thoroughly studied from the supermolecule approach. Reaction paths for the ground and first triplet electronic states have been calculated. The hypersurface is explored from a geometrical disposition of the reactants that mimics the one found in several dehydrogenases (perpendicular configuration). A hydride ion is found to be the particle transferred on the unconstrained as well as the constrained reaction pathways in the ground electronic state. In the triplet state (perpendicular configuration) the mechanism is stepwise: electron transfer followed by a hydrogen atom transfer. It has been noticed that the perpendicular geometrical disposition of the reactants plays an important role by polarizing the susceptible cyclopropene C-H bond in the sense of increasing the electronic density at the hydrogen nucleus. This provides a clue to rationalize several dehydrogenase's active site structure and mechanism. The reactant molecular complex found in the inverted potential energy curves, namely the LiH---Cp+ association has an electronic distribution which can be described as a hydride ion cementing two electron deficient centers corresponding to the cyclopropenyl and the lithium cations. Direct CI calculations confirm the overall picture obtained above.
THERMODYNAMICS OF MATERIALS: FROM AB INITIO TO PHENOMENOLOGY
Turchi, P A
2004-09-24
Quantum mechanical-based (or ab initio) methods are used to predict the stability properties of materials although their application is limited to relatively simple systems in terms of structures and number of alloy components. However thermodynamics of complex multi-component alloys requires a more versatile approach afforded within the CALPHAD formalism. Despite its success, the lack of experimental data very often prevents the design of robust thermodynamic databases. After a brief survey of ab initio methodologies and CALPHAD, it will be shown how ab initio electronic structure methods can supplement in two ways CALPHAD for subsequent applications. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics '{acute a} la CALPHAD'. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys. Finally, comments will be made on challenges and future prospects.
Bhattacharyya, Swarnendu Domcke, Wolfgang; Dai, Zuyang
2015-11-21
A diabatic three-sheeted six-dimensional potential-energy surface has been constructed for the ground state and the lowest excited state of the PH{sub 3}{sup +} cation. Coupling terms of Jahn-Teller and pseudo-Jahn-Teller origin up to eighth order had to be included to describe the pronounced anharmonicity of the surface due to multiple conical intersections. The parameters of the diabatic Hamiltonian have been optimized by fitting the eigenvalues of the potential-energy matrix to ab initio data calculated at the CASSCF/MRCI level employing the correlation-consistent triple-ζ basis. The theoretical photoelectron spectrum of phosphine and the non-adiabatic nuclear dynamics of the phosphine cation have been computed by propagating nuclear wave packets with the multiconfiguration time-dependent Hartree method. The theoretical photoelectron bands obtained by Fourier transformation of the autocorrelation function agree well with the experimental results. It is shown that the ultrafast non-radiative decay dynamics of the first excited state of PH{sub 3}{sup +} is dominated by the exceptionally strong Jahn-Teller coupling of the asymmetric bending vibrational mode together with a hyperline of conical intersections with the electronic ground state induced by the umbrella mode. Time-dependent population probabilities have been computed for the three adiabatic electronic states. The non-adiabatic Jahn-Teller dynamics within the excited state takes place within ≈5 fs. Almost 80% of the excited-state population decay to the ground state within about 10 fs. The wave packets become highly complex and delocalized after 20 fs and no further significant transfer of electronic population seems to occur up to 100 fs propagation time.
Bhattacharyya, Swarnendu; Dai, Zuyang; Domcke, Wolfgang
2015-11-21
A diabatic three-sheeted six-dimensional potential-energy surface has been constructed for the ground state and the lowest excited state of the PH3 (+) cation. Coupling terms of Jahn-Teller and pseudo-Jahn-Teller origin up to eighth order had to be included to describe the pronounced anharmonicity of the surface due to multiple conical intersections. The parameters of the diabatic Hamiltonian have been optimized by fitting the eigenvalues of the potential-energy matrix to ab initio data calculated at the CASSCF/MRCI level employing the correlation-consistent triple-ζ basis. The theoretical photoelectron spectrum of phosphine and the non-adiabatic nuclear dynamics of the phosphine cation have been computed by propagating nuclear wave packets with the multiconfiguration time-dependent Hartree method. The theoretical photoelectron bands obtained by Fourier transformation of the autocorrelation function agree well with the experimental results. It is shown that the ultrafast non-radiative decay dynamics of the first excited state of PH3 (+) is dominated by the exceptionally strong Jahn-Teller coupling of the asymmetric bending vibrational mode together with a hyperline of conical intersections with the electronic ground state induced by the umbrella mode. Time-dependent population probabilities have been computed for the three adiabatic electronic states. The non-adiabatic Jahn-Teller dynamics within the excited state takes place within ≈5 fs. Almost 80% of the excited-state population decay to the ground state within about 10 fs. The wave packets become highly complex and delocalized after 20 fs and no further significant transfer of electronic population seems to occur up to 100 fs propagation time. PMID:26590531
NASA Astrophysics Data System (ADS)
Bhattacharyya, Swarnendu; Dai, Zuyang; Domcke, Wolfgang
2015-11-01
A diabatic three-sheeted six-dimensional potential-energy surface has been constructed for the ground state and the lowest excited state of the PH3+ cation. Coupling terms of Jahn-Teller and pseudo-Jahn-Teller origin up to eighth order had to be included to describe the pronounced anharmonicity of the surface due to multiple conical intersections. The parameters of the diabatic Hamiltonian have been optimized by fitting the eigenvalues of the potential-energy matrix to ab initio data calculated at the CASSCF/MRCI level employing the correlation-consistent triple-ζ basis. The theoretical photoelectron spectrum of phosphine and the non-adiabatic nuclear dynamics of the phosphine cation have been computed by propagating nuclear wave packets with the multiconfiguration time-dependent Hartree method. The theoretical photoelectron bands obtained by Fourier transformation of the autocorrelation function agree well with the experimental results. It is shown that the ultrafast non-radiative decay dynamics of the first excited state of PH3+ is dominated by the exceptionally strong Jahn-Teller coupling of the asymmetric bending vibrational mode together with a hyperline of conical intersections with the electronic ground state induced by the umbrella mode. Time-dependent population probabilities have been computed for the three adiabatic electronic states. The non-adiabatic Jahn-Teller dynamics within the excited state takes place within ≈5 fs. Almost 80% of the excited-state population decay to the ground state within about 10 fs. The wave packets become highly complex and delocalized after 20 fs and no further significant transfer of electronic population seems to occur up to 100 fs propagation time.
Quantum Computation and Quantum Information
NASA Astrophysics Data System (ADS)
Nielsen, Michael A.; Chuang, Isaac L.
2010-12-01
Part I. Fundamental Concepts: 1. Introduction and overview; 2. Introduction to quantum mechanics; 3. Introduction to computer science; Part II. Quantum Computation: 4. Quantum circuits; 5. The quantum Fourier transform and its application; 6. Quantum search algorithms; 7. Quantum computers: physical realization; Part III. Quantum Information: 8. Quantum noise and quantum operations; 9. Distance measures for quantum information; 10. Quantum error-correction; 11. Entropy and information; 12. Quantum information theory; Appendices; References; Index.
Kłos, J; Aoiz, F J; Menéndez, M; Brouard, M; Chadwick, H; Eyles, C J
2012-07-01
Adiabatic potential energy surfaces for the ground electronic state of the Xe⋅⋅⋅NO(X(2)Π) van der Waals complex have been calculated using the spin-restricted coupled cluster method with single, double, and non-iterative triple excitations (RCCSD(T)). The scalar relativistic effects present in the Xe atom were included by an effective core potential and we extended the basis with bond functions to improve the description of the dispersion interaction. It has been found that the global minimum on the A(') adiabatic surface occurs at a T-shaped geometry with γ(e) = 94° and R(e) = 7.46 a(0), and with well depth of D(e) = 148.68 cm(-1). There is also an additional local minimum for the collinear geometry Xe-NO with a well depth of 104.5 cm(-1). The adiabat of A('') symmetry exhibits a single minimum at a distance R(e) = 7.68 a(0) and has a skewed geometry with γ(e) = 64° and a well depth of 148.23 cm(-1). Several C(nl) van der Waals dispersion coefficients are also estimated, of which C(6, 0) and C(6, 2) are in a reasonable agreement with previous theoretical results obtained by Nielson et al. [J. Chem. Phys. 64, 2055 (1976)]. The new potential energy surfaces were used to calculate bound states of the complex for total angular momentum quantum numbers up to J = 7/2. The ground state energy of Xe⋅⋅⋅NO(X(2)Π) is D(0) = 117 cm(-1), which matches the experimental value very accurately (within 3.3%). Scattering calculations of integral and differential cross sections have also been performed using fully quantum close coupling calculations and quasi-classical trajectory method at a collision energy of 63 meV. These calculations reveal the important role played by L-type rainbows in the scattering dynamics of the heavier Rg-NO(X) systems. PMID:22779653
Quantum transport calculations using periodic boundaryconditions
Wang, Lin-Wang
2004-06-15
An efficient new method is presented to calculate the quantum transports using periodic boundary conditions. This method allows the use of conventional ground state ab initio programs without big changes. The computational effort is only a few times of a normal groundstate calculations, thus is makes accurate quantum transport calculations for large systems possible.
Ab initio alpha-alpha scattering
NASA Astrophysics Data System (ADS)
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A.; Luu, Thomas; Meißner, Ulf-G.
2015-12-01
Processes such as the scattering of alpha particles (4He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei—nuclei with even and equal numbers of protons and neutrons—is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the ‘adiabatic projection method’ to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
Ab initio alpha-alpha scattering.
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A; Luu, Thomas; Meißner, Ulf-G
2015-12-01
Processes such as the scattering of alpha particles ((4)He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei--nuclei with even and equal numbers of protons and neutrons--is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the 'adiabatic projection method' to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
AB 1725: A Comprehensive Analysis.
ERIC Educational Resources Information Center
California Community Colleges, Sacramento. Board of Governors.
A summary and analysis is provided of California Assembly Bill (AB) 1725, a reform bill that provides new direction and support for the state's community colleges. The analysis addresses each of the eight sections of the bill: (1) mission, highlighting reforms related to mission statements, transfer core curriculum, remedial limits, articulation…
Recent progress in ab initio density matrix renormalization group methodology
NASA Astrophysics Data System (ADS)
Hachmann, Johannes; Dorando, Jonathan J.; Kin-Lic Chan, Garnet
2008-03-01
We present some recent developments in the ab initio density matrix renormalization group (DMRG) method for quantum chemical problems, in particular our local, quadratic scaling algorithm [1] for low dimensional systems. This method is particularly suited for the description of strong nondynamic correlation, and allows us to compute numerically exact (FCI) correlated energies for large active spaces, up to one order of magnitude larger then can be done by conventional CASCI techniques. Other features of this method are its inherent multireference nature, compactness, variational results, size-consistency and size-extensivity. In addition we will review the problems (predominantly organic electronic materials) on which we applied the ab initio DMRG: 1) metal-insulator transition in hydrogen chains [1] 2) all-trans polyacetylene [1] 3) acenes [2] 4) polydiacetylenes [3]. References [1] Hachmann, Cardoen, Chan, JCP 125 (2006), 144101. [2] Hachmann, Dorando, Avil'es, Chan, JCP 127 (2007), 134309. [3] unpublished.
Ab Initio Calculations Of Light-Ion Reactions
Navratil, P; Quaglioni, S; Roth, R; Horiuchi, W
2012-03-12
The exact treatment of nuclei starting from the constituent nucleons and the fundamental interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of nuclear forces, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD employing Hamiltonians constructed within chiral effective field theory. In this contribution, we present one of such promising techniques capable of describing simultaneously both bound and scattering states in light nuclei. By combining the resonating-group method (RGM) with the ab initio no-core shell model (NCSM), we complement a microscopic cluster approach with the use of realistic interactions and a microscopic and consistent description of the clusters. We discuss applications to light nuclei scattering, radiative capture and fusion reactions.
NASA Astrophysics Data System (ADS)
Novoselova, Anna V.; Chernavina, Mariya L.; Berezin, Kirill V.; Berezin, Valentin I.
2015-03-01
This article describes the algorithm and the creation of programs for the input process automate the scaling factors of quantum mechanical force fields calculated in the natural coordinates using ab initio methods and the density functional theory (DFT-methods).
Dispersion Interactions between Rare Gas Atoms: Testing the London Equation Using ab Initio Methods
ERIC Educational Resources Information Center
Halpern, Arthur M.
2011-01-01
A computational chemistry experiment is described in which students can use advanced ab initio quantum mechanical methods to test the ability of the London equation to account quantitatively for the attractive (dispersion) interactions between rare gas atoms. Using readily available electronic structure applications, students can calculate the…
Parallel Quantum Circuit in a Tunnel Junction.
Faizy Namarvar, Omid; Dridi, Ghassen; Joachim, Christian
2016-01-01
Spectral analysis of 1 and 2-states per line quantum bus are normally sufficient to determine the effective Vab(N) electronic coupling between the emitter and receiver states through the bus as a function of the number N of parallel lines. When Vab(N) is difficult to determine, an Heisenberg-Rabi time dependent quantum exchange process must be triggered through the bus to capture the secular oscillation frequency Ωab(N) between those states. Two different linear and regimes are demonstrated for Ωab(N) as a function of N. When the initial preparation is replaced by coupling of the quantum bus to semi-infinite electrodes, the resulting quantum transduction process is not faithfully following the Ωab(N) variations. Because of the electronic transparency normalisation to unity and of the low pass filter character of this transduction, large Ωab(N) cannot be captured by the tunnel junction. The broadly used concept of electrical contact between a metallic nanopad and a molecular device must be better described as a quantum transduction process. At small coupling and when N is small enough not to compensate for this small coupling, an N(2) power law is preserved for Ωab(N) and for Vab(N). PMID:27453262
Parallel Quantum Circuit in a Tunnel Junction
NASA Astrophysics Data System (ADS)
Faizy Namarvar, Omid; Dridi, Ghassen; Joachim, Christian
2016-07-01
Spectral analysis of 1 and 2-states per line quantum bus are normally sufficient to determine the effective Vab(N) electronic coupling between the emitter and receiver states through the bus as a function of the number N of parallel lines. When Vab(N) is difficult to determine, an Heisenberg-Rabi time dependent quantum exchange process must be triggered through the bus to capture the secular oscillation frequency Ωab(N) between those states. Two different linear and regimes are demonstrated for Ωab(N) as a function of N. When the initial preparation is replaced by coupling of the quantum bus to semi-infinite electrodes, the resulting quantum transduction process is not faithfully following the Ωab(N) variations. Because of the electronic transparency normalisation to unity and of the low pass filter character of this transduction, large Ωab(N) cannot be captured by the tunnel junction. The broadly used concept of electrical contact between a metallic nanopad and a molecular device must be better described as a quantum transduction process. At small coupling and when N is small enough not to compensate for this small coupling, an N2 power law is preserved for Ωab(N) and for Vab(N).
Parallel Quantum Circuit in a Tunnel Junction
Faizy Namarvar, Omid; Dridi, Ghassen; Joachim, Christian
2016-01-01
Spectral analysis of 1 and 2-states per line quantum bus are normally sufficient to determine the effective Vab(N) electronic coupling between the emitter and receiver states through the bus as a function of the number N of parallel lines. When Vab(N) is difficult to determine, an Heisenberg-Rabi time dependent quantum exchange process must be triggered through the bus to capture the secular oscillation frequency Ωab(N) between those states. Two different linear and regimes are demonstrated for Ωab(N) as a function of N. When the initial preparation is replaced by coupling of the quantum bus to semi-infinite electrodes, the resulting quantum transduction process is not faithfully following the Ωab(N) variations. Because of the electronic transparency normalisation to unity and of the low pass filter character of this transduction, large Ωab(N) cannot be captured by the tunnel junction. The broadly used concept of electrical contact between a metallic nanopad and a molecular device must be better described as a quantum transduction process. At small coupling and when N is small enough not to compensate for this small coupling, an N2 power law is preserved for Ωab(N) and for Vab(N). PMID:27453262
Michael E. McIlwain; Nick Thompson; Da Gao
2008-03-01
Considerable interest is given to the excellent scintillation properties of cerium doped lanthanum chloride (LaCl3) and lanthanum bromide (LaBr3). The scintillation efficiencies are much greater than other materials, even those containing cerium. This high efficiency is attributed to the high mobility of electrons and holes, unique placement of the cerium 5d states within the band gap, and energy of the band gap. To better understand the scintillation process and better define the nature of the Self Trapped Exciton (STE) within these unique scintillation materials, density functional theory (DFT), and Ab-inito (HF-MP2) calculations are reported. DFT calculations have yielded a qualitative description of the orbital composition and energy distribution of the band structure in the crystalline material. MP2 and single configuration interaction calculations have provided quantitative values for the band gap and provided energies for the possible range of excited states created following hole and electron creation. Based on this theoretical treatment, one possible description of the STE is the combination of Vk center (Br2-1) and LaBr+1 species that recombine to form a distorted geometry LaBr3* (triplet state). Depending on the distance between the LaBr and Br2, the STE emission band can be reproduced.
Michael E. McIlwain; Da Gao; Nick Thompson
2007-12-01
Considerable interest is given to the excellent scintillation properties of cerium doped lanthanum chloride (LaCl3) and lanthanum bromide (LaBr3). The scintillation efficiencies are much greater than other materials, even those containing cerium. This high efficiency is attributed to the high mobility of electrons and holes, unique placement of the cerium 5d states within the band gap, and energy of the band gap. To better understand the scintillation process and better define the nature of the Self Trapped Exciton (STE) within these unique scintillation materials, density functional theory (DFT), and Ab-inito (HF-MP2) calculations are reported. DFT calculations have yielded a qualitative description of the orbital composition and energy distribution of the band structure in the crystalline material. MP2 and single configuration interaction calculations have provided quantitative values for the band gap and provided energies for the possible range of excited states created following hole and electron creation. Based on this theoretical treatment, one possible description of the STE is the combination of Vk center (Br2-1) and LaBr+1 species that recombine to form a distorted geometry LaBr3* (triplet state). Depending on the distance between the LaBr and Br2, the STE emission band can be reproduced.
Ab initio study of neutron drops with chiral Hamiltonians
NASA Astrophysics Data System (ADS)
Potter, H. D.; Fischer, S.; Maris, P.; Vary, J. P.; Binder, S.; Calci, A.; Langhammer, J.; Roth, R.
2014-12-01
We report ab initio calculations for neutron drops in a 10 MeV external harmonic-oscillator trap using chiral nucleon-nucleon plus three-nucleon interactions. We present total binding energies, internal energies, radii and odd-even energy differences for neutron numbers N = 2- 18 using the no-core shell model with and without importance truncation. Furthermore, we present total binding energies for N = 8 , 16 , 20 , 28 , 40 , 50 obtained in a coupled-cluster approach. Comparisons with quantum Monte Carlo results, where available, using Argonne v8‧ with three-nucleon interactions reveal important dependences on the chosen Hamiltonian.
Ab Initio Studies of Halogen and Nitrogen Oxide Species of Interest in Stratospheric Chemistry
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.
Ab initio two-component Ehrenfest dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-01
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H2 and O2. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Ab initio two-component Ehrenfest dynamics
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-21
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H{sub 2} and O{sub 2}. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
The ab-initio density matrix renormalization group in practice
Olivares-Amaya, Roberto; Hu, Weifeng; Sharma, Sandeep; Yang, Jun; Chan, Garnet Kin-Lic; Nakatani, Naoki
2015-01-21
The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.
Hu, Hao; Yang, Weitao
2013-01-01
Determining the free energies and mechanisms of chemical reactions in solution and enzymes is a major challenge. For such complex reaction processes, combined quantum mechanics/molecular mechanics (QM/MM) method is the most effective simulation method to provide an accurate and efficient theoretical description of the molecular system. The computational costs of ab initio QM methods, however, have limited the application of ab initio QM/MM methods. Recent advances in ab initio QM/MM methods allowed the accurate simulation of the free energies for reactions in solution and in enzymes and thus paved the way for broader application of the ab initio QM/MM methods. We review here the theoretical developments and applications of the ab initio QM/MM methods, focusing on the determination of reaction path and the free energies of the reaction processes in solution and enzymes. PMID:24146439
Comparison of DFT and ab initio QM/MM methods for modelling reaction in chorismate synthase
NASA Astrophysics Data System (ADS)
Lawan, Narin; Ranaghan, Kara E.; Manby, Frederick R.; Mulholland, Adrian J.
2014-07-01
Quantum mechanics/molecular mechanics (QM/MM) methods are a popular tool in the investigation of enzyme reactions. Here, we compare B3LYP density functional theory (DFT) and ab initio QM/MM methods for modelling the conversion of 5-enolpyruvylshikimate-3-phosphate to chorismate in chorismate synthase. Good agreement with experimental data is only obtained at the SCS-MP2/CHARMM27 level for a reaction mechanism in which phosphate elimination precedes proton transfer. B3LYP predicts reaction energetics that are qualitatively wrong, stressing the need for ab initio QM/MM methods, and caution in interpretation of DFT results for this enzyme.
AB INITIO AND CALPHAD THERMODYNAMICS OF MATERIALS
Turchi, P A
2004-04-14
Ab initio electronic structure methods can supplement CALPHAD in two major ways for subsequent applications to stability in complex alloys. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics {acute a} la CALPHAD. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys.
Quantum correlation via quantum coherence
NASA Astrophysics Data System (ADS)
Yu, Chang-shui; Zhang, Yang; Zhao, Haiqing
2014-06-01
Quantum correlation includes quantum entanglement and quantum discord. Both entanglement and discord have a common necessary condition—quantum coherence or quantum superposition. In this paper, we attempt to give an alternative understanding of how quantum correlation is related to quantum coherence. We divide the coherence of a quantum state into several classes and find the complete coincidence between geometric (symmetric and asymmetric) quantum discords and some particular classes of quantum coherence. We propose a revised measure for total coherence and find that this measure can lead to a symmetric version of geometric quantum correlation, which is analytic for two qubits. In particular, this measure can also arrive at a monogamy equality on the distribution of quantum coherence. Finally, we also quantify a remaining type of quantum coherence and find that for two qubits, it is directly connected with quantum nonlocality.
Ultrasound Biomicroscopy Comparison of Ab Interno and Ab Externo Intraocular Lens Scleral Fixation.
Horiguchi, Lie; Garcia, Patricia Novita; Malavazzi, Gustavo Ricci; Allemann, Norma; Gomes, Rachel L R
2016-01-01
Purpose. To compare ab interno and ab externo scleral fixation of posterior chamber intraocular lenses (PCIOL) using ultrasound biomicroscopy (UBM). Methods. Randomized patients underwent ab externo or ab interno scleral fixation of a PCIOL. Ultrasound biomicroscopy was performed 3 to 6 months postoperatively, to determine PCIOL centration, IOL distance to the iris at 12, 3, 6, and 9 hours, and haptics placement in relation to the ciliary sulcus. Results. Fifteen patients were enrolled in the study. The ab externo technique was used in 7 eyes (46.6%) and the ab interno in 8 eyes (53.3%). In the ab externo technique, 14 haptics were located: 4 (28.57%) in the ciliary sulcus; 2 (14.28%) anterior to the sulcus; and 8 (57.14%) posterior to the sulcus, 6 in the ciliary body and 2 posterior to the ciliary body. In the ab interno group, 4 haptics (25.0%) were in the ciliary sulcus, 2 (12.50%) anterior to the sulcus, and 10 (75.0%) posterior to the sulcus, 4 in the ciliary body and 6 posterior to the ciliary body. Conclusions. Ab externo and ab interno scleral fixation techniques presented similar results in haptic placement. Ab externo technique presented higher vertical tilt when compared to the ab interno. PMID:27293878
Ultrasound Biomicroscopy Comparison of Ab Interno and Ab Externo Intraocular Lens Scleral Fixation
Horiguchi, Lie; Garcia, Patricia Novita; Malavazzi, Gustavo Ricci; Allemann, Norma
2016-01-01
Purpose. To compare ab interno and ab externo scleral fixation of posterior chamber intraocular lenses (PCIOL) using ultrasound biomicroscopy (UBM). Methods. Randomized patients underwent ab externo or ab interno scleral fixation of a PCIOL. Ultrasound biomicroscopy was performed 3 to 6 months postoperatively, to determine PCIOL centration, IOL distance to the iris at 12, 3, 6, and 9 hours, and haptics placement in relation to the ciliary sulcus. Results. Fifteen patients were enrolled in the study. The ab externo technique was used in 7 eyes (46.6%) and the ab interno in 8 eyes (53.3%). In the ab externo technique, 14 haptics were located: 4 (28.57%) in the ciliary sulcus; 2 (14.28%) anterior to the sulcus; and 8 (57.14%) posterior to the sulcus, 6 in the ciliary body and 2 posterior to the ciliary body. In the ab interno group, 4 haptics (25.0%) were in the ciliary sulcus, 2 (12.50%) anterior to the sulcus, and 10 (75.0%) posterior to the sulcus, 4 in the ciliary body and 6 posterior to the ciliary body. Conclusions. Ab externo and ab interno scleral fixation techniques presented similar results in haptic placement. Ab externo technique presented higher vertical tilt when compared to the ab interno. PMID:27293878
NASA Astrophysics Data System (ADS)
Le Gouët, Jean-Louis; Moiseev, Sergey
2012-06-01
Interaction of quantum radiation with multi-particle ensembles has sparked off intense research efforts during the past decade. Emblematic of this field is the quantum memory scheme, where a quantum state of light is mapped onto an ensemble of atoms and then recovered in its original shape. While opening new access to the basics of light-atom interaction, quantum memory also appears as a key element for information processing applications, such as linear optics quantum computation and long-distance quantum communication via quantum repeaters. Not surprisingly, it is far from trivial to practically recover a stored quantum state of light and, although impressive progress has already been accomplished, researchers are still struggling to reach this ambitious objective. This special issue provides an account of the state-of-the-art in a fast-moving research area that makes physicists, engineers and chemists work together at the forefront of their discipline, involving quantum fields and atoms in different media, magnetic resonance techniques and material science. Various strategies have been considered to store and retrieve quantum light. The explored designs belong to three main—while still overlapping—classes. In architectures derived from photon echo, information is mapped over the spectral components of inhomogeneously broadened absorption bands, such as those encountered in rare earth ion doped crystals and atomic gases in external gradient magnetic field. Protocols based on electromagnetic induced transparency also rely on resonant excitation and are ideally suited to the homogeneous absorption lines offered by laser cooled atomic clouds or ion Coulomb crystals. Finally off-resonance approaches are illustrated by Faraday and Raman processes. Coupling with an optical cavity may enhance the storage process, even for negligibly small atom number. Multiple scattering is also proposed as a way to enlarge the quantum interaction distance of light with matter. The
Electron Transport through Polyene Junctions in between Carbon Nanotubes: an Ab Initio Realization
NASA Astrophysics Data System (ADS)
Chen, Yiing-Rei; Chen, Kai-Yu; Dou, Kun-Peng; Tai, Jung-Shen; Lee, Hsin-Han; Kaun, Chao-Cheng
With both ab initio and tight-binding model calculations, we study a system of polyene bridged armchair carbon nanotube electrodes, considering one-polyene and two-polyene cases, to address aspects of quantum transport through junctions with multiple conjugated molecules. The ab initio results of the two-polyene cases not only show the interference effect in transmission, but also the sensitive dependence of such effect on the combination of relative contact sites, which agrees nicely with the tight-binding model. Moreover, we show that the discrepancy mainly brought by ab initio relaxation provides an insight into the influence upon transmission spectra, from the junction's geometry, bonding and effective potential. This work was supported by the Ministry of Science and Technology of the Republic of China under Grant Nos. 99-2112-M-003-012-MY2 and 103-2622-E-002-031, and the National Center for Theoretical Sciences of Taiwan.
Feet on the potential energy surface, head in the pi clouds
Smith, Quentin
2011-01-01
This work presents explorations of the potential energy surface of clusters of atoms and of the interactions between molecules. First, structures of small aluminum clusters are examined and classified as ground states, transition states, or higher-order saddle points. Subsequently, the focus shifts to dispersion-dominated π-π interactions when the potential energy surfaces of benzene, substituted benzene, and pyridine dimers are explored. Because DNA nucleotide bases can be thought of as substituted heterocycles, a natural extension of the substituted benzene and pyridine investigations is to model paired nucleotide bases. Finally, the success of the dispersion studies inspires the development of an extension to the computational method used, which will enable the dispersion energy to be modeled – and the potential energy surface explored – in additional chemical systems. The effective fragment potential (EFP) method is described, as well as various quantum mechanical methods. An ab inito quantum mechanical study of 13-atom aluminum clusters is described. EFP studies of aromatic dimers are reported in which dispersion energy makes a significant contribution to the attraction between monomers. Theory and code development toward a means of computing dispersion energy in mixed ab inito-EFP systems are described.
Ab initio nuclear structure theory
NASA Astrophysics Data System (ADS)
Negoita, Gianina Alina
Ab initio no core methods have become major tools for understanding the properties of light nuclei based on realistic nucleon-nucleon (NN) and three-nucleon (NNN) interactions. A brief description is provided for the inter-nucleon interactions that fit two-body scattering and bound state data, as well as NNN interactions. Major new progress, including the goal of applying these interactions to solve for properties of nuclei, is limited by convergence issues. That is, with the goal of obtaining high precision solutions of the nuclear many-body Hamiltonian with no core methods (all nucleons treated on the same footing), one needs to proceed to very large basis spaces to achieve a convergence pattern suitable for extrapolation to the exact result. This thesis investigates (1) the similarity renormalization group (SRG) approach to soften the interaction, while preserving its phase shift properties, and (2) adoption of a realistic basis space using Woods-Saxon (WS) single-particle wavefunctions. Both have their advantages and limitations, discussed here. For (1), SRG was demonstrated by applying it to a realistic NN interaction, JISP16, in a harmonic oscillator (HO) representation. The degree of interaction softening achieved through a regulator parameter is examined. For (2), new results are obtained with the realistic JISP16 NN interaction in ab initio calculations of light nuclei 4He, 6He and 12C, using a WS basis optimized to minimize the ground-state energy within the truncated no core shell model. These are numerically-intensive many-body calculations. Finally, to gain insight into the potential for no core investigations of heavier nuclei, an initial investigation was obtained for the odd mass A = 47 - 49 region nuclei straddling 48Ca. The motivation for selecting these nuclei stems from the aim of preparing for nuclear double beta-decay studies of 48Ca. In these heavier systems, phenomenological additions to the realistic NN interaction determined by previous
Ab initio calculations in a uniform magnetic field using periodic supercells
Cai, W; Galli, G
2003-10-21
We present a formulation of ab initio electronic structure calculations in a finite magnetic field, which retains the simplicity and efficiency of techniques widely used in first principles molecular dynamics simulations, based on plane-wave basis sets and Fourier transforms. In addition we discuss results obtained with this method for the energy spectrum of interacting electrons in quantum wells, and for the electronic properties of dense fluid deuterium in a uniform magnetic field.
Kaolin polytypes revisited ab initio.
Mercier, Patrick H J; Le Page, Yvon
2008-04-01
The well known 36 distinguishable transformations between adjacent kaolin layers are split into 20 energetically distinguishable transformations (EDT) and 16 enantiomorphic transformations, hereafter denoted EDT*. For infinitesimal energy contribution of interactions between non-adjacent layers, the lowest-energy models must result from either (a) repeated application of an EDT or (b) alternate application of an EDT and its EDT*. All modeling, quantum input preparation and interpretation was performed with Materials Toolkit, and quantum optimizations with VASP. Kaolinite and dickite are the lowest-energy models at zero temperature and pressure, whereas nacrite and HP-dickite are the lowest-enthalpy models under moderate pressures based on a rough enthalpy/pressure graph built from numbers given in the supplementary tables. Minor temperature dependence of this calculated 0 K graph would explain the bulk of the current observations regarding synthesis, diagenesis and transformation of kaolin minerals. Other stackings that we list have energies so competitive that they might crystallize at ambient pressure. A homometric pair of energetically distinguishable ideal models, one of them for nacrite, is exposed. The printed experimental structure of nacrite correctly corresponds to the stable member of the pair. In our opinion, all recent literature measurements of the free energy of bulk kaolinite are too negative by approximately 15 kJ mol(-1) for some unknown reason. PMID:18369284
NASA Astrophysics Data System (ADS)
Georgescu, I. M.; Ashhab, S.; Nori, Franco
2014-01-01
Simulating quantum mechanics is known to be a difficult computational problem, especially when dealing with large systems. However, this difficulty may be overcome by using some controllable quantum system to study another less controllable or accessible quantum system, i.e., quantum simulation. Quantum simulation promises to have applications in the study of many problems in, e.g., condensed-matter physics, high-energy physics, atomic physics, quantum chemistry, and cosmology. Quantum simulation could be implemented using quantum computers, but also with simpler, analog devices that would require less control, and therefore, would be easier to construct. A number of quantum systems such as neutral atoms, ions, polar molecules, electrons in semiconductors, superconducting circuits, nuclear spins, and photons have been proposed as quantum simulators. This review outlines the main theoretical and experimental aspects of quantum simulation and emphasizes some of the challenges and promises of this fast-growing field.
Ab initio phonon limited transport
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu
We revisit the thermoelectric (TE) transport properties of two champion materials, PbTe and SnSe, using fully first principles methods. In both cases the performance of the material is due to subtle combinations of structural effects, scattering, and phase space reduction. In PbTe anharmonic effects are completely opposite to the predicted quasiharmonic evolution of phonon frequencies and to frequently (and incorrectly) cited extrapolations of experiments. This stabilizes the material at high T, but also tends to enhance its thermal conductivity, in a non linear manner, above 600 Kelvin. This explains why PbTe is in practice limited to room temperature applications. SnSe has recently been shown to be the most efficient TE material in bulk form. This is mainly due to a strongly enhanced carrier concentration and electrical conductivity, after going through a phase transition from 600 to 800 K. We calculate the transport coefficients as well as the defect concentrations ab initio, showing excellent agreement with experiment, and elucidating the origin of the double phase transition as well as the new charge carriers. AH Romero, EKU Gross, MJ Verstraete, and O Hellman PRB 91, 214310 (2015) O. Hellman, IA Abrikosov, and SI Simak, PRB 84 180301 (2011)
NASA Astrophysics Data System (ADS)
Cruz Reyes, Greis Julieth; Franco, Roberto; Silva Valencia, Jereson; Universidad Santo Tomas Collaboration; Universidad Nacional de Colombia Collaboration
Recent research is focused on superlattices arising from optical lattices, which allow a tunable environment. Experimentally bosons present transitions from superfluid to Mott insulator by changing the energy offset in the unit cell [Nat. Commun. 5:5735 (2014)]. Many studies displayed that ground state of spinless boson systems on superlattices present superfluid, Mott insulator and an additional CDW phase created by the energy shift between the sites into the unit cell [Phys. Rev. A 83, 053621 (2011)]. The first confinement methods were magnetic traps, which freezes the spin; with optical lattices the grade of freedom of spin plays an important role. We consider bosons with spin S =1 on a superlattice made by two sites with energy offset per unit cell (AB chain). The Hamiltonian that describes the system is the Bose-Hubbard model with the superlattice potential (W) and the exchange interaction (V) parameters. This model supports CDW, Mott insulator and superfluid phases. For W near to U, with V =0, Mott phase disappears, but for V increasing, a new CDW appears due to the spin interaction, while the half-integer CDW decrease. These results are widely different from spinless boson, where the CDW phases are stables.
GINGA Observations of AB Doradus
NASA Astrophysics Data System (ADS)
Vilhu, O.; Tsuru, T.; Collier Cameron, A.
We report GINGA observations of the pre main sequence star AB Doradus (HD 36705), performed during 8 - 12 January, 1990. Some rotational modulation might be present. four X-ray flares were detected. Three of these events were similar to the EINSTEIN HRI-flare (Vilhu and Linsky, 1987), with decay times around 25 min. The last flare had long rise and decay times (100 min), resembling the EXOSAT flares observed by Collier Cameron et.al. (1988). The mean flare spectrum can be fitted by a thermal bremstrahlung with temperature 5.0 keV, or by a power-law model with photon index 2.2. The 3 upper limit of the Iron line equivalent width in the flare spectrum is 1 keV, weaker than predicted by thermal models. This Iron line anomaly was first discussed in the case of UX Ari by Tsuru et. al. (1989). However, normal equivalent widths can be derived from several EXOSAT spectra of active cool stars (Pallavicini and Tagliaferri, 1990). We discuss the possibility that the continuum from non-thermal electrons (producing also the microwave emission) could occasionally lower the apparent equivalent width. The mechanism works for reasonably low magnetic field strengths and electon power-law indexes. However, a large population of non-thermal electrons is needed (comparable to the thermal one). Stronger magnetic fields could explain the radio emission with less electrons, but then the non-thermal X-ray continuum remains small.
NASA Astrophysics Data System (ADS)
Aubert, B.; Bazan, A.; Boucham, A.; Boutigny, D.; De Bonis, I.; Favier, J.; Gaillard, J.-M.; Jeremie, A.; Karyotakis, Y.; Le Flour, T.; Lees, J. P.; Lieunard, S.; Petitpas, P.; Robbe, P.; Tisserand, V.; Zachariadou, K.; Palano, A.; Chen, G. P.; Chen, J. C.; Qi, N. D.; Rong, G.; Wang, P.; Zhu, Y. S.; Eigen, G.; Reinertsen, P. L.; Stugu, B.; Abbott, B.; Abrams, G. S.; Amerman, L.; Borgland, A. W.; Breon, A. B.; Brown, D. N.; Button-Shafer, J.; Clark, A. R.; Dardin, S.; Day, C.; Dow, S. F.; Fan, Q.; Gaponenko, I.; Gill, M. S.; Goozen, F. R.; Gowdy, S. J.; Gritsan, A.; Groysman, Y.; Hernikl, C.; Jacobsen, R. G.; Jared, R. C.; Kadel, R. W.; Kadyk, J.; Karcher, A.; Kerth, L. T.; Kipnis, I.; Kluth, S.; Kral, J. F.; Lafever, R.; LeClerc, C.; Levi, M. E.; Lewis, S. A.; Lionberger, C.; Liu, T.; Long, M.; Luo, L.; Lynch, G.; Luft, P.; Mandelli, E.; Marino, M.; Marks, K.; Matuk, C.; Meyer, A. B.; Minor, R.; Mokhtarani, A.; Momayezi, M.; Nyman, M.; Oddone, P. J.; Ohnemus, J.; Oshatz, D.; Patton, S.; Pedrali-Noy, M.; Perazzo, A.; Peters, C.; Pope, W.; Pripstein, M.; Quarrie, D. R.; Rasson, J. E.; Roe, N. A.; Romosan, A.; Ronan, M. T.; Shelkov, V. G.; Stone, R.; Strother, P. D.; Telnov, A. V.; von der Lippe, H.; Weber, T. F.; Wenzel, W. A.; Zizka, G.; Bright-Thomas, P. G.; Hawkes, C. M.; Kirk, A.; Knowles, D. J.; O'Neale, S. W.; Watson, A. T.; Watson, N. K.; Deppermann, T.; Koch, H.; Krug, J.; Kunze, M.; Lewandowski, B.; Peters, K.; Schmuecker, H.; Steinke, M.; Andress, J. C.; Barlow, N. R.; Bhimji, W.; Chevalier, N.; Clark, P. J.; Cottingham, W. N.; De Groot, N.; Dyce, N.; Foster, B.; Mass, A.; McFall, J. D.; Wallom, D.; Wilson, F. F.; Abe, K.; Hearty, C.; McKenna, J. A.; Thiessen, D.; Camanzi, B.; Harrison, T. J.; McKemey, A. K.; Tinslay, J.; Antohin, E. I.; Blinov, V. E.; Bukin, A. D.; Bukin, D. A.; Buzykaev, A. R.; Dubrovin, M. S.; Golubev, V. B.; Ivanchenko, V. N.; Kolachev, G. M.; Korol, A. A.; Kravchenko, E. A.; Mikhailov, S. F.; Onuchin, A. P.; Salnikov, A. A.; Serednyakov, S. I.; Skovpen, Yu. I.; Telnov, V. I.; Yushkov, A. N.; Booth, J.; Lankford, A. J.; Mandelkern, M.; Pier, S.; Stoker, D. P.; Zioulas, G.; Ahsan, A.; Arisaka, K.; Buchanan, C.; Chun, S.; Faccini, R.; MacFarlane, D. B.; Prell, S. A.; Rahatlou, Sh.; Raven, G.; Sharma, V.; Burke, S.; Callahan, D.; Campagnari, C.; Dahmes, B.; Hale, D.; Hart, P. A.; Kuznetsova, N.; Kyre, S.; Levy, S. L.; Long, O.; Lu, A.; May, J.; Richman, J. D.; Verkerke, W.; Witherell, M.; Yellin, S.; Beringer, J.; DeWitt, J.; Dorfan, D. E.; Eisner, A. M.; Frey, A.; Grillo, A. A.; Grothe, M.; Heusch, C. A.; Johnson, R. P.; Kroeger, W.; Lockman, W. S.; Pulliam, T.; Rowe, W.; Sadrozinski, H.; Schalk, T.; Schmitz, R. E.; Schumm, B. A.; Seiden, A.; Spencer, E. N.; Turri, M.; Walkowiak, W.; Wilder, M.; Williams, D. C.; Chen, E.; Dubois-Felsmann, G. P.; Dvoretskii, A.; Hanson, J. E.; Hitlin, D. G.; Kolomensky, Yu. G.; Metzler, S.; Oyang, J.; Porter, F. C.; Ryd, A.; Samuel, A.; Weaver, M.; Yang, S.; Zhu, R. Y.; Devmal, S.; Geld, T. L.; Jayatilleke, S.; Jayatilleke, S. M.; Mancinelli, G.; Meadows, B. T.; Sokoloff, M. D.; Bloom, P.; Broomer, B.; Erdos, E.; Fahey, S.; Ford, W. T.; Gaede, F.; van Hoek, W. C.; Johnson, D. R.; Michael, A. K.; Nauenberg, U.; Olivas, A.; Park, H.; Rankin, P.; Roy, J.; Sen, S.; Smith, J. G.; Wagner, D. L.; Blouw, J.; Harton, J. L.; Krishnamurthy, M.; Soffer, A.; Toki, W. H.; Warner, D. W.; Wilson, R. J.; Zhang, J.; Brandt, T.; Brose, J.; Dahlinger, G.; Dickopp, M.; Dubitzky, R. S.; Eckstein, P.; Futterschneider, H.; Kocian, M. L.; Krause, R.; Müller-Pfefferkorn, R.; Schubert, K. R.; Schwierz, R.; Spaan, B.; Wilden, L.; Behr, L.; Bernard, D.; Bonneaud, G. R.; Brochard, F.; Cohen-Tanugi, J.; Ferrag, S.; Fouque, G.; Gastaldi, F.; Matricon, P.; Mora de Freitas, P.; Renard, C.; Roussot, E.; T'Jampens, S.; Thiebaux, C.; Vasileiadis, G.; Verderi, M.; Anjomshoaa, A.; Bernet, R.; Di Lodovico, F.; Muheim, F.; Playfer, S.; Swain, J. E.; Falbo, M.; Bozzi, C.; Dittongo, S.; Folegani, M.; Piemontese, L.; Ramusino, A. C.; Treadwell, E.; Anulli, F.; Baldini-Ferroli, R.; Calcaterra, A.; de Sangro, R.; Falciai, D.; Finocchiaro, G.; Patteri, P.; Peruzzi, I. M.; Piccolo, M.; Xie, Y.; Zallo, A.; Bagnasco, S.; Buzzo, A.; Contri, R.; Crosetti, G.; Fabbricatore, P.; Farinon, S.; Lo Vetere, M.; Macri, M.; Minutoli, S.; Monge, M. R.; Musenich, R.; Pallavicini, M.; Parodi, R.; Passaggio, S.; Pastore, F. C.; Patrignani, C.; Pia, M. G.; Priano, C.; Robutti, E.; Santroni, A.; Bartoldus, R.; Dignan, T.; Hamilton, R.; Mallik, U.; Cochran, J.; Crawley, H. B.; Fischer, P. A.; Lamsa, J.; McKay, R.; Meyer, W. T.; Rosenberg, E. I.; Albert, J. N.; Beigbeder, C.; Benkebil, M.; Breton, D.; Cizeron, R.; Du, S.; Grosdidier, G.; Hast, C.; Höcker, A.; Lacker, H. M.; LePeltier, V.; Lutz, A. M.
2002-02-01
B AB AR, the detector for the SLAC PEP-II asymmetric e +e - B Factory operating at the ϒ(4 S) resonance, was designed to allow comprehensive studies of CP-violation in B-meson decays. Charged particle tracks are measured in a multi-layer silicon vertex tracker surrounded by a cylindrical wire drift chamber. Electromagnetic showers from electrons and photons are detected in an array of CsI crystals located just inside the solenoidal coil of a superconducting magnet. Muons and neutral hadrons are identified by arrays of resistive plate chambers inserted into gaps in the steel flux return of the magnet. Charged hadrons are identified by d E/d x measurements in the tracking detectors and by a ring-imaging Cherenkov detector surrounding the drift chamber. The trigger, data acquisition and data-monitoring systems, VME- and network-based, are controlled by custom-designed online software. Details of the layout and performance of the detector components and their associated electronics and software are presented.
Faster quantum chemistry simulation on fault-tolerant quantum computers
NASA Astrophysics Data System (ADS)
Cody Jones, N.; Whitfield, James D.; McMahon, Peter L.; Yung, Man-Hong; Van Meter, Rodney; Aspuru-Guzik, Alán; Yamamoto, Yoshihisa
2012-11-01
Quantum computers can in principle simulate quantum physics exponentially faster than their classical counterparts, but some technical hurdles remain. We propose methods which substantially improve the performance of a particular form of simulation, ab initio quantum chemistry, on fault-tolerant quantum computers; these methods generalize readily to other quantum simulation problems. Quantum teleportation plays a key role in these improvements and is used extensively as a computing resource. To improve execution time, we examine techniques for constructing arbitrary gates which perform substantially faster than circuits based on the conventional Solovay-Kitaev algorithm (Dawson and Nielsen 2006 Quantum Inform. Comput. 6 81). For a given approximation error ɛ, arbitrary single-qubit gates can be produced fault-tolerantly and using a restricted set of gates in time which is O(log ɛ) or O(log log ɛ) with sufficient parallel preparation of ancillas, constant average depth is possible using a method we call programmable ancilla rotations. Moreover, we construct and analyze efficient implementations of first- and second-quantized simulation algorithms using the fault-tolerant arbitrary gates and other techniques, such as implementing various subroutines in constant time. A specific example we analyze is the ground-state energy calculation for lithium hydride.
NASA Astrophysics Data System (ADS)
Ma, H.; Guo, S.; Hong, X.; Zhou, Y.
2015-05-01
The HJ-1A/B satellite offers free images with high spatial and temporal resolution, which are effective for dynamically monitoring cyanobacteria blooms. However, the HJ-1A/B satellite also receives distorted signals due to the influence of atmosphere. To acquire accurate information about cyanobacteria blooms, atmospheric correction is needed. HJ-1A/B images were atmosphere corrected using the FLAASH atmospheric correction model. Considering the quantum effect within a certain wavelength range, a spectral response function was included in the process. Then the model was used to process HJ-1A/B images, and the NDVI after atmospheric correction was compared with that before correction. The standard deviation improved from 0.13 to 0.158. Results indicate that atmospheric correction effectively reduces the distorted signals. Finally, NDVI was utilized to monitor the cyanobacteria bloom in Donghu Lake. The accuracy was enhanced compared with that before correction.
Ab initio water pair potential with flexible monomers.
Jankowski, Piotr; Murdachaew, Garold; Bukowski, Robert; Akin-Ojo, Omololu; Leforestier, Claude; Szalewicz, Krzysztof
2015-03-26
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects. PMID:25687650
Spectroscopic behavior of bioconjugated quantum dots
NASA Astrophysics Data System (ADS)
Chornokur, G.; Ostapenko, S.; Emirov, Yu; Korsunska, N. E.; Sellers, T.; Phelan, C.
2008-07-01
We report on a short-wavelength, 'blue' spectral shift of the photoluminescence (PL) spectrum in CdSeTe/ZnS core/shell quantum dots (QDs) caused by bioconjugation with several monoclonal cancer-related antibodies (ABs). Scanning PL spectroscopy was performed on samples dried on solid substrates at various temperatures. The influence of the AB chemical origin on the PL spectral shift was observed. The QD-AB conjugation reaction was confirmed using the agarose gel electrophoresis technique. The spectral shift was strongly increased and the process facilitated when the samples were dried above room temperature. The PL spectroscopic mapping revealed a profile of the PL spectral shift across the dried QD-AB spot. A mechanism of the blue shift is attributed to changes in the QD electronic energy levels caused by a local stress applied to the bioconjugated QD.
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, M. A.; Maris, P.; Vary, J. P.; Smith, R.
2015-08-01
Through ab initio approaches in nuclear theory, we may now seek to quantitatively understand the wealth of nuclear collective phenomena starting from the underlying internucleon interactions. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments and electromagnetic transitions. In this review, NCCI calculations of 7-9Be are used to illustrate and explore ab initio rotational structure, and the resulting predictions for rotational band properties are compared with experiment. We highlight the robustness of ab initio rotational predictions across different choices for the internucleon interaction.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
Elenewski, Justin E.; Hackett, John C
2015-01-01
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis. PMID:25681906
4He Thermophysical Properties: New Ab Initio Calculations
Hurly, John J.; Mehl, James B.
2007-01-01
Since 2000, atomic physicists have reduced the uncertainty of the helium-helium “ab initio” potential; for example, from approximately 0.6 % to 0.1 % at 4 bohr, and from 0.8 % to 0.1 % at 5.6 bohr. These results led us to: (1) construct a new inter-atomic potential ϕ07, (2) recalculate values of the second virial coefficient, the viscosity, and the thermal conductivity of 4He from 1 K to 10,000 K, and (3), analyze the uncertainties of the thermophysical properties that propagate from the uncertainty of ϕ07 and from the Born-Oppenheimer approximation of the electron-nucleon quantum mechanical system. We correct minor errors in a previous publication [J. J. Hurly and M. R. Moldover, J. Res. Nat. Inst. Standards Technol. 105, 667 (2000)] and compare our results with selected data published after 2000. The ab initio results tabulated here can serve as standards for the measurement of thermophysical properties. PMID:27110456
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
NASA Astrophysics Data System (ADS)
Elenewski, Justin E.; Hackett, John C.
2015-02-01
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
Elenewski, Justin E.; Hackett, John C
2015-02-14
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
Unified ab initio approaches to nuclear structure and reactions
NASA Astrophysics Data System (ADS)
Navrátil, Petr; Quaglioni, Sofia; Hupin, Guillaume; Romero-Redondo, Carolina; Calci, Angelo
2016-05-01
The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD-employing Hamiltonians constructed within chiral effective field theory. After a brief overview of the field, we focus on ab initio many-body approaches—built upon the no-core shell model—that are capable of simultaneously describing both bound and scattering nuclear states, and present results for resonances in light nuclei, reactions important for astrophysics and fusion research. In particular, we review recent calculations of resonances in the 6He halo nucleus, of five- and six-nucleon scattering, and an investigation of the role of chiral three-nucleon interactions in the structure of 9Be. Further, we discuss applications to the 7Be {({{p}},γ )}8{{B}} radiative capture. Finally, we highlight our efforts to describe transfer reactions including the 3H{({{d}},{{n}})}4He fusion.
Barrett, B R; Navratil, P; Vary, J P
2011-04-11
A long-standing goal of nuclear theory is to determine the properties of atomic nuclei based on the fundamental interactions among the protons and neutrons (i.e., nucleons). By adopting nucleon-nucleon (NN), three-nucleon (NNN) and higher-nucleon interactions determined from either meson-exchange theory or QCD, with couplings fixed by few-body systems, we preserve the predictive power of nuclear theory. This foundation enables tests of nature's fundamental symmetries and offers new vistas for the full range of complex nuclear phenomena. Basic questions that drive our quest for a microscopic predictive theory of nuclear phenomena include: (1) What controls nuclear saturation; (2) How the nuclear shell model emerges from the underlying theory; (3) What are the properties of nuclei with extreme neutron/proton ratios; (4) Can we predict useful cross sections that cannot be measured; (5) Can nuclei provide precision tests of the fundamental laws of nature; and (6) Under what conditions do we need QCD to describe nuclear structure, among others. Along with other ab initio nuclear theory groups, we have pursued these questions with meson-theoretical NN interactions, such as CD-Bonn and Argonne V18, that were tuned to provide high-quality descriptions of the NN scattering phase shifts and deuteron properties. We then add meson-theoretic NNN interactions such as the Tucson-Melbourne or Urbana IX interactions. More recently, we have adopted realistic NN and NNN interactions with ties to QCD. Chiral perturbation theory within effective field theory ({chi}EFT) provides us with a promising bridge between QCD and hadronic systems. In this approach one works consistently with systems of increasing nucleon number and makes use of the explicit and spontaneous breaking of chiral symmetry to expand the strong interaction in terms of a dimensionless constant, the ratio of a generic small momentum divided by the chiral symmetry breaking scale taken to be about 1 GeV/c. The resulting NN
Quantum networks reveal quantum nonlocality.
Cavalcanti, Daniel; Almeida, Mafalda L; Scarani, Valerio; Acín, Antonio
2011-01-01
The results of local measurements on some composite quantum systems cannot be reproduced classically. This impossibility, known as quantum nonlocality, represents a milestone in the foundations of quantum theory. Quantum nonlocality is also a valuable resource for information-processing tasks, for example, quantum communication, quantum key distribution, quantum state estimation or randomness extraction. Still, deciding whether a quantum state is nonlocal remains a challenging problem. Here, we introduce a novel approach to this question: we study the nonlocal properties of quantum states when distributed and measured in networks. We show, using our framework, how any one-way entanglement distillable state leads to nonlocal correlations and prove that quantum nonlocality is a non-additive resource, which can be activated. There exist states, local at the single-copy level, that become nonlocal when taking several copies of them. Our results imply that the nonlocality of quantum states strongly depends on the measurement context. PMID:21304513
Zabidi, Noriza Ahmad; Kassim, Hasan Abu; Shrivastava, Keshav N.
2008-05-20
Polonium is the only element with a simple cubic (sc) crystal structure. Atoms in solid polonium sit at the corners of a simple cubic unit cell and no where else. Polonium has a valence electron configuration 6s{sup 2}6p{sup 4} (Z = 84). The low temperature {alpha}-phase transforms into the rhombohedral (trigonal) {beta} structure at {approx}348 K. The sc {alpha}-Po unit cell constant is a = 3.345 A. The beta form of polonium ({beta}-Po) has the lattice parameters, a{sub R} = 3.359 A and a rhombohedral angle 98 deg. 13'. We have performed an ab initio electronic structure calculation by using the density functional theory. We have performed the calculation with and without spin-orbit (SO) coupling by using both the LDA and the GGA for the exchange-correlations. The k-points in a simple cubic BZ are determined by R (0.5, 0.5, 0.5), {gamma} (0, 0, 0), X (0.5, 0, 0), M (0.5, 0.5, 0) and {gamma} (0, 0, 0). Other directions of k-points are {gamma} (0, 0, 0), X (0.5, 0, 0), R (0.5, 0.5, 0.5) and {gamma} (0, 0, 0). The SO splittings of p states at the {gamma} point in the GGA+SO scheme for {alpha}-Po are 0.04 eV and 0.02 eV while for the {beta}-Po these are 0.03 eV and 0.97 eV. We have also calculated the vibrational spectra for the unit cells in both the structures. We find that exchanging of a Po atom by Pb atom produces several more bands and destabilizes the {beta} phase.
Stapp, H.P.
1988-12-01
Quantum ontologies are conceptions of the constitution of the universe that are compatible with quantum theory. The ontological orientation is contrasted to the pragmatic orientation of science, and reasons are given for considering quantum ontologies both within science, and in broader contexts. The principal quantum ontologies are described and evaluated. Invited paper at conference: Bell's Theorem, Quantum Theory, and Conceptions of the Universe, George Mason University, October 20-21, 1988. 16 refs.
Quantum Computer Games: Quantum Minesweeper
ERIC Educational Resources Information Center
Gordon, Michal; Gordon, Goren
2010-01-01
The computer game of quantum minesweeper is introduced as a quantum extension of the well-known classical minesweeper. Its main objective is to teach the unique concepts of quantum mechanics in a fun way. Quantum minesweeper demonstrates the effects of superposition, entanglement and their non-local characteristics. While in the classical…
Exact quantum cross sections for a three dimensional angle dependent model for three body reactions.
NASA Technical Reports Server (NTRS)
Baer, M.; Kouri, D. J.
1971-01-01
Exact quantum mechanical reactive cross sections are reported for a three dimensional angle dependent model surface. The surface simulates an atom-heteronuclear diatom system A + BC leading to AB + C where atom B is much heavier than A or C. The molecules BC and AB are taken to be rotating vibrators which can dissociate. Results for two angle dependent potentials are given.
NASA Astrophysics Data System (ADS)
Semenov, Andrew G.; Zaikin, Andrei D.
2010-01-01
We investigate the effect of electron-electron interactions on Aharonov-Bohm (AB) current oscillations in nanorings formed by a chain of metallic quantum dots. We demonstrate that electron-electron interactions cause electron dephasing thereby suppressing the amplitude of AB oscillations at all temperatures down to T=0. The crossover between thermal and quantum dephasing is found to be controlled by the ring perimeter. Our predictions can be directly tested in future experiments.
NASA Astrophysics Data System (ADS)
Pfeiffer, P.; Egusquiza, I. L.; di Ventra, M.; Sanz, M.; Solano, E.
2016-07-01
Technology based on memristors, resistors with memory whose resistance depends on the history of the crossing charges, has lately enhanced the classical paradigm of computation with neuromorphic architectures. However, in contrast to the known quantized models of passive circuit elements, such as inductors, capacitors or resistors, the design and realization of a quantum memristor is still missing. Here, we introduce the concept of a quantum memristor as a quantum dissipative device, whose decoherence mechanism is controlled by a continuous-measurement feedback scheme, which accounts for the memory. Indeed, we provide numerical simulations showing that memory effects actually persist in the quantum regime. Our quantization method, specifically designed for superconducting circuits, may be extended to other quantum platforms, allowing for memristor-type constructions in different quantum technologies. The proposed quantum memristor is then a building block for neuromorphic quantum computation and quantum simulations of non-Markovian systems.
Pfeiffer, P; Egusquiza, I L; Di Ventra, M; Sanz, M; Solano, E
2016-01-01
Technology based on memristors, resistors with memory whose resistance depends on the history of the crossing charges, has lately enhanced the classical paradigm of computation with neuromorphic architectures. However, in contrast to the known quantized models of passive circuit elements, such as inductors, capacitors or resistors, the design and realization of a quantum memristor is still missing. Here, we introduce the concept of a quantum memristor as a quantum dissipative device, whose decoherence mechanism is controlled by a continuous-measurement feedback scheme, which accounts for the memory. Indeed, we provide numerical simulations showing that memory effects actually persist in the quantum regime. Our quantization method, specifically designed for superconducting circuits, may be extended to other quantum platforms, allowing for memristor-type constructions in different quantum technologies. The proposed quantum memristor is then a building block for neuromorphic quantum computation and quantum simulations of non-Markovian systems. PMID:27381511
Pfeiffer, P.; Egusquiza, I. L.; Di Ventra, M.; Sanz, M.; Solano, E.
2016-01-01
Technology based on memristors, resistors with memory whose resistance depends on the history of the crossing charges, has lately enhanced the classical paradigm of computation with neuromorphic architectures. However, in contrast to the known quantized models of passive circuit elements, such as inductors, capacitors or resistors, the design and realization of a quantum memristor is still missing. Here, we introduce the concept of a quantum memristor as a quantum dissipative device, whose decoherence mechanism is controlled by a continuous-measurement feedback scheme, which accounts for the memory. Indeed, we provide numerical simulations showing that memory effects actually persist in the quantum regime. Our quantization method, specifically designed for superconducting circuits, may be extended to other quantum platforms, allowing for memristor-type constructions in different quantum technologies. The proposed quantum memristor is then a building block for neuromorphic quantum computation and quantum simulations of non-Markovian systems. PMID:27381511
Free Quantum Field Theory from Quantum Cellular Automata
NASA Astrophysics Data System (ADS)
Bisio, Alessandro; D'Ariano, Giacomo Mauro; Perinotti, Paolo; Tosini, Alessandro
2015-10-01
After leading to a new axiomatic derivation of quantum theory (see D'Ariano et al. in Found Phys, 2015), the new informational paradigm is entering the domain of quantum field theory, suggesting a quantum automata framework that can be regarded as an extension of quantum field theory to including an hypothetical Planck scale, and with the usual quantum field theory recovered in the relativistic limit of small wave-vectors. Being derived from simple principles (linearity, unitarity, locality, homogeneity, isotropy, and minimality of dimension), the automata theory is quantum ab-initio, and does not assume Lorentz covariance and mechanical notions. Being discrete it can describe localized states and measurements (unmanageable by quantum field theory), solving all the issues plaguing field theory originated from the continuum. These features make the theory an ideal framework for quantum gravity, with relativistic covariance and space-time emergent solely from the interactions, and not assumed a priori. The paper presents a synthetic derivation of the automata theory, showing how the principles lead to a description in terms of a quantum automaton over a Cayley graph of a group. Restricting to Abelian groups we show how the automata recover the Weyl, Dirac and Maxwell dynamics in the relativistic limit. We conclude with some new routes about the more general scenario of non-Abelian Cayley graphs. The phenomenology arising from the automata theory in the ultra-relativistic domain and the analysis of corresponding distorted Lorentz covariance is reviewed in Bisio et al. (Found Phys 2015, in this same issue).
On the Origin of Quantum Diffusion Coefficient and Quantum Potential
NASA Astrophysics Data System (ADS)
Gupta, Aseem
2016-03-01
Synchronizability of space and time experiences between different inhabitants of a spacetime is abstracted as a fundamental premise of Classical physics. Absence thereof i.e. desynchronization between space and time experiences of a system under study and the observer is then studied for a single dimension single particle system. Desynchronization fundamentally makes probability concepts enter physics ab-initio and not as secondary tools to deal with situations wherein incomplete information in situation following perfectly deterministic dynamics demands its introduction. Desynchronization model based on Poisson distribution of events vis-à-vis an observer, leads to expectation of particle's motion as a Brownian motion deriving Nelson's quantum diffusion coefficient naturally, without needing to postulate it. This model also incorporates physical effects akin to those of Bohm's Quantum Potential, again without needing any sub-quantum medium. Schrodinger's equation is shown to be derivable incorporating desynchronization only of space while Quantum Field Theory is shown to model desynchronization of time as well. Fundamental suggestion of the study is that it is desynchronization that is at the root of quantum phenomena rather than sub-micro scales of spacetime. Absence of possibility of synchronization between system's space and time and those of observer is studied. Mathematical modeling of desynchronized evolution explains some intriguing aspects of Quantum Mechanical theory.
Accelerating Ab Initio Path Integral Simulations via Imaginary Multiple-Timestepping.
Cheng, Xiaolu; Herr, Jonathan D; Steele, Ryan P
2016-04-12
This work investigates the use of multiple-timestep schemes in imaginary time for computationally efficient ab initio equilibrium path integral simulations of quantum molecular motion. In the simplest formulation, only every n(th) path integral replica is computed at the target level of electronic structure theory, whereas the remaining low-level replicas still account for nuclear motion quantum effects with a more computationally economical theory. Motivated by recent developments for multiple-timestep techniques in real-time classical molecular dynamics, both 1-electron (atomic-orbital basis set) and 2-electron (electron correlation) truncations are shown to be effective. Structural distributions and thermodynamic averages are tested for representative analytic potentials and ab initio molecular examples. Target quantum chemistry methods include density functional theory and second-order Møller-Plesset perturbation theory, although any level of theory is formally amenable to this framework. For a standard two-level splitting, computational speedups of 1.6-4.0x are observed when using a 4-fold reduction in time slices; an 8-fold reduction is feasible in some cases. Multitiered options further reduce computational requirements and suggest that quantum mechanical motion could potentially be obtained at a cost not significantly different from the cost of classical simulations. PMID:26966920
Magnetic edge states in Aharonov-Bohm graphene quantum rings
NASA Astrophysics Data System (ADS)
Farghadan, R.; Saffarzadeh, A.; Heidari Semiromi, E.
2013-12-01
The effect of electron-electron interaction on the electronic structure of Aharonov-Bohm (AB) graphene quantum rings (GQRs) is explored theoretically using the single-band tight-binding Hamiltonian and the mean-field Hubbard model. The electronic states and magnetic properties of hexagonal, triangular, and circular GQRs with different sizes and zigzag edge terminations are studied. The results show that, although the AB oscillations in the all types of nanoring are affected by the interaction, the spin splitting in the AB oscillations strongly depends on the geometry and the size of graphene nanorings. We found that the total spin of hexagonal and circular rings is zero and therefore, no spin splitting can be observed in the AB oscillations. However, the non-zero magnetization of the triangular rings breaks the degeneracy between spin-up and spin-down electrons, which produces spin-polarized AB oscillations.
NASA Astrophysics Data System (ADS)
Khaliullin, Rustam Z.; Eshet, Hagai; Kühne, Thomas D.; Behler, Jörg; Parrinello, Michele
2010-03-01
An interatomic potential for the diamond and graphite phases of carbon has been created using a neural-network (NN) representation of the ab initio potential energy surface. The NN potential combines the accuracy of a first-principles description of both phases with the efficiency of empirical force fields and allows one to perform a molecular-dynamics study, of ab initio quality, of the thermodynamics of graphite-diamond coexistence. Good agreement between the experimental and calculated coexistence curves is achieved if nuclear quantum effects are included in the simulation.
Quantum robots and quantum computers
Benioff, P.
1998-07-01
Validation of a presumably universal theory, such as quantum mechanics, requires a quantum mechanical description of systems that carry out theoretical calculations and systems that carry out experiments. The description of quantum computers is under active development. No description of systems to carry out experiments has been given. A small step in this direction is taken here by giving a description of quantum robots as mobile systems with on board quantum computers that interact with different environments. Some properties of these systems are discussed. A specific model based on the literature descriptions of quantum Turing machines is presented.
Quantum hair and quantum gravity
Coleman, S. ); Krauss, L.M. ); Preskill, J. ); Wilczek, F. )
1992-01-01
A black hole may carry quantum numbers that are not associated with massless gauge fields, contrary to the spirit of the 'no-hair' theorems. The 'quantum hair' is invisible in the classical limit, but measurable via quantum interference experiments. Quantum hair alters the temperature of the radiation emitted by a black hole. It also induces non-zero expectation values for fields outside the event horizon; these expectation values are non-perturbative in [Dirac h], and decay exponentially far from the hole. The existence of quantum hair demonstrates that a black hole can have an intricate quantum-mechanical structure that is completely missed by standard semiclassical theory.
Zurek, Wojciech H
2008-01-01
Quantum Darwinism - proliferation, in the environment, of multiple records of selected states of the system (its information-theoretic progeny) - explains how quantum fragility of individual state can lead to classical robustness of their multitude.
NASA Astrophysics Data System (ADS)
Harju, Antti J.
2016-06-01
This is a study of orbifold-quotients of quantum groups (quantum orbifolds {Θ } rightrightarrows Gq). These structures have been studied extensively in the case of the quantum S U 2 group. A generalized theory of quantum orbifolds over compact simple and simply connected quantum groups is developed. Associated with a quantum orbifold there is an invariant subalgebra and a crossed product algebra. For each spin quantum orbifold, there is a unitary equivalence class of Dirac spectral triples over the invariant subalgebra, and for each effective spin quantum orbifold associated with a finite group action, there is a unitary equivalence class of Dirac spectral triples over the crossed product algebra. A Hopf-equivariant Fredholm index problem is studied as an application.
Pfeiffer, P.; Egusquiza, I. L.; Di Ventra, M.; Sanz, M.; Solano, E.
2016-07-06
Technology based on memristors, resistors with memory whose resistance depends on the history of the crossing charges, has lately enhanced the classical paradigm of computation with neuromorphic architectures. However, in contrast to the known quantized models of passive circuit elements, such as inductors, capacitors or resistors, the design and realization of a quantum memristor is still missing. Here, we introduce the concept of a quantum memristor as a quantum dissipative device, whose decoherence mechanism is controlled by a continuous-measurement feedback scheme, which accounts for the memory. Indeed, we provide numerical simulations showing that memory effects actually persist in the quantummore » regime. Our quantization method, specifically designed for superconducting circuits, may be extended to other quantum platforms, allowing for memristor-type constructions in different quantum technologies. As a result, the proposed quantum memristor is then a building block for neuromorphic quantum computation and quantum simulations of non-Markovian systems.« less
Ab interno trabeculectomy: patient selection and perspectives.
Vinod, Kateki; Gedde, Steven J
2016-01-01
Ab interno trabeculectomy is one among several recently introduced minimally invasive glaucoma surgeries that avoid a conjunctival incision and full-thickness sclerostomy involved in traditional glaucoma surgery. Ablation of the trabecular meshwork and inner wall of Schlemm's canal is performed in an arcuate fashion via a clear corneal incision, alone or in combination with phacoemulsification cataract surgery. Intraocular pressure reduction following ab interno trabeculectomy is limited by resistance in distal outflow pathways and generally stabilizes in the mid-to-high teens. Relief of medication burden has been demonstrated by some studies. A very low rate of complications, most commonly transient hyphema and intraocular pressure elevations in the immediate postoperative period, have been reported. However, available data are derived from small retrospective and prospective case series. Randomized, controlled trials are needed to better elucidate the potential merits of ab interno trabeculectomy in the combined setting versus phacoemulsification cataract surgery alone and to compare it with other minimally invasive glaucoma surgeries. PMID:27574396
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of chlorine oxide and nitrogen oxide species will be demonstrated by presentation of some example studies. In particular the geometrical structures, vibrational spectra, and heats of formation Of ClNO2, CisClONO, and trans-ClONO are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the ab initio results are shown to fill in the gaps and to resolve the experimental controversy. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of ClONO2, HONO2, ClOOC17 ClOOH, and HOOH will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of the experimental studies.
CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule.
Barletta, Paolo; Shirin, Sergei V; Zobov, Nikolai F; Polyansky, Oleg L; Tennyson, Jonathan; Valeev, Edward F; Császár, Attila G
2006-11-28
The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed. PMID:17144700
MELCOR 1.8.2 assessment: Aerosol experiments ABCOVE AB5, AB6, AB7, and LACE LA2
Souto, F.J.; Haskin, F.E.; Kmetyk, L.N.
1994-10-01
The MELCOR computer code has been used to model four of the large-scale aerosol behavior experiments conducted in the Containment System Test Facility (CSTF) vessel. Tests AB5, AB6 and AB7 of the ABCOVE program simulate the dry aerosol conditions during a hypothetical severe accident in an LMFBR. Test LA2 of the LACE program simulates aerosol behavior in a condensing steam environment during a postulated severe accident in an LWR with failure to isolate the containment. The comparison of code results to experimental data show that MELCOR is able to correctly predict most of the thermal-hydraulic results in the four tests. MELCOR predicts reasonably well the dry aerosol behavior of the ABCOVE tests, but significant disagreements are found in the aerosol behavior modelling for the LA2 experiment. These results tend to support some of the concerns about the MELCOR modelling of steam condensation onto aerosols expressed in previous works. During these analyses, a limitation in the MELCOR input was detected for the specification of the aerosol parameters for more than one component. A Latin Hypercube Sampling (LHS) sensitivity study of the aerosol dynamic constants is presented for test AB6. The study shows the importance of the aerosol shape factors in the aerosol deposition behavior, and reveals that MELCOR input/output processing is highly labor intensive for uncertainty and sensitivity analyses based on LHS.
NASA Astrophysics Data System (ADS)
Moulick, Subhayan Roy; Panigrahi, Prasanta K.
2016-06-01
We propose the idea of a quantum cheque scheme, a cryptographic protocol in which any legitimate client of a trusted bank can issue a cheque, that cannot be counterfeited or altered in anyway, and can be verified by a bank or any of its branches. We formally define a quantum cheque and present the first unconditionally secure quantum cheque scheme and show it to be secure against any no-signalling adversary. The proposed quantum cheque scheme can been perceived as the quantum analog of Electronic Data Interchange, as an alternate for current e-Payment Gateways.
NASA Astrophysics Data System (ADS)
Moulick, Subhayan Roy; Panigrahi, Prasanta K.
2016-03-01
We propose the idea of a quantum cheque scheme, a cryptographic protocol in which any legitimate client of a trusted bank can issue a cheque, that cannot be counterfeited or altered in anyway, and can be verified by a bank or any of its branches. We formally define a quantum cheque and present the first unconditionally secure quantum cheque scheme and show it to be secure against any no-signalling adversary. The proposed quantum cheque scheme can been perceived as the quantum analog of Electronic Data Interchange, as an alternate for current e-Payment Gateways.
NASA Astrophysics Data System (ADS)
Brown, Matthew J.
2014-02-01
The framework of quantum frames can help unravel some of the interpretive difficulties i the foundation of quantum mechanics. In this paper, I begin by tracing the origins of this concept in Bohr's discussion of quantum theory and his theory of complementarity. Engaging with various interpreters and followers of Bohr, I argue that the correct account of quantum frames must be extended beyond literal space-time reference frames to frames defined by relations between a quantum system and the exosystem or external physical frame, of which measurement contexts are a particularly important example. This approach provides superior solutions to key EPR-type measurement and locality paradoxes.
Characterization of mAb dimers reveals predominant dimer forms common in therapeutic mAbs.
Plath, Friederike; Ringler, Philippe; Graff-Meyer, Alexandra; Stahlberg, Henning; Lauer, Matthias E; Rufer, Arne C; Graewert, Melissa A; Svergun, Dmitri; Gellermann, Gerald; Finkler, Christof; Stracke, Jan O; Koulov, Atanas; Schnaible, Volker
2016-07-01
The formation of undesired high molecular weight species such as dimers is an important quality attribute for therapeutic monoclonal antibody formulations. Therefore, the thorough understanding of mAb dimerization and the detailed characterization mAb dimers is of great interest for future pharmaceutical development of therapeutic antibodies. In this work, we focused on the analyses of different mAb dimers regarding size, surface properties, chemical identity, overall structure and localization of possible dimerization sites. Dimer fractions of different mAbs were isolated to a satisfactory purity from bulk material and revealed 2 predominant overall structures, namely elongated and compact dimer forms. The elongated dimers displayed one dimerization site involving the tip of the Fab domain. Depending on the stress applied, these elongated dimers are connected either covalently or non-covalently. In contrast, the compact dimers exhibited non-covalent association. Several interaction points were detected for the compact dimers involving the hinge region or the base of the Fab domain. These results indicate that mAb dimer fractions are rather complex and may contain more than one kind of dimer. Nevertheless, the overall appearance of mAb dimers suggests the existence of 2 predominant dimeric structures, elongated and compact, which are commonly present in preparations of therapeutic mAbs. PMID:27031922
Ab initio methods for nuclear properties - a computational physics approach
NASA Astrophysics Data System (ADS)
Maris, Pieter
2011-04-01
A microscopic theory for the structure and reactions of light nuclei poses formidable challenges for high-performance computing. Several ab-initio methods have now emerged that provide nearly exact solutions for some nuclear properties. The ab-initio no-core full configuration (NCFC) approach is based on basis space expansion methods and uses Slater determinants of single-nucleon basis functions to express the nuclear wave function. In this approach, the quantum many-particle problem becomes a large sparse matrix eigenvalue problem. The eigenvalues of this matrix give us the binding energies, and the corresponding eigenvectors the nuclear wave functions. These wave functions can be employed to evaluate experimental quantities. In order to reach numerical convergence for fundamental problems of interest, the matrix dimension often exceeds 1 billion, and the number of nonzero matrix elements may saturate available storage on present-day leadership class facilities. I discuss different strategies for distributing and solving this large sparse matrix on current multicore computer architectures, including methods to deal with with memory bottleneck. Several of these strategies have been implemented in the code MFDn, which is a parallel fortran code for nuclear structure calculations. I will show scaling behavior and compare the performance of the pure MPI version with the hybrid MPI/OpenMP code on Cray XT4 and XT5 platforms. For large core counts (typically 5,000 and above), the hybrid version is more efficient than pure MPI. With this code, we have been able to predict properties of the unstable nucleus 14F, which have since been confirmed by experiments. I will also give an overview of other recent results for nuclei in the A = 6 to 16 range with 2- and 3-body interactions. Supported in part by US DOE Grant DE-FC02-09ER41582.
Ab initio calculation of valley splitting in monolayer δ-doped phosphorus in silicon.
Drumm, Daniel W; Budi, Akin; Per, Manolo C; Russo, Salvy P; L Hollenberg, Lloyd C
2013-01-01
: The differences in energy between electronic bands due to valley splitting are of paramount importance in interpreting transport spectroscopy experiments on state-of-the-art quantum devices defined by scanning tunnelling microscope lithography. Using vasp, we develop a plane-wave density functional theory description of systems which is size limited due to computational tractability. Nonetheless, we provide valuable data for the benchmarking of empirical modelling techniques more capable of extending this discussion to confined disordered systems or actual devices. We then develop a less resource-intensive alternative via localised basis functions in siesta, retaining the physics of the plane-wave description, and extend this model beyond the capability of plane-wave methods to determine the ab initio valley splitting of well-isolated δ-layers. In obtaining an agreement between plane-wave and localised methods, we show that valley splitting has been overestimated in previous ab initio calculations by more than 50%. PMID:23445785
Ab initio calculation of valley splitting in monolayer δ-doped phosphorus in silicon
2013-01-01
The differences in energy between electronic bands due to valley splitting are of paramount importance in interpreting transport spectroscopy experiments on state-of-the-art quantum devices defined by scanning tunnelling microscope lithography. Using vasp, we develop a plane-wave density functional theory description of systems which is size limited due to computational tractability. Nonetheless, we provide valuable data for the benchmarking of empirical modelling techniques more capable of extending this discussion to confined disordered systems or actual devices. We then develop a less resource-intensive alternative via localised basis functions in siesta, retaining the physics of the plane-wave description, and extend this model beyond the capability of plane-wave methods to determine the ab initio valley splitting of well-isolated δ-layers. In obtaining an agreement between plane-wave and localised methods, we show that valley splitting has been overestimated in previous ab initio calculations by more than 50%. PMID:23445785
Guterding, Daniel; Jeschke, Harald O; Valentí, Roser
2016-01-01
Electronic states with non-trivial topology host a number of novel phenomena with potential for revolutionizing information technology. The quantum anomalous Hall effect provides spin-polarized dissipation-free transport of electrons, while the quantum spin Hall effect in combination with superconductivity has been proposed as the basis for realizing decoherence-free quantum computing. We introduce a new strategy for realizing these effects, namely by hole and electron doping kagome lattice Mott insulators through, for instance, chemical substitution. As an example, we apply this new approach to the natural mineral herbertsmithite. We prove the feasibility of the proposed modifications by performing ab-initio density functional theory calculations and demonstrate the occurrence of the predicted effects using realistic models. Our results herald a new family of quantum anomalous Hall and quantum spin Hall insulators at affordable energy/temperature scales based on kagome lattices of transition metal ions. PMID:27185665
Guterding, Daniel; Jeschke, Harald O.; Valentí, Roser
2016-01-01
Electronic states with non-trivial topology host a number of novel phenomena with potential for revolutionizing information technology. The quantum anomalous Hall effect provides spin-polarized dissipation-free transport of electrons, while the quantum spin Hall effect in combination with superconductivity has been proposed as the basis for realizing decoherence-free quantum computing. We introduce a new strategy for realizing these effects, namely by hole and electron doping kagome lattice Mott insulators through, for instance, chemical substitution. As an example, we apply this new approach to the natural mineral herbertsmithite. We prove the feasibility of the proposed modifications by performing ab-initio density functional theory calculations and demonstrate the occurrence of the predicted effects using realistic models. Our results herald a new family of quantum anomalous Hall and quantum spin Hall insulators at affordable energy/temperature scales based on kagome lattices of transition metal ions. PMID:27185665
NASA Astrophysics Data System (ADS)
Guterding, Daniel; Jeschke, Harald O.; Valentí, Roser
2016-05-01
Electronic states with non-trivial topology host a number of novel phenomena with potential for revolutionizing information technology. The quantum anomalous Hall effect provides spin-polarized dissipation-free transport of electrons, while the quantum spin Hall effect in combination with superconductivity has been proposed as the basis for realizing decoherence-free quantum computing. We introduce a new strategy for realizing these effects, namely by hole and electron doping kagome lattice Mott insulators through, for instance, chemical substitution. As an example, we apply this new approach to the natural mineral herbertsmithite. We prove the feasibility of the proposed modifications by performing ab-initio density functional theory calculations and demonstrate the occurrence of the predicted effects using realistic models. Our results herald a new family of quantum anomalous Hall and quantum spin Hall insulators at affordable energy/temperature scales based on kagome lattices of transition metal ions.
NASA Astrophysics Data System (ADS)
Liu, Jian-Heng; Tu, Matisse Wei-Yuan; Zhang, Wei-Min
2016-07-01
By considering a nanoscale Aharonov-Bohm (AB) interferometer consisting of a laterally coupled double dot coupled to the source and drain electrodes, we investigate the AB phase dependence of the bonding and antibonding states and the transport currents via the bonding and antibonding state channels. The relations of the AB phase dependence between the quantum states and the associated transport current components are analyzed, which provides useful information for the reconstruction of quantum states through the measurement of the transport current in such systems. We also obtain the validity of the experimental analysis [given in T. Hatano et al., Phys. Rev. Lett. 106, 076801 (2011), 10.1103/PhysRevLett.106.076801] that bonding state currents in different energy configurations are almost the same. With the coherent properties in the quantum dot states as well as in the transport currents, we also provide a way to manipulate the bonding and antibonding states through the AB magnetic flux.
NASA Astrophysics Data System (ADS)
Steffen, Matthias
2013-03-01
Quantum mechanics plays a crucial role in many day-to-day products, and has been successfully used to explain a wide variety of observations in Physics. While some quantum effects such as tunneling limit the degree to which modern CMOS devices can be scaled to ever reducing dimensions, others may potentially be exploited to build an entirely new computing architecture: The quantum computer. In this talk I will review several basic concepts of a quantum computer. Why quantum computing and how do we do it? What is the status of several (but not all) approaches towards building a quantum computer, including IBM's approach using superconducting qubits? And what will it take to build a functional machine? The promise is that a quantum computer could solve certain interesting computational problems such as factoring using exponentially fewer computational steps than classical systems. Although the most sophisticated modern quantum computing experiments to date do not outperform simple classical computations, it is increasingly becoming clear that small scale demonstrations with as many as 100 qubits are beginning to be within reach over the next several years. Such a demonstration would undoubtedly be a thrilling feat, and usher in a new era of controllably testing quantum mechanics or quantum computing aspects. At the minimum, future demonstrations will shed much light on what lies ahead.
NASA Astrophysics Data System (ADS)
Ryabov, V. A.
2015-08-01
Quantum systems in a mechanical embedding, the breathing mode of a small particles, optomechanical system, etc. are far not the full list of examples in which the volume exhibits quantum behavior. Traditional consideration suggests strain in small systems as a result of a collective movement of particles, rather than the dynamics of the volume as an independent variable. The aim of this work is to show that some problem here might be essentially simplified by introducing periodic boundary conditions. At this case, the volume is considered as the independent dynamical variable driven by the internal pressure. For this purpose, the concept of quantum volume based on Schrödinger’s equation in 𝕋3 manifold is proposed. It is used to explore several 1D model systems: An ensemble of free particles under external pressure, quantum manometer and a quantum breathing mode. In particular, the influence of the pressure of free particle on quantum oscillator is determined. It is shown also that correction to the spectrum of the breathing mode due to internal degrees of freedom is determined by the off-diagonal matrix elements of the quantum stress. The new treatment not using the “force” theorem is proposed for the quantum stress tensor. In the general case of flexible quantum 3D dynamics, quantum deformations of different type might be introduced similarly to monopole mode.
Realizing Controllable Quantum States
NASA Astrophysics Data System (ADS)
Takayanagi, Hideaki; Nitta, Junsaku
-- 4. Mesoscopic superconductivity with unconventional superconductor or ferromagnet. Ultraefficient microrefrigerators realized with ferromagnet-superconductor junctions / F. Giazotto et al. Anomalous charge transport in triplet superconductor junctions by the synergy effect of the proximity effect and the mid gap Andreev resonant states / Y. Tanaka and S. Kashiwaya. Paramagnetic and glass states in superconductive YBa[symbol]Cu[symbol]O[symbol] ceramics of sub-micron scale grains / H. Deguchi et al. Quantum properties of single-domain triplet superconductors / A. M. Gulian and K. S. Wood. A numerical study of Josephson current in p wave superconducting junctions / Y. Asano et al. Tilted bi-crystal sapphire substrates improve properties of grain boundary YBa[symbol]Cu[symbol]O[symbol] junctions and extend their Josephson response to THZ frequencies / E. Stepantsov et al. Circuit theory analysis of AB-plane tunnel junctions of unconventional superconductor Bi[symbol]Sr[symbol]Ca[symbol]Cu[symbol]O[symbol] / I. Shigeta et al. Transport properties of normal metal/anisotropic superconductor junctions in the eutectic system Sr[symbol]RuO[symbol]Ru / M. Kawamura et al. Macroscopic quantum tunneling in d-wave superconductor Josephson / S. Kawabata et al. Quasiparticle states of high-T[symbol] oxides observed by a Zeeman magnetic field response / S. Kashiwaya et al. Experimentally realizable devices for controlling the motion of magnetic flux quanta in anisotropic superconductors: vortex lenses, vortex diodes and vortex pumps / S. Savel'ev and F. Nori. Stability of vortex-antivortex "molecules" in mesoscopic superconducting triangles / V. R. Misko et al. Superconducting network with magnetic decoration - Hofstadter butterfly in spatially modulated magnetic field / Y. Iye et al. Observation of paramagnetic supercurrent in mesoscopic superconducting rings and disks using multiple-small-tunnel-junction method / A. Kanda et al. Guidance of vortices in high
Quantum games as quantum types
NASA Astrophysics Data System (ADS)
Delbecque, Yannick
In this thesis, we present a new model for higher-order quantum programming languages. The proposed model is an adaptation of the probabilistic game semantics developed by Danos and Harmer [DH02]: we expand it with quantum strategies which enable one to represent quantum states and quantum operations. Some of the basic properties of these strategies are established and then used to construct denotational semantics for three quantum programming languages. The first of these languages is a formalisation of the measurement calculus proposed by Danos et al. [DKP07]. The other two are new: they are higher-order quantum programming languages. Previous attempts to define a denotational semantics for higher-order quantum programming languages have failed. We identify some of the key reasons for this and base the design of our higher-order languages on these observations. The game semantics proposed in this thesis is the first denotational semantics for a lambda-calculus equipped with quantum types and with extra operations which allow one to program quantum algorithms. The results presented validate the two different approaches used in the design of these two new higher-order languages: a first one where quantum states are used through references and a second one where they are introduced as constants in the language. The quantum strategies presented in this thesis allow one to understand the constraints that must be imposed on quantum type systems with higher-order types. The most significant constraint is the fact that abstraction over part of the tensor product of many unknown quantum states must not be allowed. Quantum strategies are a new mathematical model which describes the interaction between classical and quantum data using system-environment dialogues. The interactions between the different parts of a quantum system are described using the rich structure generated by composition of strategies. This approach has enough generality to be put in relation with other
Quantum Monte Carlo Methods for First Principles Simulation of Liquid Water
ERIC Educational Resources Information Center
Gergely, John Robert
2009-01-01
Obtaining an accurate microscopic description of water structure and dynamics is of great interest to molecular biology researchers and in the physics and quantum chemistry simulation communities. This dissertation describes efforts to apply quantum Monte Carlo methods to this problem with the goal of making progress toward a fully "ab initio"…
On the Strain Rate Sensitivity of Abs and Abs Plus Fused Deposition Modeling Parts
NASA Astrophysics Data System (ADS)
Vairis, A.; Petousis, M.; Vidakis, N.; Savvakis, K.
2016-06-01
In this work the effect of strain rate on the tensile strength of fused deposition modeling parts built with Acrylonitrile-butadiene-styrene (ABS) and ABS plus material is presented. ASTM D638-02a specimens were built with ABS and ABS plus and they were tested on a Schenck Trebel Co. tensile test machine at three different test speeds, equal, lower, and higher to the test speed required by the ASTM D638-02a standard. The experimental tensile strength results were compared and evaluated. The fracture surfaces of selected specimens were examined with a scanning electron microscope, to determine failure mode of the filament strands. It was found that, as the test speed increases, specimens develop higher tensile strength and have higher elastic modulus. Specimens tested in the highest speed of the experiment had on average about 10% higher elastic modulus and developed on average about 11% higher tensile strength.
Rapid digestion of Cry34Ab1 and Cry35Ab1 in simulated gastric fluid.
Herman, Rod A; Schafer, Barry W; Korjagin, Valerie A; Ernest, April D
2003-11-01
Two genes were identified in Bacillus thuringiensis Berliner (Bt) that code for the proteins that comprise a Cry34Ab1/Cry35Ab1 binary insecticidal crystal protein. Maize, Zea mays L., plants have been transformed to express the Cry34Ab1/Cry35Ab1 proteins, and as a result, these plants are resistant to attack by western corn rootworm, Diabrotica virgifera virgifera LeConte, a major pest in the Midwestern corn-growing area of the U.S.A. As part of the safety assessment for the proteins, digestibility studies were conducted. Digestion experiments with both proteins demonstrated rapid degradation in simulated gastric fluid, comparable to other registered plant-incorporated protectants. Quantitative and qualitative approaches for determining digestibility are illustrated. PMID:14582981
NASA Astrophysics Data System (ADS)
Chibotaru, L. F.; Ungur, L.
2012-08-01
A methodology for the rigorous nonperturbative derivation of magnetic pseudospin Hamiltonians of mononuclear complexes and fragments based on ab initio calculations of their electronic structure is described. It is supposed that the spin-orbit coupling and other relativistic effects are already taken fully into account at the stage of quantum chemistry calculations of complexes. The methodology is based on the establishment of the correspondence between the ab initio wave functions of the chosen manifold of multielectronic states and the pseudospin eigenfunctions, which allows to define the pseudospin Hamiltonians in the unique way. Working expressions are derived for the pseudospin Zeeman and zero-field splitting Hamiltonian corresponding to arbitrary pseudospins. The proposed calculation methodology, already implemented in the SINGLE_ANISO module of the MOLCAS-7.6 quantum chemistry package, is applied for a first-principles evaluation of pseudospin Hamiltonians of several complexes exhibiting weak, moderate, and very strong spin-orbit coupling effects.
NASA Astrophysics Data System (ADS)
Levy, Amikam; Diósi, Lajos; Kosloff, Ronnie
2016-05-01
In this work we present the concept of a quantum flywheel coupled to a quantum heat engine. The flywheel stores useful work in its energy levels, while additional power is extracted continuously from the device. Generally, the energy exchange between a quantum engine and a quantized work repository is accompanied by heat, which degrades the charging efficiency. Specifically when the quantum harmonic oscillator acts as a work repository, quantum and thermal fluctuations dominate the dynamics. Quantum monitoring and feedback control are applied to the flywheel in order to reach steady state and regulate its operation. To maximize the charging efficiency one needs a balance between the information gained by measuring the system and the information fed back to the system. The dynamics of the flywheel are described by a stochastic master equation that accounts for the engine, the external driving, the measurement, and the feedback operations.
NASA Astrophysics Data System (ADS)
Xu, Ping
We introduce a general notion of quantum universal enveloping algebroids (QUE algebroids), or quantum groupoids, as a unification of quantum groups and star-products. Some basic properties are studied including the twist construction and the classical limits. In particular, we show that a quantum groupoid naturally gives rise to a Lie bialgebroid as a classical limit. Conversely, we formulate a conjecture on the existence of a quantization for any Lie bialgebroid, and prove this conjecture for the special case of regular triangular Lie bialgebroids. As an application of this theory, we study the dynamical quantum groupoid , which gives an interpretation of the quantum dynamical Yang-Baxter equation in terms of Hopf algebroids.
NASA Astrophysics Data System (ADS)
Braun, Daniel; Giraud, Olivier; Braun, Peter A.
2010-03-01
We introduce and study a measure of ``quantumness'' of a quantum state based on its Hilbert-Schmidt distance from the set of classical states. ``Classical states'' were defined earlier as states for which a positive P-function exists, i.e. they are mixtures of coherent states [1]. We study invariance properties of the measure, upper bounds, and its relation to entanglement measures. We evaluate the quantumness of a number of physically interesting states and show that for any physical system in thermal equilibrium there is a finite critical temperature above which quantumness vanishes. We then use the measure for identifying the ``most quantum'' states. Such states are expected to be potentially most useful for quantum information theoretical applications. We find these states explicitly for low-dimensional spin-systems, and show that they possess beautiful, highly symmetric Majorana representations. [4pt] [1] Classicality of spin states, Olivier Giraud, Petr Braun, and Daniel Braun, Phys. Rev. A 78, 042112 (2008)
Mismatched quantum filtering and entropic information
NASA Astrophysics Data System (ADS)
Tsang, Mankei
2014-03-01
Quantum filtering is a signal processing technique that estimates the posterior state of a quantum system under continuous measurements and has become a standard tool in quantum information processing, with applications in quantum state preparation, quantum metrology, and quantum control. If the filter assumes a wrong model due to assumptions or approximations, however, the estimation accuracy is bound to suffer. In this talk I shall present formulas that relate the error penalty caused by quantum filter mismatch to the relative entropy between the true model and the nominal model, with one formula for Gaussian measurements, such as homodyne detection, and another for Poissonian measurements, such as photon counting. These formulas generalize recent seminal results in classical information theory and provide new operational meanings to relative entropy, mutual information, and channel capacity in the context of quantum experiments. See http://arxiv.org/abs/1310.0291 for details. This work is supported by the Singapore National Research Foundation under NRF Grant No. NRF-NRFF2011-07.
ERIC Educational Resources Information Center
Halpern, Arthur M.; Glendening, Eric D.
2013-01-01
A project for students in an upper-level course in quantum or computational chemistry is described in which they are introduced to the concepts and applications of a high quality, ab initio treatment of the ground-state potential energy curve (PEC) for H[subscript 2] and D[subscript 2]. Using a commercial computational chemistry application and a…
Coleman, Piers; Schofield, Andrew J
2005-01-20
As we mark the centenary of Albert Einstein's seminal contribution to both quantum mechanics and special relativity, we approach another anniversary--that of Einstein's foundation of the quantum theory of solids. But 100 years on, the same experimental measurement that puzzled Einstein and his contemporaries is forcing us to question our understanding of how quantum matter transforms at ultra-low temperatures. PMID:15662409