2,498 research outputs found

    Electrostatic considerations affecting the calculated HOMO-LUMO gap in protein molecules.

    Get PDF
    A detailed study of energy differences between the highest occupied and lowest unoccupied molecular orbitals (HOMO-LUMO gaps) in protein systems and water clusters is presented. Recent work questioning the applicability of Kohn-Sham density-functional theory to proteins and large water clusters (E. Rudberg, J. Phys.: Condens. Mat. 2012, 24, 072202) has demonstrated vanishing HOMO-LUMO gaps for these systems, which is generally attributed to the treatment of exchange in the functional used. The present work shows that the vanishing gap is, in fact, an electrostatic artefact of the method used to prepare the system. Practical solutions for ensuring the gap is maintained when the system size is increased are demonstrated. This work has important implications for the use of large-scale density-functional theory in biomolecular systems, particularly in the simulation of photoemission, optical absorption and electronic transport, all of which depend critically on differences between energies of molecular orbitals.Comment: 13 pages, 4 figure

    Linear-scaling quantum Monte Carlo technique with non-orthogonal localized orbitals

    Get PDF
    We have reformulated the quantum Monte Carlo (QMC) technique so that a large part of the calculation scales linearly with the number of atoms. The reformulation is related to a recent alternative proposal for achieving linear-scaling QMC, based on maximally localized Wannier orbitals (MLWO), but has the advantage of greater simplicity. The technique we propose draws on methods recently developed for linear-scaling density functional theory. We report tests of the new technique on the insulator MgO, and show that its linear-scaling performance is somewhat better than that achieved by the MLWO approach. Implications for the application of QMC to large complex systems are pointed out

    O(N) methods in electronic structure calculations

    Full text link
    Linear scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high performance computers. The linear scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas is then discussed. The applications of linear scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear scaling methods are discussed.Comment: 85 pages, 15 figures, 488 references. Resubmitted to Rep. Prog. Phys (small changes

    Second Stage String Fragmentation Model

    Get PDF
    A string model, advocated by Bowler, provides a physical and intuitive picture of heavy quark fragmentation. When supplemented by an ad hoc factor of (1-z), to suppress fragmentation near z=1, it supplies an excellent fit to the data. We extend Bowler's model by accounting for the further decay of the massive mesonic states produced by the initial string breaking. We find that each subsequent string break and cascade decay beyond the first, introduces a factor of (1-z). Furthermore we find that including a finite mass for the quarks, which pop out of the vacuum and split the string, forces the first string breaking to produce massive states requiring further decay. This sequence terminates at the second stage of fragmentation where only relatively "light" heavy meson systems are formed. Thus we naturally account for the phenomenologically required factor of (1-z). We also predict that the ratio of (primary) fragments-vector/(vector plus scalar) should be .61. Our second stage string fragmentation model provides an appealing picture of heavy quark fragmentation.Comment: 15 page
    • …
    corecore