Ab-Initio calculations including Van der Waals interactions: the SnS2 layered material

Tin disulfide SnS2 was recently proposed as a high efficiency solar cell precursor [1]. The aim of this work is a deep study of the structural disposition of the most important polytipes of this layered material, not only describing the electronic correlation but also the interatomic Van der Waals interactions that is present between the layers. The two recent implementations to take Van der Waals interactions into account in the VASP code are the self-consistent Dion et al. [2] functional optimized for solids by Michaelides et al [3] and the Grimme [4] dispersion correction that is applied after each autoconsistent PBE electronic calculation. In this work these two methods are compared with DFT PBE functional. The results we will presented at this Conference, demonstrates the enhancement of the geometric parameters by the use of the Van der Waals interactions in agreement with the experimental values.

Van der Waals-corrected density functional theory: benchmarking for hydrogen-nanotube and nanotube-nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H-2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H-2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

Van der Waals interactions between two parallel infinitely long single-walled nanotubes

The potential energies of van der Waals (VDW) interactions between two parallel, infinitely long and perfect SWNTs with identical, and different sizes were studied based on the continuum Lennard-Jones model. The conclusion of Girifalco's work on (n, n) SWNTs that the potentials of SWNT-SWNT fell on a single curve, is also applicable to SWNTs with different sizes. We further obtained the corresponding constants of the well depth phi(0) and equilibrium VDW gap g(0) for SWNTs with a radius from 2 to 25 Angstrom. (C) 2005 Elsevier B.V. All rights reserved.

London-van der Waals adhesiveness of rough particles

Van der Waals forces often dominate interactions and adhesion between fine particles and, in turn, decisively influence the bulk behaviour of powders. However, so far there is no effective means to characterize the adhesive behaviour of such particles. A complication is that most powder particles have rough surfaces, and it is the asperities on the surfaces that touch, confounding the actual surface that is in contact. Conventional approaches using surface energy provide limited information regarding adhesion, and pull-off forces measured through atomic force microscope (AFM) are highly variable and difficult to interpret. In this paper we develop a model which combines the Rumpf-Rabinovich and the JKR-DMT theories to account simultaneously for the effects of surface roughness and deformation on adhesion. This is applied to a 'characteristic asperity' which may be easily obtained from AFM measurements. The concept of adhesiveness, a material property reflecting the influences of elastic deformability, surface roughness, and interfacial surface energy, is introduced as an efficient and quantitative measure of the adhering tendency of a powder. Furthermore, a novel concept of specific adhesiveness is proposed as a convenient tool for characterizing and benchmarking solid materials. This paper provides an example to illustrate the use of the proposed theories. (c) 2005 Elsevier B.V. All rights reserved.

van der Waals coefficients for positronium interactions with atoms

The random-phase approximation with exchange (RPAE) is used with a B-spline basis to compute dynamic dipole polarizabilities of noble-gas atoms and several other closed-shell atoms (Be, Mg, Ca, Zn, Sr, Cd, and Ba). From these, values of the van der Waals C₆ constants for positronium interactions with these atoms are determined and compared with existing data. After correcting the RPAE polarizabilities to fit the most accurate static polarizability data, our best predictions of C₆ for Ps–noble-gas pairs are expected to be accurate to within 1%, and to within a few percent for the alkaline-earth metals. We also used accurate dynamic dipole polarizabilities from the literature to compute the C₆ coefficients for the alkali-metal atoms. Implications of increased C₆ values for Ps scattering from more polarizable atoms are discussed.

Epitaxial growth of molecular crystals on van der Waals substrates for high-performance organic electronics

Epitaxial van der Waals (vdW) heterostructures of organic and layered materials are demonstrated to create high-performance organic electronic devices. High-quality rubrene films with large single-crystalline domains are grown on h-BN dielectric layers via vdW epitaxy. In addition, high carrier mobility comparable to free-standing single-crystal counterparts is achieved by forming interfacial electrical contacts with graphene electrodes.

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).

Bestimmung der Hamaker-Konstanten von polymeren Nanopartikeln in organischen Lösungsmitteln und wässrigen Elektrolytlösungen mittels asymmetrischer Fluss-Feld-Fluss-Fraktionierung

Die zwischen allen Objekten vorhandenen Wechselwirkungen können repulsiver und attraktiver Natur sein. Bei den attraktiven Kräften kommt der Bestimmung von Dispersionskräften eine besondere Bedeutung zu, da sie in allen kolloidalen Systemen vorhanden sind und entscheidenden Einfluss auf die Eigenschaften und Prozesse dieser Systeme nehmen. Eine der Möglichkeiten, Theorie und Experiment zu verbinden, ist die Beschreibung der London-Van der Waals-Wechselwirkung durch die Hamaker-Konstante, welche durch Berechnungen der Wechselwirkungsenergie zwischen Objekten erhalten werden kann. Für die Beschreibung von Oberflächenphänomenen wie Adhäsion, die in Termen der totalen potentiellen Energie zwischen Partikeln und Substrat beschrieben werden, benötigt man exakt bestimmte Hamaker-Konstanten. In der vorliegenden Arbeit wurde die asymmetrische Fluss Feld-Fluss Fraktionierung in Kombination mit einem auf dem Newton-Algorithmus basierenden Iterationsverfahren zur Bestimmung der effektiven Hamaker-Konstanten verschiedener Nanopartikeln sowie Polystyrollatex-Partikel in Toluol bzw. Wasser verwendet. Der Einfluss verschiedener Systemparameter und Partikeleigenschaften wurde im Rahmen der klassischen DLVO-Theorie untersucht.

Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Estudo teórico QTAIM e DFT dos compostos de coordenação: efeito quelato, titanocenos e ligação química

