Differences in interaction of graphene/graphene oxide with bacterial and mammalian cell membranes.

Graphene, a single layer, hexagonally packed two-dimensional carbon sheet is an attractive candidate for diverse applications including antibacterial potential and drug delivery. One of the knowledge gaps in biomedical application of graphene is the interaction of these materials with the cells. To address this, we investigated the interaction between graphene materials (graphene and graphene oxide) and plasma membranes of cells (bacterial and mammalian cells). The interactions of four of the most abundant phospholipids in bacteria and mammalian plasma membranes with graphene materials were studied using density functional theory (DFT) at the atomic level. The calculations showed that the mammalian phospholipids have stronger bonding to each other compared to bacterial phospholipids. When the graphene/graphene oxide sheet is approaching the phospholipid pairs, the bacterial pairs exhibit less repulsive interactions, thereby a more stable system with the sheets was found. We also assembled bacterial and mammalian phospholipids into liposomes. We further observed that the bacterial liposomes and cells let the graphene flakes penetrate the membrane. The differential scanning calorimetry measurements of liposomes revealed that the bacterial liposomes have the lowest heat capacity; this strengthens the theoretical predictions of weaker interaction between the bacterial phospholipids compared to the mammalian phospholipids. We further demonstrated that graphene oxide could be internalized into the mammalian liposomes without disrupting the membrane integrity. The results suggest that the weak bonding among bacteria phospholipids and less repulsive force when graphene materials approach, result in graphene materials interacting differently with the bacteria compared to mammalian cells.

Filter function of graphene oxide: Trapping perfluorinated molecules.

We need clean drinking water, but current water purification methods are not always sufficient. This study examines the binding and binding mechanisms when graphene oxide is used as a filter material for removing perfluorinated substances and trihalomethanes. We use density functional theory calculations to examine the binding of the harmful molecules on graphene oxide. Our results indicate that the binding energies between graphene oxide and the investigated molecules are in the range of 370-1450 meV per molecule, similar to the binding energies obtained in other studies, where adsorption of similar size molecules onto graphene oxide has been investigated. This indicates that graphene oxide has the potential to separate the molecules of interest from the water. Significant contribution to the binding energies comes from the van der Waals (dispersion) interaction between the molecule and graphene oxide, while the hydrogen bonding between the functional groups of graphene oxide and the hydrogen atoms in functional groups on the molecules also plays a role in the binding.

Extent of Fock-exchange mixing for a hybrid van der Waals density functional?

The vdW-DF-cx0 exchange-correlation hybrid design [K. Berland et al., J. Chem. Phys. 146, 234106 (2017)] has a truly nonlocal correlation component and aims to facilitate concurrent descriptions of both covalent and non-covalent molecular interactions. The vdW-DF-cx0 design mixes a fixed ratio, a, of the Fock exchange into the consistent-exchange van der Waals density functional, vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)]. The mixing value a is sometimes taken as a semi-empirical parameter in hybrid formulations. Here, instead, we assert a plausible optimum average a value for the vdW-DF-cx0 design from a formal analysis; A new, independent determination of the mixing a is necessary since the Becke fit [A. D. Becke, J. Chem. Phys. 98, 5648 (1993)], yielding a' = 0.2, is restricted to semilocal correlation and does not reflect non-covalent interactions. To proceed, we adapt the so-called two-legged hybrid construction [K. Burke et al., Chem. Phys. Lett. 265, 115 (1997)] to a starting point in the vdW-DF-cx functional. For our approach, termed vdW-DF-tlh, we estimate the properties of the adiabatic-connection specification of the exact exchange-correlation functional, by combining calculations of the Fock exchange and of the coupling-constant variation in vdW-DF-cx. We find that such vdW-DF-tlh hybrid constructions yield accurate characterizations of molecular interactions (even if they lack self-consistency). The accuracy motivates trust in the vdW-DF-tlh determination of system-specific values of the Fock-exchange mixing. We find that an average value a' = 0.2 best characterizes the vdW-DF-tlh description of covalent and non-covalent interactions, although there exists some scatter. This finding suggests that the original Becke value, a' = 0.2, also represents an optimal average Fock-exchange mixing for the new, truly nonlocal-correlation hybrids. To enable self-consistent calculations, we furthermore define and test a zero-parameter hybrid functional vdW-DF-cx0p (having fixed mixing a' = 0.2) and document that this truly nonlocal correlation hybrid works for general molecular interactions (at reference and at relaxed geometries). It is encouraging that the vdW-DF-cx0p functional remains useful also for descriptions of some extended systems.

