Direct Prediction of Calcite Surface Wettability with First-Principles

Oct 6, 2017 - Prediction of intrinsic surface wettability from first-principles offers great opportunities in probing new physics of natural phenomena...
0 downloads 16 Views 794KB Size
Subscriber access provided by University of Virginia Libraries & VIVA (Virtual Library of Virginia)

Letter

Direct Prediction of Calcite Surface Wettability with First-Principles Quantum Simulation Jin-You Lu, Qiaoyu Ge, Hongxia Li, Aikifa Raza, and TieJun Zhang J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b02270 • Publication Date (Web): 06 Oct 2017 Downloaded from http://pubs.acs.org on October 7, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Direct Prediction of Calcite Surface Wettability with First-principles Quantum Simulation Jin You Lu, Qiaoyu Ge, Hongxia Li, Aikifa Raza, TieJun Zhang* Department of Mechanical and Materials Engineering, Masdar Institute, Khalifa University of Science and Technology, P.O. Box 54224, Abu Dhabi, UAE * Address correspondence to: [email protected]

Abstract: Prediction of intrinsic surface wettability from first principles offers great opportunities in probing new physics of natural phenomena and enhancing energy production or transport efficiency. In this paper, we propose a general quantum mechanical approach to predict the macroscopic wettability of any solid crystal surfaces for different liquids directly through atomic-level density functional simulation.

As a benchmark, the wetting characteristics of

calcite crystal (10.4) under different types of fluids (water, hexane, and mercury), including either contact angle or spreading coefficient, are predicted and further validated with experimental measurements. A unique feature of our approach lies in its capability of capturing the interactions among various polar fluid molecules and solid surface ions, particularly their charge density difference distributions. Moreover, this approach provides insightful and quantitative predictions of complicated surface wettability alternation problems and wetting behaviors of liquid/liquid/solid tri-phase systems.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 23

TOC GRAPHICS

Surface wettability, the tendency of how a liquid wets a solid surface is an indispensable and intriguing topic in nature and many practical applications, such as enhanced oil recovery1–3, water treatment4, nanomaterials synthesis and fabrication5,6, heat transfer7–9, and self-cleaning applications10. When an immiscible liquid droplet comes into contact with an impermeable solid surface, the surface adhesion force at the solid/liquid interface competes with the intermolecular cohesion force inside the liquid droplet. Consequently, two regimes of wetting behaviors could be observed: total wetting with the droplet spreading over the solid surface11 or partial wetting with a spherical liquid cap forming on the surface12. Theoretical predictions of surface wettability have been developed, and the most widely used approach is based on the semi-empirical Van Oss-Chaudhury-Good theory13. This approach works well for low surface energy materials, such as polymers14,15, however, it is not applicable to high energy surfaces like oxides and hydroxides. Later on, numerical molecular dynamics (MD) simulation studies16–18 have been carried out to predict the surface wettability for various problems, such as a water droplet on graphene19–21. However, it is still challenging for classical MD simulation to choose tens of suitable parameters in the force potential model so as to

ACS Paragon Plus Environment

2

Page 3 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

characterize interatomic interactions at the solid/liquid interface and predict wetting behaviors of a liquid/liquid/solid tri-phase system. For the past thirty years, density functional theory (DFT) has proven to be a dominant method for the electronic structure calculations of periodic systems22. With DFT simulation, insights into the intermolecular interactions at the solid/liquid interface can be obtained. A great advantage of DFT simulation is that it mainly relies on the intrinsic electronic configurations of atoms. When compared to MD simulation, much fewer input parameters are required in DFT simulations. Recently, systematic DFT studies23–27 have effectively predicted the wettability changes of calcite/aragonite surfaces under different ion substitutions. However, an effective method is still unresolved to directly predict the absolute (or intrinsic) surface wettability of a solid surface by different liquids from DFT simulations. In the present study, we propose a new approach to directly predict the intrinsic wettability of a solid surface based on a DFT slab model. The scheme is designed as follows: the wetting behavior of a liquid droplet on a solid surface, in equilibrium with the liquid vapor as illustrated in Fig. 1(a), can be described by the Young- Dupré equation28,29 WSL   SV + LV - SL = LV (1  cos( ))

(1)

where WSL is the work of adhesion at the solid/liquid interface,  SV ,  SL , and  LV are the interfacial tensions at the solid/vapor, solid/liquid, and liquid/vapor interfaces, respectively, and  is the contact angle. The Young-Dupré equation connects the vapor-liquid-solid contact angle

