Abstract:Electronic-structure calculations play a fundamental role in predicting important physical (optical, electrical, etc) properties of condensed matter. Based on density functional theory (DFT), Kohn-Sham (KS) equation replaces the interacting many-body (electrons) problem by an equivalent set of self-consistent single-particle equations. Different from well-adopted basis set approach, we present a real-space finite-difference method to discretize the KS equation where the Laplacian operator is represented by high-order differences. After spatial discretization, a derived nonlinear eigenvalue problem is solved by a self-consistent field (SCF) iteration scheme with the update of electron density and potentials. In spite of the simplicity and versatility of the finite-difference method, one has to perform large-scale calculations of both eigenvalue and electrostatic problems at each iteration step. For eigenvalue problems, we minimize the Rayleigh quotients by the conjugate gradient method; and find minimum eigenvalues of interest through subspace diagonalization without Gram-Schmidt procedure. For obtaining electrostatic potentials, we use Dirichlet boundary condition to truncate the computational domain after inserting electron density and Gaussian compensating charge density together into Poisson's equation. Then a preconditioned Krylov subspace solver is employed to solve the Poisson's equation. A three-dimensional quantum dot with few spin-polarized electrons (each electron has spin up) is simulated by the real-space DFT incorporating local spin density (LSD) approximation. Wave function, electron density, eigen-energy, and total energy are calculated and compared to published results. The real-space finite-difference method produces reliable simulation results with a high computational efficiency. This work is fundamentally important to quantum-mechanical ab-initio calculation.

A Discrete Time-Dependent Method for Metastable Atoms and Molecules in Intense Fields

A Unified Hamiltonian Solution to Maxwell–Schrödinger Equations for Modeling Electromagnetic Field–particle Interaction

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Application of Coulomb Wave Function Discrete Variable Representation to Atomic Systems in Strong Laser Fields.

Modeling Rabi Oscillation By Rigorously Solving Maxwell-Schrodinger Equation

Application of a Numerical-Basis-state Method to Strong-Field Excitation and Ionization of Hydrogen Atoms

Dynamic Tunnelling Ionization of H_2^+ in Intense Fields

A Hybrid Method for Integrated Atomistic-Continuum Simulation of Ultrashort-Pulsed Laser Interaction with Semiconductors

A pr 2 00 6 Coulomb Wave Function DVR : Application to Atomic Systems in Strong Laser Fields

A combined R-matrix eigenstate basis set and finite-differences propagation method for the time-dependent Schr\"{od}dinger equation: the one-electron case

Above-Threshold Ionization and Laser-Induced Electron Diffraction in Diatomic Molecules

Few-Electron Dynamics of Atoms and Molecules in Intense Laser Fields

Real-time solutions of coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations: fundamentals, implementation, and nano-optical applications

Efficient Time-Sampling Method in Coulomb-corrected Strong-Field Approximation

Application of Discrete Variable Representation to Planar ${\rm H}_2^+$H2+ in Strong Xuv Laser Fields

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Atomic hydrogen in circularly polarized, high-intensity and high frequency laser fields

Accurate Computation of above Threshold Ionization Spectra for Stretched H-2(+) in Strong Laser Fields

Quantitative theory for electron-nuclear energy sharing in molecular ionization