DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation−dissipation (ACFD) theorem (namely the Rutgers−Chalmers vdW-DF, Vydrov−Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko−Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PYSCF as a development environment. We then summarize the capabilities of PYSCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PYSCF across the domains of quantum chemistry, materials science, machine learning and quantum information science.

No abstract

Non-covalent π−π interactions are central to chemical and biological processes, yet the full understanding of their origin that would unite the simplicity of empirical approaches with the accuracy of quantum calculations is still missing. Here we employ a quantum-mechanical Hamiltonian model for van der Waals interactions, to demonstrate that intermolecular electron correlation in large supramolecular complexes at equilibrium distances is appropriately described by collective charge fluctuations. We visualize these fluctuations and provide connections both to orbital-based approaches to electron correlation, as well as to the simple London pairwise picture. The reported binding energies of ten supramolecular complexes obtained from the quantum-mechanical fluctuation model joined with density functional calculations are within 5% of the reference energies calculated with the diffusion quantum Monte-Carlo method. Our analysis suggests that π−π stacking in supramolecular complexes can be characterized by strong contributions to the binding energy from delocalized, collective charge fluctuations—in contrast to complexes with other types of bonding.

Noncovalent van der Waals (vdW) interactions are responsible for a wide range of phenomena in matter. Popular density-functional methods that treat vdW interactions use disparate physical models for these intricate forces, and as a result the applicability of existing methods is often restricted to a subset of relevant molecules and materials. Aiming towards a general-purpose density functional model of vdW interactions, here we unify two complementary approaches: nonlocal vdW functionals for polarization and interatomic methods for many-body interactions. The developed nonlocal many-body dispersion method (MBD-NL) increases the accuracy and efficiency of existing vdW functionals and is shown to be broadly applicable to molecules, soft and hard materials including ionic and metallic compounds, as well as organic/inorganic interfaces.arXiv:1910.03073v2 [cond-mat.mtrl-sci]

In recent years, we have been witnessing a paradigm shift in computational materials science. In fact, traditional methods, mostly developed in the second half of the XXth century, are being complemented, extended, and sometimes even completely replaced by faster, simpler, and often more accurate approaches. The new approaches, that we collectively label by machine learning, have their origins in the fields of informatics and artificial intelligence, but are making rapid inroads in all other branches of science. With this in mind, this Roadmap article, consisting of multiple contributions from experts across the field, discusses the use of machine learning in materials science, and share perspectives on current and future challenges in problems as diverse as the prediction of materials properties, the construction of force-fields, the development of exchange correlation functionals for density-functional theory, the solution of the many-body problem, and more. In spite of the already numerous and exciting success stories, we are just at the beginning of a long path that will reshape materials science for the many challenges of the XXIth century.

We present an approach for computing long-range van der Waals (vdW) interactions between complex molecular systems and arbitrarily shaped macroscopic bodies, melding atomistic treatments of electronic fluctuations based on density functional theory in the former with continuum descriptions of strongly shapedependent electromagnetic fields in the latter, thus capturing many-body and multiple scattering effects to all orders. Such a theory is especially important when considering vdW interactions at mesoscopic scales, i.e., between molecules and structured surfaces with features on the scale of molecular sizes, in which case the finite sizes, complex shapes, and resulting nonlocal electronic excitations of molecules are strongly influenced by electromagnetic retardation and wave effects that depend crucially on the shapes of surrounding macroscopic bodies. We show that these effects together can modify vdW interaction energies and forces, as well as molecular shapes deformed by vdW interactions, by orders of magnitude compared to previous treatments based on Casimir-Polder, nonretarded, or pairwise approximations, which are valid only at macroscopically large or atomic-scale separations or in dilute insulating media, respectively. DOI: 10.1103/PhysRevLett.118.266802 Van der Waals (vdW) interactions play an essential role in noncovalent phenomena throughout biology, chemistry, and condensed-matter physics [1][2][3]. It has long been known that vdW interactions among a system of polarizable atoms are not pairwise additive but instead strongly depend on geometric and material properties [2,4,5]. However, only recently developed theoretical methods have made it possible to account for short-range quantum interactions in addition to long-range many-body screening in molecular ensembles [3,[6][7][8][9][10][11][12][13][14][15], demonstrating that nonlocal many-body effects cannot be captured by simple, pairwise-additive descriptions; these calculations typically neglect electromagnetic retardation effects in molecular systems. Simultaneously, recent theoretical and experimental work has characterized dipolar Casimir-Polder (CP) interactions between macroscopic metallic or dielectric objects and atoms, molecules, or Bose-Einstein condensates, further extending to nonzero temperatures, dynamical situations, and fluctuations in excited states (as in so-called Rydberg atoms) [16][17][18][19][20][21][22][23][24][25]. Yet, while theoretical treatments have thus far accounted for the full electrodynamic response of macroscopic bodies (including retardation), they often treat molecules as point dipoles of some effective bulk permittivities or as collections of noninteracting atomic dipoles, ignoring finite size and other many-body electromagnetic effects.In this Letter, motivated by the aforementioned theoretical developments [1,[16][17][18][24][25][26][27][28], we describe an approach that seamlessly connects atomistic descriptions of large molecules to continuum descriptions of arbitrary macroscopic bodies, characterizing t...