to the work of adhesion at the liquid/solid interface WSL and the surface tension of the liquid  LV . In order to predict the WSL and  LV , we need two related physical quantities, the adhesion energy at the liquid/solid interface U SL and the cohesion energy between the liquid layers U LV / 2 , which can be obtained directly from our DFT simulation. Then, the WSL and  LV are calculated by dividing the adhesion and cohesion energy by the surface area of a liquid molecule, i.e.,

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

WSL  U SL / a 2 and  LV  U LV / (2a 2 ) , respectively, where ܽ is the molecular length. Here, the cross-

section area of a liquid molecule at the liquid/solid interface is assumed to be same as that at the liquid/liquid interface, as indicated by our studies. Therefore, the value of WSL /  LV is obtained by calculating USL / (U LV / 2) from DFT simulations, which solves the contact angle of Eq.(1).

(a)

Liquid(L)

 LV Vapor(V)  SV

U adm 3~   

θ Liquid(L) SL U adm 1  U SL

WSL   SV   LV   SL Solid (S)

(b)

U LV 2

Liquid (L)

Solid (S)

(c)

(d)

m Ucoh

UL

U mML

Um Top liquid layer remains

Solid (S)

m U coh  U m  nL  U L

Top view

Figure 1. (a) Schematic diagram to predict the surface wettability of a solid by a liquid. The work of adhesion WSL at the solid/liquid interface and the interfacial tension  LV at liquid/vapor interface can be predicted with the DFT solid/liquid slab model.

m U ad

is the adhesion energy of

the mth liquid layer on the solid surface covered with (m-1) layers of liquids. (b-d) the procedures to calculate the cohesion energy

m U coh

within the individual liquid (mth) monolayer. (b) Lowest

energy configuration of the solid surface with m liquid layers, (c) remove all atoms below top layer of liquid molecules and calculate Um for this layer, and (d) Calculate

m U coh

for the top layer (a

top view is shown) using Eq. (3).

ACS Paragon Plus Environment

4

Page 5 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

The following procedures are used to obtain USL and U LV / 2 from the DFT total energy studies on the solid/liquid system. The total energy includes electron kinetic energy, the ionelectron/ion-ion/electron-electron Coulomb potentials, exchange-correlation energy, and van der Waals dispersion correction. Layer-by-layer films of liquid molecules are added on the solid surface with the lowest energy configuration. In order to ensure the configuration reaches the overall minimum energy state, the total energies of a number of possible initial configurations were calculated by performing the DFT geometry optimization30. The detailed DFT simulation parameters and how initial configurations were generated by ab initio MD simulations31 are included in the Supporting Information (S.I). Then, the adhesion energy of mth liquid layer on the solid surface covered with (m-1) layers of liquids is calculated using Eq. (2), U adm 

U S  mL  U S  ( m 1) L  U mL nL

(2)

U where U S  mL is the total energy of a solid surface and m liquid layers above the surface, S ( m1) L is

the total energy of a solid surface and (m-1) overlaying liquid layers, U mL is the total energy of the mth liquid monolayer, and nL is the number of molecules in one liquid monolayer. U mL

is calculated as follows: by removing the atoms below the mth liquid layer from the

lowest energy solid/liquid configuration, the remaining top monolayer is calculated through a self-consistent calculation without geometry optimization. In this way, the total energy of mth m 1 liquid monolayer is obtained. When m is equal to one, the adhesion energy U ad is U SL , which is

the adhesion energy at the solid/liquid interface. As m is larger than one, the adhesion energy m U ad

of the mth liquid layer on the solid surface covered with (m-1) layers of liquids is the half

cohesion energy ~ U LV / 2 between two liquid monolayers, as shown in Fig. 1(a) and Eq. (3). When calculating cohesion energy between liquid layers, we take the averaged value of the

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 23

adhesion energy from the third liquid layer to top layer, to avoid the influence of the solid surface, as indicated in this study. U m  3~   LV U adm 1  U SL U ad 2 ;

(3)

Then, a governing equation is given to predict the vapor/liquid/solid contact angle based on conventional DFT simulations, as shown in Eq, (4). cos( ) 

U SL  U LV / 2 U adm 1  U adm 3~   U LV / 2 U adm 3~ 

(4)

In addition to the adhesion energy of the respective liquid layer on the solid or liquid covering solid surfaces, the cohesion energy

m U coh

within the individual mth liquid monolayer can be

calculated using Eq. (5), as illustrated in Fig. 1(b-d). m U coh  U mL  nL  U L

(5)

where U mL is the total energy of the mth liquid monolayer, nL is the number of molecules in one liquid monolayer, and the U L is the total energy of a single liquid molecule. By calculating this cohesion energy within the liquid layer, we can investigate the effect of solid surface on the intermolecular interaction within individual liquid layers. The proposed method is applied to predict the surface wettability of the calcite (10.4) surface for three different liquids, including water, hexane, and mercury. The plane (10.4) in calcite crystals is the most thermodynamically stable. We chose calcite because of its wide presence on earth, serving as a pH buffer and concentration of dissolved calcium in geological energy and environmental systems2,32,33.

