Nanoscale van der Waals interactions

Exploring HCl-HCl interactions: QZVPP calculations, improved Lennard-Jones potential, and second virial coefficient analysis for thermodynamics and industrial applications

In this paper, we present a comprehensive analysis of HCl-HCl interactions, including QZVPP calculations, energy fitting, conformation validation, and the determination of the second virial coefficient B using improved Lennard-Jones (ILJ) potential parameters. To acquire accurate interaction energies, initial QZVPP calculations are performed on approximately 1851 randomly generated HCl-HCl conformations. Then, these energies are used to fit an improved Lennard-Jones potential energy surface, allowing for a robust description of HCl-HCl interactions. The ILJ potential parameters are then used to validate particular HCl dimer conformations, ensuring their stability and consistency with experimental observations. The correlation between calculated and experimental conformations strengthens the validity of the ILJ potential parameters. In addition, the second viral coefficient B is calculated at various temperatures using the ILJ potential. The obtained B values are compared to experimental data, demonstrating close agreement, and validating the ILJ potential's ability to accurately capture the intermolecular interactions and gas-phase behavior of the HCl-HCl system. The results of this study demonstrate the effective implementation of QZVPP calculations, energy fitting, and ILJ potential parameters in validating HCl-HCl conformations and accurately determining the second virial coefficient B. The high degree of concordance between calculated B values and experimental data demonstrates the validity of the ILJ potential and its suitability for modeling HCl-HCl interactions. This research contributes to a greater comprehension of HCl-HCl interactions and their implications for numerous chemical and atmospheric processes. The validated conformations, energy fitting method, and calculated second virial coefficients provide valuable instruments for future research and pave the way for more accurate modeling and simulations of HCl-HCl systems.

Including dispersion in density functional theory for adsorption on flat oxide surfaces, in metal-organic frameworks and in acidic zeolites