Involving High School Students in Computational Physics University Research: Theory Calculations of Toluene Adsorbed on Graphene.

To increase public awareness of theoretical materials physics, a small group of high school students is invited to participate actively in a current research projects at Chalmers University of Technology. The Chalmers research group explores methods for filtrating hazardous and otherwise unwanted molecules from drinking water, for example by adsorption in active carbon filters. In this project, the students use graphene as an idealized model for active carbon, and estimate the energy of adsorption of the methylbenzene toluene on graphene with the help of the atomic-scale calculational method density functional theory. In this process the students develop an insight into applied quantum physics, a topic usually not taught at this educational level, and gain some experience with a couple of state-of-the-art calculational tools in materials research.

Graphene oxide and adsorption of chloroform: A density functional study.

Chlorinated hydrocarbon compounds are of environmental concerns, since they are toxic to humans and other mammals, and are widespread, and exposure is hard to avoid. Understanding and improving methods to reduce the amount of the substances are important. We present an atomic-scale calculational study of the adsorption of chlorine-based substance chloroform (CHCl3) on graphene oxide, as a step in estimating the capacity of graphene oxide for filtering out such substances, e.g., from drinking water. The calculations are based on density functional theory, and the recently developed consistent-exchange functional for the van der Waals density-functional method is employed. We obtain values of the chloroform adsorption energy varying from roughly 0.2 to 0.4 eV per molecule. This is comparable to previously found results for chloroform adsorbed directly on clean graphene, using similar calculations. In a wet environment, like filters for drinking water, the graphene will not stay clean and will likely oxidize, and thus adsorption onto graphene oxide, rather than clean graphene, is a more relevant process to study.

van der Waals forces in density functional theory: a review of 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 emphasizing 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.

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. 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 (2014)] 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 PbTiO3, the adsorption of small molecules within metal-organic frameworks, 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.

A van der Waals density functional study of chloroform and other trihalomethanes on graphene.

A computational study of chloroform (CHCl(3)) and other trihalomethanes (THMs) adsorbed on graphene is presented. The study uses the van der Waals density functional method to obtain adsorption energies and adsorption structures for these molecules of environmental concern. In this study, chloroform is found to adsorb with the H atom pointing away from graphene, with adsorption energy 357 meV (34.4 kJ/mol). For the other THMs studied the calculated adsorption energy values vary from 206 meV (19.9 kJ/mol) for fluoroform (CHF(3)) to 404 meV (39.0 kJ/mol) for bromoform (CHBr(3)). The corrugation of graphene as seen by the THMs is small, the difference in adsorption energy along the graphene plane is less than 6 meV for chloroform.

Physisorption of nucleobases on graphene: a comparative van der Waals study.

The physisorption of the nucleobases adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) on graphene is studied using several variants of the density functional theory (DFT): the generalized gradient approximation with the inclusion of van der Waals interaction (vdW) based on the TS approach (Tkatchenko and Scheffer 2009 Phys. Rev. Lett. 102 073005) and our simplified version of this approach (here called sTS), the van der Waals density functional vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401) and vdW-DF2 (Lee et al 2010 Phys. Rev. B 82 081101), and DFT-D2 (Grimme 2006 J. Comput. Chem. 27 1787) and DFT-D3 (Grimme et al 2010 J. Chem. Phys. 132 154104) methods. The binding energies of nucleobases on graphene are found to be in the following order: G > A > T > C > U within TS, sTS, vdW-DF, and DFT-D2, and in the following order: G > A > T ~ C > U within DFT-D3 and vdW-DF2. The binding separations are found to be different within different methods and in the following order: DFT-D2 < TS < DFT-D3 ~ vdW-DF2 < vdW-DF. We also comment on the efficiency of combining the DFT-D approach and vdW-DF to study systems with van der Waals interactions.

Desorption of n-alkanes from graphene: a van der Waals density functional study.

