429 research outputs found

    Alternative sampling for variational quantum Monte Carlo

    Full text link
    Expectation values of physical quantities may accurately be obtained by the evaluation of integrals within Many-Body Quantum mechanics, and these multi-dimensional integrals may be estimated using Monte Carlo methods. In a previous publication it has been shown that for the simplest, most commonly applied strategy in continuum Quantum Monte Carlo, the random error in the resulting estimates is not well controlled. At best the Central Limit theorem is valid in its weakest form, and at worst it is invalid and replaced by an alternative Generalised Central Limit theorem and non-Normal random error. In both cases the random error is not controlled. Here we consider a new `residual sampling strategy' that reintroduces the Central Limit Theorem in its strongest form, and provides full control of the random error in estimates. Estimates of the total energy and the variance of the local energy within Variational Monte Carlo are considered in detail, and the approach presented may be generalised to expectation values of other operators, and to other variants of the Quantum Monte Carlo method.Comment: 14 pages, 9 figure

    Accurate structure factors from pseudopotential methods

    Full text link
    Highly accurate experimental structure factors of silicon are available in the literature, and these provide the ideal test for any \emph{ab initio} method for the construction of the all-electron charge density. In a recent paper [J. R. Trail and D. M. Bird, Phys. Rev. B {\bf 60}, 7863 (1999)] a method has been developed for obtaining an accurate all-electron charge density from a first principles pseudopotential calculation by reconstructing the core region of an atom of choice. Here this method is applied to bulk silicon, and structure factors are derived and compared with experimental and Full-potential Linear Augmented Plane Wave results (FLAPW). We also compare with the result of assuming the core region is spherically symmetric, and with the result of constructing a charge density from the pseudo-valence density + frozen core electrons. Neither of these approximations provide accurate charge densities. The aspherical reconstruction is found to be as accurate as FLAPW results, and reproduces the residual error between the FLAPW and experimental results.Comment: 6 Pages, 3 figure

    Papers in Southeast Asian Linguistics No. 7

    Get PDF

    Core reconstruction in pseudopotential calculations

    Full text link
    A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {\bf 14}, 3795 (1981)]. In this method the core region is \emph{reconstructed}, and none of the simplifying approximations (such as spherical symmetry of the charge density/potential or frozen core electrons) that previous solutions to this problem have required are made. The embedding method requires an accurate real space Green function, and an analysis of the errors introduced in constructing this from a set of numerical eigenstates is given. Results are presented for an all-electron reconstruction of bulk aluminium, for both the charge density and the density of states.Comment: 14 pages, 5 figure

    Hydrogen incorporation and charge balance in natural zircon

    Get PDF
    The water and trace element contents of natural igneous zircons were determined to constrain the mechanism of hydrogen incorporation. The low radiation-damage zircons were derived from Fe–Ti oxide gabbros from the Vema Fracture Zone (11°N, Mid-Atlantic Ridge). They contain up to 1212 ppmw H2O, 1.9 wt.% Y2O3 and 0.6 wt.% P2O5 and are generally strongly zoned. REE + Y are partially charge-balanced by P (Y, REE3+ + P5+ = Zr4+ + Si4+), but a large REE excess is present. On an atomic basis, this excess is closely approximated by the amount of H present in the zircons. We therefore conclude that H is incorporated by a charge-balance mechanism (H+ + REE3+ = Zr4+). This interpretation is consistent with FTIR data of the Vema zircons, which shows a strongly polarised main absorption band at ca. 3100 cm−1, similar to experimentally grown Lu-doped hydrous zircon. The size of this 3100 cm−1 band scales with H and REE contents. Apart from a small overlapping band at 3200 cm−1, no other absorption bands are visible, indicating that a hydrogrossular-type exchange mechanism does not appear to be operating in these zircons. Because of charge-balanced uptake of H, P and REE in zircon, the partitioning of these elements into zircon is dependent on each of their concentrations. For instance, DREEzrc/melt increases with increasing H and P contents of the melt, whereas DHzrc/melt increases with increasing REE content but decreases with increasing P content. In addition, H–P–REE systematics of sector zoning indicate kinetic effects may play an important role. Hence, using H in zircon to determine the water content of melts is problematic, and REE partitioning studies need to take into account P and H2O contents of the melt

    An analysis of arthroscopic knot techniques – looking for the perfect knot

    Full text link

    Density-functional embedding using a plane-wave basis

    Full text link
    The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {\bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {\bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded electron densities, the sum of which is the electron density of the system of interest. Coupling between the substrate and embedded systems is achieved via approximate kinetic energy functionals. Bulk aluminium is examined as a test case for which there is a strong interaction between the substrate and embedded systems. A number of approximations to the kinetic-energy functional, both semi-local and non-local, are investigated. It is found that Kohn-Sham results can be well reproduced using a non-local kinetic energy functional, with the total energy accurate to better than 0.1 eV per atom and good agreement between the electron densities.Comment: 11 pages, 4 figure
    • …
    corecore