Screening and Antiscreening Effects in Endohedral Nanotubes.

Recently we investigated from first-principles screening properties in systems where small molecules, characterized by a finite electronic dipole moment, are encapsulated in different nanocages. The most relevant result was the observation of an antiscreening effect in alkali-halide nanocages characterized by ionic bonds: in fact, due to the relative displacement of positive and negative ions, induced by the dipole moment of the encapsulated molecule, these cages act as dipole-field amplifiers, different from what is observed in carbon fullerene nanocages, which exhibit instead a pronounced screening effect. Here we extend the study to another class of nanostructures: the nanotubes. Using first-principles techniques based on density functional theory, we studied the properties of endohedral nanotubes obtained by encapsulation of a water molecule or a linear HF molecule. A detailed analysis of the effective dipole moment of the complexes and of the electronic charge distribution suggests that screening effects crucially depend not only on the nature of the intramolecular bonds but also on the size and the shape of the nanotubes and on the specific encapsulated molecule. As observed in endohedral nanocages, screening is maximum in covalent-bond carbon nanotubes, while it is reduced in partially ionic nanotubes, and an antiscreening effect is observed in some ionic nanotubes. However, in this case, the scenario is more complex than in corresponding ionic nanocages. In fact the specific geometric structure of alkali-halide nanotubes turns out to be crucial for determining the screening/antiscreening behavior: while nanotubes with octagonal transversal section can exhibit an antiscreening effect, which quantitatively depends on the number of layers in the longitudinal direction, instead nanotubes with dodecagonal section are always characterized by a reduction of the total dipole moment so that a screening behavior is observed. Our results show that, even in nanotube structures, in principle one can tune the dipole moment and generate electrostatic fields at the nanoscale without the aid of external potentials.

Superfluidity Meets the Solid State: Frictionless Mass Transport through a (5,5) Carbon Nanotube.

Superfluidity is a well-characterized quantum phenomenon which entails frictionless motion of mesoscopic particles through a superfluid, such as ^{4}He or dilute atomic gases at very low temperatures. As shown by Landau, the incompatibility between energy and momentum conservation, which ultimately stems from the spectrum of the elementary excitations of the superfluid, forbids quantum scattering between the superfluid and the moving mesoscopic particle, below a critical speed threshold. Here, we predict that frictionless motion can also occur in the absence of a standard superfluid, i.e., when a He atom travels through a narrow (5,5) carbon nanotube (CNT). Because of the quasilinear dispersion of the plasmon and phonon modes that could interact with He, the (5,5) CNT embodies a solid-state analog of the superfluid, thereby enabling straightforward transfer of Landau's criterion of superfluidity. As a result, Landau's equations acquire broader generality and may be applicable to other nanoscale friction phenomena, whose description has been so far purely classical.

Transport properties in liquids from first-principles: The case of liquid water and liquid argon.

Shear and bulk viscosities of liquid water and argon are evaluated from first-principles in the density functional theory (DFT) framework, by performing molecular dynamics simulations in the NVE ensemble and using the Kubo-Greenwood equilibrium approach. The standard DFT functional is corrected in such a way to allow for a reasonable description of van der Waals effects. For liquid argon, the thermal conductivity has been also calculated. Concerning liquid water, to our knowledge, this is the first estimate of the bulk viscosity and of the shear-viscosity/bulk-viscosity ratio from first-principles. By analyzing our results, we can conclude that our first-principles simulations, performed at a nominal average temperature of 366 to guarantee that the systems are liquid-like, actually describe the basic dynamical properties of liquid water at about 330 K. In comparison with liquid water, the normal, monatomic liquid Ar is characterized by a much smaller bulk-viscosity/shear-viscosity ratio (close to unity) and this feature is well reproduced by our first-principles approach, which predicts a value of the ratio in better agreement with experimental reference data than that obtained using the empirical Lennard-Jones potential. The computed thermal conductivity of liquid argon is also in good agreement with the experimental value.

Quantum-mechanical water-flow enhancement through a sub-nanometer carbon nanotube.

