Introduction to Configuration Path Integral Monte Carlo

In low-temperature high-density plasmas quantum effects of the electrons are becoming increasingly important. This requires the development of new theoretical and computational tools. Quantum Monte Carlo methods are among the most successful approaches to first-principle simulations of many-body quantum systems. In this chapter we present a recently developed method---the configuration path integral Monte Carlo (CPIMC) method for moderately coupled, highly degenerate fermions at finite temperatures. It is based on the second quantization representation of the $N$-particle density operator in a basis of (anti-)symmetrized $N$-particle states (configurations of occupation numbers) and allows to tread arbitrary pair interactions in a continuous space. We give a detailed description of the method and discuss the application to electrons or, more generally, Coulomb-interacting fermions. As a test case we consider a few quantum particles in a one-dimensional harmonic trap. Depending on the coupling parameter (ratio of the interaction energy to kinetic energy), the method strongly reduces the sign problem as compared to direct path integral Monte Carlo (DPIMC) simulations in the regime of strong degeneracy which is of particular importance for dense matter in laser plasmas or compact stars. In order to provide a self-contained introduction, the chapter includes a short introduction to Metropolis Monte Carlo methods and the second quantization of quantum mechanics.Comment: chapter in book "Introduction to Complex Plasmas: Scientific Challenges and Technological Opportunities", Michael Bonitz, K. Becker, J. Lopez and H. Thomsen (Eds.) Springer Series "Atomic, Optical and Plasma Physics", vol. 82, Springer 2014, pp. 153-194 ISBN: 978-3-319-05436-0 (Print) 978-3-319-05437-7 (Online

Coupled coarse graining and Markov Chain Monte Carlo for lattice systems

We propose an efficient Markov Chain Monte Carlo method for sampling equilibrium distributions for stochastic lattice models, capable of handling correctly long and short-range particle interactions. The proposed method is a Metropolis-type algorithm with the proposal probability transition matrix based on the coarse-grained approximating measures introduced in a series of works of M. Katsoulakis, A. Majda, D. Vlachos and P. Plechac, L. Rey-Bellet and D.Tsagkarogiannis,. We prove that the proposed algorithm reduces the computational cost due to energy differences and has comparable mixing properties with the classical microscopic Metropolis algorithm, controlled by the level of coarsening and reconstruction procedure. The properties and effectiveness of the algorithm are demonstrated with an exactly solvable example of a one dimensional Ising-type model, comparing efficiency of the single spin-flip Metropolis dynamics and the proposed coupled Metropolis algorithm.Comment: 20 pages, 4 figure

Non-Fermi-liquid d-wave metal phase of strongly interacting electrons

Developing a theoretical framework for conducting electronic fluids qualitatively distinct from those described by Landau's Fermi-liquid theory is of central importance to many outstanding problems in condensed matter physics. One such problem is that, above the transition temperature and near optimal doping, high-transition-temperature copper-oxide superconductors exhibit `strange metal' behaviour that is inconsistent with being a traditional Landau Fermi liquid. Indeed, a microscopic theory of a strange-metal quantum phase could shed new light on the interesting low-temperature behaviour in the pseudogap regime and on the d-wave superconductor itself. Here we present a theory for a specific example of a strange metal---the 'd-wave metal'. Using variational wavefunctions, gauge theoretic arguments, and ultimately large-scale density matrix renormalization group calculations, we show that this remarkable quantum phase is the ground state of a reasonable microscopic Hamiltonian---the usual t-J model with electron kinetic energy $t$ and two-spin exchange $J$ supplemented with a frustrated electron `ring-exchange' term, which we here examine extensively on the square lattice two-leg ladder. These findings constitute an explicit theoretical example of a genuine non-Fermi-liquid metal existing as the ground state of a realistic model.Comment: 22 pages, 12 figures: 6 pages, 7 figures of main text + 16 pages, 5 figures of Supplementary Information; this is approximately the version published in Nature, minus various subedits in the main tex

Theoretical study of the insulating oxides and nitrides: SiO2, GeO2, Al2O3, Si3N4, and Ge3N4

An extensive theoretical study is performed for wide bandgap crystalline oxides and nitrides, namely, SiO_{2}, GeO_{2}, Al_{2}O_{3}, Si_{3}N_{4}, and Ge_{3}N_{4}. Their important polymorphs are considered which are for SiO_{2}: $\alpha$-quartz, $\alpha$- and $\beta$-cristobalite and stishovite, for GeO_{2}: $\alpha$-quartz, and rutile, for Al_{2}O_{3}: $\alpha$-phase, for Si_{3}N_{4} and Ge_{3}N_{4}: $\alpha$- and $\beta$-phases. This work constitutes a comprehensive account of both electronic structure and the elastic properties of these important insulating oxides and nitrides obtained with high accuracy based on density functional theory within the local density approximation. Two different norm-conserving \textit{ab initio} pseudopotentials have been tested which agree in all respects with the only exception arising for the elastic properties of rutile GeO_{2}. The agreement with experimental values, when available, are seen to be highly satisfactory. The uniformity and the well convergence of this approach enables an unbiased assessment of important physical parameters within each material and among different insulating oxide and nitrides. The computed static electric susceptibilities are observed to display a strong correlation with their mass densities. There is a marked discrepancy between the considered oxides and nitrides with the latter having sudden increase of density of states away from the respective band edges. This is expected to give rise to excessive carrier scattering which can practically preclude bulk impact ionization process in Si_{3}N_{4} and Ge_{3}N_{4}.Comment: Published version, 10 pages, 8 figure

Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Recent experiments have confirmed the importance of nuclear quantum effects even in large biomolecules at physiological temperature. Here we describe how the path integral formalism can be used to describe rigorously the nuclear quantum effects on equilibrium and kinetic properties of molecules. Specifically, we explain how path integrals can be employed to evaluate the equilibrium (EIE) and kinetic (KIE) isotope effects, and the temperature dependence of the rate constant. The methodology is applied to the [1,5] sigmatropic hydrogen shift in pentadiene. Both the KIE and the temperature dependence of the rate constant confirm the importance of tunneling and other nuclear quantum effects as well as of the anharmonicity of the potential energy surface. Moreover, previous results on the KIE were improved by using a combination of a high level electronic structure calculation within the harmonic approximation with a path integral anharmonicity correction using a lower level method.Comment: 9 pages, 4 figure

Liquid metal nanodroplet dynamics inside nanocontainers

Here we report direct observations of spatial movements of nanodroplets of Pb metal trapped inside sealed carbon nanocontainers. We find drastic changes in the mobility of the liquid droplets as the particle size increases from a few to a few ten nanometers. In open containers the droplet becomes immobile and readily evaporates to the vacuum environment. The particle mobility strongly depends on confinement, particle size, and wetting on the enclosed surface. The collisions between droplets increase mobility but the tendency is reversed if collisions lead to droplet coalescence. The dynamics of confined nanodroplets could provide new insights into the activity of nanostructures in spatially constrained geometries

Uniform electron gases

We show that the traditional concept of the uniform electron gas (UEG) --- a homogeneous system of finite density, consisting of an infinite number of electrons in an infinite volume --- is inadequate to model the UEGs that arise in finite systems. We argue that, in general, a UEG is characterized by at least two parameters, \textit{viz.} the usual one-electron density parameter $\rho$ and a new two-electron parameter $\eta$. We outline a systematic strategy to determine a new density functional $E(\rho,\eta)$ across the spectrum of possible $\rho$ and $\eta$ values.Comment: 8 pages, 2 figures, 5 table

Static and dynamic ionization levels of transition metal-doped zinc chalcogenides

Transition metal (TM) impurities in semiconductors have a considerable effect on the electronic properties and on the lattice vibrations. The unfilled d shell permits the impurity atoms to exist in a variety of charge states. In this work, the static donor and acceptor ionization energies of ZnX:M, with X = S, Se, Te and M:Sc, Ti, V, Fe, Co, Ni are obtained from first principles total energy calculations and compared with experimental results in the literature where they exist. From these results, many of the TM-doped zinc chalogenides have an amphoteric behavior. To analyze the rule of the deep gap levels in both the radiative and non-radiative processes, the dynamic ionization energies are obtained as a function of the inward and outward M–X displacements. In many cases, the changes in the mass and the force constants resulting from the substitution of an impurity center for a lattice atom are small. When the charge or the environment of the impurity changes, the electron population tend to remain compensated. As consequence, the changes in the lattice vibrational modes are small

Oxidation behavior of graphene-coated copper at intrinsic graphene defects of different origins

The development of ultrathin barrier films is vital to the advanced semiconductor industry. Graphene appears to hold promise as a protective coating; however, the polycrystalline and defective nature of engineered graphene hinders its practical applications. Here, we investigate the oxidation behavior of graphene-coated Cu foils at intrinsic graphene defects of different origins. Macro-scale information regarding the spatial distribution and oxidation resistance of various graphene defects is readily obtained using optical and electron microscopies after the hot-plate annealing. The controlled oxidation experiments reveal that the degree of structural deficiency is strongly dependent on the origins of the structural defects, the crystallographic orientations of the underlying Cu grains, the growth conditions of graphene, and the kinetics of the graphene growth. The obtained experimental and theoretical results show that oxygen radicals, decomposed from water molecules in ambient air, are effectively inverted at Stone-Wales defects into the graphene/Cu interface with the assistance of facilitators

Gamification and Simulation

Gamification and simulation methods are two of the most important components of serious games. In order to create an effective training tool, it is imperative to understand these methods and their relationship to each other. If designed correctly, gamification techniques can build upon simulations to provide an effective training medium, which enhances learning, engagement and motivation in users. This chapter discusses their uses, strengths and weaknesses whilst identifying how to most effectively utilise them in developing serious games