An efficient variational many-body wave function for electrons in a quantum dot is introduced. We show that the proposed bosonic two-body correlation factor is able to capture nearly exactly the missing correlations of the lowest-Landau-level many-body wave function. We compare the results with exact diagonalizations. The comparison shows that quantum Monte Carlo techniques combined with the proposed wave function are a powerful tool for studying quantum dots. As further evidence for this, we apply the scheme to a maximum density droplet state for up to 40 electrons.
The ground states of parabolically confined electrons in a quantum dot are studied by both direct numerical diagonalization and quantum Monte Carlo (QMC) methods. We present a simple but accurate variational many-body wave function for the dot in the limit of a weak magnetic field. The wave function has the center-of-mass motion restricted to the lowest-energy state and the electron-electron interaction is taken into account by a Jastrow two-body correlation factor. The optimized wave function has an accuracy very close to the state-of-the-art numerical diagonalization calculations. The results and the computational efficiency indicate that the presented wave function combined with the QMC method suits ideally for studies of large quantum dots.
This thesis investigates the use of wave-function methods for the study of quantum-dot systems. It investigates single dots, using quantum Monte Carlo techniques in a wide range of magnetic field values, and a double-dot system, using the exact-diagonalization method.\n\nThe thesis proposes simple yet accurate many-particle wave functions for various angular-momentum and spin states, for both weak and strong magnetic fields. Using these trial wave functions, it evaluates various properties of dots and studies Wigner crystallization and spin polarization for the weak-field limit. For strong magnetic fields, the thesis investigates ground states of different spin polarizations as a function of the magnetic field and tests the commonly used lowest-Landau-level approximation. The results are compared to calculations from the density-functional theory.\n\nFinally, the thesis presents a method that combines the accuracy of the exact-diagonalization method and the scalability of Monte Carlo methods. The Monte Carlo-based diagonalization is a promising tool for use in situations that cannot be handled with a simple trial wave function, and have too many particles for an exact-diagonalization treatment. For quantum dots, methods for efficiently evaluating the high-magnetic-field basis functions and their gradients are demonstrated.
Recent works indicate that coupled quantum dots (QD) are promising nanostructures for realizing quantum gates that are the basic elements of a quantum computer. In QDs, fundamental quantum logic operations can be achieved in a quantum gate by controlling the entanglement between two electron spins by means of external electric and magnetic fields. This chapter describes the computational approach used to solve the many-body Schrodinger equation. It provides results of the simulations for both the strongly coupled and weakly coupled QDs and analyzes the important differences between the two cases. The utilization of simple trial wavefunctions in VMC have lead to good agreement with the ED method, when electrons are decoupled, which validates two-level approaches such as the Heitler-London approximation in weakly coupled QDs.
Variational two-electron reduced-density-matrix (2RDM) theory is applied to computing energy spectra and properties of few-electron quantum dots. The model system is a two-dimensional electron gas with a central confinement potential. For each orbital angular momentum
J, the energy and 2RDM are computed by the variational 2RDM method in which the energy is minimized as a functional of the 2RDM. In the minimization, which is performed by semidefinite programming, the 2RDM is constrained to represent a
N-electron wave function with angular momentum
J by
N- and
J-representability conditions [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Advantages of the variational 2RDM method include (i) lower bounds on the energies for all
J values, (ii) calculation of approximate 2RDMs in polynomial time without many-electron wave functions, (iii) exploitation of angular symmetry in the sparse block-diagonal structure of the 2RDM, (iv) accurate description of multireference correlation (entanglement) effects, and (v) direct calculation of one- and two-electron properties from the 2RDM. With the 2RDM we directly compute pair-correlation functions, radial charge densities, and average radial electron displacements in the quantum dot. Energies and properties are compared to those from solving the
N-electron Schr\"odinger equation by large-scale, exact diagonalization. It is found that the accuracy of the variational 2RDM approach is sensitive to the total orbital angular momentum and the symmetry of the final wave function. For quantum dots of high symmetry, the variational algorithm isolates a highly accurate solution that recovers the correlation energy within a few percent.
We develop a variational Monte Carlo (VMC) method for electron-phonon coupled systems. The VMC method has been extensively used for investigating strongly correlated electrons over the last decades. However, its applications to electron-phonon coupled systems have been severely restricted because of its large Hilbert space. Here, we propose a variational wave function with a large number of variational parameters, which is suitable and tractable for systems with electron-phonon coupling. In the proposed wave function, we implement an unexplored electron-phonon correlation factor, which takes into account the effect of the entanglement between electrons and phonons. The method is applied to systems with diagonal electron-phonon interactions, i.e., interactions between charge densities and lattice displacements (phonons). As benchmarks, we compare VMC results with previous results obtained by the exact diagonalization, the Green function Monte Carlo method and the density matrix renormalization group for the Holstein and Holstein-Hubbard model. From these benchmarks, we show that the present method offers an efficient way to treat strongly coupled electron-phonon systems.