ACS Paragon Plus Environment

6

Page 7 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 2(a-b). Two different side views for the lowest energy configuration of a four-layered water on a calcite (10.4) surface. Each water layer consists of four H2O molecules, and each calcite layer has four CaCO3. Blue spheres are the Ca2+ ions, and the trigonal planar CO32- ions are shown with golden carbon and red oxygen. A water molecule consists of two white hydrogen atoms and one red oxygen atom. (c) Charge density difference at the water/calcite interface. (d) Adhesion energy of the mth water layer on the calcite (10.4) surface covered with (m-1) layers of liquids calculated using Eq. (2). (e) Cohesion energy within individual mth water layers on the calcite (10.4) surface calculated using Eq. (3). We start with the study of the calcite surface (10.4) for water molecules. The details of how to get the lowest energy configuration of calcite surface (10.4) in simulation domain are included in the Supporting Information. By comparing total energies among different configurations, we have found that four water molecules behaving as one monolayer have the lowest total energy and individual water layer arranges in a crisscross pattern (Fig. 2(a-b)). The

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 23

lowest energy configuration and the charge density difference distribution (Fig. 2(c)) shows that a strong binding at the water/calcite interface is formed because the negatively polarized oxygen atoms of water molecules are coupled to the positively charged calcium ion. The top view at the water/calcite interface and the lowest energy configuration of eight-layered water on the calcite surface are provided in Fig. S1(a-c) of Supporting Information. m The adhesion energy U ad predicted by Eq. (2) changes dramatically as the number of

water layers increases, as shown in Fig. 2(c). As the number of water layers increases, the adhesion energy of mth water layer on the calcite surface covered with (m-1) layers of water slightly oscillates around a constant value. The adhesion energy of the water layer on the wet calcite surface is essentially the cohesion energy between two water monolayers. The oscillation of the adhesion energy for higher liquid layers comes from the strong hydrogen bonding among water molecules, as indicated by the dashed lines in Fig. 2(a-b). In addition, the cohesion energy m U coh

within individual mth water layers is predicted by using Eq. (3). Figure 2(e) shows extremely

weak intermolecular interactions within the first layer on the calcite surface owing to the strong dipole binding at the calcite/water interface. Moreover, the DFT simulation results are consistent with each other under different pseudo-potentials, namely, norm-conserving (NC), ultrasoft (US), and projector augmented waves (PAW) sets. This consistency indicates that different approximations of electrons at the inner shells give essentially identical results. According to the Eq. (5), the predicted contact angles of a water droplet on the calcite surface are 43.1o, 36.3o, and 36.6o, for DFT simulations with US, NC, and PAW pseudopotentials, respectively. The contact angle predictions match our contact angle measurement of sessile water drop on the calcite (10.4) surface with an average result of 40.8 ±3.97° (Fig. S4(a) in Supporting Information), also in agreement with the ref.34, as summarized in Table 1.

ACS Paragon Plus Environment

8

Page 9 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Table 1. Prediction results and experiments of the calcite (10.4) surface for water, hexane, and mercury. Species Water Hexane Mercury Hexane

ULV/2 (kJ/mol) USL (kJ/mol) Contact Angle (o) ULV/2 (kJ/mol) USL (kJ/mol) Contact Angle (o) ULV/2 (kJ/mol) USL (kJ/mol) Contact Angle (o) S/γLV*

PBE-US 57.7 99.8 43.1 19.2 53.2 0 89.2 34.3 128.1 0.80

PBE-NC 52.6 95.0 36.3 25.7 36.3 0 31.7 25.4 101.6 0.93

PBE-PAW 53.5 96.5 36.6 18.4 36.6 0 89.1 34.6 127.8 0.86

Experiment 40.8 ±3.97 0 136+3/133+434 >0

*a dimensionless quantity obtained from spreading coefficient divided by the surface tension of hexane

Figure 3(a-b). Two different side views for the lowest energy configuration for a four-layered hexane on the calcite (10.4) surface. Each layer consists of two hexane molecules, as indicated by MD generated configurations, and each hexane molecule has fourteen white hydrogen atoms

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

and six golden carbon atoms. (c) Charge density difference at the hexane/calcite interface. (d) Adhesion energy of the mth hexane layer on the calcite surface covered with (m-1) hexane layers. m m (e) Cohesion energy within the mth hexane layer on calcite surface. (f-g) U ad and U coh for mercury

