1
|
Carlson S, Schullian O, Becker MR, Netz RR. Modeling Water Interactions with Graphene and Graphite via Force Fields Consistent with Experimental Contact Angles. J Phys Chem Lett 2024; 15:6325-6333. [PMID: 38856977 PMCID: PMC11194815 DOI: 10.1021/acs.jpclett.4c01143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 05/31/2024] [Accepted: 06/04/2024] [Indexed: 06/11/2024]
Abstract
Accurate simulation models for water interactions with graphene and graphite are important for nanofluidic applications, but existing force fields produce widely varying contact angles. Our extensive review of the experimental literature reveals extreme variation among reported values of graphene-water contact angles and a clustering of graphite-water contact angles into groups of freshly exfoliated (60° ± 13°) and not-freshly exfoliated graphite surfaces. The carbon-oxygen dispersion energy for a classical force field is optimized with respect to this 60° graphite-water contact angle in the infinite-force-cutoff limit, which in turn yields a contact angle for unsupported graphene of 80°, in agreement with the mean of the experimental results. Interaction force fields for finite cutoffs are also derived. A method for calculating contact angles from pressure tensors of planar equilibrium simulations that is ideally suited to graphite and graphene surfaces is introduced. Our methodology is widely applicable to any liquid-surface combination.
Collapse
Affiliation(s)
- Shane
R. Carlson
- Fachbereich
Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
| | - Otto Schullian
- Fachbereich
Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
- Department
of Biomaterials, Max Planck Institute of
Colloids and Interfaces, D-14424 Potsdam, Germany
| | - Maximilian R. Becker
- Fachbereich
Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
| |
Collapse
|
2
|
Saito T, Kubo M, Tsukada T, Shoji E, Kikugawa G, Surblys D, Kubo M. Molecular dynamics simulations for interfacial structure and affinity between carboxylic acid-modified Al2O3 and polymer melts. J Chem Phys 2023; 159:164708. [PMID: 37888762 DOI: 10.1063/5.0169721] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 10/10/2023] [Indexed: 10/28/2023] Open
Abstract
Controlling the dispersion state of nanoparticles in a polymer matrix is necessary to produce polymer nanocomposites. The surface modification of nanoparticles is used to enable their dispersion in polymers. Moreover, molecular dynamics (MD) simulations are useful for revealing the interfacial properties between nanoparticles and polymers to aid in the design of materials. In this study, the effect of surface coverage, modifier length, and polymer species on the interfacial structure and affinity between surface-modified Al2O3 and polymer melts were investigated using all-atom MD simulations. Hexanoic, decanoic, and tetradecanoic acids were used as surface modifiers, and polypropylene (PP), polystyrene (PS), and poly (methyl methacrylate) (PMMA) were used as polymers. The work of adhesion Wadh and the work of immersion Wimm were selected as quantitative measures of affinity. Wadh was calculated using the phantom-wall approach, and Wimm was calculated by simply subtracting the surface tension of polymers γL from Wadh. The results showed that Wadh and Wimm were improved by surface modification with low coverage, owing to a good penetration of the polymer. The effect of modifier length on Wadh and Wimm was small. Whereas Wadh increased in the following order: PP < PS < PMMA, Wimm increased as follows: PMMA < PS < PP. Finally, the trend of Wadh and Wimm was organized using the Flory-Huggins interaction parameter χ between the modifier and the polymer. This study demonstrates that the interfacial affinity can be improved by tuning the surface coverage and modifier species depending on the polymer matrix.
Collapse
Affiliation(s)
- Takamasa Saito
- Department of Chemical Engineering, Tohoku University, 6-6-07 Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Masaki Kubo
- Department of Chemical Engineering, Tohoku University, 6-6-07 Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Takao Tsukada
- New Industry Creation Hatchery Center, Tohoku University, 6-6-10 Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Eita Shoji
- Department of Mechanical Systems Engineering, Tohoku University, 6-6-01 Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan
| | - Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan
| | - Momoji Kubo
- Institute for Materials Research, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan
| |
Collapse
|
3
|
Marks SM, Vicars Z, Thosar AU, Patel AJ. Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling. J Phys Chem B 2023. [PMID: 37378637 DOI: 10.1021/acs.jpcb.3c01627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.
Collapse
Affiliation(s)
- Sean M Marks
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zachariah Vicars
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Aniket U Thosar
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amish J Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
4
|
Alizadeh Sahraei A, Azizi D, Mokarizadeh AH, Boffito DC, Larachi F. Emerging Trends of Computational Chemistry and Molecular Modeling in Froth Flotation: A Review. ACS ENGINEERING AU 2023; 3:128-164. [PMID: 37362006 PMCID: PMC10288516 DOI: 10.1021/acsengineeringau.2c00053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2022] [Revised: 04/04/2023] [Accepted: 04/06/2023] [Indexed: 06/28/2023]
Abstract
Froth flotation is the most versatile process in mineral beneficiation, extensively used to concentrate a wide range of minerals. This process comprises mixtures of more or less liberated minerals, water, air, and various chemical reagents, involving a series of intermingled multiphase physical and chemical phenomena in the aqueous environment. Today's main challenge facing the froth flotation process is to gain atomic-level insights into the properties of its inherent phenomena governing the process performance. While it is often challenging to determine these phenomena via trial-and-error experimentations, molecular modeling approaches not only elicit a deeper understanding of froth flotation but can also assist experimental studies in saving time and budget. Thanks to the rapid development of computer science and advances in high-performance computing (HPC) infrastructures, theoretical/computational chemistry has now matured enough to successfully and gainfully apply to tackle the challenges of complex systems. In mineral processing, however, advanced applications of computational chemistry are increasingly gaining ground and demonstrating merit in addressing these challenges. Accordingly, this contribution aims to encourage mineral scientists, especially those interested in rational reagent design, to become familiarized with the necessary concepts of molecular modeling and to apply similar strategies when studying and tailoring properties at the molecular level. This review also strives to deliver the state-of-the-art integration and application of molecular modeling in froth flotation studies to assist either active researchers in this field to disclose new directions for future research or newcomers to the field to initiate innovative works.
Collapse
Affiliation(s)
- Abolfazl Alizadeh Sahraei
- Department
of Chemical Engineering, Université
Laval, 1065 Avenue de la Médecine, Québec, Québec G1V 0A6, Canada
| | - Dariush Azizi
- Department
of Chemical Engineering, École Polytechnique
de Montréal, 2900 Boulevard Édouard-Montpetit, Montréal H3T 1J4, Canada
| | - Abdol Hadi Mokarizadeh
- School
of Polymer Science and Polymer Engineering, University of Akron, Akron, Ohio 44325, United States
| | - Daria Camilla Boffito
- Department
of Chemical Engineering, École Polytechnique
de Montréal, 2900 Boulevard Édouard-Montpetit, Montréal H3T 1J4, Canada
| | - Faïçal Larachi
- Department
of Chemical Engineering, Université
Laval, 1065 Avenue de la Médecine, Québec, Québec G1V 0A6, Canada
| |
Collapse
|
5
|
Hong X, Yu H, Xu H, Wang X, Jin X, Wu H, Wang F. Competitive adsorption of asphaltene and n-heptane on quartz surfaces and its effect on crude oil transport through nanopores. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119312] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
6
|
Lbadaoui-Darvas M, Garberoglio G, Karadima KS, Cordeiro MNDS, Nenes A, Takahama S. Molecular simulations of interfacial systems: challenges, applications and future perspectives. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1980215] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Affiliation(s)
- Mária Lbadaoui-Darvas
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
| | - Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (FBK-ECT*), Trento, Italy
- Trento Institute for Fundamental Physics and Applications (TIFPA-INFN), Trento, Italy
| | - Katerina S. Karadima
- Department of Chemical Engineering, University of Patras, Patras, Greece
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas(FORTH-ICE/HT), Patras, Greece
| | | | - Athanasios Nenes
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas(FORTH-ICE/HT), Patras, Greece
| | - Satoshi Takahama
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
7
|
Saito T, Shoji E, Kubo M, Tsukada T, Kikugawa G, Surblys D. Evaluation of the work of adhesion at the interface between a surface-modified metal oxide and an organic solvent using molecular dynamics simulations. J Chem Phys 2021; 154:114703. [PMID: 33752377 DOI: 10.1063/5.0040900] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Advancing the practical applications of surface-modified nanoparticles requires that their dispersion in solvents can be controlled. The degree of dispersion depends on the affinity between surface-modified nanoparticles and solvents, which can be quantified using the work of adhesion at the interface. Herein, the affinity between a surface-modified inorganic solid and an organic solvent was evaluated by calculating the work of adhesion at the interface. The phantom-wall method, which is a thermodynamic route for evaluating the work of adhesion at an interface using molecular dynamics simulations, was applied to the decanoic acid-modified Al2O3/hexane interface. Molecular dynamics simulations were performed for flat interface systems to focus on the interactions between substances that affect the affinity on the surface. As a result, the surface coverage of decanoic acid was found to affect the work of adhesion, with a maximum value of 45.66 ± 0.75 mJ/m2 at a surface coverage of 75%. An analysis of the mass density profiles of Al2O3, decanoic acid, and hexane in the vicinity of the interface showed that the increase in the work of adhesion with the surface coverage was due to the penetration of hexane molecules into the decanoic acid layer on the Al2O3 surface. At a surface coverage of 75%, some hexane molecules were trapped in the layer of oriented decanoic acid molecules. These results suggested that the interfacial affinity can be enhanced by controlling the surface modification so that the solvent can penetrate the layer of the modifier.
Collapse
Affiliation(s)
- Takamasa Saito
- Department of Chemical Engineering, Tohoku University, 6-6-07 Aramaki, Aoba-ku, Sendai 980-8579, Japan
| | - Eita Shoji
- Department of Mechanical Systems Engineering, Tohoku University, 6-6-01 Aramaki, Aoba-ku, Sendai 980-8579, Japan
| | - Masaki Kubo
- Department of Chemical Engineering, Tohoku University, 6-6-07 Aramaki, Aoba-ku, Sendai 980-8579, Japan
| | - Takao Tsukada
- Department of Chemical Engineering, Tohoku University, 6-6-07 Aramaki, Aoba-ku, Sendai 980-8579, Japan
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| |
Collapse
|
8
|
Playing the long game wins the cohesion-adhesion rivalry. Proc Natl Acad Sci U S A 2019; 116:23874-23876. [PMID: 31659047 DOI: 10.1073/pnas.1916911116] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
9
|
Recent advances in estimating contact angles using molecular simulations and enhanced sampling methods. Curr Opin Chem Eng 2019. [DOI: 10.1016/j.coche.2019.03.012] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|