Chiral organic molecular structures supported by planar surfaces

We employ the molecular dynamics simulations to study the dynamics of acetanilide (ACN) molecules placed on a flat surface of planar multilayer hexagonal boron nitride. We demonstrate that the ACN molecules, known to be achiral in the three-dimensional space, become chiral after being placed on the substrate. Homochirality of the ACN molecules leads to stable secondary structures stabilized by hydrogen bonds between peptide groups of the molecules. By employing molecular dynamics simulations, we reveal that the structure of the resulting hydrogen-bond chains depends on the isomeric composition of the molecules. If all molecules are homochiral (i.e., with only one isomer being present), they form secondary structures (chains of hydrogen bonds in the shapes of arcs, circles, and spirals). If the molecules at the substrate form a racemic mixture, then no regular secondary structures appear, and only curvilinear chains of hydrogen bonds of random shapes emerge. A hydrogen-bond chain can form a zigzag array only if it has an alternation of isomers. Such chains can create two-dimensional (2D) regular lattices or 2D crystals. The melting scenarios of such 2D crystals depend on density of its coverage of the substrate. At 25% coverage, melting occurs continuously in the temperature interval 295-365 K. For a complete coverage, melting occurs at 415-470 K due to a shift of 11% of all molecules into the second layer of the substrate.

Review of research of nanocomposites based on graphene quantum dots

Graphene quantum dots (GQDs) belong to the vast and versatile family of carbon nanomaterials. Their unique position amongst versatile carbon nanoparticles (NPs) originates from the properties of quantum confinement and edge effects. GQDs are similar to conventional semiconductor QDs due to their tunable band gaps and high photoluminescence activity. However, GQDs have superior characteristics due to their excellent biocompatibility, low toxicity, good water dispersibility, large optical absorptivity, high fluorescence activity and photostability. These properties have generated significant interest in GQDs applications in various fields: nanosensor fabrication, drug delivery, photocatalysis, photovoltaics, and photodynamic therapy. Numerous GQD-based nanocomposites/nanohybrides have been synthesized and/or studied computationally. This review focuses on recent computational studies of various GQD-based nanocomposites/nanohybrides and systems which can be related to them.

Combined DFT and Molecular Mechanics Modeling of the Adsorption of Semiconducting Molecules on an Ionic Substrate: PTCDA and CuPc on NaCl

Experimental results suggest that molecular geometry and energies can be influenced by the presence of thin film substrates as well as surrounding molecules. It is imperative that computational models take this influence into account. The accurate computational modeling of these molecules is an efficient way of carrying out chemistry calculations and reinforcing experimental findings. In our study, density functional theory (DFT) and molecular mechanics (MM) are used to model the configurations of the organic semiconducting materials, 3,4,9,10-perylene tetracarboxylic dianhydride, C₂₄H₈O₆ (PTCDA), and copper(II) phthalocyanine, C₃₂H₁₆CuN₈ (CuPc), as adsorbed on single- and double-layer NaCl substrates of various dimensions and charge settings. After a geometry and charge optimization of the molecules using DFT, the molecular geometries are optimized under different environments using computational calculations with specific force-field settings in HyperChem Professional 8.0(TM) software using MM. Energies and geometries of the molecules are then recorded, and our data are compared to experimental results of similar systems. We find that, with the appropriate choice of substrate properties, the calculated molecular configurations directly reflect those found experimentally. Our results support the idea that this method of simulation can produce reliable models in the field of physical chemistry.

Development of a schwarzite-based moving bed 3D printed water treatment system for nanoplastic remediation

The impact of micro and nanoplastic debris on our aquatic ecosystem is among the most prominent environmental challenges we face today. In addition, nanoplastics create significant concern for environmentalists because of their toxicity and difficulty in separation and removal. Here we report the development of a 3D printed moving bed water filter (M-3DPWF), which can perform as an efficient nanoplastic scavenger. The enhanced separation of the nanoplastics happens due to the creation of a charged filter material that traps the more surface charged nanoparticles selectively. Synthetic contaminated water from polycarbonate waste has been tested with the filter, and enhanced nanoplastic removal has been achieved. The proposed filtration mechanism of surface-charge based water cleaning is further validated using density function theory (semi-empirical) based simulation. The filter has also shown good structural and mechanical stability in both static and dynamic water conditions. The field suitability of the novel treatment system has also been confirmed using water from various sources, such as sea, river, and pond. Our results suggest that the newly developed water filter can be used for the removal of floating nanoparticles in water as a robust advanced treatment system.