liquid on the calcite surface, respectively. As a non-polar simple hydrocarbon liquid, hexane (C6H14)

is a suitable candidate

representing crude oils for studying the work of adhesion at the liquid/calcite interface, which provides the insights into subsurface wettability alteration for enhancing oil recovery35. By comparing a number of MD generated configurations, the lowest energy configuration shows two hexane molecules act as one monolayer on the calcite surface, as shown in Fig. 3(a-b). The charge density difference distribution at the hexane/calcite interface in Fig. 3(c) is negligible as compared to water molecules on calcite surface (Fig 2.(c)) owing to the non-polar nature of hexane molecules. It should be noted that, for easy comparison, we use the same scale when plotting electron difference density distributions for different liquids on the calcite surface. Figure 3(d) shows that the adhesion strength at the calcite/hexane interface is larger than that at the liquid/liquid interface by three folds owing to the low interfacial tension of the hexane liquid. The cohesion energy within individual mth hexane liquid layers does not vary for the different layers, as shown in Fig. 3(e), owing to its non-polar nature. By following Eq. (4), the drastic difference between USL and ULV/2 ( WSL  2 L ) leads to the positive spreading coefficients, which corresponds to the total wetting condition (contact angle = 0o), consistent with experimental observation. For another type of non-polar liquid, the metallic mercury liquid has the strongest cohesive forces as compared to other fluids, such as water and hydrocarbon. Owing to its nonwetting property, mercury has been widely used in characterizing pore structures in rocks, such

ACS Paragon Plus Environment

10

Page 11 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

as the permeability and pore-throat distribution, through mercury injection capillary pressure techniques3637. As discussed in Supporting Information (S. VI), the DFT-predicted contact angles of a mercury droplet on the calcite (10.4) surface is larger than 90o, which indicates that ULV/2 is larger than USL (Fig. 3(f)), owing to strong cohesion between mercury layers. The contact angle predictions from our DFT simulation are around 128o , which is close to the static contact angle measurement of 133~136o34, except for that with NC pseudo-potential, as shown in Fig. 3(f-g). This deviation comes from the difference of van der Waals corrections for DFT-NC, DFT-US, and DFT-PAW simulations. Unlike the DFT-US and DFT-PAW simulations using the DFT-D238 van der Waals correction, the DFT-NC simulations in this study uses TkatchenkoScheffler (TS) van der Waals correction39, which is problematic to deal with the highly delocalized nature of the electrons40of the metallic mercury atom. Therefore, the DFT simulation based on DFT-TS van der Waals corrections is not applicable to the metallic liquids. Beside the finite contact angle prediction, it is also interesting to gain more insight into the total wetting condition from DFT simulation, such as hexane on the calcite surface. Here, we introduce a dimensionless physical quantity S /  LV , which denotes the spreading coefficient S41,42 divided by the interfacial tension  LV . S /  LV can be used as an input parameter to characterize the spreading capability of liquid film on the solid surface, as discussed in ref43,44. The definition of spreading coefficient S of a vapor/liquid/solid tri-phase system is given by Eq. (6), same as in Eq. (2.5-2.6) of ref.45 S  WSL  2 LV

(6)

where WSL is the work of adhesion at the solid/liquid interface, and  LV is the surface tension of the liquid. From DFT predicted quantities of USL and ULV/2, we are able to obtain S /  LV , as shown in Eq. (7).

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

S /  LV 

WSL

 LV

2 

U SL 2 U LV / 2

Page 12 of 23

(7)

The DFT prediction S /  LV as summarized in Table 1 can be used to estimate the minimal liquid film thickness e at equilibrium, as proposed by Joanny and De Gennes43 ea

3 LV 2S

(8)

where a is molecular length. The predicted thickness shows that as the hexane liquid spreads over the calcite surface, the final truncated thickness e of liquid film on the edge is ~ 1.4a, slightly more than the hexane monolayer. This prediction result agrees with our adhesion energy studies at different hexane layers, where only the first liquid layer is strongly binding to the solid surface. The calculated U SL and U LV / 2 of the solid surface for different liquids can be extended to predict the wettability of the same solid surface by two different liquids. A calcite-water-hexane system is shown in Fig. S4 (Supporting Information), and the contact angle of a water droplet in hexane fluid on the calcite surface can be determined by the Young equation  SO   SW   OW Cos( )

(9)

where  SO ,  SW , and  OW are the interfacial tensions at calcite/hexane, calcite/water and hexane/water interfaces, respectively. By introducing the Dupré equations into  SO and  SW , a governing equation is given to predict the contact angle of a water droplet on the calcite surface in hexane ( OV   WV )  (WOS  WWS )   OW Cos( )