Experimental observations unambiguously reveal quasi-frictionless water flow through nanometer-scale carbon nanotubes (CNTs). Classical fluid mechanics is deemed unfit to describe this enhanced flow, and recent investigations indicated that quantum mechanics is required to interpret the extremely weak water-CNT friction. In fact, by quantum scattering, water can only release discrete energy upon excitation of electronic and phononic modes in the CNT. Here, we analyze in detail how a traveling water molecule couples to both plasmon and phonon excitations within a sub-nanometer, periodic CNT. We find that the water molecule needs to exceed a minimum speed threshold of ∼50 m/s in order to scatter against CNT electronic and vibrational modes. Below this threshold, scattering is suppressed, as in standard superfluidity mechanisms. The scattering rates, relevant for faster water molecules, are also estimated.

Many-body van der Waals interactions beyond the dipole approximation.

Long-ranged van der Waals (vdW) interactions are most often treated via Lennard-Jones approaches based on the combination of two-body and dipolar approximations. While beyond-dipole interactions and many-body contributions were separately addressed, little is known about their combined effect, especially in large molecules and relevant nanoscale systems. Here, we provide a full many-body description of vdW interactions beyond the dipole approximation, efficiently applicable to large-scale systems. Dipole-quadrupole interactions consistently exhibit large magnitude up to nm-scale separations, while many-body effects lead to system-dependent screening effects, which can reduce vdW interactions by a large fraction. Combined many-body and multipolar terms emerge as an essential ingredient for the reliable description of vdW interactions in molecular and nanoscale systems.

Tunable van der Waals interactions in low-dimensional nanostructures.

Non-covalent van der Waals interactions play a major role at the nanoscale, and even a slight change in their asymptotic decay could produce a major impact on surface phenomena, self-assembly of nanomaterials, and biological systems. By a full many-body description of vdW interactions in coupled carbyne-like chains and graphenic structures, here, we demonstrate that both modulus and a range of interfragment forces can be effectively tuned, introducing mechanical strain and doping (or polarizability change). This result contrasts with conventional pairwise vdW predictions, where the two-body approximation essentially fixes the asymptotic decay of interfragment forces. The present results provide viable pathways for detailed experimental control of nanoscale systems that could be exploited both in static geometrical configurations and in dynamical processes.

van der Waals interactions in DFT using Wannier functions without empirical parameters.

A new implementation is proposed for including van der Waals (vdW) interactions in Density Functional Theory (DFT) using the Maximally Localized Wannier Functions (MLWFs), which is free from empirical parameters. With respect to the previous DFT/vdW-WF2 method, in the present DFT/vdW-WF2-x approach, the empirical, short-range, damping function is replaced by an estimate of the Pauli exchange repulsion, also obtained by the MLWF properties. Applications to systems contained in the popular S22 molecular database and to the case of an Ar atom interacting with graphite and comparison with reference data indicate that the new method, besides being more physically founded, also leads to a systematic improvement in the description of vdW-bonded systems.

Anomalous van der Waals-Casimir interactions on graphene: A concerted effect of temperature, retardation, and non-locality.

Dispersion forces play a major role in graphene, largely influencing adhesion of adsorbate moieties and stabilization of functional multilayered structures. However, the reliable prediction of dispersion interactions on graphene up to the relevant ∼10 nm scale is an extremely challenging task: in fact, electromagnetic retardation effects and the highly non-local character of π electrons can imply sizeable qualitative variations of the interaction with respect to known pairwise approaches. Here we address both issues, determining the finite-temperature van der Waals (vdW)-Casimir interaction for point-like and extended adsorbates on graphene, explicitly accounting for the non-local dielectric permittivity. We find that temperature, retardation, and non-locality play a crucial role in determining the actual vdW scaling laws and the stability of both atomic and larger molecular adsorbates. Our results highlight the importance of these effects for a proper description of systems of current high interest, such as graphene interacting with biomolecules, and self-assembly of complex nanoscale structures. Due to the generality of our approach and the observed non-locality of other 2D materials, our results suggest non-trivial vdW interactions from hexagonal mono-layered materials from group 14 of the periodic table, to transition metal dichalcogenides.