The impact of micro and nanoplastic debris on our aquatic ecosystem is among the most prominent environmental challenges we face today.

Adsorption-induced clustering of CO2 on graphene

Utilization of graphene-based materials for selective carbon dioxide capture has been demonstrated recently as a promising technological approach. In this study we report results from density functional theory calculations and molecular dynamics simulations on the adsorption of CO2, N2, and CH4 gases on a graphene sheet. We calculate adsorption isotherms of ternary and binary mixtures of these gases and reproduce the larger selectivity of CO2 to graphene relative to the other two gases. Furthermore it is shown that the confinement to two-dimensions, associated with adsorbing the CO2 gas molecules on the plane of graphene, increases their propensity to form clusters on the surface. Above a critical surface coverage (or partial pressure) of the gas, these CO2-CO2 interactions augment the effective adsorption energy to graphene, and, in part, contribute to the high selectivity of carbon dioxide with respect to nitrogen and methane. The origin of the attractive interaction between the CO2 molecules adsorbed on the surface is of electric quadrupole-quadrupole nature, in which the positively-charged carbon of one molecule interacts with the negatively-charged oxygen of another molecule. The energy of attraction of forming a CO2 dimer is predicted to be around 5-6 kJ mol-1, much higher than the corresponding values calculated for N2 and CH4. We also evaluated the adsorption energies of these gases to a graphene sheet and found that the attractions obtained using the classical force-fields might be over-exaggerated. Nevertheless, even when the magnitudes of these (classical force-field) graphene-gas interactions are scaled-down sufficiently, the tendency of CO2 molecules to cluster on the surface is still observed.

Exact matrix treatment of the statistical mechanics of adsorption of large aromatic molecules on graphene

Experimental studies of adsorption from solution of the large aromatic molecules 1,2-dihydroxybenzene (catechol) and phenyl hydroquinone on graphene nanoplatelets show that at low coverage adsorption is followed by a transition which occurs from adsorbed molecules in flat to more vertically oriented states. Catechol adsorption isotherms exhibit 2 plateaus while phenyl hydroquinone shows 3 plateaus indicating 2 and 3 active conformers respectively participating in the adsorption process. Modelling such adsorption isotherms presents a challenge. Here, an exact matrix treatment of the statistical mechanics of a one-dimensional model of adsorption of catechol and dihydroquinone on graphene nanoplatelets is presented. The theoretical adsorption isotherms successfully reproduce all the features of both the catechol and dihydroquinone experimental adsorption isotherms. As suggested by the experimentalists, our theoretical model demonstrates that adsorbed phenyl hydroquinone molecules adopt a flat orientation at low concentrations and an edge orientation at higher coverage before eventually adopting a vertical configuration. Both catechol and phenyl hydroquinone can be described by our interconvertible monomer-dimer-trimer model. The theoretical adsorption isotherms obtained show several plateaus reflecting the types of conformer on the graphene surface.

Density Functional Study on the Adsorption of 5-Membered N-Heterocycles on B/N/BN-Doped Graphene: Coronene as a Model System

Adsorption of seven 5-membered N-heterocycles on B/N/BN-doped graphene (with coronene as a model system) has been studied using density functional theory (DFT). The geometry of the complexes validated the involvement of both π···π stacking and N-H···π interaction in the adsorption process. The stability of the complexes is measured in terms of stabilization energy, and the results suggested that the complexes are stable enough (stabilization energies are in the range of 7.61-14.77 kcal mol^-1). Studies confirmed the stability of complexes in the solvent phase too irrespective of the dielectric of the solvent. Dispersive force is the major mode of interaction in stabilizing the complexes. Natural bond orbital analysis indicated a small contribution from electrostatic and covalent interactions. Thermochemical analysis revealed that the complexation is exothermic in nature and favorable at a lower temperature. Adsorption of N-heterocycles exerts a nominal impact on the electronic properties of the undoped/doped graphene. The study presents a simple approach to introduce an arbitrary functionality to undoped/doped graphene by preserving its electronic properties.