(10)

where  OV and  WV are the surface tensions of oil and water liquids. The WOS and WWS are the work of adhesion at the hexane/calcite and water/calcite interfaces, respectively. To calculate the contact angle θ, the surface tensions of liquids and interfacial tensions between these liquids are

ACS Paragon Plus Environment

12

Page 13 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

required to be as input parameters for predicting surface wettability. These parameters can be predicted by DFT simulations, as pointed by ref.

4647

, However, we took the measured values

from the literature4849 for convenience.. Three physical quantities:  OV ,  WV , and  OW are 18.4 , 72.8, and 47.1 mN/m, respectively, at room temperature are taken from the ref.4849. WOS and WWS are estimated from our DFT simulation with three different pseudo-potentials. The predicted contact angle of a water droplet on calcite in hexane liquid at room temperature is around 60.1~64.8o, as summarized in Table 2. Comparing with our experimental contact angle of 63.4±3.45°, (Fig. S4(b) in Supporting Information), the extension of the proposed approach to the calcite-water-hexane system gives reasonable agreement. Moreover, our calcite-water-hexane contact angle measurements are consistent with the previous experimental results of 65o ~80o in the ref50. Table 2. A comparison between the predicted and experimentally measured calcite-water-hexane contact angle

*

DFT

 ov   wv (mN/m)

 ow (mN/m)

DFT θ (o)

US

-54.4

-74.4

47.1

64.8

NC

-54.4

-77.6

47.1

60.1

PAW

-54.4

-77.8

47.1

60.1

Wso  Wsw

(mN/m)

Exp. θ (o) 63.4±3.45* 65~8050

Experimental measurement of this work.

In addition to the calcite surface, we can apply this DFT approach to study the wettability of other types of solid crystal surfaces, such as metallic Ag surface, which further confirm the generality of our methodology, as discussed in Supporting Information. Moreover, the proposed method is not only able to predict intrinsic surface wettability but it is also capable of quantifying wettability alteration of solid surfaces under different ion substitutions, such as alkaline earth metal ions, i.e. Be2+, Mg2+, Sr2+, and Ba2+, which is valuable for enhanced oil recovery and

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

geological carbon sequestration applications. In the DFT slab model, we replace one of the four calcium ions in the calcite surface layer with the other metal ion, i.e., 25% ion substitution on the calcite surface. The simulation results show only adhesion energy USL at the calcite/water interface is altered under different ion substitutions while the cohesion energy ULV/2 between water layers remains the same, as shown in Fig. 4(a). According to Eq. (4), in opposition to Sr2+ and Ba2+ ion substitutions, only Mg2+ and Be2+ ion substitutions result in more water-wet phenomena owing to stronger adhesion strength at the calcite/water interface, which agrees with the measurements from previous lab experiments35. Moreover, the charge density difference distributions indicate how the charge density changes upon adsorption of the water molecules to the ion-substituted surfaces. These distributions further reveal the adhesion energy changes at the calcite/water interface are related to the strength of interaction between the substituted cation and water dipole moment, as shown in Fig. 4(c-f). The influence of Mg2+ ion substation on calcite surface wettability is further studied for higher ion substitution rates on the calcite surface. The simulation results show that the adhesion energy at the first and second layers of water is increased as the number of Mg2+ substituted ions increases. Interestingly, as shown in Figure 4(b), the adhesion strength at the first layer reaches the maximum value when the ion substitution rate is 75%. As for that of the second layer, the adhesion strength continuously increases with the substitution rate increases. However, the adhesion energy of higher water layers (ULV/2) is not affected by different ion substitutions and different concentrations of ion substitutions at the liquid/solid interface. The quantitative prediction of wettability changes for the calcite surface are summarized in Table 3 and 4 under different ion substitutions and different substitution rates, respectively. These studies on wettability alternation further demonstrate the correctness of our methodology in Fig. 1(a), where

ACS Paragon Plus Environment

14

Page 15 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

the adhesion energy of higher liquid layers on the liquid covering calcite surface is essentially the half cohesion energy between liquid layers U LV / 2 , as indicated in Fig. 4(a-b).

Figure 4. (a)Adhesion energy at respective water layers on calcite (10.4) surface under different ion substitutions (b) and under Mg2+ ion substitutions with different substitution rates. (c-f) Charge density difference of the water on calcite surface substituted with Be2+, Mg2+, Sr2+, and Ba2+ ions, respectively. Table 3. Predicted water contact angles on the calcite (10.4) surface under different ion substitutions (25% ion substitution rate) Be2+

Mg2+

Ca2+

Sr2+

Ba2+

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23

ULV/2 [kJ/mol]

57.57

57.62

57.71

