multiscale-sampling-heterogeneous-watermetal-catalyst-interface-using

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics

A significant number of heterogeneously-catalyzed chemical processes occur under liquid conditions, but simulating catalyst function under such conditions ...

Clemson University

Many proposed applications using dendrimers, such as drug delivery and environmental remediation, involve dendrimer interactions with small molecules. Understanding the details of these interactions is important for designing dendrimers with tunable association with guest molecules. In this work, we investigate dendrimer interactions with small aromatic hydrocarbons using all-atom molecular dynamics simulations. We study the association of naphthalene (NPH)-the smallest polycyclic aromatic hydrocarbon-with 3rd-6th generation (G3-G6) polyamidoamine (PAMAM) dendrimers. Our work emphasizes that the association of small aromatic molecules with PAMAM dendrimers involves the formation of dynamic pocket-like association sites through interactions between flexible dendrimer branches and NPH molecules. The association sites are primarily formed by branches from the two outermost dendrimer subgenerations, and often involve the tertiary amine groups. Irrespective of their location on the dendrimer-whether buried or near the outer surface-these pocket-like structures lower the hydration of the associated NPH molecules. We show that on average NPH molecules with a lower hydration have a greater tendency to remain associated with the dendrimer for longer times. In general, the association sites are similar for the G3-G6 PAMAM dendrimers, indicating similarities in the association mechanisms across different dendrimer generations. 

Association of small aromatic molecules with PAMAM dendrimers.

We compare the free energies of adsorption (ΔAads) and the structural preferences of amino acids on graphene obtained using the non-polarizable force fields-Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid-graphene interactions are favorable irrespective of the force field. While the magnitudes of ΔAads differ between the force fields, the relative free energy of adsorption across amino acids is similar for the studied force fields. ΔAads positively correlates with amino acid-graphene and negatively correlates with graphene-water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that the force fields studied will give qualitatively consistent relative strength of adsorption across proteins but different structural preferences in the graphene adsorbed state.

Adsorption of amino acids on graphene: assessment of current force fields.

Solvent plays an important role in liquid phase heterogeneous catalysis; however, methods for calculating the free energies of catalytic phenomena at the solid-liquid interface are not well-established. For example, solvent molecules alter the energies of catalytic species and participate in catalytic reactions and can thus significantly influence catalytic performance. In this work, we begin to establish methods for calculating the free energies of such phenomena, specifically, by employing an explicit solvation method using a multiscale sampling (MSS) approach. This MSS approach combines classical molecular dynamics with density functional theory. We use it to calculate the free energies of solvation of catalytic species, specifically adsorbed NH*, NH*, CO*, COH*, CHOH*, and CHO* on Pt(111) surfaces under aqueous phase and under a mixed HO/CHOH solvent. We compare our calculated values with analogous values from implicit solvation for validation and to identify situations where implicit solvation is sufficient versus where explicit solvent is needed to compute adsorbate free energies. Our results indicate that explicit quantum-based methods are needed when adsorbates form chemical bonds and/or strong hydrogen bonds with HO solvent. Using MSS, we further separate the calculated free energies into energetic and entropic contributions in order to understand how each influences the free energy. We find that adsorbates that exhibit strong energies also exhibit strong and negative entropies, and we attribute this relationship to hydrogen bonding between the adsorbates and the solvent molecules, which provides a large energetic contribution but reduces the overall mobility of the solvent.

Sapna Sarupria

Department of Chemical and Biomolecular Engineering

Clemson University

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics