Source author record

Per Hyldgaard

Per Hyldgaard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Stability of and conduction in single-walled Si$_2$BN nanotubes

We explore the possibility and potential benefit of rolling a Si2BN sheet into single-walled nanotubes (NTs). Using density functional theory (DFT), we consider both structural stability and the impact on the nature of chemical bonding and conduction. The structure is similar to carbon NTs and hexagonal boron-nitride (hBN) NTs and we consider both armchair and zigzag Si2BN configurations with varying diameters. The stability of these Si$_2$BN NTs is confirmed by first-principles molecular dynamics calculations, by an exothermal formation, an absence of imaginary modes in the phonon spectra. Also, we find the nature of conduction varies semiconducting, from semi-metallic to metallic, reflecting differences in armchair/zigzag-type structures, curvature effects, and the effect of quantum confinement. We present the detailed characterization of how these properties lead to differences in both the bonding nature and electronic structures

preprint2022arXiv

vdW-DF-ahcx: a range-separated van der Waals density functional hybrid

Hybrid density functionals replace a fraction of an underlying generalized-gradient approximation (GGA) exchange description with a Fock-exchange component. Range-separated hybrids (RSHs) also effectively screen the Fock-exchange component and thus open the door for characterizations of metals and adsorption at metal surfaces. The RSHs are traditionally based on a robust GGA, such as PBE [PRL $\textbf{77}$, 3865 (1996)], for example, as implemented in the HSE design [JPC $\textbf{118}$, 8207 (2003)]. Here we define a RSH extension to the van der Waals density functional (vdW-DF) method [ROPP $\textbf{78}$, 066501 (2015)], launching vdW-DF-ahcx. We use an analytical-hole (AH) framework [JCP $\textbf{128}$, 194105 (2008)] to characterize the GGA-type exchange in the vdW-DF-cx version [PRB $\textbf{89}$, 075148 (2014)], isolate the short-ranged exchange component, and define the new RSH. We find that the performance vdW-DF-ahcx compares favorably to (dispersion-corrected) HSE for descriptions of bulk (broad molecular) properties. We also find that it provides accurate descriptions of noble-metal surface properties, including CO adsorption.

preprint2017arXiv

libvdwxc: A library for exchange-correlation functionals in the vdW-DF family

We present libvdwxc, a general library for evaluating the energy and potential for the family of vdW-DF exchange--correlation functionals. libvdwxc provides an efficient implementation of the vdW-DF method and can be interfaced with various general-purpose DFT codes. Currently, the GPAW and Octopus codes implement interfaces to libvdwxc. The present implementation emphasizes scalability and parallel performance, and thereby enables \textit{ab initio} calculations of nanometer-scale complexes. The numerical accuracy is benchmarked on the S22 test set whereas parallel performance is benchmarked on ligand-protected gold nanoparticles ($\text{Au}_{144}(\text{SC}_{11}\text{NH}_{25})_{60}$) up to 9696 atoms.

preprint2016arXiv

Finite-temperature properties of non-magnetic transition metals: Comparison of the performance of constraint-based semi and nonlocal functionals

We assess the performance of nonempirical, truly nonlocal and semi-local functionals with regard to structural and thermal properties of $3d$, $4d$, and $5d$ non-magnetic transition metals. We focus on constraint-based functionals and consider the new consistent-exchange van der Waals density functional version vdW-DF-cx [Phys. Rev. B 89, 035412 (2014)], the semi-local PBE [Phys. Rev. Lett. 77, 3865 (1996)] and PBEsol functionals [Phys. Rev. Lett. 100, 136406 (2008)] as well as the AM05 meta-functional [Phys. Rev. B 72, 085108 (2005)]. Using the quasi-harmonic approximation structural parameters, elastic response, and thermal expansion at finite temperatures are computed and compared to experimental data. We also compute cohesive energies explicitly including zero-point vibrations. It is shown that overall vdW-DF-cx provides an accurate description of thermal properties and retains a level of transferability and accuracy that is comparable to or better than some of the best constraint-based semi-local functionals. Especially, with regard to the cohesive energies the consistent inclusion of spin polarization effects in the atoms turns out to be crucial and it is important to use the rigorous spin-vdW-DF-cx formulation [Phys. Rev. Lett. 115, 136402 (2015)]. This demonstrates that vdW-DF-cx has general-purpose character and can be used to study systems that have both sparse and dense electron distributions.