We examine the performance of nine commonly used methods for including dispersion interactions in density functional theory (DFT): three different parametrizations of damped 1/Rn terms (n = 6, 8, …) added to the DFT energy (Grimme's D2 and D3 parameterizations as well as that of Tkatchenko and Scheffler), three different implementations of the many-body dispersion approach (MBD, MBD/HI and MBD/FI), the density-dependent energy correction, called dDsC, and two "first generation" van der Waals density functionals, revPBE-vdW and optB86b-vdW. As test set we use eight molecule-surface systems for which agreement has been reached between experiment and hybrid QM:QM calculations within chemical accuracy limits (±4.2 kJ mol-1). It includes adsorption of carbon monoxide and dioxide in the Mg2(2,5-dioxido-1,4-benzenedicarboxylate) metal-organic framework (Mg-MOF-74, CPO-27-Mg), adsorption of carbon monoxide as well as of monolayers of methane and ethane on the MgO(001) surface, as well as adsorption of methane, ethane and propane in H-chabazite (H-CHA). D2 with Ne parameters for Mg2+, D2(Ne), MBD/HI and MBD/FI perform best. With the PBE functional, the mean unsigned errors are 6.1, 5.6 and 5.4 kJ mol-1, respectively.

Van der Waals forces in free and wetting liquid films

Van der Waals interactions induced by fluctuations of electromagnetic field bear universal nature and act between individual atoms, condensed particles or bodies of any type. Continuously growing interest to theoretical understanding as well as to precise evaluation of van der Waals forces is caused by their fundamental role in many physical, chemical, and biological processes. In this paper, we scrutinize progress in the studies of van der Waals forces, related to recent active development of Coupled Dipole Method (CDM) for the analysis of the behavior and properties of nanosized systems. The application of CDM for the analysis of thin liquid films allowed achieving substantial progress in understanding the behavior of free and wetting films. It was shown that both the macroscopic properties, such as excess free energy and Hamaker constants and the local microscopic parameters, such as polarizabilities, can be successfully calculated based only on properties of individual molecules. The impact of lateral film confinement on the specific excess free energy and the film stability was elucidated, and effect of spatial constraints on the spectrum of vibrational states for liquid film and the underlying substrate was analyzed. It was shown that van der Waals interactions between molecules represent the universal mechanism for dynamic structuring and formation of boundary layers and that the CDM allows self-consistently calculating the properties of these layers in both solid and liquid phases.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

Tailoring van der Waals dispersion interactions with external electric charges

van der Waals (vdW) dispersion interactions strongly impact the properties of molecules and materials. Often, the description of vdW interactions should account for the coupling with pervasive electric fields, stemming from membranes, ionic channels, liquids, or nearby charged functional groups. However, this quantum-mechanical effect has been omitted in atomistic simulations, even in widely employed electronic-structure methods. Here, we develop a model and study the effects of an external charge on long-range vdW correlations. We show that a positive external charge stabilizes dispersion interactions, whereas a negative charge has an opposite effect. Our analytical results are benchmarked on a series of (bio)molecular dimers and supported by calculations with high-level correlated quantum-chemical methods, which estimate the induced dispersion to reach up to 35% of intermolecular binding energy (4 kT for amino-acid dimers at room temperature). Our analysis bridges electrostatic and electrodynamic descriptions of intermolecular interactions and may have implications for non-covalent reactions, exfoliation, dissolution, and permeation through biological membranes.

Phonon-Polariton Mediated Thermal Radiation and Heat Transfer among Molecules and Macroscopic Bodies: Nonlocal Electromagnetic Response at Mesoscopic Scales

Thermal radiative phenomena can be strongly influenced by the coupling of phonons and long-range electromagnetic fields at infrared frequencies. Typically employed macroscopic descriptions of thermal fluctuations often ignore atomistic effects that become relevant at nanometric scales, whereas purely microscopic treatments ignore long-range, geometry-dependent electromagnetic effects. We describe a mesoscopic framework for modeling thermal fluctuation phenomena among molecules near macroscopic bodies, conjoining atomistic treatments of electronic and vibrational fluctuations obtained from density functional theory in the former with continuum descriptions of electromagnetic scattering in the latter. The interplay of these effects becomes particularly important at mesoscopic scales, where phonon polaritons can be strongly influenced by the objects' finite sizes, shapes, and nonlocal or many-body response to electromagnetic fluctuations. We show that, even in small but especially in elongated low-dimensional molecules, such effects can modify thermal emission and heat transfer by orders of magnitude and produce qualitatively different behavior compared to predictions based on local, dipolar, or pairwise approximations.

Unifying Microscopic and Continuum Treatments of van der Waals and Casimir Interactions

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 shape-dependent 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.

Many-body dispersion effects in the binding of adsorbates on metal surfaces

A correct description of electronic exchange and correlation effects for molecules in contact with extended (metal) surfaces is a challenging task for first-principles modeling. In this work, we demonstrate the importance of collective van der Waals dispersion effects beyond the pairwise approximation for organic-inorganic systems on the example of atoms, molecules, and nanostructures adsorbed on metals. We use the recently developed many-body dispersion (MBD) approach in the context of density-functional theory [Tkatchenko et al., Phys. Rev. Lett. 108, 236402 (2012) and Ambrosetti et al., J. Chem. Phys. 140, 18A508 (2014)] and assess its ability to correctly describe the binding of adsorbates on metal surfaces. We briefly review the MBD method and highlight its similarities to quantum-chemical approaches to electron correlation in a quasiparticle picture. In particular, we study the binding properties of xenon, 3,4,9,10-perylene-tetracarboxylic acid, and a graphene sheet adsorbed on the Ag(111) surface. Accounting for MBD effects, we are able to describe changes in the anisotropic polarizability tensor, improve the description of adsorbate vibrations, and correctly capture the adsorbate-surface interaction screening. Comparison to other methods and experiment reveals that inclusion of MBD effects improves adsorption energies and geometries, by reducing the overbinding typically found in pairwise additive dispersion-correction approaches.

van der Waals dispersion interactions in molecular materials: beyond pairwise additivity

van der Waals (vdW) dispersion interactions are a key ingredient in the structure, stability, and response properties of many molecular materials and essential for us to be able to understand and design novel intricate molecular systems. Pairwise-additive models of vdW interactions are ubiquitous, but neglect their true quantum-mechanical many-body nature. In this perspective we focus on recent developments and applications of methods that can capture collective and many-body effects in vdW interactions. Highlighting a number of recent studies in this area, we demonstrate both the need for and usefulness of explicit many-body treatments for obtaining qualitative and quantitative accuracy for modelling molecular materials, with applications presented for small-molecule dimers, supramolecular host-guest complexes, and finally stability and polymorphism in molecular crystals.

Image-charge forces in thin interlayers due to surface charges in electrolyte

The surface forces arising in wetting films of nonpolar liquids or in thin air interlayers between an electrolyte and a nonpolar medium in the case of discrete charging of the dielectric-electrolyte interface are considered. The contributions of polarization effects to the distribution of the electrostatic potential in the three contacting media were calculated. Within the Debye-Hückel approximation, the analytical solutions were derived for the disjoining pressure in thin films, for the case of either dilute or relatively concentrated electrolyte solutions in the aforementioned systems. Analysis of the analytical and numerical results demonstrated that for dilute solutions the contribution of image forces to the disjoining pressure may significantly exceed the van der Waals forces for films from a few to tens of nanometers thick.

Non-pairwise additivity of the leading-order dispersion energy

The leading-order (i.e., dipole-dipole) dispersion energy is calculated for one-dimensional (1D) and two-dimensional (2D) infinite lattices, and an infinite 1D array of infinitely long lines, of doubly occupied locally harmonic wells. The dispersion energy is decomposed into pairwise and non-pairwise additive components. By varying the force constant and separation of the wells, the non-pairwise additive contribution to the dispersion energy is shown to depend on the overlap of density between neighboring wells. As well separation is increased, the non-pairwise additivity of the dispersion energy decays. The different rates of decay for 1D and 2D lattices of wells is explained in terms of a Jacobian effect that influences the number of nearest neighbors. For an array of infinitely long lines of wells spaced 5 bohrs apart, and an inter-well spacing of 3 bohrs within a line, the non-pairwise additive component of the leading-order dispersion energy is -0.11 kJ mol(-1) well(-1), which is 7% of the total. The polarizability of the wells and the density overlap between them are small in comparison to that of the atomic densities that arise from the molecular density partitioning used in post-density-functional theory (DFT) damped dispersion corrections, or DFT-D methods. Therefore, the nonadditivity of the leading-order dispersion observed here is a conservative estimate of that in molecular clusters.

Many-body van der Waals interactions in molecules and condensed matter

This work reviews the increasing evidence that many-body van der Waals (vdW) or dispersion interactions play a crucial role in the structure, stability and function of a wide variety of systems in biology, chemistry and physics. Starting with the exact expression for the electron correlation energy provided by the adiabatic connection fluctuation-dissipation theorem, we derive both pairwise and many-body interatomic methods for computing the long-range dispersion energy by considering a model system of coupled quantum harmonic oscillators within the random-phase approximation. By coupling this approach to density functional theory, the resulting many-body dispersion (MBD) method provides an accurate and efficient scheme for computing the frequency-dependent polarizability and many-body vdW energy in molecules and materials with a finite electronic gap. A select collection of applications are presented that ascertain the fundamental importance of these non-bonded interactions across the spectrum of intermolecular (the S22 and S66 benchmark databases), intramolecular (conformational energies of alanine tetrapeptide) and supramolecular (binding energy of the 'buckyball catcher') complexes, as well as molecular crystals (cohesive energies in oligoacenes). These applications demonstrate that electrodynamic response screening and beyond-pairwise many-body vdW interactions--both captured at the MBD level of theory--play a quantitative, and sometimes even qualitative, role in describing the properties considered herein. This work is then concluded with an in-depth discussion of the challenges that remain in the future development of reliable (accurate and efficient) methods for treating many-body vdW interactions in complex materials and provides a roadmap for navigating many of the research avenues that are yet to be explored.

Scaling laws for van der Waals interactions in nanostructured materials

Van der Waals interactions have a fundamental role in biology, physics and chemistry, in particular in the self-assembly and the ensuing function of nanostructured materials. Here we utilize an efficient microscopic method to demonstrate that van der Waals interactions in nanomaterials act at distances greater than typically assumed, and can be characterized by different scaling laws depending on the dimensionality and size of the system. Specifically, we study the behaviour of van der Waals interactions in single-layer and multilayer graphene, fullerenes of varying size, single-wall carbon nanotubes and graphene nanoribbons. As a function of nanostructure size, the van der Waals coefficients follow unusual trends for all of the considered systems, and deviate significantly from the conventionally employed pairwise-additive picture. We propose that the peculiar van der Waals interactions in nanostructured materials could be exploited to control their self-assembly.

Van der Waals interactions have a large influence on phenomena that occur at short-length scales. Gobre et al. demonstrate that van der Waals interactions in low-dimensional materials act at very large distances, and can significantly influence the self-assembly of nanostructured systems.

Radiative Exchange between Graphitic Nanostructures: A Microscopic Perspective

Electromagnetic radiative heat exchange involving graphene nanostrucrures is studied using an atomistic approach based on the coupled dipole method modified by the fluctuation dissipation theorem. This method includes taking into account many-particle electromagnetic contributions and enables treating two or more nanostructures with nontrivial boundary conditions at different temperatures. We present a microscopic picture of the heat exchange process in graphene nanostructured based systems in terms of a transmission coefficient, characteristic temperature function, and atomic morphology. Our studies provide general pathways of near-field radiation control at the nanoscale.

Interatomic methods for the dispersion energy derived from the adiabatic connection fluctuation-dissipation theorem

Interatomic pairwise methods are currently among the most popular and accurate ways to include dispersion energy in density functional theory calculations. However, when applied to more than two atoms, these methods are still frequently perceived to be based on ad hoc assumptions, rather than a rigorous derivation from quantum mechanics. Starting from the adiabatic connection fluctuation-dissipation (ACFD) theorem, an exact expression for the electronic exchange-correlation energy, we demonstrate that the pairwise interatomic dispersion energy for an arbitrary collection of isotropic polarizable dipoles emerges from the second-order expansion of the ACFD formula upon invoking the random-phase approximation (RPA) or the full-potential approximation. Moreover, for a system of quantum harmonic oscillators coupled through a dipole-dipole potential, we prove the equivalence between the full interaction energy obtained from the Hamiltonian diagonalization and the ACFD-RPA correlation energy. This property makes the Hamiltonian diagonalization an efficient method for the calculation of the many-body dispersion energy. In addition, we show that the switching function used to damp the dispersion interaction at short distances arises from a short-range screened Coulomb potential, whose role is to account for the spatial spread of the individual atomic dipole moments. By using the ACFD formula, we gain a deeper understanding of the approximations made in the interatomic pairwise approaches, providing a powerful formalism for further development of accurate and efficient methods for the calculation of the dispersion energy.

Electronically coarse-grained model for water

We introduce an electronically coarse-grained description of water representing all long range, many-body electronic responses via an embedded quantum oscillator. Leading-order response coefficients and gas phase electrostatic moments are exactly reproduced. Molecular dynamics, using electronic path integral sampling, shows that this framework is sufficient for a realistic liquid to emerge naturally with transferability extending further to nonambient state points and to the free water surface. The model allows the strength of many-body dispersion and polarization to be adjusted independently and these are found to have significant effects on the condensed phase.

Seamless and Accurate Modeling of Organic Molecular Materials

The near endless possibilities for assembling molecular materials has long posed a difficult challenge for theory. All crystal-structure prediction methods acknowledge the crucial contribution of van der Waals or dispersion interactions, but few go beyond a pairwise additive description of dispersion, ignoring its many-body nature. Here we use two databases to show how a many-body approach to dispersion can seamlessly model both solid and gas-phase interactions within the coveted "chemical accuracy" benchmark, while the underlying pairwise approach fails for solid-state interactions due to the absence of many-body polarization and energy contributions. Our results show that recently developed methods that treat the truly collective nature of dispersion interactions are able to reach the accuracy required for predicting molecular materials, when coupled with nonempirical density functionals.

Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory

Electron dispersion forces play a crucial role in determining the structure and properties of biomolecules, molecular crystals, and many other systems. However, an accurate description of dispersion is highly challenging, with the most widely used electronic structure technique, density functional theory (DFT), failing to describe them with standard approximations. Therefore, applications of DFT to systems where dispersion is important have traditionally been of questionable accuracy. However, the last decade has seen a surge of enthusiasm in the DFT community to tackle this problem and in so-doing to extend the applicability of DFT-based methods. Here we discuss, classify, and evaluate some of the promising schemes to emerge in recent years. A brief perspective on the outstanding issues that remain to be resolved and some directions for future research are also provided.

Van der waals coefficients for nanostructures: fullerenes defy conventional wisdom

The van der Waals coefficients between quasispherical nanostructures can be modeled accurately and analytically by those of classical solid spheres (for nanoclusters) or spherical shells (for fullerenes) of uniform valence electron density, with the true static dipole polarizability. Here, we derive analytically and confirm numerically from this model the size dependencies of the van der Waals coefficients of all orders, showing, for example, that the asymptotic dependence for C(6) is the expected n(2) for pairs of nanoclusters A(n)-A(n), each containing n atoms, but n(2.75) for pairs of single-walled fullerenes C(n)-C(n). Large fullerenes are argued to have much larger polarizabilities and dispersion coefficients than those predicted by either the standard atom pair-potential model or widely used nonlocal van der Waals correlation energy functionals.

A Density Functional with Spherical Atom Dispersion Terms

A new hybrid density functional, APF, is introduced, which avoids the spurious long-range attractive or repulsive interactions that are found in most density functional theory (DFT) models. It therefore provides a sound baseline for the addition of an empirical dispersion correction term, which is developed from a spherical atom model (SAM). The APF-D empirical dispersion model contains nine adjustable parameters that were selected based on a very small training set (15 noble gas dimers and 4 small hydrocarbon dimers), along with two computed atomic properties (ionization potential and effective atomic polarizability) for each element. APF-D accurately describes a large portion of the potential energy surfaces of complexes of noble gas atoms with various diatomic molecules involving a wide range of elements and of dimers of small hydrocarbons, and it reproduces the relative conformational energies of organic molecules. The accuracy for these weak interactions is comparable to that of CCSD(T)/aug-cc-pVTZ calculations. The accuracy in predicting the geometry of hydrogen bond complexes is competitive with other models involving DFT and empirical dispersion.

Collective many-body van der Waals interactions in molecular systems

Van der Waals (vdW) interactions are ubiquitous in molecules and condensed matter, and play a crucial role in determining the structure, stability, and function for a wide variety of systems. The accurate prediction of these interactions from first principles is a substantial challenge because they are inherently quantum mechanical phenomena that arise from correlations between many electrons within a given molecular system. We introduce an efficient method that accurately describes the nonadditive many-body vdW energy contributions arising from interactions that cannot be modeled by an effective pairwise approach, and demonstrate that such contributions can significantly exceed the energy of thermal fluctuations--a critical accuracy threshold highly coveted during molecular simulations--in the prediction of several relevant properties. Cases studied include the binding affinity of ellipticine, a DNA-intercalating anticancer agent, the relative energetics between the A- and B-conformations of DNA, and the thermodynamic stability among competing paracetamol molecular crystal polymorphs. Our findings suggest that inclusion of the many-body vdW energy is essential for achieving chemical accuracy and therefore must be accounted for in molecular simulations.

Accurate and efficient method for many-body van der Waals interactions

An efficient method is developed for the microscopic description of the frequency-dependent polarizability of finite-gap molecules and solids. This is achieved by combining the Tkatchenko-Scheffler van der Waals (vdW) method [Phys. Rev. Lett. 102, 073005 (2009)] with the self-consistent screening equation of classical electrodynamics. This leads to a seamless description of polarization and depolarization for the polarizability tensor of molecules and solids. The screened long-range many-body vdW energy is obtained from the solution of the Schrödinger equation for a system of coupled oscillators. We show that the screening and the many-body vdW energy play a significant role even for rather small molecules, becoming crucial for an accurate treatment of conformational energies for biomolecules and binding of molecular crystals. The computational cost of the developed theory is negligible compared to the underlying electronic structure calculation.