Hydrogen bonding characterization in water and small molecules.

The prototypical hydrogen bond in water dimer and hydrogen bonds in the protonated water dimer, in other small molecules, in water cyclic clusters, and in ice, covering a wide range of bond strengths, are theoretically investigated by first-principles calculations based on density functional theory, considering not only a standard generalized gradient approximation functional but also, for the water dimer, hybrid and van der Waals corrected functionals. We compute structural, energetic, and electrostatic (induced molecular dipole moments) properties. In particular, hydrogen bonds are characterized in terms of differential electron density distributions and profiles, and of the shifts of the centres of maximally localized Wannier functions. The information from the latter quantities can be conveyed to a single geometric bonding parameter that appears to be correlated with the Mayer bond order parameter and can be taken as an estimate of the covalent contribution to the hydrogen bond. By considering the water trimer, the cyclic water hexamer, and the hexagonal phase of ice, we also elucidate the importance of cooperative/anticooperative effects in hydrogen-bonding formation.

Communication: Enhanced chemical reactivity of graphene on a Ni(111) substrate.

Due to the unique combination of structural, mechanical, and transport properties, graphene has emerged as an exceptional candidate for catalysis applications. The low chemical reactivity caused by sp(2) hybridization and strongly delocalized π electrons, however, represents a main challenge for straightforward use of graphene in its pristine, free-standing form. Following recent experimental indications, we show that due to charge hybridization, a Ni(111) substrate can enhance the chemical reactivity of graphene, as exemplified by the interaction with the CO molecule. While CO only physisorbs on free-standing graphene, chemisorption of CO involving formation of ethylene dione complexes is predicted in Ni(111)-graphene. Higher chemical reactivity is also suggested in the case of oxidized graphene, opening the way to a simple and efficient control of graphene chemical properties, devoid of complex defect patterning or active metallic structures deposition.

Including screening in van der Waals corrected density functional theory calculations: the case of atoms and small molecules physisorbed on graphene.

The Density Functional Theory (DFT)/van der Waals-Quantum Harmonic Oscillator-Wannier function (vdW-QHO-WF) method, recently developed to include the vdW interactions in approximated DFT by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique, is applied to the cases of atoms and small molecules (X=Ar, CO, H2, H2O) weakly interacting with benzene and with the ideal planar graphene surface. Comparison is also presented with the results obtained by other DFT vdW-corrected schemes, including PBE+D, vdW-DF, vdW-DF2, rVV10, and by the simpler Local Density Approximation (LDA) and semilocal generalized gradient approximation approaches. While for the X-benzene systems all the considered vdW-corrected schemes perform reasonably well, it turns out that an accurate description of the X-graphene interaction requires a proper treatment of many-body contributions and of short-range screening effects, as demonstrated by adopting an improved version of the DFT/vdW-QHO-WF method. We also comment on the widespread attitude of relying on LDA to get a rough description of weakly interacting systems.

Van der Waals interactions in density functional theory by combining the quantum harmonic oscillator-model with localized Wannier functions.

We present a new scheme to include the van der Waals (vdW) interactions in approximated Density Functional Theory (DFT) by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique. With respect to the recently developed DFT/vdW-WF2 method, also based on Wannier Functions, the new approach is more general, being no longer restricted to the case of well separated interacting fragments. Moreover, it includes higher than pairwise energy contributions, coming from the dipole-dipole coupling among quantum oscillators. The method is successfully applied to the popular S22 molecular database, and also to extended systems, namely graphite and H2 adsorbed on the Cu(111) metal surface (in this case metal screening effects are taken into account). The results are also compared with those obtained by other vdW-corrected DFT schemes.

On the polarity of buckminsterfullerene with a water molecule inside.