This work is a study of coordination compounds by quantum theory of atoms in molecules (QTAIM), based on the topological analysis of the electron density of molecular systems, both theoretically and experimentally obtained. The coordination chemistry topics which were studied are the chelate effect, bent titanocene and chemical bond in coordination complexes. The chelate effect was investigated according to topological and thermodynamic parameters. The exchange of monodentate ligands on polydentate ligands from same transition metal increases the stability of the complex both from entropy and enthalpy contributions. In some cases, the latter had a higher contribution to the stability of the complex in comparison with entropy. This enthalpic contribution is explained according to topological analysis of the M-ligand bonds where polidentate complex had higher values of electron density of bond critical point, Laplacian of electron density of bond critical point and delocalization index (number of shared electrons between two atoms). In the second chapter, was studied bent titanocenes with bulky cyclopentadienyl derivative π-ligand. The topological study showed the presence of secondary interactions between the atoms of π-ligands or between atoms of π-ligand and -ligand. It was found that, in the case of titanocenes with small difference in point group symmetry and with bulky ligands, there was an nearly linear relationship between stability and delocalization index involving the ring carbon atoms (Cp) and the titanium. However, the titanocene stability is not only related to the interaction between Ti and C atoms of Cp ring, but secondary interactions also play important role on the stability of voluminous titanocenes. The third chapter deals with the chemical bond in coordination compounds by means of QTAIM. The quantum theory of atoms in molecules so far classifies bonds and chemical interactions in two categories: closed shell interaction (ionic bond, hydrogen bond, van der Waals interaction, etc) and shared interaction (covalent bond). Based on topological parameters such as electron density, Laplacian of electron density, delocalization index, among others, was classified the chemical bond in coordination compounds as an intermediate between closed shell and shared interactions

Structural properties and energetics of diffuse Rb-87 clusters in three-dimension

A correlated two-body basis function is used to describe the three-dimensional bosonic clusters interacting via two-body van der Waals potential. We calculate the ground state and the zero orbital angular momentum excited states for Rb-N clusters with up to N = 40. We solve the many-particle Schrodinger equation by potential harmonics expansion method, which keeps all possible two-body correlations in the calculation and determines the lowest effective many-body potential. We study energetics and structural properties for such diffuse clusters both at dimer and tuned scattering length. The motivation of the present study is to investigate the possibility of formation of N-body clusters interacting through the van der Waals interaction. We also compare the system with the well studied He, Ne, and Ar clusters. We also calculate correlation properties and observe the generalised Tjon line for large cluster. We test the validity of the shape-independent potential in the calculation of the ground state energy of such diffuse cluster. These are the first such calculations reported for Rb clusters. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730972]

Modeling and simulation of fibrinogen and its adsorption behavior

In this thesis different approaches for the modeling and simulation of the blood protein fibrinogen are presented. The approaches are meant to systematically connect the multiple time and length scales involved in the dynamics of fibrinogen in solution and at inorganic surfaces. The first part of the thesis will cover simulations of fibrinogen on an all atom level. Simulations of the fibrinogen protomer and dimer are performed in explicit solvent to characterize the dynamics of fibrinogen in solution. These simulations reveal an unexpectedly large and fast bending motion that is facilitated by molecular hinges located in the coiled-coil region of fibrinogen. This behavior is characterized by a bending and a dihedral angle and the distribution of these angles is measured. As a consequence of the atomistic detail of the simulations it is possible to illuminate small scale behavior in the binding pockets of fibrinogen that hints at a previously unknown allosteric effect. In a second step atomistic simulations of the fibrinogen protomer are performed at graphite and mica surfaces to investigate initial adsorption stages. These simulations highlight the different adsorption mechanisms at the hydrophobic graphite surface and the charged, hydrophilic mica surface. It is found that the initial adsorption happens in a preferred orientation on mica. Many effects of practical interest involve aggregates of many fibrinogen molecules. To investigate such systems, time and length scales need to be simulated that are not attainable in atomistic simulations. It is therefore necessary to develop lower resolution models of fibrinogen. This is done in the second part of the thesis. First a systematically coarse grained model is derived and parametrized based on the atomistic simulations of the first part. In this model the fibrinogen molecule is represented by 45 beads instead of nearly 31,000 atoms. The intra-molecular interactions of the beads are modeled as a heterogeneous elastic network while inter-molecular interactions are assumed to be a combination of electrostatic and van der Waals interaction. A method is presented that determines the charges assigned to beads by matching the electrostatic potential in the atomistic simulation. Lastly a phenomenological model is developed that represents fibrinogen by five beads connected by rigid rods with two hinges. This model only captures the large scale dynamics in the atomistic simulations but can shed light on experimental observations of fibrinogen conformations at inorganic surfaces.

Electrostatically mediated specific adsorption of small molecules in metallo-organic frameworks

We investigate the interaction of ethylene and ethane with a Cu-tricarboxylate complex and show that at low loadings the lighter molecule has a higher binding energy as a result of an increased interaction with the framework Cu and stronger hydrogen bonding with the basic framework oxygens. This leads to selective adsorption of ethylene by a factor of about 2 at low pressure, which is overcome by the stronger van der Waals interaction of ethane at high loadings, explaining recent literature data. The results suggest the possibility of separation of light hydrocarbons at low pressures or in trace amounts.