preprint2016arXiv

Structural and excited-state properties of oligoacene crystals from first principles

Molecular crystals are a prototypical class of van der Waals (vdW) bound organic materials with excited state properties relevant for optoelectronics applications. Predicting the structure and excited state properties of molecular crystals presents a challenge for electronic structure theory, as standard approximations to density functional theory (DFT) do not capture long range vdW dispersion interactions and do not yield excited state properties. In this work, we use a combination of DFT including vdW forces) using both non local correlation functionals and pair wise correction methods (together with many body perturbation theory (MBPT) to study the geometry and excited states, respectively, of the entire series of oligoacene crystals, from benzene to hexacene. We find that vdW methods can predict lattice constants within 1 percent of the experimental measurements, on par with the previously reported accuracy of pairwise approximations for the same systems. We further find that excitation energies are sensitive to geometry, but if optimized geometries are used MBPT can yield excited state properties within a few tenths of an eV from experiment. We elucidate trends in MBPT computed charged and neutral excitation energies across the acene series and discuss the role of common approximations used in MBPT.

preprint2014arXiv

Interpretation of van der Waals density functionals

The nonlocal correlation energy in the van der Waals density functional (vdW-DF) method [Phys. Rev. Lett. 92, 246401 (2004); Phys. Rev. B 76, 125112 (2007); Phys. Rev. B 89, 035412 (2014)] can be interpreted in terms of a coupling of zero-point energies of characteristic modes of semilocal exchange-correlation (xc) holes. These xc holes reflect the internal functional in the framework of the vdW-DF method [Phys. Rev. B 82, 081101(2010)]. We explore the internal xc hole components, showing that they share properties with those of the generalized-gradient approximation. We use these results to illustrate the nonlocality in the vdW-DF description and analyze the vdW-DF formulation of nonlocal correlation.

preprint2014arXiv

van der Waals density functionals built upon the electron-gas tradition: Facing the challenge of competing interactions

The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory (DFT). One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO$_3$, the adsorption of small molecules within metal-organic frameworks (MOFs), the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general purpose functional that could be applied to a range of materials problems with a variety of competing interactions.

preprint2014arXiv

van der Waals forces in density functional theory: The vdW-DF method

A density functional theory (DFT) that accounts for van der Waals (vdW) interactions in condensed matter, materials physics, chemistry, and biology is reviewed. The insights that led to the construction of the Rutgers-Chalmers van der Waals Density Functional (vdW-DF) are presented with the aim of giving a historical perspective, while also emphasising more recent efforts which have sought to improve its accuracy. In addition to technical details, we discuss a range of recent applications that illustrate the necessity of including dispersion interactions in DFT. This review highlights the value of the vdW-DF method as a general-purpose method, not only for dispersion bound systems, but also in densely packed systems where these types of interactions are traditionally thought to be negligible.

preprint2013arXiv

A van der Waals density functional mapping of attraction in DNA dimers

The dispersion interaction between a pair of parallel DNA double-helix structures is investigated by means of the van der Waals density functional (vdW-DF) method. Each double-helix structure consists of an infinite repetition of one B-DNA coil with 10 base pairs. This parameter-free density functional theory (DFT) study illustrates the initial step in a proposed vdW-DF computational strategy for large biomolecular problems. The strategy is to first perform a survey of interaction geometries, based on the evaluation of the van der Waals (vdW) attraction, and then limit the evaluation of the remaining DFT parts (specifically the expensive study of the kinetic-energy repulsion) to the thus identified interesting geometries. Possibilities for accelerating this second step is detailed in a separate study. For the B-DNA dimer, the variation in van der Waals attraction is explored at relatively short distances (although beyond the region of density overlap) for a 360 degrees rotation. This study highlights the role of the structural motifs, like the grooves, in enhancing or reducing the vdW interaction strength. We find that to a first approximation, it is possible to compare the DNA double strand at large wall-to-wall separations to the cylindrical shape of a carbon nanotube (which is almost isotropic under rotation). We compare our first-principles results with the atom-based dispersive interaction predicted by DFT-D2 [J. Comp. Chem. 27, 1787 (2006)] and find agreement in the asymptotic region. However, we also find that the differences in the enhancement that occur at shorter distances reveal characteristic features that result from the fact that the vdW-DF method is an electron-based (as opposed to atom-based) description.

preprint2013arXiv

An analysis of van der Waals density functional components: Binding and corrugation of benzene and C60 on boron nitride and graphene

The adsorption of benzene and C60 on graphene and boron nitride (BN) is studied using density functional theory with the non-local correlation functional vdW-DF. By comparing these systems we can systematically investigate their adsorption nature and differences between the two functional versions vdW-DF1 and vdW-DF2. The bigger size of the C60 molecule makes it bind stronger to the surface than benzene, yet the interface between the molecules and the sheets are similar in nature. The binding separation is more sensitive to the exchange variant used in vdW-DF than to the correlation version. This result is related to the exchange and correlation components of the potential energy curve (PEC). We show that a moderate dipole forms for C60 on graphene, unlike for the other adsorption systems. We find that the corrugation is very sensitive to the variant or version of vdW-DF used, in particular the exchange. Further, we show that this sensitivity arise indirectly through the shift in binding separation caused by changing vdW-DF variant. Based on our results, we suggest a concerted theory-experiment approach to assess the exchange and correlation contributions to physisorption. Using DFT calculations, the corrugation can be linked to the optimal separation, allowing us to extract the exchange-correlation part of the adsorption energy. Molecules with same interfaces to the surface, but different geometries, can in turn cast light on the role of van der Waals forces.

preprint2013arXiv

An exchange functional that tests the robustness of the plasmon description of the van der Waals density functional

Is the plasmon description within the non-local correlation of the van der Waals density functional by Dion and coworkers (vdW-DF1) robust enough to describe all exchange-correlation components? To address this question, we design an exchange functional, termed LV-PW86r based on this plasmon description as well as recent analysis on exchange in the large $s$-regime. In the regime with reduced gradients $s=|\nabla n|/2n k_{\rm F}(n)$ smaller than $\approx 2.5$, dominating the non-local correlation part of the binding energy, the enhancement factor $F_x(s)$ closely resembles the Langreth-Vosko screened exchange. In the $s$-regime beyond, dominated by exchange, $F_x(s)$ passes smoothly over to the revised Perdew-Wang-86 form. Our tests indicate that vdW-DF1(LV-PW86r) produces accurate separations and binding energies of the S22 data set of molecular dimers as well as accurate lattice constants of layered materials and tightly-bound solids. These results suggest that vdW-DF1 has a good plasmon description in the low-to-moderate $s$-regime.

preprint2012arXiv

Benchmarking van der Waals Density Functionals with Experimental Data: Potential Energy Curves for H2 Molecules on Cu(111), (100), and (110) Surfaces

Detailed physisorption data from experiment for the H_2 molecule on low-index Cu surfaces challenge theory. Recently, density-functional theory (DFT) has been developed to account for nonlocal correlation effects, including van der Waals (dispersion) forces. We show that the functional vdW-DF2 gives a potential-energy curve, potential-well energy levels, and difference in lateral corrugation promisingly close to the results obtained by resonant elastic backscattering-diffraction experiments. The backscattering barrier is found selective for choice of exchange-functional approximation. Further, the DFT-D3 and TS-vdW corrections to traditional DFT formulations are also benchmarked, and deviations are analyzed.

preprint2012arXiv

Nonequilibrium thermodynamics of interacting tunneling transport: variational grand potential, density-functional formulation, and nature of steady-state forces

The standard formulation of tunneling transport rests on an open-boundary modeling. There, conserving approximations to nonequilibrium Green function or quantum-statistical mechanics provide consistent but computational costly approaches; alternatively, use of density-dependent ballistic-transport calculations [e.g., Phys. Rev. B 52, 5335 (1995)], here denoted `DBT', provide computationally efficient (approximate) atomistic characterizations of the electron behavior but has until now lacked a formal justification. This paper presents an exact, variational nonequilibrium thermodynamic theory for fully interacting tunneling and provides a rigorous foundation for frozen-nuclei DBT calculations as a lowest order approximation to an exact nonequilibrium thermodynamics density functional evaluation. The theory starts from the complete electron nonequilibrium quantum statistical mechanics and I identify the operator for the nonequilibrium Gibbs free energy. I demonstrate a minimal property of a functional for the nonequilibrium thermodynamic grand potential which thus uniquely identifies the solution as the exact nonequilibrium density matrix. I also show that a uniqueness-of-density proof from a closely related study [Phys. Rev. B 78, 165109 (2008)] makes it possible to provide a single-particle formulation based on universal electron-density functionals. I illustrate a formal evaluation of the thermodynamics grand potential value which is closely related to the variation in scattering phase shifts and hence to Friedel density oscillations. This paper also discusses the difference between the here-presented exact thermodynamics forces and the often-used electrostatic forces. Finally the paper documents an inherent adiabatic nature of the thermodynamics forces and observes that these are suited for a nonequilibrium implementation of the Born-Oppenheimer approximation.

preprint2011arXiv

Computational scheme for ab-initio predictions of chemical compositions interfaces realized by deposition growth

We present a novel computational scheme to predict chemical compositions at interfaces as they emerge in a growth process. The scheme uses the Gibbs free energy of reaction associated with the formation of interfaces with a specific composition as predictor for their prevalence. It explicitly accounts for the growth conditions by rate-equation modeling of the deposition environment. The Bell-Evans-Polanyi principle motivates our emphasis on an effective nonequilibrium thermodynamic description inspired by chemical reaction theory. We illustrate the scheme by characterizing the interface between TiC and alumina. Equilibrium thermodynamics favors a nonbinding interface, being in conflict with the wear-resistant nature of TiC/alumina multilayer coatings. Our novel scheme predicts that deposition of a strongly adhering interface is favored under realistic conditions.

preprint2011arXiv

Evaluation of New Density Functional with Account of van der Waals Forces by Use of Experimental H2 Physisorption Data on Cu(111)

Detailed experimental data for physisorption potential-energy curves of H2 on low-indexed faces of Cu challenge theory. Recently, density-functional theory has been developed to also account for nonlocal correlation effects, including van der Waals forces. We show that one functional, denoted vdW-DF2, gives a potential-energy curve promisingly close to the experiment-derived physisorptionenergy curve. The comparison also gives indications for further improvements of the functionals.

preprint2011arXiv

Polarization-balanced design of AlN/GaN heterostructures: Application to double-barrier structures

Inversion- and depletion-regions generally form at the interfaces between doped leads (cladding layers) and the active region of polar heterostructures like AlN/GaN and other nitride compounds. The band bending in the depletion region sets up a barrier which may seriously impede perpendicular electronic transport. This may ruin the performance of devices such as quantum-cascade lasers and resonant-tunneling diodes. Here we introduce the concepts of polarization balance and polarization-balanced designs: A structure is polarization balanced when the applied bias match the voltage drop arising from spontaneous and piezeolectric fields. Devices designed to operate at this bias have polarization-balanced designs. These concepts offer a systematic approach to avoid the formation of depletion regions. As a test case, we consider the design of AlN/GaN double barrier structures with Al$_{\tilde{x}}$Ga$_{1-\tilde{x}}$N leads. To guide our efforts, we derive a simple relation between the intrinsic voltage drop arising from polar effects, average alloy composition of the active region, and the alloy concentration of the leads. Polarization-balanced designs secure good filling of the ground state for unbiased structures, while for biased structures with efficient emptying of the active structure it removes the depletion barriers.

preprint2011arXiv

Relative stability of $6H$-SiC$\{0001\}$ surface terminations and formation of graphene overlayers by Si evaporation

We present density functional theory (DFT) calculations for 6H-SiC$\{0001\}$ surfaces with different surface stackings and terminations. We compare the relative stability of different $(0001)$ and $(000\bar1)$ surfaces in terms of their surface free energies. Removing surface and subsurface Si atoms, we simulate the formation of graphene and graphene-like overlayers by Si evaporation. We find that overlayers with a different nature of bonding are preferred at the two non-equivalent surface orientations. At $(0001)$, a chemically bonded, highly strained and buckled film is predicted. At $(000\bar1)$, a van der Waals (vdW) bonded overlayer is preferred. We quantify the vdW binding and show that it can have a doping effect on electron behavior in the overlayer.

preprint2011arXiv

Stacking and band structure of van der Waals bonded graphane multilayers

We use density functional theory and the van der Waals density functional (vdW-DF) method to determine the binding separation in bilayer and bulk graphane and study the changes in electronic band structure that arise with the multilayer formation. The calculated binding separation (distance between center-of-mass planes) and binding energy are 4.5-5.0 Å (4.5-4.8 Å) and 75-102 meV/cell (93-127 meV/cell) in the bilayer (bulk), depending on the choice of vdW-DF version. We obtain the corresponding band diagrams using calculations in the ordinary generalized gradient approximation for the geometries specified by our vdW-DF results, so probing the indirect effect of vdW forces on electron behavior. We find significant band-gap modifications by up to -1.2 eV (+4.0 eV) in various regions of the Brillouin zone, produced by the bilayer (bulk) formation.

preprint2011arXiv

Surface assembly and ultrafast operation of all-nanoscale resonant-tunneling transistors

Realization of a robust nanotube-heterostructure tunneling transistors [Solid State Comm. 116, p. 569 (2000)] requires the difficult formation [Science 293, p. 76 (2001)] of a central nanoscale barrier separating a pair of outside metallic leads. Here I suggest an alternative surface-based assembly based on self-organization of one-dimensional metallic states on oxides and trapping well-resolved resonant orbitals in a central island. I present and explain the (universal) transistor characteristics and robustness. In addition, I calculate typical the island/level-to-gate capacitance to predict ultrafast (beyond-THz) switching but also document a limited importance of Coulomb blockade effects in the (nanotube) resonant-tunneling transistors.

preprint2010arXiv

Ab-initio thermodynamics of deposition growth: surface terminations of CVD titanium carbide and nitride

We present a calculational method to predict terminations of growing or as-deposited surfaces as a function of the deposition conditions. Such characterizations are valuable for understanding catalysis and growth phenomena. The method combines ab-initio density functional theory (DFT) calculations and experimental thermodynamical data with a rate-equations description of partial pressures in the reaction chamber. The use of rate equations enables a complete description of a complex gas environment in terms of a few, (experimentally accessible) parameters. The predictions are based on comparisons between free energies of reaction associated with the formation of surfaces with different terminations. The method has an intrinsic non-equilibrium character. In the limit of dynamic equilibrium (with equal chemical potential in the surface and the gas phase) we find that the predictions of the method coincide with those of standard ab-initio based equilibrium thermodynamics. We illustrate the method for chemical vapor deposition (CVD) of TiC(111) and TiN(111), and find that the emerging termination can be controlled both by the environment and the growth rate.

preprint2010arXiv

Structure and binding in crystals of cage-like molecules: hexamine and platonic hydrocarbons

In this paper, we show that first-principle calculations using a van der Waals density functional (vdW-DF), [Phys. Rev. Lett. $\mathbf{92}$, 246401 (2004)] permits determination of molecular crystal structure. We study the crystal structures of hexamine and the platonic hydrocarbons (cubane and dodecahedrane). The calculated lattice parameters and cohesion energy agree well with experiments. Further, we examine the asymptotic accounts of the van der Waals forces by comparing full vdW-DF with asymptotic atom-based pair potentials extracted from vdW-DF. The character of the binding differ in the two cases, with vdW-DF giving a significant enhancement at intermediate and relevant binding separations. We analyze consequences of this result for methods such as DFT-D, and question DFT-D's transferability over the full range of separations.

preprint2010arXiv

Temperature stability of intersubband transitions in AlN/GaN quantum wells

Temperature dependence of intersubband transitions in AlN/GaN multiple quantum wells grown with molecular beam epitaxy is investigated both by absorption studies at different temperatures and modeling of conduction-band electrons. For the absorption study, the sample is heated in increments up to $400^\circ$C. The self-consistent Schrödinger-Poisson modeling includes temperature effects of the band-gap and the influence of thermal expansion on the piezoelectric field. We find that the intersubband absorption energy decreases only by $\sim 6$ meV at $400^\circ$C relative to its room temperature value.

preprint2010arXiv

Understanding adhesion at as-deposited interfaces from ab initio thermodynamics of deposition growth: thin-film alumina on titanium carbide

We investigate the chemical composition and adhesion of chemical vapour deposited thin-film alumina on TiC using and extending a recently proposed nonequilibrium method of ab initio thermodynamics of deposition growth (AIT-DG) [Rohrer J and Hyldgaard P 2010 Phys. Rev. B 82 045415]. A previous study of this system [Rohrer J, Ruberto C and Hyldgaard P 2010 J. Phys.: Condens. Matter 22 015004] found that use of equilibrium thermodynamics leads to predictions of a non-binding TiC/alumina interface, despite the industrial use as a wear-resistant coating. This discrepancy between equilibrium theory and experiment is resolved by the AIT-DG method which predicts interfaces with strong adhesion. The AIT-DG method combines density functional theory calculations, rate-equation modelling of the pressure evolution of the deposition environment and thermochemical data. The AIT-DG method was previously used to predict prevalent terminations of growing or as-deposited surfaces of binary materials. Here we extent the method to predict surface and interface compositions of growing or as-deposited thin films on a substrate and find that inclusion of the nonequilibrium deposition environment has important implications for the nature of buried interfaces.

preprint2009arXiv

Ab initio structure modeling of complex thin-film oxides: thermodynamical stability of TiC/thin-film alumina

We present an efficient and general method to identify promising candidate configurations for thin-film oxides and to determine structural characteristics of (metastable) thin-film structures using ab initio calculations. At the heart of this method is the complexity of the oxide bulk structure, from which a large number of thin films with structural building blocks, that is motifs, from metastable bulk oxide systems can be extracted. These span a coarse but well-defined network of initial configurations for which density functional theory (DFT) calculations predict and implement dramatic atomic relaxations in the corresponding, resulting thin-film candidates. The network of thin-film candidates (for various film thicknesses and stoichiometries) can be ordered according to their variation in ab initio total energy or in ab initio equilibrium Gibbs free energy. Analysis of the relaxed atomic structures for the most favored structures gives insight into the nature of stable and metastable thin-film oxides. We investigate ultrathin alumina nucleated on TiC as a model system to illustrate this method.

preprint2007arXiv

Van der Waals density functional: Self-consistent potential and the nature of the van der Waals bond

We derive the exchange-correlation potential corresponding to the nonlocal van der Waals density functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)]. We use this potential for a self-consistent calculation of the ground state properties of a number of van der Waals complexes as well as crystalline silicon. For the latter, where little or no van der Waals interaction is expected, we find that the results are mostly determined by semilocal exchange and correlation as in standard generalized gradient approximations (GGA), with the fully nonlocal term giving little effect. On the other hand, our results for the van der Waals complexes show that the self-consistency has little effect at equilibrium separations. This finding validates previous calculations with the same functional that treated the fully nonlocal term as a post GGA perturbation. A comparison of our results with wave-function calculations demonstrates the usefulness of our approach. The exchange-correlation potential also allows us to calculate Hellmann-Feynman forces, hence providing the means for efficient geometry relaxations as well as unleashing the potential use of other standard techniques that depend on the self-consistent charge distribution. The nature of the van der Waals bond is discussed in terms of the self-consistent bonding charge.