zbMATH — the first resource for mathematics

Real-time adaptive finite element solution of time-dependent Kohn-Sham equation. (English) Zbl 1354.65187
Summary: In our previous paper [J. Comput. Phys. 231, No. 14, 4967–4979 (2012; Zbl 1245.65125)], a general framework of using adaptive finite element methods to solve the Kohn-Sham equation has been presented. This work is concerned with solving the time-dependent Kohn-Sham equations. The numerical methods are studied in the time domain, which can be employed to explain both the linear and the nonlinear effects. A Crank-Nicolson scheme and linear finite element space are employed for the temporal and spatial discretizations, respectively. To resolve the trouble regions in the time-dependent simulations, a heuristic error indicator is introduced for the mesh adaptive methods. An algebraic multigrid solver is developed to efficiently solve the complex-valued system derived from the semi-implicit scheme. A mask function is employed to remove or reduce the boundary reflection of the wavefunction. The effectiveness of our method is verified by numerical simulations for both linear and nonlinear phenomena, in which the effectiveness of the mesh adaptive methods is clearly demonstrated.

65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
Full Text: DOI
[1] Bao, G.; Hu, G. H.; Liu, D., An h-adaptive finite element solver for the calculations of the electronic structures, J. Comput. Phys., 231, 14, 4967-4979, (2012) · Zbl 1245.65125
[2] Bao, G.; Hu, G. H.; Liu, D., Numerical solution of the Kohn-Sham equation by finite element methods with an adaptive mesh redistribution technique, J. Sci. Comput., 55, 2, 372-391, (2013) · Zbl 1272.82003
[3] Bao, G.; Hu, G. H.; Liu, D.; Luo, S. T., Multi-physical modeling and multi-scale computation of nano-optical responses, (Recent Advances in Scientific Computing and Applications, Contemporary Mathematics, vol. 586, (2013)), 43-55 · Zbl 1330.78005
[4] Brandt, Achi; Livne, Oren E., Multigrid techniques, (2011), Society for Industrial and Applied Mathematics · Zbl 1227.65121
[5] Castro, A.; Appel, H.; Oliveira, M.; Rozzi, C. A.; Andrade, X.; Lorenzen, F.; Marques, M. A.L.; Gross, E. K.U.; Rubio, A., Octopus: a tool for the application of time-dependent density functional theory, Phys. Status Solidi B, 243, 2465-2488, (2006)
[6] Chen, Z. J., Efficient modeling techniques for time-dependent quantum system with applications to carbon nanotubes, (2010), University of Massachusetts Amherst, Master’s thesis
[7] Cleary, A.; Falgout, R.; Henson, V.; Jones, J.; Manteuffel, T.; McCormick, S.; Miranda, G.; Ruge, J., Robustness and scalability of algebraic multigrid, SIAM J. Sci. Comput., 21, 5, 1886-1908, (2000) · Zbl 0959.65049
[8] Cohen, O.; Kronik, L.; Brandt, A., Locally refined multigrid solution of the all-electron Kohn-Sham equation, J. Chem. Theory Comput., 9, 11, 4744-4760, (2013)
[9] Corkum, P. B., Plasma perspective on strong field multiphoton ionization, Phys. Rev. Lett., 71, 1994-1997, (1993)
[10] Fattebert, J. L.; Hornung, R. D.; Wissink, A. M., Finite elements approach for density functional theory calculations on locally refined meshes, J. Comput. Phys., 223, 759-773, (2007) · Zbl 1132.65106
[11] Fiolhais, C.; Nogueira, F.; Marques, M. A.L., A primer in density functional theory, Lecture Notes in Physics, (2003), Springer
[12] Fowe, E. P.; Bandrauk, A. D., Nonlinear time-dependent density-functional-theory study of ionization and harmonic generation in CO2 by ultrashort intense laser pulses: orientational effects, Phys. Rev. A, 81, 023411, (Feb 2010)
[13] Gross, E. K.U.; Dobson, J. F.; Petersilka, M., Topics in current chemistry, vol. 181, (1996), Springer-Verlag Heidelberg
[14] Gygi, F., Architecture of qbox: a scalable first-principles molecular dynamics code, IBM J. Res. Dev., 52, 1-2, 137-144, (2008)
[15] Hong, X. H.; Wang, F., TDDFT calculation for photoabsorption spectra of li_n clusters, Phys. Lett. A, 375, 18, 1883-1888, (2011)
[16] Kohn, W.; Sham, L. J., Self-consistent equations including exchange and correlation effects, Phys. Rev., 140, A1133-A1138, (1965)
[17] Lehtovaara, L.; Havu, V.; Puska, M., All-electron density functional theory and time-dependent density functional theory with high-order finite elements, J. Chem. Phys., 131, 054103, (2009)
[18] Lehtovaara, L.; Havu, V.; Puska, M., All-electron time-dependent density functional theory with finite elements: time-propagation approach, J. Chem. Phys., 135, 154104, (2011)
[19] Marques, M. A.L.; Castro, A.; Bertsch, G. F.; Rubio, A., Octopus: a first-principles tool for excited electron-ion dynamics, Comput. Phys. Commun., 151, 60-78, (2003)
[20] Marques, M. A.L.; Gross, E. K.U., Time-dependent density functional theory, Annu. Rev. Phys. Chem., 55, 427-455, (2004)
[21] Mortensen, J. J.; Hansen, L. B.; Jacobsen, K. W., Real-space grid implementation of the projector augmented wave method, Phys. Rev. B, 71, 3, 035109, (2005)
[22] Motamarri, P.; Nowak, M. R.; Leiter, K.; Knap, J.; Gavini, V., Higher-order adaptive finite-element methods for Kohn-Sham density functional theory, J. Comput. Phys., 253, 0, 308-343, (2013) · Zbl 1349.74331
[23] Pask, J. E.; Sterne, P. A., Finite element methods in ab initio electronic structure calculations, Model. Simul. Mater. Sci. Eng., 13, 3, R71, (2005)
[24] Runge, E.; Gross, E. K.U., Density-functional theory for time-dependent systems, Phys. Rev. Lett., 52, 997-1000, (1984)
[25] Schauer, V.; Linder, C., All-electron Kohn-Sham density functional theory on hierarchic finite element spaces, J. Comput. Phys., 250, 644-664, (2013)
[26] Schleife, André; Draeger, Erik W.; Kanai, Yosuke; Correa, Alfredo A., Plane-wave pseudopotential implementation of explicit integrators for time-dependent Kohn-Sham equations in large-scale simulations, J. Chem. Phys., 137, 22A546, (2012)
[27] Soba, A., A finite element method solver for time-dependent and stationary Schrödinger equations with a generic potential, Commun. Comput. Phys., 5, 928-941, (2009) · Zbl 1364.65201
[28] Soler, J. M.; Artacho, E.; Gale, J. D.; Garcia, A.; Junquera, J.; Ordejon, P.; Sanchez-Portal, D., The SIESTA method for ab initio order-N materials simulation, J. Phys. Condens. Matter, 14, 11, 2745, (2002)
[29] Suryanarayana, P.; Gavini, V.; Blesgen, T.; Bhattacharya, K.; Ortiz, M., Non-periodic finite-element formulation of Kohn-Sham density functional theory, J. Mech. Phys. Solids, 58, 256-280, (2010) · Zbl 1193.81006
[30] Ullrich, C. A., Time-dependent density-functional theory: concepts and applications, Oxford Graduate Texts, (2012), OUP Oxford · Zbl 1237.81004
[31] van Leeuwen, R., Causality and symmetry in time-dependent density-functional theory, Phys. Rev. Lett., 80, 1280-1283, (1998)
[32] Verfürth, R., Robust a posteriori error estimates for nonstationary convection-diffusion equations, SIAM J. Numer. Anal., 43, 4, 1783-1802, (2005) · Zbl 1099.65077
[33] Van Der Vorst, H. A., Bi-CGSTAB: a fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems, SIAM J. Sci. Comput., 13, 631-644, (1992) · Zbl 0761.65023
[34] Walter, M.; Häkkinen, H.; Lehtovaara, L.; Puska, M.; Enkovaara, J.; Rostgaard, C.; Mortensen, J. J., Time-dependent density-functional theory in the projector augmented-wave method, J. Chem. Phys., 128, 24, 244101, (2008)
[35] Zhang, D.; Shen, L.; Zhou, A.; Gong, X., Finite element method for solving Kohn-Sham equations based on self-adaptive tetrahedral mesh, Phys. Lett. A, 372, 5071-5076, (2008) · Zbl 1221.81225
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.