Long-Range Repulsion Between Spatially Confined van der Waals Dimers

libMBD: A general-purpose package for scalable quantum many-body dispersion calculations

Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.

Confinement-induced nonlocality and casimir force in transdimensional systems

We study within the framework of the Lifshitz theory the long-range Casimir force for in-plane isotropic and anisotropic free-standing transdimensional material slabs. In the former case, we show that the confinement-induced nonlocality not only weakens the attraction of ultrathin slabs but also changes the distance dependence of the material-dependent correction to the Casimir force to go as contrary to the ∼1/l dependence of that of the local Lifshitz force. In the latter case, we use closely packed array of parallel aligned single-wall carbon nanotubes in a dielectric layer of finite thickness to demonstrate strong orientational anisotropy and crossover behavior for the inter-slab attractive force in addition to its reduction with decreasing slab thickness. We give physical insight as to why such a pair of ultrathin slabs prefers to stick together in the perpendicularly oriented manner, rather than in the parallel relative orientation as one would customarily expect.

Optimized Quantum Drude Oscillators for Atomic and Molecular Response Properties

The quantum Drude oscillator (QDO) is an efficient yet accurate coarse-grained approach that has been widely used to model electronic and optical response properties of atoms and molecules as well as polarization and dispersion interactions between them. Three effective parameters (frequency, mass, and charge) fully characterize the QDO Hamiltonian and are adjusted to reproduce response properties. However, the soaring success of coupled QDOs for many-atom systems remains fundamentally unexplained, and the optimal mapping between atoms/molecules and oscillators has not been established. Here we present an optimized parametrization (OQDO) where the parameters are fixed by using only dipolar properties. For the periodic table of elements as well as small molecules, our model accurately reproduces atomic (spatial) polarization potentials and multipolar dispersion coefficients, elucidating the high promise of the presented model in the development of next-generation quantum-mechanical force fields for (bio)molecular simulations.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

Molecular Interactions Induced by a Static Electric Field in Quantum Mechanics and Quantum Electrodynamics

By means of quantum mechanics and quantum electrodynamics applied to coupled harmonic Drude oscillators, we study the interaction between two neutral atoms or molecules subject to a uniform static electric field. Our focus is to understand the interplay between leading contributions to field-induced electrostatics/polarization and dispersion interactions, as considered within the employed Drude model for both non-retarded and retarded regimes. For the first case, we present an exact solution for two coupled oscillators obtained by diagonalizing the corresponding quantum-mechanical Hamiltonian and demonstrate that the external field can control the strength of different intermolecular interactions and relative orientations of the molecules. In the retarded regime described by quantum electrodynamics, our analysis shows that field-induced electrostatic and polarization energies remain unchanged (in isotropic and homogeneous vacuum) compared to the non-retarded case. For interacting species modeled by quantum Drude oscillators, the developed framework based on quantum mechanics and quantum electrodynamics yields the leading contributions to molecular interactions under the combined action of external and vacuum fields.

Unprecedently large 37Cl/35Cl equilibrium isotopic fractionation on nano-confinement of chloride anion

Confinement can result in unusual properties leading to new, exciting discoveries in the nano-realm. One such consequence of confinement at the nanoscale is extremally large isotopic fractionation, especially at sub-van der Waals distances. Herein, on the example of chlorine isotope effects, we show that at conditions of nanoencapsulation these effects may reach values by far larger than observed for the bulk environment, which in the case of nanotubes can lead to practical applications (e.g., in isotopic enrichment) and needs to be considered in analytical procedures that employ nanomaterials.

Interactions between large molecules pose a puzzle for reference quantum mechanical methods

Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol^-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.

Motion of water monomers reveals a kinetic barrier to ice nucleation on graphene

The interfacial behaviour of water remains a central question to fields as diverse as protein folding, friction and ice formation. While the properties of water at interfaces differ from those in the bulk, major gaps in our knowledge limit our understanding at the molecular level. Information concerning the microscopic motion of water comes mostly from computation and, on an atomic scale, is largely unexplored by experiment. Here, we provide a detailed insight into the behaviour of water monomers on a graphene surface. The motion displays remarkably strong signatures of cooperative behaviour due to repulsive forces between the monomers, enhancing the monomer lifetime ( ≈ 3 s at 125 K) in a free-gas phase that precedes the nucleation of ice islands and, in turn, provides the opportunity for our experiments to be performed. Our results give a molecular perspective on a kinetic barrier to ice nucleation, providing routes to understand and control the processes involved in ice formation.

