Arxiv:1609.00252V1

Challenges in Large Scale Quantum Mechanical Calculations Laura E. Ratcliff,1 Stephan Mohr,2 Georg Huhs,2 Thierry Deutsch,3, 4 Michel Masella,5 and Luigi Genovese3,4, ∗ 1Argonne Leadership Computing Facility, Argonne National Laboratory, Illinois 60439, USA 2Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC-CNS), Barcelona, Spain 3Univ. Grenoble Alpes, INAC-MEM, L Sim, F-38000 Grenoble, France 4CEA, INAC-MEM, L Sim, F-38000 Grenoble, France 5Laboratoire de Biologie Structurale et Radiologie, Service de Bioénergétique, Biologie Structurale et Mécanisme, Institut de Biologie et de Technologie de Saclay, CEA Saclay, F-91191 Gif-sur-Yvette Cedex, France (Dated: September 2, 2016) During the past decades, quantum mechanical methods have undergone an amazing transition from pioneering investigations of experts into a wide range of practical applications, made by a vast community of researchers. First principles calculations of systems containing up to a few hundred atoms have become a standard in many branches of science. The sizes of the systems which can be simulated have increased even further during recent years, and quantum-mechanical calculations of systems up to many thousands of atoms are nowadays possible. This opens up new appealing possibilities, in particular for interdisciplinary work, bridging together communities of different needs and sensibilities. In this review we will present the current status of this topic, and will also give an outlook on the vast multitude of applications, challenges and opportunities stimulated by electronic structure calculations, making this field an important working tool and bringing together researchers of many different domains. I. INTRODUCTION non-interacting fermions that provide the same distribu- tion of the density. These are the fundamental ideas of Density Functional Theory (DFT). The fundamental laws for a Quantum Mechani- cal (QM) description of atomistic systems up to the DFT has been, for more than twenty years, the nanoscale are known and have been well established for workhorse method for simulations within the solid state a little less than a century. Yet, there are many chal- community. Moreover, in spite of the fact that DFT dras- lenges related to the Quantum Mechanical treatment of tically reduces the complexity with respect to the ab ini- large systems. In the vast majority of cases, we are still tio methods of Quantum Chemistry, the success of such unable to solve the fundamental Schrödinger equation a treatment in the latter community is undeniable. This for systems of realistic sizes in such a way that the re- is mainly due to the fact that, on one hand, the qual- sults satisfy “universal” requirements of accuracy, pre- ity of the exchange and correlation functionals available cision and especially predictability. Unfortunately, this permits the calculation of certain properties with almost also implies that we are still far from being able to quan- chemical accuracy, and, on the other hand, there are nu- titatively predict experimental results at the nanoscale. merous software packages, which are relatively easy to use, that have contributed to the diffusion of the compu- The problems are not only related to the computa- tational approach. tional complexity needed to solve the equations of QM, there are also intrinsic obstacles. To give an example, let During the past years, there has been a multiplication us remind the so-called “Coulson’s challenge”. In 1960, of DFT software packages that are able to treat systems Coulson1 noticed that the most compact object needed of increasingly large size. This has, on one hand, been to characterize quantum mechanically an N-electron sys- enabled by both the advances in supercomputing archi- arXiv:1609.00252v1 [physics.chem-ph] 1 Sep 2016 tem (at least in its ground state) is the two-body reduced tectures and the code developers continuously improving density matrix (2RDM). However, it turns out that we do their codes to exploit the steadily increasing performance not know all the necessary conditions for the 2RDM to provided, but it is at the same time also motivated by be N-representable, i.e. coming from an anti-symmetric various scientific needs. This fact clearly extends the wavefunction of an N-electron system. Thus, even if a range of possible applications to new fields, and to com- compact (and, in principle, computationally accessible) munities traditionally focused on larger systems. In a object exists, theoretical and algorithmic bottlenecks hin- similar manner to the uptake of DFT in the Quantum der its practical usage. In 19642 and 19653, Kohn, Ho- Chemistry community, things are progressing as if we henberg and Sham further reduced the complexity by are entering a “second era” of DFT calculations, where showing that the electronic density is in a one-to-one DFT and, more generally, large-scale quantum mechani- correspondence with the ground state energy of a sys- cal treatments, are susceptible to wide diffusion in other tem of interacting electrons, and that such an interact- communities. ing system can be replaced by a mean-field problem of N In this review paper we will present some of the mo- 2 tivations that led computational physicists and quantum typical QM approaches are developed and improved is chemists into this second era. The different aspects will of the order of few atoms, even though some of these be separated into various subcategories, while trying to concepts are then also applied to larger systems. John give a general overview, and will be completed by no- Perdew introduced the renowned metaphor of Jacob’s table examples in the literature. This inspection of the Ladder 4, where the computational complexity of the im- state-of-the-art will provide the reader with an outlook plementation of the DFT exchange and correlation func- on the present capabilities of QM approaches. We will tional is (in principle) directly related to the accuracy then continue our discussion by presenting the key con- of the description, aiming at the “heaven” of chemical cepts that have emerged in the last decade. These con- accuracy. cepts are often specific to QM calculations at large scale Likewise, for more than ten years, a lot of work has and are rather different from those which are typical of been done to extend the range of applicability of QM traditional calculations, where the systems’ sizes are lim- methods to larger systems. The so-called nearsighted- ited to a few hundreds of orbitals. These concepts are ness property5 suggested that, at least in principle, one therefore of high importance for potential users of such could exploit locality to build linear scaling methods that advanced DFT methods. are able to reach larger scales. Initially, the develop- ment of such computational methods was driven by the “academic” purpose of verifying the computational con- A. The need for large-scale QM sequences of nearsightedness. In Sec. II we will overview the most important advancements in this topic and some Given the unbiased predictive power of Quantum Me- of the established computational approaches in large- chanics, there is obviously no need to explain why sys- scale QM. This is by no means new and there are a num- tems containing only a few atoms should be modelled ber of valid review papers on the topic, to which we will using this approach. However, for simulations of sys- also refer. Our aim is not to be fully exhaustive on this tems at the nanoscale, composed of many thousands of topic as there has been many research studies in this di- atoms, the question of the need for a Quantum Mechan- rection. However we would like to put the emphasis on ical treatment might appear legitimate. For systems of the fact that nowadays the panorama is so rich and there these sizes, the electronic degrees of freedom are seldom is enough diversity in the computational approaches to of interest, and the interatomic potential might be de- claim that such a discipline is now mature enough to be scribed by more compact approaches like Force Fields, largely diffused also among non-specialists. possibly fine-tuned for describing experimentally known The reason for this diffusion is related to the opportu- structural and dynamical or polarizability properties. In nities that a QM approach opens for systems composed other terms, the intimate nature of the problem changes: of thousands, if not hundreds of thousands, atoms. On instead of focussing solely on the correct estimation of in- the one hand, there are quantities which are intrinsi- teractions and correlation between electrons, one rather cally only accessible using QM, for example all investiga- has to concentrate on the exhaustive sampling of the con- tions dealing with electronic excitations6; we will present figuration space, thereby losing the need for an intrinsic some more examples in Sec. III and IV. On the other quantum mechanical description. hand, large QM calculation are also needed to access In addition we know that, even if a QM approach were error bars and statistics of the results. An example feasible, this would not necessarily lead to a better de- is the need to get good statistics among different con- scription. Although the complexity of the model is cer- stituents in a morphology—a task which is not possible tainly higher and the description is less biased, there are by implicit, classical modelling of the environment. In still many approximations which are hidden in a Quan- addition, another aspect where first-principles QM ap- tum Mechanical calculation; therefore we are—even for proaches are important is the need for validation of non- a system containing only a handful of atoms—in general QM approaches. still far from chemical accuracy. The situation for large For all of these tasks, there is a typical length scale, systems will be the same or even worse, since more se- ranging from a few hundred to many thousands of atoms, vere approximations have to be adopted.

Load more