Modeling the adsorption equilibrium of small-molecule gases on graphene: effect of the volume to surface ratio

The adsorption probability of molecules at surfaces depends on the available volume-to-surface ratio with important implications for 2D self-assembly.

Interplay Between π-Stacking and Hydrogen Bonding in the Self-Association of Different Isomers of Naphthalenedicarboxylic Acid

Using proton and carbon chemical shifts, we investigated the self-association of three isomers of naphthalenedicarboxylic acid, a model for the aggregation of asphaltenes. Experimental proton chemical shifts of each isomer were measured as a function of concentration in an aprotic solvent. Several potential structures of the monomer and dimer of each naphthalenedicarboxylic acid were considered, and calculated proton chemical shifts for the potential monomer and dimer structures were compared to the experimental chemical shifts to find the weighted average structure that best fit the experimental shifts. Calculated carbon chemical shifts were also compared to experimental values. The chemical shift comparison and calculated energies indicate that π-stacked dimers are not likely to contribute significantly to the dimer structure of any of the three naphthalenedicarboxylic acid isomers studied.

Approximate quantum chemical methods for modelling carbohydrate conformation and aromatic interactions: β-cyclodextrin and its adsorption on a single-layer graphene sheet

Adsorption of carbohydrates on graphene has the potential to improve graphene dispersibility in water. Here we assess the ability of DFTB-based and NDDO-based quantum chemical methods to model β-cyclodextrin conformations and interactions with graphene.

Shielding the chemical reactivity using graphene layers for controlling the surface properties of carbon materials

Graphene can efficiently shield chemical interactions and gradually decrease the binding to reactive defect areas. In the present study, we have used the observed graphene shielding effect to control the reactivity patterns on the carbon surface. The experimental findings show that a surface coating with a tiny carbon layer of 1-2 nm thickness is sufficient to shield the defect-mediated reactivity and create a surface with uniform binding ability. The shielding effect was directly observed using a combination of microscopy techniques and evaluated with computational modeling. The theoretical calculations indicate that a few graphene layers can drastically reduce the binding energy of the metal centers to the surface defects by 40-50 kcal mol(-1). The construction of large carbon areas with controlled surface reactivity is extremely difficult, which is a key limitation in many practical applications. Indeed, the developed approach provides a flexible and simple tool to change the reactivity patterns on large surface areas within a few minutes.

Recent Progress in Treating Protein-Ligand Interactions with Quantum-Mechanical Methods

We review the first successes and failures of a “new wave” of quantum chemistry-based approaches to the treatment of protein/ligand interactions. These approaches share the use of “enhanced”, dispersion (D), and/or hydrogen-bond (H) corrected density functional theory (DFT) or semi-empirical quantum mechanical (SQM) methods, in combination with ensemble weighting techniques of some form to capture entropic effects. Benchmark and model system calculations in comparison to high-level theoretical as well as experimental references have shown that both DFT-D (dispersion-corrected density functional theory) and SQM-DH (dispersion and hydrogen bond-corrected semi-empirical quantum mechanical) perform much more accurately than older DFT and SQM approaches and also standard docking methods. In addition, DFT-D might soon become and SQM-DH already is fast enough to compute a large number of binding modes of comparably large protein/ligand complexes, thus allowing for a more accurate assessment of entropic effects.

Noncovalent Interaction of Graphene with Heterocyclic Compounds: Benzene, Imidazole, Tetracene, and Imidazophenazines