Coulomb interactions between dipolar quantum fluctuations in van der Waals bound molecules and materials

Mutual Coulomb interactions between electrons lead to a plethora of interesting physical and chemical effects, especially if those interactions involve many fluctuating electrons over large spatial scales. Here, we identify and study in detail the Coulomb interaction between dipolar quantum fluctuations in the context of van der Waals complexes and materials. Up to now, the interaction arising from the modification of the electron density due to quantum van der Waals interactions was considered to be vanishingly small. We demonstrate that in supramolecular systems and for molecules embedded in nanostructures, such contributions can amount to up to 6 kJ/mol and can even lead to qualitative changes in the long-range van der Waals interaction. Taking into account these broad implications, we advocate for the systematic assessment of so-called Dipole-Correlated Coulomb Singles in large molecular systems and discuss their relevance for explaining several recent puzzling experimental observations of collective behavior in nanostructured materials.

Impact of nanoconfinement on acetylacetone Equilibria in Ordered Mesoporous Silicates

Nanoconfinement is one of the most intriguing nanoscale effects and affects several physical and chemical properties of molecules and materials, including viscosity, reaction kinetics, and glass transition temperature. In this work, liquid nuclear magnetic resonance (NMR) was used to analyze the behavior of 2,4-pentadienone in ordered mesoporous materials with a pore diameter of between 3 and 10 nm. The liquid NMR results showed meaningful changes in the hydrogen chemical shift and the keto-enol chemical equilibrium, which were associated with the pore diameter, allowing the authors to observe the effects of nanoconfinement. An interesting phenomenon was observed where the chemical equilibria of 2,4-pentadienone confined in a mesoporous material with a pore diameter of 3.5 nm was similar to that obtained with free (bulk) 2,4-pentadienone in larger pore materials. Another interesting result was observed for the enthalpy and entropy of the tautomeric equilibria of 2,4-pentadienone confined in mesoporous materials with a 5.5 nm pore diameter being -7.9 kJ mol^-1 and -15.9 J mol^-1.K. These values are similar to those obtained by dimethyl sulfoxide. This phenomenon indicates the possible use of ordered mesoporous materials as a reaction substitute in organic solvents. It was further observed that while the values of enthalpy (ΔH) and entropy (ΔS) had been modified by confinement, the Gibbs free energy (ΔG) value remained closer to that observed in free (bulk) 2,4-pentadienone. It is expected that this study will help in understanding the effects of nanoconfinement and provide a simple method to employ NMR techniques to analyze these phenomena.

Bridging transitions and capillary forces for colloids in a slit

Capillary bridges can form between colloids immersed in a two-phase fluid, e.g., in a binary liquid mixture, if the surface of the colloids prefers the species other than the one favored in the bulk liquid. Here, we study the formation of liquid bridges induced by confining colloids to a slit, with the slit walls having a preference opposite to the one of the colloid surface. Using mean field theory, we show that there is a line of first-order phase transitions between the bridge and the no-bridge states, which ends at a critical point. By decreasing the slit width, this critical point is shifted toward smaller separations between the colloids. However, at very small separations and far from criticality, we observe only a minor influence of the slit width on the location of the transition. Monte Carlo simulations of the Ising model, which mimics incompressible binary liquid mixtures, confirm the occurrence of the bridging transitions, as manifested by the appearance of "spinodal" regions where both bridge and no-bridge configurations are stable or metastable. Interestingly, we find that there is no such spinodal region in the case of small colloids, but we observe a sharpening of the transition when the colloid size increases. In addition, we demonstrate that the capillary force acting between the colloids can depend sensitively on the slit width and varies drastically with temperature, thus achieving strengths orders of magnitude higher than at criticality of the fluid.

Ab initio nanofluidics: disentangling the role of the energy landscape and of density correlations on liquid/solid friction