57.90

58.02

USL [kJ/mol]

107.05

103.32

99.84

94.97

91.87

Contact Angle (o)

30.8

37.5

43.1

50.2

54.3

Table 4. Predicted water contact angles on the calcite (10.4) surface under Mg2+ ion substitutions with different substitution rates 100%-Ca2+

25%-Mg2+

50%-Mg2+

75%-Mg2+

100%-Mg2+

ULV/2 [kJ/mol]

57.71

57.62

57.47

57.36

57.42

USL [kJ/mol]

99.84

103.32

107.86

108.32

108.34

Contact Angle (o)

43.1

37.5

28.7

27.3

27.5

In summary, we have developed an approach to predict the intrinsic surface wettability of a solid by a liquid based on a simple DFT slab model. As compared to MD simulations or Van Oss-Chaudhury-Good theory, the proposed first-principle prediction method is much more capable, as it requires less input parameters and can even deal with high-energy oxide surfaces under the total wetting condition. The predicted physical quantities can be potentially used in large-scale computational fluid simulation for practical oil and gas reservoir applications. Although we consider calcite surface as a benchmark in this study, the proposed method is very general and applicable to any crystal surfaces and surrounding fluids. It offers a powerful tool for predicting intrinsic wettability and wettability alteration characteristics of multi-phase systems in nature and various energy-water applications. ACKNOWLEDGMENT

ACS Paragon Plus Environment

16

Page 17 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

This work was funded by Abu Dhabi National Oil Company and mainly undertaken on the High Performance Computing facilities at Masdar Institute. The authors would like to acknowledge the technical support of Masdar Institute Research Computing team (S. Benna, D. El Moghraby, S.J. Illikkal, S. Sanchez) and generous GPU cloud computing support from Alibaba Cloud. The authors thank Prof. Kean Wang at Petroleum Institute for helping with XRD characterization. Supporting Information Available: computational details; calcite (10.4) surface cleaved from bulk calcite crystal; configurations of eight-layered water on the calcite (10.4) surface; schematic of a water droplet on calcite surface in hexane; experimental measurements for contact angles of water on calcite surface in air and in hexane; non-polar mercury liquid on the calcite (10.4) surface; surface wettability of a metal surface for water liquid; list of all used pseudo-potentials in QuantumEspresso package. REFERENCES (1)

Wang, Y.; Xu, H.; Yu, W.; Bai, B.; Song, X.; Zhang, J.Surfactant Induced Reservoir Wettability Alteration: Recent Theoretical and Experimental Advances in Enhanced Oil Recovery. Pet. Sci. 2011, 8 , 463–476.

(2)

Hirasaki, G.; Zhang, D. L.Surface Chemistry of Oil Recovery from Fractured, Oil-Wet, Carbonate Formations. Spe J. 2004, 9 , 151–162.

(3)

Alvarado, V.; Manrique, E.Enhanced Oil Recovery: An Update Review. Energies 2010, 3 , 1529–1575.

(4)

Liu, H.; Raza, A.; Aili, A.; Lu, J.; AlGhaferi, A.; Zhang, T.Sunlight-Sensitive AntiFouling Nanostructured TiO2 Coated Cu Meshes for Ultrafast Oily Water Treatment. Sci. Rep. 2016, 6.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(5)

Page 18 of 23

Pastero, L.; Costa, E.; Bruno, M.; Rubbo, M.; Sgualdino, G.; Aquilano, D.Morphology of Calcite (CaCO3) Crystals Growing from Aqueous Solutions in the Presence of Li+ Ions. Surface Behavior of the {0001} Form. Cryst. Growth Des. 2004, 4 , 485–490.

(6)

Aquilano, D.; Otálora, F.; Pastero, L.; García-Ruiz, J. M.Three Study Cases of Growth Morphology in Minerals: Halite, Calcite and Gypsum. Prog. Cryst. Growth Charact. Mater. 2016, 62 , 227–251.

(7)

Li, G.; Alhosani, M. H.; Yuan, S.; Liu, H.; Ghaferi, A.Al; Zhang, T.Microscopic Droplet Formation and Energy Transport Analysis of Condensation on Scalable Superhydrophobic Nanostructured Copper Oxide Surfaces. Langmuir 2014, 30 , 14498–14511.

(8)

Aili, A.; Li, H.; Alhosani, M. H.; Zhang, T.Unidirectional Fast Growth and Forced Jumping of Stretched Droplets on Nanostructured Microporous Surfaces. ACS Appl. Mater. Interfaces 2016, 8 , 21776–21786.

(9)

Cho, H. J.; Preston, D. J.; Zhu, Y.; Wang, E. N.Nanoengineered Materials for Liquid– vapour Phase-Change Heat Transfer. Nat. Rev. Mater. 2016, 2, 16092.