Since the recent achievement of Kurotobi and Murata to capture a water molecule in a C(60) fullerene (Science 2011, 333, 613), there has been a debate about the properties of this H(2)O@C(60) complex. In particular, the polarity of the complex, which is thought to be underlying the easy separation of H(2)O@C(60) from the empty fullerene by HPLC, was calculated and found to be almost equal to that of an isolated water molecule. Here we present our detailed analysis of the charge distribution of the water-encapsulated C(60) complex, which shows that the polarity of the complex is, with 0.5 ± 0.1 D, indeed substantial, but significantly smaller than that of H(2)O. This may have important implications for the aim to design water-soluble and biocompatible fullerenes.

Optical van-der-Waals forces in molecules: from electronic Bethe-Salpeter calculations to the many-body dispersion model.

Molecular forces induced by optical excitations are connected to a wide range of phenomena, from chemical bond dissociation to intricate biological processes that underpin vision. Commonly, the description of optical excitations requires the solution of computationally demanding electronic Bethe-Salpeter equation (BSE). However, when studying non-covalent interactions in large-scale systems, more efficient methods are desirable. Here we introduce an effective approach based on coupled quantum Drude oscillators (cQDO) as represented by the many-body dispersion model. We find that the cQDO Hamiltonian yields semi-quantitative agreement with BSE calculations and that both attractive and repulsive optical van der Waals (vdW) forces can be induced by light. These optical-vdW interactions dominate over vdW dispersion in the long-distance regime, showing a complexity that grows with system size. Evidence of highly non-local forces in the human formaldehyde dehydrogenase 1MC5 protein suggests the ability to selectively activate collective molecular vibrations by photoabsorption, in agreement with recent experiments.

Exact Sum-Rule Approach to Polarizability and Asymptotic van der Waals Functionals─Derivation of Exact Single-Particle Benchmarks.

Using a sum-rule approach, we develop an exact theoretical framework for polarizability and asymptotic van der Waals correlation energy functionals of small isolated objects. The functionals require only monomer ground-state properties as input. Functional evaluation proceeds via solution of a single position-space differential equation, without the usual summations over excited states or frequency integrations. Explicit functional forms are reported for reference physical systems, including atomic hydrogen and single electrons subject to harmonic confinement, and immersed in a spherical-well potential. A direct comparison to the popular Vydrov-van Voorhis density functional shows that the best performance is obtained when density decay occurs at atomic scales. The adopted sum-rule approach implies general validity of our theory, enabling exact benchmarking of van der Waals density functionals and direct inspection of the subtle long-range correlation effects that constitute a major challenge for approximate (semi)local density functionals.

Experimental and Theoretical Studies toward Superior Anti-corrosive Nanocomposite Coatings of Aminosilane Wrapped Layer-by-Layer Graphene Oxide@MXene/Waterborne Epoxy.

Herein, layer-by-layer MXene/graphene oxide nanosheets wrapped with 3-aminopropyltriethoxy silane (abbreviated as F-GO@MXene) are proposed as an anti-corrosion promoter for waterborne epoxies. The GO@MXene nanohybrid is synthesized by a solvothermal reaction to produce a multi-layered 2D structure without defects. Then, the GO@MXene is modified by silane wrapping under a reflux reaction, in order to achieve chemical stability and to create active sites on the nanohybrid surface for reaction with the polymer matrix of the coating. The organic coating modified with 0.1 wt % F-GO@MXene has revealed superior corrosion protection efficiency than the organic coatings modified with either F-GO or F-MXene nanosheets. The impedance modulus at low frequency for the pure epoxy, epoxy/F-MXene, epoxy/F-GO, and epoxy/F-GO@MXene coatings is 4.17 × 105, 5.5 × 108, 4.46 × 108, and 1.14 × 1010 Ω·cm2 after 30 days of immersion in the corrosive media, respectively. The remarkable anti-corrosion property is assigned to the intense effect of the nanohybrid on the barrier performance, surface roughness, and adhesion strength of the epoxy coating. The complemental analysis based on first-principles density functional theory reveals that the adhesion strength related to the silane functional groups in its complexes follows the order F-GO@MXene > F-MXene > F-GO. The enhanced stabilization predicted on the GO@MXene nanohybrid ultimately stems from the combined role of the electrostatic and van der Waals forces, suggesting an increase in the penetration path of the corrosive media.