Despite relevance to water purification and renewable energy conversion membranes, the molecular mechanisms underlying water slip are poorly understood. We disentangle the static and dynamical origin of water slippage on graphene, hBN and MoS2 by means of large-scale ab initio molecular dynamics. Accounting for the role of the electronic structure of the interface is essential to determine that water slips five and eleven times faster on graphene compared to hBN and to MoS2, respectively. Intricate changes in the water energy landscape as well as in the density correlations of the fluid provide, respectively, the main static and dynamical origin of water slippage. Surprisingly, the timescales of the density correlations are the same on graphene and hBN, whereas they are longer on MoS2 and yield a 100% slowdown in the flow of water on this material. Our results pave the way for an in silico first principles design of materials with enhanced water slip, through the modification of properties connected not only to the structure, but also to the dynamics of the interface.

Understanding non-covalent interactions in larger molecular complexes from first principles

Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.

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.

Unexpectedly High Salt Accumulation inside Carbon Nanotubes Soaked in Dilute Salt Solutions

We experimentally demonstrate the formation of salt aggregations with unexpectedly high concentration inside multiwalled carbon nanotubes (CNTs) soaked only in dilute salt solution sand even in solutions containing only traces of salts. This finding suggests the blocking of fluid across CNTs by the salt aggregations when CNTs are soaked in a dilute salt solution with the concentration of seawater or even lower, which may open new avenues for the development of novel CNT-based desalination techniques. The high salt accumulation of CNTs also provides a new CNT-based strategy for the collection or extraction of noble metal salts in solutions containing traces of noble metal salts. Theoretical analyses reveal that this high salt accumulation inside CNTs can be mainly attributed to the strong hydrated cation-π interactions of hydrated cations and π electrons in the aromatic rings of CNTs.

Quantum-Mechanical Relation between Atomic Dipole Polarizability and the van der Waals Radius

The atomic dipole polarizability α and the van der Waals (vdW) radius R_{vdW} are two key quantities to describe vdW interactions between atoms in molecules and materials. Until now, they have been determined independently and separately from each other. Here, we derive the quantum-mechanical relation R_{vdW}=const×α^{1/7}, which is markedly different from the common assumption R_{vdW}∝α^{1/3} based on a classical picture of hard-sphere atoms. As shown for 72 chemical elements between hydrogen and uranium, the obtained formula can be used as a unified definition of the vdW radius solely in terms of the atomic polarizability. For vdW-bonded heteronuclear dimers consisting of atoms A and B, the combination rule α=(α_{A}+α_{B})/2 provides a remarkably accurate way to calculate their equilibrium interatomic distance. The revealed scaling law allows us to reduce the empiricism and improve the accuracy of interatomic vdW potentials, at the same time suggesting the existence of a nontrivial relation between length and volume in quantum systems.

Probing interface localization-delocalization transitions by colloids

Interface localization-delocalization transitions (ILDT) occur in two-phase fluids confined in a slit with competing preferences of the walls for the two fluid phases. At low temperatures the interface between the two phases is localized at one of the walls. Upon increasing temperature it unbinds. Although intensively studied theoretically and computationally, such transitions have not yet been observed experimentally due to severe challenges in resolving fine details of the fluid structure. Here, using mean field theory and Monte Carlo simulations of the Ising model, we propose to detect these ILDT by using colloids. We show that the finite-size and fluctuation induced force acting on a colloid confined in such a system experiences a vivid change if, upon lowering the temperature, the interface localizes at one of the walls. This change can serve as a more easily accessible experimental indicator of the transition.

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.

Nonadditive interactions and phase transitions in strongly confined colloidal systems

The behaviour of colloids can be controlled effectively by tuning the solvent-mediated interactions among them. An extensively studied example is the temperature-induced aggregation of suspended colloids close to the consolute point of their binary solvent. Here, using mean field theory and Monte Carlo simulations, we study the behaviour of colloids confined to a narrow slit containing a nearly-critical binary liquid mixture. We found that the effective interactions in this system are highly non-additive. In particular, the effective interactions among the colloids can be a few times stronger than the corresponding sum of the effective pair potentials. Inter alia, this non-additivity manifests itself in the phase behaviour of confined colloids, which depends sensitively on the slit width and temperature. In addition, we demonstrate the possibility of a first-order bridging transition between colloids confined to a slit and suspended in a phase-separated fluid well below the critical point of the solvent and at its critical composition in the bulk. This transition is accompanied by a remarkably large hysteresis loop, in which the force between the colloids varies by two orders of magnitude.