(10) Zhang, M.; Feng, S.; Wang, L.; Zheng, Y.Lotus Effect in Wetting and Self-Cleaning. Biotribology 2016, 5, 31–43. (11) Xu, Y.; Dibble, C. J.; Petrik, N. G.; Smith, R. S.; Kay, B. D.; Kimmel, G. A.Complete Wetting of Pt (111) by Nanoscale Liquid Water Films. J. Phys. Chem. Lett. 2016, 7 , 541– 547. (12) Deng, Y.; Chen, L.; Liu, Q.; Yu, J.; Wang, H.Nanoscale View of Dewetting and Coating on Partially Wetted Solids. J. Phys. Chem. Lett. 2016, 7 , 1763–1768. (13) VanOss, C. J.Interfacial Forces in Aqueous Media; CRC press, 2006.

ACS Paragon Plus Environment

18

Page 19 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(14) DellaVolpe, C.; Siboni, S.; Maniglio, D.; Morra, M.; Cassinelli, C.; Anderle, M.; Speranza, G.; Canteri, R.; Pederzolli, C.; Gottardi, G.Recent Theoretical and Experimental Advancements in the Application of the van Oss–Chaudhury–Good Acid–base Theory to the Analysis of Polymer Surfaces II. Some Peculiar Cases. J. Adhes. Sci. Technol. 2003, 17 , 1425–1456. (15) DellaVolpe, C.; Maniglio, D.; Siboni, S.; Morra, M.Recent Theoretical and Experimental Advancements in the Application of van Oss–Chaudury–Good Acid–base Theory to the Analysis of Polymer Surfaces I. General Aspects. J. Adhes. Sci. Technol. 2003, 17 , 1477– 1505. (16) Saville, G.Computer Simulation of the Liquid–solid–vapour Contact Angle. J. Chem. Soc. Faraday Trans. 2 Mol. Chem. Phys. 1977, 73 , 1122–1132. (17) Matsumoto, S. M. T. K. S.; Kimura, Y. Y. T.Liquid Droplet in Contact with a Solid Surface. Microscale Thermophys. Eng. 1998, 2 , 49–62. (18) Pandey, P. R.; Roy, S.Is It Possible to Change Wettability of Hydrophilic Surface by Changing Its Roughness? J. Phys. Chem. Lett. 2013, 4 , 3692–3697. (19) Raj, R.; Maroo, S. C.; Wang, E. N.Wettability of Graphene. Nano Lett. 2013, 13 , 1509– 1515. (20) Adams, P.; Henderson, J. R.Molecular Dynamics Simulations of Wetting and Drying in LJ Models of Solid-Fluid Interfaces in the Presence of Liquid-Vapour Coexistence. Mol. Phys. 1991, 73 , 1383–1399. (21) Wei, N.; Lv, C.; Xu, Z.Wetting of Graphene Oxide: A Molecular Dynamics Study. Langmuir 2014, 30 , 3572–3578.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

(22) Harrison, N. M.An Introduction to Density Functional Theory. NATO Sci. Ser. SUB Ser. III Comput. Syst. Sci. 2003, 187, 45–70. (23) Sakuma, H.; Andersson, M. P.; Bechgaard, K.; Stipp, S. L. S.Surface Tension Alteration on Calcite, Induced by Ion Substitution. J. Phys. Chem. C 2014, 118 , 3078–3087. (24) Andersson, M. P.; Sakuma, H.; Stipp, S. L. S.Strontium, Nickel, Cadmium, and Lead Substitution into Calcite, Studied by Density Functional Theory. Langmuir 2014, 30 , 6129–6133. (25) Kawano, J.; Sakuma, H.; Nagai, T.Incorporation of Mg2+ in Surface Ca2+ Sites of Aragonite: An Ab Initio Study. Prog. earth Planet. Sci. 2015, 2 , 7. (26) Andersson, M. P.; Dideriksen, K.; Sakuma, H.; Stipp, S. L. S.Modelling How Incorporation of Divalent Cations Affects Calcite Wettability–implications for Biomineralisation and Oil Recovery. Sci. Rep. 2016, 6. (27) Andersson, M. P.; Stipp, S. L. S.Sensitivity Analysis of Cluster Models for Calculating Adsorption Energies for Organic Molecules on Mineral Surfaces. J. Phys. Chem. C 2011, 115 (20), 10044–10055. (28) Decker, E. L.; Frank, B.; Suo, Y.; Garoff, S.Physics of Contact Angle Measurement. Colloids Surfaces A Physicochem. Eng. Asp. 1999, 156 , 177–189. (29) Schrader, M. E.Young-Dupre Revisited. Langmuir 1995, 11 , 3585–3589. (30) Sheppard, D.; Terrell, R.; Henkelman, G.Optimization Methods for Finding Minimum Energy Paths. J. Chem. Phys. 2008, 128 , 134106. (31) See Http://quantumwise.com/ for QuantumWise Atomistix ToolKit (ATK) with Virtual NanoLab.

