Manolo C. Per

Anisotropic intracule densities and electron correlation in H2: A quantum Monte Carlo study

We derive efficient quantum Monte Carlo estimators for the anisotropic intracule and extracule densities. These estimators are used in conjunction with an accurate explicitly correlated wave function to investigate the bond-length dependence of electron correlation effects in the ground-state H(2) molecule. It is shown that the localized increase in the magnitude of the correlation energy as the bond is stretched is accompanied by highly anisotropic correlation effects. In addition, we find a small long-range part of the Coulomb hole, which is present even at the equilibrium bond length.

Electron-nucleus cusp correction and forces in quantum Monte Carlo

A simple method is presented which ensures the electron-nucleus cusp condition is satisfied by the Slater-Jastrow wavefunctions commonly employed in quantum Monte Carlo simulations. The method is applied in variational energy calculations of the neon atom and a selection of molecules using both Gaussian and Slater basis sets. In addition, we discuss the relationship between the electron-nucleus cusps and the variance of forces, and investigate the sensitivity of forces to the quality of the cusps for various diatomic molecules.

How Important is Orbital Choice in Single-Determinant Diffusion Quantum Monte Carlo Calculations?

The accuracy of total electronic energies obtained using the fixed-node diffusion quantum Monte Carlo (FN-DMC) method is determined by the choice of the many-body nodal surface. Here, we perform a systematic comparison of the quality of FN-DMC energies for a selection of atoms and diatomic molecules using nodal surfaces defined by single determinants of Hartree-Fock, B3LYP, and LDA orbitals. Through comparison with experimental results, we show that the use of Kohn-Sham orbitals results in significantly improved FN-DMC atomization energies over those obtained using Hartree-Fock orbitals. We also discuss the effect of spin contamination in the orbitals.

Ab initio calculation of valley splitting in monolayer δ-doped phosphorus in silicon

AbstractThe 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%.

Quantum Monte Carlo calculations of the potential energy curve of the helium dimer.

We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.

Delocalized oxygen as the origin of two-level defects in Josephson junctions.

One of the key problems facing superconducting qubits and other Josephson junction devices is the decohering effects of bistable material defects. Although a variety of phenomenological models exist, the true microscopic origin of these defects remains elusive. For the first time we show that these defects may arise from delocalization of the atomic position of the oxygen in the oxide forming the Josephson junction barrier. Using a microscopic model, we compute experimentally observable parameters for phase qubits. Such defects are charge neutral but have nonzero response to both applied electric field and strain. This may explain the observed long coherence time of two-level defects in the presence of charge noise, while still coupling to the junction electric field and substrate phonons.

Thermodynamic stability of neutral Xe defects in diamond

Optically active defect centers in diamond are of considerable interest, and ab initio calculations have provided valuable insight into the physics of these systems. Candidate structures for the Xe center in diamond, for which little structural information is known, are modeled using density functional theory. The relative thermodynamic stabilities were calculated for two likely structural arrangements. The split-vacancy structure is found to be the most stable for all temperatures up to 1500 K. A vibrational analysis was also carried out, predicting Raman- and IR-active modes which may aid in distinguishing between center structures.

An analysis of the correlation energy contribution to the interaction energy of inert gas dimers

An accurate description of electron correlation is essential for the calculation of interaction energies in cases where dispersion energy is a major component, for example, for the rare gas atoms, physisorption on graphite, and graphene-graphene interactions. Such calculations are computationally demanding using supermolecule methods and the energies calculated lack a simple, physical interpretation. Alternatively density functional theories (DFTs) may be used to give an approximate estimate of the correlation energy. However, the physical nature of this DFT estimate of electron correlation energy is not well understood and, in fact, most current DFT methods do not describe dispersion energy at all. Hence, an analysis of the correlation energy contribution to interaction energies where dispersion energy is important is needed. In order to do this we provide an analysis of the correlation energy contribution to the potential energy curves of He(2), Ne(2), and Ar(2) in terms of the Hartree-Fock (HF) interaction term DeltaE(int) (HF), a dispersion energy term E(disp) and an electron correlation term DeltaE(int) (C). DeltaE(int) (C) includes all other correlation energy effects besides E(disp) and is shown to be repulsive, of a similar short range character to, but of smaller magnitude than DeltaE(int) (HF). This analysis was used to develop a theoretical model which gives a very good estimate of the potential energy wells for He(2), Ne(2), Ar(2), HeNe, HeAr, and NeAr.

Ab initio electronic properties of dual phosphorus monolayers in silicon

In the midst of the epitaxial circuitry revolution in silicon technology, we look ahead to the next paradigm shift: effective use of the third dimension - in particular, its combination with epitaxial technology. We perform ab initio calculations of atomically thin epitaxial bilayers in silicon, investigating the fundamental electronic properties of monolayer pairs. Quantitative band splittings and the electronic density are presented, along with effects of the layers’ relative alignment and comments on disordered systems, and for the first time, the effective electronic widths of such device components are calculated.

Zero-variance zero-bias quantum Monte Carlo estimators for the electron density at a nucleus

We derive new quantum Monte Carlo (QMC) estimators for the electronic density at the position of a point nucleus using the zero-variance and zero-bias principles. The resulting estimators are highly efficient, and are significantly simpler to implement and use than alternative methods, as they contain no adjustable parameters. In addition, they can be used in both variational and diffusion QMC calculations. Our best estimator is used to calculate the most accurate available estimates of the total electron density at the nucleus for the first-row atoms Li-Ne, the Ar atom, and the diatomic molecules B(2), N(2), and F(2).