The initiation mechanisms for surface hydrosilylation with 1-alkenes.

Hydrosilylation provides a route to form substituted silanes in solution. A similar reaction has been observed in the formation of covalent organic monolayers on a hydrogen-terminated silicon surface and is called thermal hydrosilylation. In solution, the mechanism requires a catalyst to add the basal silicon and saturating hydrogen to the C=C double bond. On the silicon surface, however, the reaction proceeds efficiently at 200 °C, initiated by visible light, and more slowly at room temperature in the dark. Such low activation energy barriers for the reactions on a surface relative to that required for solution hydrosilylation are remarkable, and although many explanations have been suggested, controversy still exists. In this work using a constrained molecular dynamics approach within the density functional theory framework, we show that the free energy activation barrier for abstraction of a hydrogen from silicon by an alkene molecule can be overcome by visible light or thermal excitation. Furthermore, we show that by concerted transfer of a hydrogen from the α-carbon to the ß-carbon, a 1-alkene can insert its α-carbon into a surface Si-H bond to accomplish hydrosilylation.

Trends in the Change in Graphene Conductivity upon Gas Adsorption: The Relevance of Orbital Distortion.

The experimental ability to alter graphene (G) conductivity by adsorption of a single gas molecule is promoting the development of ultra-high-sensitivity gas detectors and could ultimately provide a novel playground for future nanoelectronics devices. At present, the underpinning effect is broadly attributed to a variation of G carrier concentration, caused by an adsorption-induced Fermi-level shift. By means of first-principle Kubo-Greenwood calculations, here we demonstrate that adsorbate-induced orbital distortion could also lead to small but finite G conductivity changes, even in the absence of Fermi-level shifts. This mechanism enables a sound physical interpretation of the observed variable sensitivity of G devices to different chemical moieties, and it can be strongly enhanced by using a suitable Ni substrate, thereby opening new pathways for the optimal design of operational nanoscale detectors.

Atomic Structure of Glassy GeTe4 as a Playground to Assess the Performances of Density Functional Schemes Accounting for Dispersion Forces.

The atomic structure of glassy GeTe4 is obtained in the framework of first-principles molecular dynamics (FPMD) by considering five different approaches for the description of the electronic structure within density functional theory (DFT). Among these schemes, one is not corrected by accounting for the dispersion forces and it is based on the BLYP exchange-correlation (XC) functional, while all of the others consider the dispersion forces according to different theoretical strategies. In particular, by maintaining the BLYP expression for the XC functional, two of them (BLYP-D2 and BLYP-D3) exploit the Grimme expressions for the dispersion forces, while the fourth scheme is based on the maximally localized Wannier functions (MLWFs). Finally, we also considered the rVV10 functional constructed to include seamlessly the dispersion part. Our results point out the better performances of BLYP-D3 and MLWF in terms of comparison with experimental data for the total pair correlation functions, with BLYP-D2 and rVV10 being closer to the uncorrected BLYP data. The implications of such findings are discussed by considering the overall limited impact of dispersion forces on the atomic structure of glassy GeTe4.

Are there immobilized water molecules around hydrophobic groups? Aqueous solvation of methanol from first principles.

Structural, dynamical, bonding, and electronic properties of water molecules around a soluted methanol molecule are studied from first principles. The results are compatible with experiments and qualitatively support the conclusions of recent classical molecular dynamics simulations concerning the controversial issue on the presence of "immobilized" water molecules around hydrophobic groups: the hydrophobic solute slightly reduces the mobility of many surrounding water molecules rather than immobilizing just the few ones which are closest to the methyl group. By generating maximally localized Wannier functions, a detailed description of the polarization effects in both solute and solvent molecules is obtained, which better elucidates the solvation process.