ACS Paragon Plus Environment

20

Page 21 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(32) Morrow, N. R.Wettability and Its Effect on Oil Recovery. J. Pet. Technol. 1990, 42 , 1– 476. (33) Johnson, C. A.; Furrer, G.Influence of Biodegradation Processes on the Duration of CaCO3 as a pH Buffer in Municipal Solid Waste Incinerator Bottom Ash. Environ. Sci. Technol. 2002, 36 , 215–220. (34) Ethington, E. F.Interfacial Contact Angle Measurements of Water, Mercury, and 20 Organic Liquids on Quartz, Calcite, Biotite, and Ca-Montmorillonite Substrates; US Geological Survey, 1990. (35) Zhang, P.; Tweheyo, M. T.; Austad, T.Wettability Alteration and Improved Oil Recovery by Spontaneous Imbibition of Seawater into Chalk: Impact of the Potential Determining Ions Ca 2+, Mg 2+, and SO 4 2−. Colloids Surfaces A Physicochem. Eng. Asp. 2007, 301 , 199–208. (36) Pittman, E. D.Relationship of Porosity and Permeability to Various Parameters Derived from Mercury Injection-Capillary Pressure Curves for Sandstone (1). Am. Assoc. Pet. Geol. Bull. 1992, 76 , 191–198. (37) Schmitt, M.; Fernandes, C. P.; daCunha Neto, J. A. B.; Wolf, F. G.; dosSantos, V. S. S.Characterization of Pore Systems in Seal Rocks Using Nitrogen Gas Adsorption Combined with Mercury Injection Capillary Pressure Techniques. Mar. Pet. Geol. 2013, 39 , 138–149. (38) Grimme,

S.Semiempirical

GGA type

Density Functional Constructed with a

Long range Dispersion Correction. J. Comput. Chem. 2006, 27 , 1787–1799.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 23

(39) Tkatchenko, A.; Scheffler, M.Accurate Molecular van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. Phys. Rev. Lett. 2009, 102 , 73005. (40) Bučko, T.; Lebègue, S.; Hafner, J.; Angyan, J. G.Tkatchenko-Scheffler van Der Waals Correction Method with and without Self-Consistent Screening Applied to Solids. Phys. Rev. B 2013, 87 , 64110. (41) Bonn, D.; Eggers, J.; Indekeu, J.; Meunier, J.; Rolley, E.Wetting and Spreading. Rev. Mod. Phys. 2009, 81 , 739. (42) Cazabat, A.-M.How Does a Droplet Spread? Contemp. Phys. 1987, 28 , 347–364. (43) Joanny, J. F.; DeGennes, P.-G.Upward Creep of a Wetting Fluid: A Scaling Analysis. J. Phys. 1986, 47 , 121–127. (44) Joanny, J.-F.; Andelman, D.Steady-State Motion of a Liquid/liquid/solid Contact Line. J. Colloid Interface Sci. 1987, 119 , 451–458. (45) DeGennes, P.-G.Wetting: Statics and Dynamics. Rev. Mod. Phys. 1985, 57 , 827. (46) Gloor, G. J.; Jackson, G.; Blas, F. J.; DelRio, E. M.; DeMiguel, E.Prediction of the Vapor− Liquid Interfacial Tension of Nonassociating and Associating Fluids with the SAFT-VR Density Functional Theory. J. Phys. Chem. C 2007, 111 , 15513–15522. (47) Andersson, M. P.; Bennetzen, M.V; Klamt, A.; Stipp, S. L. S.First-Principles Prediction of Liquid / Liquid Interfacial Tension. J. Chem. Theory Comput 2014, 10, 3401−3408. (48) Rulison, C.Two-Component Surface Energy Characterization as a Predictor of Wettability and Dispersability. Krus. Appl. note AN213 2000, 1–22. (49) Zeppieri, S.; Rodríguez, J.; López de Ramos, A. L.Interfacial Tension of Alkane+ Water Systems. J. Chem. Eng. Data 2001, 46 , 1086–1088.

ACS Paragon Plus Environment

22

Page 23 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(50) Okayama, T.; Keller, D. S.; Luner, P.The Wetting of Calcite Surfaces. J. Adhes. 1997, 63 , 231–252.

ACS Paragon Plus Environment

23