1
|
Shintaku M, Oga H, Kusudo H, Smith ER, Omori T, Yamaguchi Y. Measuring line tension: Thermodynamic integration during detachment of a molecular dynamics droplet. J Chem Phys 2024; 160:224502. [PMID: 38856068 DOI: 10.1063/5.0201973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Accepted: 05/23/2024] [Indexed: 06/11/2024] Open
Abstract
The contact line (CL) is where solid, liquid, and vapor phases meet, and Young's equation describes the macroscopic force balance of the interfacial tensions between these three phases. These interfacial tensions are related to the nanoscale stress inhomogeneity appearing around the interface, and for curved CLs, e.g., a three-dimensional droplet, another force known as the line tension must be included in Young's equation. The line tension has units of force, acting parallel to the CL, and is required to incorporate the extra stress inhomogeneity around the CL into the force balance. Considering this feature, Bey et al. [J. Chem. Phys. 152, 094707 (2020)] reported a mechanical approach to extract the value of line tension τℓ from molecular dynamics (MD) simulations. In this study, we show a novel thermodynamics interpretation of the line tension as the free energy per CL length, and based on this interpretation, through MD simulations of a quasi-static detachment process of a quasi-two-dimensional droplet from a solid surface, we obtained the value τℓ as a function of the contact angle. The simulation scheme is considered to be an extension of a thermodynamic integration method, previously used to calculate the solid-liquid and solid-vapor interfacial tensions through a detachment process, extended here to the three-phase system. The obtained value agreed well with the result by Bey et al. and showed the validity of thermodynamic integration at the three-phase interface.
Collapse
Affiliation(s)
- Minori Shintaku
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Haruki Oga
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Hiroki Kusudo
- Department of Mechanical Systems Engineering, Tohoku University, 6-6-01 Aramaki, Aoba-ku, Sendai 980-8579, Japan
| | - Edward R Smith
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge UB8 3PH, United Kingdom
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Research Center (WaTUS), Research Institute for Science and Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| |
Collapse
|
2
|
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
|
3
|
Surblys D, Müller-Plathe F, Ohara T. Computing the Work of Solid-Liquid Adhesion in Systems with Damped Coulomb Interactions via Molecular Dynamics: Approaches and Insights. J Phys Chem A 2022; 126:5506-5516. [PMID: 35929812 PMCID: PMC9393893 DOI: 10.1021/acs.jpca.2c03934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Recently, the dry-surface method [Langmuir2016, 31, 8335−8345] has been developed
to compute the work of adhesion of solid–liquid and other interfaces
using molecular dynamics via thermodynamic integration. Unfortunately,
when long-range Coulombic interactions are present in the interface,
a special treatment is required, such as solving additional Poisson
equations, which is usually not implemented in generic molecular dynamics
software, or as fixing some groups of atoms in place, which is undesirable
most of the time. In this work, we replace the long-range Coulombic
interactions with damped Coulomb interactions, and explore several
thermal integration paths. We demonstrate that regardless of the integration
path, the same work of adhesion values are obtained as long as the
path is reversible, but the numerical efficiency differs vastly. Simple
scaling of the interactions is most efficient, requiring as little
as 8 sampling points, followed by changing the Coulomb damping parameter,
while modifying the Coulomb interaction cutoff length performs worst.
We also demonstrate that switching long-range Coulombic interactions
to damped ones results in a higher work of adhesion by about 10 mJ/m2 because of slightly different liquid molecule orientation
at the solid–liquid interface, and this value is mostly unchanged
for surfaces with substantially different Coulombic interactions at
the solid–liquid interface. Finally, even though it is possible
to split the work of adhesion into van der Waals and Coulomb components,
it is known that the specific per-component values are highly dependent
on the integration path. We obtain an extreme case, which demonstrates
that caution should be taken even when restricting to qualitative
comparison.
Collapse
Affiliation(s)
- Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, 980-8577, Japan
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, D-64287, Germany
| | - Taku Ohara
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, 980-8577, Japan
| |
Collapse
|
4
|
Watanabe K, Kusudo H, Bistafa C, Omori T, Yamaguchi Y. Quantifying the solid–fluid interfacial tensions depending on the substrate curvature: Young’s equation holds for wetting around nanoscale cylinder. J Chem Phys 2022; 156:054701. [DOI: 10.1063/5.0079816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Keitaro Watanabe
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Carlos Bistafa
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Takeshi Omori
- Deptartment of Mechanical Engineering, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Research Center (WaTUS), Research Institute for Science and Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| |
Collapse
|
5
|
Martinez MN, Smith AG, Nowack LM, Lin B, Rice SA. Interaction between dilute water vapor and dodecane thiol ligated Au nanoparticles: Hydrated structure and pair potential of mean force. J Chem Phys 2021; 155:144902. [PMID: 34654291 DOI: 10.1063/5.0065718] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The interaction between two ligated nanoparticles depends on whether they are isolated or immersed in a liquid solvent. However, very little is known about the influence of solvent vapor on the interaction between two ligated nanoparticles. Recent experiments yield the surprising result that the cyclic exposure of solvent free suspended monolayers of dodecane thiol ligated gold nanoparticles (AuNPs) to water vapor and dry nitrogen generates reversible cyclic decreases and increases in Young's modulus of the monolayer, implying corresponding cyclic changes in the AuNP-AuNP interaction. We examine how water vapor interacts with an isolated dodecane thiol dressed AuNP and how water vapor affects the interaction between a pair of nanoparticles, using all-atom molecular-dynamics simulations. We find that there is condensation of water molecules onto the ligand shell of an AuNP in the form of clusters of 100-2000 molecules that partially cover the shell, with most of the water in a few large clusters. A water cluster bridges the AuNPs, with a sensibly constant number of water molecules for AuNP-AuNP separations from the edge-to-edge contact up to center-to-center separations of 100 Å. The wet AuNP-AuNP interaction has a slightly deeper and wider asymmetric well than does the dry interaction, a change that is qualitatively consistent with that implied by the observed water vapor induced change in Young's modulus of a monolayer of these AuNPs. We find that macroscopic analyses of water drop-deformable surface interactions and dynamics provide both guidance to understanding and qualitatively correct predictions of the phenomena observed in our simulations.
Collapse
Affiliation(s)
- Michael N Martinez
- James Franck Institute, University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Alex G Smith
- James Franck Institute, University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Linsey M Nowack
- James Franck Institute, University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Binhua Lin
- James Franck Institute, University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Stuart A Rice
- James Franck Institute, University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| |
Collapse
|
6
|
Trentin LN, Pereira CS, Silveira RL, Hill S, Sorieul M, Skaf MS. Nanoscale Wetting of Crystalline Cellulose. Biomacromolecules 2021; 22:4251-4261. [PMID: 34515474 DOI: 10.1021/acs.biomac.1c00801] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cellulose possesses considerable potential for a wide range of sustainable applications. Nanocellulose-based material properties are primarily dependent on the structural surface characteristics of its crystalline planes. Experimental measurements of the affinity of crystalline nanocellulose surfaces with water are scarce and challenging to obtain. Therefore, the relative hydrophilicity of different cellulose allomorphs crystalline planes is often inferred from qualitative assessments of their surface and the exposition of polar groups to the solvent. This work investigates the relative hydrophilicity of cellulose surfaces using molecular dynamics simulations. The behavior of a water droplet laid on different crystal planes was used to determine their relative hydrophilicity. The water molecules fully spread onto highly hydrophilic surfaces. However, a water droplet placed on less hydrophilic surfaces equilibrates as an oblate spheroidal cap allowing the measurement of a contact angle. The results indicate that the Iα (010), Iα (11̅0), Iβ (010), and Iβ (110) faces, as well as the faces of human-made celluloses II and III_I (100), (11̅0), (010), and (110) are all highly hydrophilic. They all have a contact angle value inferior to 11°. Not unexpectedly, the Iα (001) and Iβ (100) surfaces are less hydrophilic with contact angles of 48 and 34°, respectively. However, the Iβ (11̅0) plane, often referred to as a hydrophilic surface, forms a contact angle of about 32°. The results are rationalized in terms of structure, exposure of hydroxyl groups to the solvent, and degree of cellulose-cellulose versus cellulose-water hydrogen bonds on each face. The simulations also show that the surface oxidation degree tunes the surface hydrophilicity in a nonlinear manner due to cooperative effects involving water-cellulose interactions. Our study helps us to understand how the degree of hydrophilicity of cellulose emerges from specific structural features of each crystalline surface.
Collapse
Affiliation(s)
- Lucas N Trentin
- Institute of Chemistry and Center for Computing in Engineering & Sciences, University of Campinas, Campinas, São Paulo 13084-862, Brazil
| | - Caroline S Pereira
- Institute of Chemistry and Center for Computing in Engineering & Sciences, University of Campinas, Campinas, São Paulo 13084-862, Brazil
| | - Rodrigo L Silveira
- Institute of Chemistry and Center for Computing in Engineering & Sciences, University of Campinas, Campinas, São Paulo 13084-862, Brazil.,Institute of Chemistry, Federal University of Rio de Janeiro, Rio de Janeiro, RJ 21941-909, Brazil
| | - Stefan Hill
- Scion, Private Bag 3020, Rotorua 3046, New Zealand
| | | | - Munir S Skaf
- Institute of Chemistry and Center for Computing in Engineering & Sciences, University of Campinas, Campinas, São Paulo 13084-862, Brazil
| |
Collapse
|
7
|
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
|
8
|
Granados-Bazán EL, Quiñones-Cisneros SE, Deiters UK. Structure and Contact Angle in Sessile Droplets of Binary Mixtures of Lennard-Jones Chains: A Molecular Dynamics Study. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:10945-10957. [PMID: 34478317 DOI: 10.1021/acs.langmuir.1c01354] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Molecular dynamics simulations were carried out to investigate cylindrical droplets consisting of binary mixtures of Lennard-Jones (LJ) fluids in contact with a solid substrate. The droplets are composed of mixtures of the monomeric LJ fluid plus linear-tangent chains of 2, 10, 20, and 30 segments per chain that interact through a harmonic potential and the spherically truncated and shifted potential Lennard-Jones. The solid surface was modeled as a semi-infinite platinum substrate with an FCC structure that interacts with the fluid by means of a LJ 9-3 potential. We place emphasis on the effect of mixing a monomeric LJ fluid with heavy components on the contact angle and on the droplet structure, especially in the liquid-solid region. The density profiles of the droplets reveal a strong discrete layering of the fluid in the vicinity of the solid-liquid interface. The layering is more pronounced at low temperatures and for mixtures of short chains (symmetric mixtures). The ordering of the fluid was much less intense for fluids of long chains (asymmetric mixtures), and some cases even show gas enrichment at the solid-liquid interface. Enrichment at the vapor-liquid interfaces and density inversion can also be observed. However, these effects are not as marked as in planar interfaces. The contact angle between the droplet and the substrate is calculated by fitting an ellipse to the vapor-liquid interface defined by the Gibbs dividing surface. In general, an increment in the concentration of the heavy component and a reduction of the temperature resulted in an increase of the contact angle, which in turn disfavored the wetting of the droplet.
Collapse
Affiliation(s)
- Eder L Granados-Bazán
- Institute of Physical Chemistry, University of Cologne, Luxemburger Str. 116, 50939 Köln, Germany
| | - Sergio E Quiñones-Cisneros
- Institute of Thermo- and Fluid Dynamics, Ruhr-University Bochum, Universitätsstr. 150, 44801 Bochum, Germany
| | - Ulrich K Deiters
- Institute of Physical Chemistry, University of Cologne, Luxemburger Str. 116, 50939 Köln, Germany
| |
Collapse
|
9
|
Bistafa C, Surblys D, Kusudo H, Yamaguchi Y. Water on hydroxylated silica surfaces: Work of adhesion, interfacial entropy, and droplet wetting. J Chem Phys 2021; 155:064703. [PMID: 34391348 DOI: 10.1063/5.0056718] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the last few years, much attention has been devoted to the control of the wettability properties of surfaces modified with functional groups. Molecular dynamics (MD) simulation is one of the powerful tools for microscopic analysis providing visual images and mean geometrical shapes of the contact line, e.g., of nanoscale droplets on solid surfaces, while profound understanding of wetting demands quantitative evaluation of the solid-liquid (SL) interfacial tension. In the present work, we examined the wetting of water on neutral and regular hydroxylated silica surfaces with five different area densities of OH groups ρA OH, ranging from a non-hydroxylated surface to a fully hydroxylated one through two theoretical methods: thermodynamic integration (TI) and MD simulations of quasi-two-dimensional equilibrium droplets. For the former, the work of adhesion needed to quasi-statically strip the water film off the solid surface was computed by the phantom wall TI scheme to evaluate the SL interfacial free energy, whereas for the latter, the apparent contact angle θapp was calculated from the droplet density distribution. The theoretical contact angle θYD and the apparent one θapp, both indicating the enhancement of wettability by an increase in ρA OH, presented good quantitative agreement, especially for non-hydroxylated and highly hydroxylated surfaces. On partially hydroxylated surfaces, in which θYD and θapp slightly deviated, the Brownian motion of the droplet was suppressed, possibly due to the pinning of the contact line around the hydroxyl groups. Relations between work of adhesion, interfacial energy, and entropy loss were also analyzed, and their influence on the wettability was discussed.
Collapse
Affiliation(s)
- Carlos Bistafa
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| |
Collapse
|
10
|
Dewetting transition of water on nanostructured and wettability patterned surfaces: A molecular dynamics study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116869] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
11
|
Abstract
Fluid interfaces with nanoscale radii of curvature are generating great interest, both for their applications and as tools to probe our fundamental understanding. One important question is what is the smallest radius of curvature at which the three main thermodynamic combined equilibrium equations are valid: the Kelvin equation for the effect of curvature on vapor pressure, the Gibbs-Thomson equation for the curvature-induced freezing point depression, and the Ostwald-Freundlich equation for the curvature-induced increase in solubility. The objective of this Perspective is to provide conceptual, molecular modeling, and experimental support for the validity of these thermodynamic combined equilibrium equations down to the smallest interfacial radii of curvature. Important concepts underpinning thermodynamics, including ensemble averaging and Gibbs's treatment of bulk phase heterogeneities in the region of an interface, give reason to believe that these equations might be valid to smaller scales than was previously thought. There is significant molecular modeling and experimental support for all three of the Kelvin equation, the Gibbs-Thomson equation, and the Ostwald-Freundlich equation for interfacial radii of curvature from 1 to 4 nm. There is even evidence of sub-nanometer quantitative accuracy for the Kelvin equation and the Gibbs-Thomson equation.
Collapse
Affiliation(s)
- Janet A W Elliott
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 1H9, Canada
| |
Collapse
|
12
|
Wang P, He L, Wang Z. The effect of surface structure and arrangement on wettability of substrate surface. Colloids Surf A Physicochem Eng Asp 2021. [DOI: 10.1016/j.colsurfa.2021.126165] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
13
|
Ding W, Han D, Zhang J, Ma Q, Li X, Zhang J, Wang X. Molecular dynamics study of anisotropic behaviours of water droplet on textured surfaces with various energies. Mol Phys 2021. [DOI: 10.1080/00268976.2020.1785028] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Wenyang Ding
- Institute of Thermal Science and Technology, Shandong University, Jinan, People’s Republic of China
| | - Dan Han
- Institute of Thermal Science and Technology, Shandong University, Jinan, People’s Republic of China
| | - Jingzhi Zhang
- School of Energy and Power Engineering, Shandong University, Jinan, People’s Republic of China
| | - Qingming Ma
- Department of Pharmaceutics, School of Pharmacy, Qingdao University, Qingdao, People’s Republic of China
| | - Xiaoyan Li
- Shandong Provincial Key Laboratory of Oral Tissue Regeneration, School of Stomatology, Shandong University, Jinan, People’s Republic of China
| | - Jingchao Zhang
- Holland Computing Center, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Xinyu Wang
- Institute of Thermal Science and Technology, Shandong University, Jinan, People’s Republic of China
| |
Collapse
|
14
|
Imaizumi Y, Omori T, Kusudo H, Bistafa C, Yamaguchi Y. Wilhelmy equation revisited: A lightweight method to measure liquid–vapor, solid–liquid, and solid–vapor interfacial tensions from a single molecular dynamics simulation. J Chem Phys 2020; 153:034701. [DOI: 10.1063/5.0011979] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuta Imaizumi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Carlos Bistafa
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Science & Technology Research Center (W-FST), Research Institute for Science & Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| |
Collapse
|
15
|
Kusudo H, Omori T, Yamaguchi Y. Extraction of the equilibrium pinning force on a contact line exerted from a wettability boundary of a solid surface through the connection between mechanical and thermodynamic routes. J Chem Phys 2019; 151:154501. [DOI: 10.1063/1.5124014] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
16
|
Jiménez-Serratos G, Cárdenas H, Müller EA. Extension of the effective solid-fluid Steele potential for Mie force fields. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1669836] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
| | - Harry Cárdenas
- Department of Chemical Engineering, Imperial College London, London, UK
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, London, UK
| |
Collapse
|
17
|
Abramov A, Iglauer S. Analysis of individual molecular dynamics snapshots simulating wetting of surfaces using spheroidal geometric constructions. J Chem Phys 2019. [DOI: 10.1063/1.5113852] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Aleksandr Abramov
- School of Engineering, Edith Cowan University, 270 Joondalup Drive, Joondalup, WA 6027 Western Australia, Australia
| | - Stefan Iglauer
- School of Engineering, Edith Cowan University, 270 Joondalup Drive, Joondalup, WA 6027 Western Australia, Australia
| |
Collapse
|
18
|
Omori T, Kobayashi Y, Yamaguchi Y, Kajishima T. Understanding the asymmetry between advancing and receding microscopic contact angles. SOFT MATTER 2019; 15:3923-3928. [PMID: 31011723 DOI: 10.1039/c9sm00521h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
By means of molecular dynamics simulation, the advancing and receding microscopic contact angles were analyzed for a shear flow of two mono-atomic fluids confined between parallel non-polar solid walls. We defined the microscopic dynamic contact angle based on the coarse-grained microscopic density distribution of the fluids (the instantaneous interface method [Willard and Chandler, J. Phys. Chem. B, 2010, 114, 1954-1958]) near the moving contact line. We have found that the asymmetric change of fluid density near the wall with respect to the moving contact line results in a different dependence between the advancing and receding contact angles on the contact line velocity in a system where the two fluids across the interface have unequal wettability to the solid wall. This difference between the advancing and receding contact angles leads to different flow resistance caused by the advancing and receding contact lines, which should have impact on the industrial applications of the fine fluid transportation with contact lines.
Collapse
Affiliation(s)
- T Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan.
| | | | | | | |
Collapse
|
19
|
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]
|
20
|
Jiang H, Fialoke S, Vicars Z, Patel AJ. Characterizing surface wetting and interfacial properties using enhanced sampling (SWIPES). SOFT MATTER 2019; 15:860-869. [PMID: 30644500 DOI: 10.1039/c8sm02317d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We introduce an accurate and efficient method for characterizing surface wetting and interfacial properties, such as the contact angle made by a liquid droplet on a solid surface, and the vapor-liquid surface tension of a fluid. The method makes use of molecular simulations in conjunction with the indirect umbrella sampling technique to systematically wet the surface and estimate the corresponding free energy. To illustrate the method, we study the wetting of a family of Lennard-Jones surfaces by water. For surfaces with a wide range of attractions for water, we estimate contact angles using our method, and compare them with contact angles obtained using droplet shapes. Notably, our method is able to capture the transition from partial to complete wetting as surface-water attractions are increased. Moreover, the method is straightforward to implement and is computationally efficient, providing accurate contact angle estimates in roughly 5 nanoseconds of simulation time.
Collapse
Affiliation(s)
- Hao Jiang
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA.
| | | | | | | |
Collapse
|
21
|
Yamaguchi Y, Kusudo H, Surblys D, Omori T, Kikugawa G. Interpretation of Young’s equation for a liquid droplet on a flat and smooth solid surface: Mechanical and thermodynamic routes with a simple Lennard-Jones liquid. J Chem Phys 2019; 150:044701. [DOI: 10.1063/1.5053881] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Science and Technology Research Center (W-FST), Research Institute for Science and Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| |
Collapse
|
22
|
Essafri I, Le breton JC, Saint-Jalmes A, Soldera A, Szymczyk A, Malfreyt P, Ghoufi A. Contact angle and surface tension of water on a hexagonal boron nitride monolayer: a methodological investigation. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1502427] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Ilham Essafri
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), Rennes, France
| | | | | | - Armand Soldera
- Laboratory of Physical Chemistry of Matter (LPCM), Department of Chemistry, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Anthony Szymczyk
- Univ Rennes, CNRS, ISCR (Institut des Sciences Chimiques de Rennes), Rennes, France
| | - Patrice Malfreyt
- Universite Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), Clermont-Ferrand, France
| | - Aziz Ghoufi
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), Rennes, France
| |
Collapse
|