A recent study of temperature-programmed desorption (TPD) measurements of small linear alkane molecules (n-alkanes, with formula C(N)H(2N+2)) from C(0001) deposited on Pt(111) shows a linear relationship of the desorption energy with increasing n-alkane chain length N. We here present a van der Waals density functional study of the desorption barrier energy of the ten smallest n-alkanes (of carbon chain length N = 1-10) from graphene. We find linear scaling with N, including a non-zero intercept with the energy axis, i.e. an offset at the extrapolation to N = 0. This calculated offset is quantitatively similar to the results of the TPD measurements. From further calculations of the polyethylene polymer we offer a suggestion for the origin of the offset.

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 sensitive to the choice of exchange functional approximation. Further, the DFT-D3 and TS-vdW corrections to traditional DFT formulations are also benchmarked, and deviations are analyzed.

A van der Waals density functional study of adenine on graphene: single-molecular adsorption and overlayer binding.

The adsorption of an adenine molecule on graphene is studied using a first-principles van der Waals functional, vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401). The cohesive energy of an ordered adenine overlayer is also estimated. For the adsorption of a single molecule, we determine the optimal binding configuration and adsorption energy by translating and rotating the molecule. The adsorption energy for a single molecule of adenine is found to be 711 meV, which is close to the calculated adsorption energy of the similarly sized naphthalene. On the basis of the single-molecular binding configuration, we estimate the cohesive energy of a two-dimensional ordered overlayer. We find a significantly stronger binding energy for the ordered overlayer than for single-molecule adsorption.

Stacking interactions and the twist of DNA.

The importance of stacking interactions for the Twist and stability of DNA is investigated using the fully ab initio van der Waals density functional (vdW-DF). Our results highlight the role that binary interactions between adjacent sets of base pairs play in defining the sequence-dependent Twists observed in high-resolution experiments. Furthermore, they demonstrate that additional stability gained by the presence of thymine is due to methyl interactions with neighboring bases, thus adding to our understanding of the mechanisms that contribute to the relative stability of DNA and RNA. Our mapping of the energy required to twist each of the 10 unique base pair steps should provide valuable information for future studies of nucleic acid stability and dynamics. The method introduced will enable the nonempirical theoretical study of significantly larger pieces of DNA or DNA/amino acid complexes than previously possible.

Application of van der Waals density functional to an extended system: adsorption of benzene and naphthalene on graphite.

It is shown that it is now possible to include van der Waals (vdW) interactions via a nonempirical implementation of density functional (DF) theory to describe the correlation energy in electronic structure calculations on infinite systems of no particular symmetry. The vdW-DF theory [Phys. Rev. Lett. 92, 246401 (2004)] is applied to the adsorption of benzene and naphthalene on an infinite sheet of graphite, as well as the binding between two graphite sheets. A comparison with recent thermal-desorption data [Phys. Rev. B 69, 155406 (2004)] shows great promise for the vdW-DF method.

van der Waals interaction of simple, parallel polymers.

We study the mutual interactions of simple parallel polymers within the framework of density-functional theory (DFT). As the conventional implementations of DFT do not treat the long-range dispersion [van der Waals (vdW)] interactions, we develop a systematic correction scheme for the nonlocal energy contribution of the polymer interaction at the intermediate to the asymptotic separations. We primarily focus on the three polymers, polyethylene, isotactic polypropylene, and isotactic polyvinylchloride, but the scheme presented applies also more generally to other simple polymers. From first-principle calculations we extract the geometrical and electronic structures of the polymers and the local part of their interaction energy, as well as the static electric response. The dynamic electrodynamic response is modeled on the basis of these static calculations, from which the nonlocal vdW interaction of the polymers is extracted.

van der Waals interactions of polycyclic aromatic hydrocarbon dimers.

Density functional theory is in principle exact and includes also long-range interactions, such as the van der Waals interactions. These are, however, part of the exchange-correlation energy functional that needs to be approximated, and are absent in the local and semilocal standard implementations. Recently a density functional which includes van der Waals interactions for planar systems has been developed, which we show can be extended to provide a treatment of planar molecules. We use this functional to calculate binding distances and energies for dimers of three of the smallest polycyclic aromatic hydrocarbons (PAHs)--naphthalene, anthracene, and pyrene.