Noncovalent functionalization of graphene with organic molecules offers a direct route to multifunctional modification of this nanomaterial, leading to its various possible practical applications. In this work, the structures of hybrids formed by linear heterocyclic compounds such as imidazophenazine (F1) and its derivatives (F2-F4) with graphene and the corresponding interaction energies are studied by using the DFT method. Special attention is paid to the hybrids where the attached molecule is located along the graphene zigzag (GZZ ) and armchair (GAC ) directions. The interaction energies corresponding to the graphene hybrids of the F1-F4 compounds for the two directions are found to be distinct, while tetracene (being a symmetrical molecule) shows a small difference between these binding energies. It is found that the back-side CH3 and CF3 groups have an important influence on the arrangements of F1 derivatives on graphene and on their binding energies. The contribution of the CF3 group to the total binding energy of the F3 molecule with graphene is the largest (3.4 kcal mol(-1) ) (the GZZ direction) while the CH3 group increases this energy of F2 only by 2.0 kcal mol(-1) (the GAC direction). It is shown that replacing the carbons with other atoms or adding a back-side group enables one to vary the polarizability of graphene.

Theoretical study of adsorption of amino acids on graphene and BN sheet in gas and aqueous phase with empirical DFT dispersion correction

Suitability of BN nanosheet/graphene towards the adsorption of amino acids established by dispersion corrected DFT calculations.

Modelling of graphene functionalization

This perspective describes the available theoretical methods and models for simulating graphene functionalization based on quantum and classical mechanics.

Perchlorination of Coronene Enhances its Propensity for Self-Assembly on Graphene

Providing a quantitative understanding of the thermodynamics involved in molecular adsorption and self-assembly at a nanostructured carbon material is of fundamental importance and finds outstanding applications in the graphene era. Here, we study the effect of edge perchlorination of coronene, which is a prototypical polyaromatic hydrocarbon, on the binding affinity for the basal planes of graphite. First, by comparing the desorption barrier of hydrogenated versus perchlorinated coronene measured by temperature-programmed desorption, we quantify the enhancement of the strength of physisorption at the single-molecule level though chlorine substitution. Then, by a thermodynamic analysis of the corresponding monolayers based on force-field calculations and statistical mechanics, we show that perchlorination decreases the free energy of self-assembly, not only enthalpically (by enhancing the strength of surface binding), but also entropically (by decreasing the surface concentration). The functional advantage of a chemically modulated 2D self-assembly is demonstrated in the context of the molecule-assisted liquid-phase exfoliation of graphite into graphene.

Benchmark calculations of the adsorption of aromatic molecules on graphene

Selecting the saturated graphene fragment as a model of graphene, we have investigated seven popular density functionals, including ωB97X-D, B97-D, B-LYP-D3, M05-2X, M06-2X, M11-L, and N12, for their performance in describing the adsorption of aromatic molecules on graphene. The best performing functionals are B97-D, B-LYP-D3, and ωB97X-D. M05-2X, M06-2X, and M11-L significantly underestimate the adsorption strengths, while N12 fails completely in this respect. The effects of the basis sets and size of the saturated graphene fragments on the geometries, energies, and properties for the adsorption of aromatic molecules on graphene have also been studied. It was found that the small basis sets such as 6-31G(d) and jun-cc-pVDZ are not suitable for the accurate description of the adsorption of aromatic molecules on graphene. The size of selected graphene fragments has a little effect on both the ωB97X-D and SCS-SAPT0 interaction energies, but the effects of the size of selected graphene fragments on the energy components are significant in some cases of the adsorption of aromatic molecules on graphene. The surprising weakness of electrostatic interactions by F substitution for the adsorption of F-substituted benzenes on graphene was explained using the energy component analysis.

Enhanced semiempirical QM methods for biomolecular interactions

Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.

Noncovalent π⋅⋅⋅π interaction between graphene and aromatic molecule: structure, energy, and nature

Noncovalent π⋅⋅⋅π interactions between graphene and aromatic molecules have been studied by using density functional theory with empirical dispersion correction (ωB97X-D) combined with zeroth-order symmetry adapted perturbation theory (SAPT0). Excellent agreement of the interaction energies computed by means of ωB97X-D and spin component scaled (SCS) SAPT0 methods, respectively, shows great promise for the two methods in the study of the adsorption of aromatic molecules on graphene. The other important finding in this study is that, according to SCS-SAPT0 analyses, π⋅⋅⋅π interactions between graphene and aromatic molecules are largely dependent on both dispersion and electrostatic type interactions. It is also noticed that π⋅⋅⋅π interactions become stronger and more dispersive (less electrostatic) upon substitution of the very electronegative fluorine atoms onto the aromatic molecules.