Capillary Condensation in 8 nm Deep Channels - The Journal of

Jan 11, 2018 - ... Development Grant, as well as ongoing research funding through the Discovery Grants program and the Canada Research Chairs program...
0 downloads 12 Views 7MB Size
Subscriber access provided by READING UNIV

Letter

Capillary Condensation in 8-nm Deep Channels Junjie Zhong, Jason Riordon, Seyed Hadi Zandavi, Yi Xu, Aaron Harrinarine Persad, Farshid Mostowfi, and David Sinton J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b03003 • Publication Date (Web): 11 Jan 2018 Downloaded from http://pubs.acs.org on January 12, 2018

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 20 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

Capillary Condensation in 8-nm Deep Channels Junjie Zhong,† Jason Riordon,† Seyed Hadi Zandavi,‡ Yi Xu,† Aaron H. Persad,†, ‡ Farshid Mostowfi,§ and David Sinton†,* †

Department of Mechanical and Industrial Engineering, University of Toronto, Toronto, Ontario

M5S 3G8, Canada ‡

Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge,

Massachusetts 02139, United States §

Schlumberger-Doll Research, Cambridge, Massachusetts 02139, United States

Corresponding Author *

E-mail: [email protected]

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 20

ABSTRACT

Condensation at the nanoscale is essential to understanding many natural and synthetic systems relevant to water, air and energy. Despite its importance, the underlying physics of condensation initiation and propagation remain largely unknown at sub-10 nm, mainly due to the challenges of controlling and probing such small systems. Here we study the condensation of n-propane down to 8-nm confinement in a nanofluidic system, distinct from previous studies at ~100 nm. The condensation initiates significantly earlier in the 8-nm channels and it initiates from the entrance, in contrast to channels just 10-times larger. The condensate propagation is observed to be governed by two liquid-vapor interfaces with an interplay between film and bridging effects. We model the experimental results using classical theories and find good agreement, demonstrating that this 8-nm non-polar fluid system can be treated as a continuum from a thermodynamic perspective, despite having only 10-20 molecular layers.

TOC GRAPHICS

Liquid

Film

Vapor

8 nm

ACS Paragon Plus Environment

2

Page 3 of 20 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

Condensation in sub-10 nm systems has garnered increased interest recently due the important role of phase change in climate engineering,1-4 water purification,5-6 energy transport and storage,6-7 shale gas production,8-9 surface treatments,10-12 heat transfer13-14 and nanomaterials1517

. The underlying physics of the condensation process at the sub-10 nm scale, however, remains

poorly understood. For example, condensation initiation within sub-100 nm confinement can be described through the Kelvin equation, where the highly curved liquid-vapor interface introduces a shift in saturation conditions.18 However, there is much debate as to the validity of Kelvin equation below 10 nm since such small systems span a small, discrete number of molecules. While certain studies claim to have validated the Kelvin equation down to sub-10 nm systems,1920

other works indicate otherwise pointing to, for example, the breakdown of the ideal gas and

incompressible liquid assumptions inherent in Kelvin equation.21-22 In addition, the condensate propagation kinetics within a confined sub-10 nm space is difficult to directly observe. Numerical simulations of nanoscale systems23 indicate that condensation propagates by the liquid bridging of films that grow on opposite channel walls, while experiments in nanoporous media suggest that the condensate propagates by capillary flow24 and vapor diffusion.25 The validity of the theoretical approaches at the sub-10 nm scale can only be examined after overcoming the challenges of obtaining accurate experimental data, specifically the difficulties of precisely controlling and probing sub-10 nm systems. Currently, the most established experimental methods for studying nanoscale condensation require measuring adsorption and desorption isotherms in nanoporous media.26 However, these techniques are often affected by nonuniformity in pore size and geometry, and generally do not allow direct observation of condensation dynamics in pores. Other approaches such as electron scanning microscopy,27 atomic force microscopy28 and electrical conductivity29 have been applied to study nanoscale

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 20

condensation of water,30 electrolytes31 and hydrocarbons26 in cylindrical pores,32-33 nanoslits32, 34 or carbon nanotubes.35 However, accurately controlling, changing and probing the thermodynamic conditions within large testing ranges in such systems remains a challenge. Recent advances in nanofluidics have led to excellent control over system geometry,36-37 thermodynamic conditions38-39 and optical accessibility40, providing an opportunity to directly study condensation at the nanoscale. Studies using nanofluidic devices have examined the condensation of water18 and hydrocarbons41 at ~100 nm scales. However, condensation within precisely manufactured, controlled and directly visualized sub-10 nm systems – a scale relevant to a broad range of natural and synthetic systems– has not yet been achieved. In this article, we visualize the condensation of n-propane in precisely fabricated 8-nm silicon nanochannels. A silicon nitride layer is deposited at the base of the nanochannels to form FabryPerot optical resonance cavities,42 enabling the imaging of discrete fluid phases at this scale. We simultaneously observe the phase behavior in microchannels (i.e the bulk phase) and find that the initiation of condensation in the 8-nm channels agrees closely with the Kelvin equation. Condensate propagation in the 8-nm channels is found to be governed by two liquid-vapor interfaces: one that is pinned at the inlet due to capillary condensation and another that grows through capillary filling. The condensation initiation location, film and bridging effects are discussed in detail and compared with previous findings in sub-100 nm systems. In providing a clear picture of condensation phenomena within a well-defined nanosystem, this work provides fundamental insight into the physics behind condensation down to 8-nm confinement, which we find to deviate significantly from the physics at larger scales. The validity of the Kelvin equation down to 8-nm based on close agreement between classical theoretical modeling and the

ACS Paragon Plus Environment

4

Page 5 of 20

experimental results supports the continuous assumption for a fluid system with only 10-20 layers of molecules.

b

a

A

A 10

A A

Height (nm)

Liquid in Nano

8

Vapor in Nano

Liquid in Bulk

Vapor in Bulk

c

6 4 2

Reservoir

0

0

2 4 Position (µ µm)

6

VBulk

(i) 340

VNano VBulk (ii)

LNano VBulk

VNano LBulk

Temperature (K)

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

320

300 Condensation in Nanochannels Condensation in Bulk Kelvin Equation (R = 4 nm) Saturation Condition (Bulk)43

(iii)

LNano

VNano 280 0.5

1.0

1.5 Pressure (MPa)

2.0

2.5

Figure 1. Thermodynamic conditions for condensation in 8-nm channels and in the bulk scale. (a) Schematic of the nanofluidic device with both nanochannels and bulk microchannels under the same thermodynamic conditions for comparison. (b) AFM profile of a single nanochannel (8.0 ± 0.2 nm). (c) Measured condensation conditions at 8 nm and in the bulk, along with corresponding saturation curves obtained using the Kelvin equation and NIST values.43 Each point for the 8-nm case represents the average of a total of 33 measurements (3 runs with 11 nanochannels). Each point for the bulk case represents the average of 6 measurements (3 runs

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 20

with 2 microchannels). Error bars are smaller than the symbols. Images show system behavior at a pressure of 1.31 MPa at one of three temperatures. When the temperature is above the Kelvin equation prediction, both micro- and nanochannels remain in the vapor phase (light red and light yellow, respectively), (i). When the temperature approaches the Kelvin equation prediction, capillary condensation occurs, with the condensate growing into the 8-nm nanochannels (the dark red phase, Lnano denoting the liquid, replacing the light red phase, Vnano denoting the vapor, while microchannels remain in the vapor phase (light yellow, Vbulk), (ii). When temperature reaches the bulk prediction, condensation occurs and the condensate grows in both the nanochannels and microchannels (dark red and dark yellow, Lnano and Lbulk, respectively), (iii). The scale bar corresponds to a distance of 100 µm and applies to all microscope images. As shown in Figure 1a, parallel dead-end nanochannels (500 µm × 4 µm × 8 nm) were etched (Figure 1b), and flanked by two deep microchannels (500 µm × 20 µm × 20 µm) to allow direct comparison of condensation behavior between sub-10 nm confinement and bulk under the same thermodynamic conditions. All channels connect to the same vapor reservoir and were sealed above via a glass layer, with a 200-nm thick silicon nitride base layer to enhance optical contrast between liquid and vapor phases (see Supporting Information Section 1). During all experiments, the pressure was controlled by a syringe pump and monitored using a pressure transducer. Temperature was controlled using a water bath via a copper chiller fixed at the bottom of the nanofluidic chip and monitored using thermocouples inserted within the device. Before the experiment, the devices were vacuumed for 6 hrs (see Supporting Information Section 2). We measured the highest condensation temperature for different pressures, in both the bulk (i.e., microchannels) and 8-nm nanochannels. In each test, propane was first set to the desired pressure, and then pumped into the nanochannel after the target temperature was established on

ACS Paragon Plus Environment

6

Page 7 of 20 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 chip. For each temperature, the system was monitored for up to 1 hr to determine whether condensation initiates, and was then was fully vacuumed for 30 min between measurements (see details in Supporting Information Section 2). Results are shown in Figure 1c, with error bars representing both the experimental uncertainty and the resolution between runs for a given pressure (a zoom-in of the T-P data is provided in Supporting Information Section 3). The condensation initiation temperatures in the 8-nm channels are plotted with those predicted by the Kelvin equation:  −   −

 



ln

  



=

(1)

where  is the equilibrium pressure for capillary condensation with a curved meniscus (i.e.,

the Kelvin pressure),   is the saturation pressure in bulk,  is the universal gas constant,  is

temperature of the vapor phase,  is the liquid-vapor interface tension and  is the main radius of curvature of the liquid-vapor interface. We also assumed that propane fully wets the nanochannel surfaces during condensation42.

We found that within a wide pressure and temperature range (~0.6-2.3 MPa, 286.15-339.15 K), predictions from the Kelvin equation based on bulk fluid properties matched the experimental initiation of condensation results well, where the measured temperature shift between nanoscale and bulk condensation (gap between red and blue points, Figure 1c, ~1-3 K) was at most 8 % different from that predicted by classical theory (see Supporting Information Section 4). It is important to note that the disjoining pressure might influence the equilibrium contact angle,44 and further affect the Kelvin pressure. However, we estimated the disjoining pressure in our system and found to be less than 0.01 MPa. The capillary pressure can thus change up to 0.4% and is negligible in the context of this experiment (see Supporting Information Section 5).

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 20

This result directly supports the applicability of the Kelvin equation under extreme nanoconfinement, which is an important indication of the validity of classical thermodynamics down to the sub-10 nm scale – a matter of debate in the field. For example, experiments using an interferometer with an ~8-nm gap suggest the applicability of the Kelvin equation for modeling the capillary condensation of cyclohexane.20, 45 In addition, recent experimental results for both water and hydrocarbons adsorption in porous media support the Kelvin equation at this scale.19 In contrast, another adsorption study in porous material46 indicates that the Kelvin equation cannot predict capillary condensation within the mesopore size range because of the instability of the adsorbed liquid film. In addition, the assumptions of an ideal gas phase and an incompressible liquid phase inherent to the Kelvin equation have been questioned, and more complicated models based on equations-of-state have been introduced to explain capillary condensation of both non-polar and polar molecules in sub-10 nm nanoporous media.21-22 In contrast to these works, our method provides direct visual distinction between fluid phases within a well-defined geometry and strong evidence of the applicability of the Kelvin equation on non-polar molecules down to 8-nm.

ACS Paragon Plus Environment

8

Page 9 of 20

a

d

300

282.15 K 282.55 K 282.95 K 283.35 K 283.75 K 284.15 K

200 100

300 200 100 0

0

b

Model Prediction Experimental Results (282.15 K)

500 400

400

Length (µm)

Condensate Length (µ m)

500

0 1

2 3 Time (s)

4

0

5

Position (m) 250

0

500

e

0.1

0.2

0.3 Time (s)

0.4

0.5

0.6

Model Prediction Experimental Results (284.15 K)

500

0 400 Length (µ m)

0.2 Time (s)

0.4 0.6 Top view

300 200 100 0 0

1

2

8 nm

c

Liquid

θ

PIL

Film

PDL

RD

VL RI

PIV Condensation Interface

Vapor

f 8 nm

Monolayer adsorption

PDV Filling Interface

500 µm

T

3 Time (s)

(°)

Cross-sectional view

Contact Angle of Condensation Interface,

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

4

0 µm 1 µm 2 µm 3 µm 4 µm

90 75

5

5 µm 10 µm 15 µm 20 µm



60 45 30 15 0 281.5

282.0

Bulk Saturation Temperature

282.5

283.0

283.5

Temperature (K)

284.0

284.5 Nano Saturation Temperature

Figure 2. Measurement of condensate growth as a function of superheat at isobaric conditions. (a) Total condensate length versus time for different temperatures, at a pressure of 0.61 MPa. Each point represents the average of a total of 33 measurements (3 runs with 11 nanochannels). For clarity, representative error bars are shown. (b) Images depicting condensate growth over time within a group of 5 nanochannels at 282.15 K. Images were stretched horizontally by a factor of 2.5 and contrast-enhanced in black and white for clarity. Red lines separate each frame (see Supporting Information Movie). (c) Schematic of the condensate growth model. The enlarged image of the interface was stretched vertically by a factor of 4 with respect to the

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 20

original image to highlight the presence of an interfacial film. Liquid (black), film (gray) and vapor (light gray) are distinguished. (d and e) Model predictions compared to experimental results at 282.15 K and 284.15 K, respectively. (f) The change of contact angle at the condensation interface for different temperatures, as predicted using the model. The ability to differentiate phases within the sub-10 nm channels uniquely enables the observation of condensation dynamics. Condensate propagation was studied under isobaric conditions (0.61 MPa), with temperature varying between 282.15 K and 284.15 K at 0.4 K intervals. Condensation initiated at the inlet of the nanochannel and grew into the nanochannel with the interface pinned at the inlet. The total condensate length changing over time was measured, and results for different conditions are shown in Figure 2a. Figure 2b shows the typical filling behavior for a group of 5 nanochannels at 282.15 K. The initiation of condensation within the 8-nm nanochannels is caused by the growth of two adsorption layers on opposite channel walls, until they eventually meet and coalesce (Figure 3); when the pressure is above the Kelvin pressure, vapor molecules enter the nanochannel in a locally-supersaturated condition, and are adsorbed quickly near the inlet to form the liquid phase. This phenomenon is in marked contrast to condensation in channels only ~10 times larger where the condensation has been observed to initiate at the dead end.41 A key difference in the initiation of condensation between sub-10 nm and 100 nm (and larger) channels is due to the role of vapor transport resistance, which scales inversely with the square of the channel dimension.47 Vapor adsorption in larger channels (i.e., sub-100 nm) is thus less limited by the molecular transport and the adsorption rate becomes similar along the entire channel, with the dead end providing the internal corner structure preferential for adsorbed layers to collect. It is worth noting that previous simulations indicate that condensation at the sub-10 nm scale may also initiate from the

ACS Paragon Plus Environment

10

Page 11 of 20 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

dead end when the channel length is close to the channel depth.23 The result observed here is thus limited to the channel geometry where the length is orders of magnitude larger than the depth dimension. In addition to clearly identifying the vapor-liquid interface location, the imaging method used here allowed for sub-8 nm resolution at the interface. Specifically, a sharp high-contrast transition at the leading edge of the liquid film to the vapor phase was observed (Figure 2b and c). This sharp onset indicates a non-zero, finite advancing contact angle during the condensate growth, with the radius of curvature of the interface on the order of ~1 nm (the inferred filling interface shape is illustrated in Figure 2c). Similar observations have been made in other systems featuring long, highly confined geometries, most notably in water condensation in carbon nanotubes with diameters of

~200 nm.48 This stepwise transition of the light intensity at the

interface lies in contrast to the more gradual signal decay expected for a near-zero film contact angle assumed in previous research, notably for capillary filling tests in nanochannels.49-50 As is shown in Figure 2c, once the initial condensate forms at the inlet, there are two effects that combine to drive condensate growth in the nanochannel: (i) the accumulation of vapor molecules at the condensation interface as a result of capillary condensation, and (ii) the liquid flow as a result of the capillary pressure at the filling interface. We modeled (i) and (ii) with the HertzKnudsen equation and Poiseuille law, respectively, and linked these with the conservation of mass (details in Supporting Information Section 6). The resulting model expresses condensate length () as a function of time (): 



  !"#$% & ) '(

=

!

-

/0

, *+ . 

− 1

23

 ( 4(

5

(2)

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

Page 12 of 20

where 0 , the main radius of curvature of the condensation interface, can be expressed as: 7





 7:

= 89 − 

' 

+

2'3

 8  < =&2>?( 4 

(3)

where ℎ is the depth of the nanochannel, A is the condensation coefficient (assumed to be unity for non-polar molecules51), B is the liquid phase density, C is the molar mass of the fluid,  is

the temperature for the system, 0 is the vapor pressure at the condensation liquid-vapor

interface, D is the viscosity of the liquid phase, and E is the thickness of the liquid film in front of

the filling interface (Figure 2b and c). The film thickness was measured through imaging and assumed to be uniform (see details in Supporting Information Section 6). Solving this system of equations (eq 2 and 3) numerically for the highest and lowest temperature cases provides the model predictions in Figure 2d and e (results for other conditions are shown in Supporting Information Section 6). In all cases, the model predicts the total condensate growth well. For example, the predicted time for full condensate filling of the nanochannel was within 6% of the experimental time (See Supporting Information Section 6). The change of the contact angle at the condensation interface (F, Figure 2c) is essential in condensate propagation kinetics, as it affects both capillary condensation and capillary filling and can be determined directly from the main radius of curvature of the condensation interface (0 ) (Supporting Information Section 6). Within the first ~20 µm of condensate length, the contact angle at the condensation interface (F) increases rapidly from zero to ~97 % of the equilibrium value (F VL. An initial liquid bridge forms and a nanobubble (with the main radius of curvature of a few nanometer) is trapped, (ii). Secondary liquid bridge forms as a result of adsorption at a pressure above the Kelvin pressure, (iii). Nanobubble collapses as a result of liquid bridging at multiple locations within the nanochannel, (iv). For low subcooling conditions, liquid bridging effects were observed in addition to continuous condensate growth (Figure 3a).

Liquid bridging during condensation has been previously

reported at the ~100 nm scale, where the liquid bridge formed at relatively high subcooling conditions as a result of capillary pressure-driven flow.41 However, within an 8-nm deep channel, multi-liquid-bridging was observed at low subcooling. This result cannot be explained through similar condensate growth theory developed for ~100 nm confinement or larger systems.41 Figure 3b shows a series of schematics illustrating the mechanisms behind the liquid bridging phenomenon was observed here. At any given condition above the Kelvin pressure, two liquid formation mechanisms act within the nanochannel: (i) the propagation of the liquid column as a result of both capillary condensation and capillary filling (VL, shown in Figure 3b (i)), and (ii) the molecular adsorption of vapor onto the solid channel surfaces (Va, shown in Figure 3b (i)). When subcooling is high (e.g., 2.3 K at 282.15 K/0.61 MPa), continuous condensate growth is observed because VL is larger than Va (Figure 2b). When subcooling is low (e.g., 0.3 K at 284.15 K/0.61 MPa), VL is reduced significantly (here by a factor of ~10 based on both experimental observations (Figure 2a) and model predictions), while Va will remain largely unchanged since it

ACS Paragon Plus Environment

14

Page 15 of 20 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

is only proportional to the saturation ratio, here from 1.032 to 1.008.52 In addition, VL decreases as the condensate length increases (Figure 2a). Also, surface defect may exist from the nanochannel fabrication process, and adsorption onto such defects may further increase the sensitivity of Va to subcooling (similar to the sensitivity observed during capillary evaporation in larger nanochannel).53 As a result, liquid bridging is most likely to occur at lower subcooling conditions for longer liquid columns, as observed here (Figure 3a). Once the initial liquid bridge forms (Figure 3b (ii)), the vapor pressure in the resulting nanobubble will be above the Kelvin pressure due to the pinned interface,54 resulting in internal adsorption and condensation. Fast adsorption and condensation within the nanobubble can lead to secondary liquid bridging (Figure 3b (iii)), which continues to repeat until the nanobubble collapses (Figure 3b (iv)). Such a multiliquid-bridging effect is also expected from previous unstable liquid film theories, where the sudden coalescence of two wetter surfaces can happen when the channel depth is below a certain threshold (~10 nm).55 Our experimental results shed new light on liquid bridging phenomena at the sub-10 nm scale and provide direct support for the predictions of previous simulations.23, 56-57 It is also worth noting that, for theoretical analysis, sub-100 nm and larger fluidic systems are typically assumed to act as continuous systems.39, 41 However, the validity of such an assumption in sub-10 nm systems is not immediately clear. Here, the theoretical explanations for both the initiation of condensation and condensation dynamics inside 8-nm high channels are based entirely on classical thermodynamics and continuum fluid mechanics. Since the experimental results agree well with classical and continuum theories, the assumption of treating a non-polar fluid system with only 10-20 molecular layers as a continuous system appears valid from a condensation perspective. Furthermore, the condensate initiation, the film formation, and the liquid bridging phenomena were rendered observable in real time by our nanofluidic system and

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 20

indicate fundamental deviations of the condensation phenomena at the sub-10 nm scale compared to the sub-100 nm scale, even though both systems may be treated as continuous. In summary, we directly measured condensation initiation and dynamics within 8-nm deep silicon nanochannels. The initiation of condensation deviates from the bulk condition and was found to agree closely with predictions from the Kelvin equation. The condensation dynamics at the sub-10 nm scale are distinct from those of similar systems that are just 10 times larger (~100 nm). Specifically, adsorption dominates the condensation initiation in sub-10 nm channels, resulting in liquid formation starting from the channel's entrance. A model combining capillary condensation and capillary filling has been developed and was found to give good predictions of the condensation propagation observed experimentally. Additional bridging dynamics were observed at low subcooling conditions – also in contrast to sub-100 nm and larger systems. Collectively, these results characterize condensation at the sub-10 nm scale. With direct, realtime imaging of condensation phenomena within a well-controlled environment, we also provide strong support of the applicability of classical theory in these systems down to the 8-nm scale.

ACS Paragon Plus Environment

16

Page 17 of 20 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

ASSOCIATED CONTENT Supporting Information Nanofluidic chip fabrication procedure/Experimental setup and procedures for measuring the initiation of condensation and condensation propagation/Zoom-in of the testing saturation conditions for condensation at 8 nm and bulk scale/Deviation of condensation temperature shift of 8-nm confinement from bulk between theoretical predictions and experimental results/Estimation on disjoining pressure and its effect on Kelvin pressure/Condensate growth model and model predictions/Condensation propagation and liquid bridging video. Notes The authors declare no competing financial interests. ACKNOWLEDGMENT The authors gratefully acknowledge support from Schlumberger Canada Ltd., Alberta InnovatesEnergy and Environment Solutions, and from the Natural Sciences and Engineering Council of Canada through a Collaborative Research and Development Grant, as well as on-going research funding through the Discovery Grants program and the Canada Research Chairs program. In addition, infrastructure funding provided by the Canada Foundation for Innovation and Ontario Research Fund is gratefully acknowledged. REFERENCES (1). Tröstl, J.; Chuang, W. K.; Gordon, H.; Heinritzi, M.; Yan, C.; Molteni, U.; Ahlm, L.; Frege, C.; Bianchi, F.; Wagner, R. The Role of Low-volatility Organic Compounds in Initial Particle Growth in the Atmosphere. Nature 2016, 533, 527-531. (2). Ehn, M.; Thornton, J. A.; Kleist, E.; Sipilä, M.; Junninen, H.; Pullinen, I.; Springer, M.; Rubach, F.; Tillmann, R.; Lee, B. A Large Source of Low-volatility Secondary Organic Aerosol. Nature 2014, 506, 476-479.

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

Page 18 of 20

(3). Kovács, T.; Meldrum, F.; Christenson, H. Crystal Nucleation without Supersaturation. J. Phys. Chem. Lett. 2012, 3, 1602-1606. (4). Li, X.; Hede, T.; Tu, Y.; Leck, C.; Ågren, H. Surface-Active cis-Pinonic Acid in Atmospheric Droplets: A Molecular Dynamics Study. J. Phys. Chem. Lett. 2010, 1, 769-773. (5). Lee, J.; Laoui, T.; Karnik, R. Nanofluidic Transport Governed by the Liquid/Vapour Interface. Nat. Nanotechnol. 2014, 9, 317-323. (6). Kim, H.; Yang, S.; Rao, S. R.; Narayanan, S.; Kapustin, E. A.; Furukawa, H.; Umans, A. S.; Yaghi, O. M.; Wang, E. N. Water Harvesting from Air with Metal-organic Frameworks Powered by Natural Sunlight. Science 2017, 356, 430-434. (7). Park, C. H.; Lee, S. Y.; Hwang, D. S.; Shin, D. W.; Cho, D. H.; Lee, K. H.; Kim, T.-W.; Kim, T.-W.; Lee, M.; Kim, D.-S. Nanocrack-regulated Self-humidifying Membranes. Nature 2016, 532, 480-483. (8). Lee, T.; Bocquet, L.; Coasne, B. Activated Desorption at Heterogeneous Interfaces and Long-Time Kinetics of Hydrocarbon Recovery from Nanoporous Media. Nat. Commun. 2016, 7, 11890. (9). Wu, K.; Chen, Z.; Li, X.; Guo, C.; Wei, M. A Model for Multiple Transport Mechanisms through Nanopores of Shale Gas Reservoirs with Real Gas Effect–adsorption-mechanic Coupling. Int. J. Heat Mass Transfer 2016, 93, 408-426. (10). Yuan, J.; Liu, X.; Akbulut, O.; Hu, J.; Suib, S. L.; Kong, J.; Stellacci, F. Superwetting Nanowire Membranes for Selective Absorption. Nat. Nanotechnol. 2008, 3, 332-6. (11). Lee, M.; Kim, B.; Kim, J.; Jhe, W. Noncontact Friction Via Capillary Shear Interaction at Nanoscale. Nat. Commun. 2015, 6, 7359. (12). Lin, C.-H.; Liu, X.; Wu, S.-H.; Liu, K.-H.; Mou, C.-Y. Corking and Uncorking a Catalytic Yolk-Shell Nanoreactor: Stable Gold Catalyst in Hollow Silica Nanosphere. J. Phys. Chem. Lett. 2011, 2, 2984-2988. (13). Miljkovic, N.; Enright, R.; Nam, Y.; Lopez, K.; Dou, N.; Sack, J.; Wang, E. N. Jumpingdroplet-enhanced Condensation on Scalable Superhydrophobic Nanostructured Surfaces. Nano Lett. 2012, 13, 179-187. (14). Tian, J.; Zhu, J.; Guo, H.-Y.; Li, J.; Feng, X.-Q.; Gao, X. Efficient Self-propelling of Small-scale Condensed Microdrops by Closely Packed ZnO Nanoneedles. J. Phys. Chem. Lett. 2014, 5, 2084-2088. (15). Nambo, A.; He, J.; Nguyen, T. Q.; Atla, V.; Druffel, T. L.; Sunkara, M. K. Ultrafast Carbon Dioxide Sorption Kinetics Using Lithium Silicate Nanowires. Nano Lett. 2017, 17, 3327–3333. (16). Wang, J.; Li, T.; Xia, B.; Jin, X.; Wei, H.; Wu, W.; Wei, Y.; Wang, J.; Liu, P.; Zhang, L. Vapor-condensation-assisted Optical Microscopy for Ultralong Carbon Nanotubes and Other Nanostructures. Nano Lett. 2014, 14, 3527-3533. (17). Wu, Q.; Ridge, C. J.; Zhao, S.; Zakharov, D.; Cen, J.; Tong, X.; Connors, E.; Su, D.; Stach, E. A.; Lindsay, C. M. Development of a New Generation of Stable, Tunable, and Catalytically Active Nanoparticles Produced by the Helium Nanodroplet Deposition Method. J. Phys. Chem. Lett. 2016, 7, 2910-2914. (18). Eijkel, J. C.; Van Den Berg, A. Water in Micro-and Nanofluidics Systems Described Using the Water Potential. Lab Chip 2005, 5, 1202-1209. (19). Zandavi, S. H.; Ward, C. A. Nucleation and Growth of Condensate in Nanoporous Materials. Phys. Chem. Chem. Phys. 2015, 17, 9828-34.

ACS Paragon Plus Environment

18

Page 19 of 20 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

(20). Fisher, L.; Israelachvili, J. Direct Experimental Verification of the Kelvin Equation for Capillary Condensation. Nature 1979, 277, 548-549. (21). Tan, S. P.; Piri, M. Equation-of-state Modeling of Confined-fluid Phase Equilibria in Nanopores. Fluid Phase Equilib. 2015, 393, 48-63. (22). Tan, S. P.; Piri, M. Equation-of-State Modeling of Associating-fluids Phase Equilibria in Nanopores. Fluid Phase Equilib. 2015, 405, 157-166. (23). Schneider, D.; Valiullin, R.; Monson, P. A. Filling Dynamics of Closed End Nanocapillaries. Langmuir 2014, 30, 1290-1294. (24). Morishige, K.; Kittaka, S. Kinetics of Capillary Condensation of Water in Mesoporous Carbon: Nucleation and Meniscus Growth. J. Phys. Chem. C 2015, 119, 18287-18292. (25). Vincent, O.; Marguet, B.; Stroock, A. D. Imbibition Triggered by Capillary Condensation in Nanopores. Langmuir 2017, 33, 1655-1661. (26). Barsotti, E.; Tan, S. P.; Saraji, S.; Piri, M.; Chen, J.-H. A Review on Capillary Condensation in Nanoporous Media: Implications for Hydrocarbon Recovery from Tight Reservoirs. Fuel 2016, 184, 344-361. (27). Nune, S. K.; Lao, D. B.; Heldebrant, D. J.; Liu, J.; Olszta, M. J.; Kukkadapu, R. K.; Gordon, L. M.; Nandasiri, M. I.; Whyatt, G.; Clayton, C. Anomalous Water Expulsion from Carbon-based Rods at High Humidity. Nat. Nanotechnol. 2016, 11, 791-797. (28). Xu, K.; Cao, P.; Heath, J. R. Graphene Visualizes the First Water Adlayers on Mica at Ambient Conditions. Science 2010, 329, 1188-1191. (29). Dukhin, S. S.; Zimmermann, R.; Werner, C. Surface Conductivity Reveals Counterion Condensation within Grafted Polyelectrolyte Layers. J. Phys. Chem. B 2007, 111, 979-981. (30). Park, M. J.; Downing, K. H.; Jackson, A.; Gomez, E. D.; Minor, A. M.; Cookson, D.; Weber, A. Z.; Balsara, N. P. Increased Water Retention in Polymer Electrolyte Membranes at Elevated Temperatures Assisted by Capillary Condensation. Nano Lett. 2007, 7, 3547-3552. (31). Aswal, V.; Kohlbrecher, J.; Goyal, P.; Amenitsch, H.; Bernstorff, S. Counterion Condensation on Charged Micelles in an Aqueous Electrolyte Solution as Studied with Combined Small-angle Neutron Scattering and Small-angle X-ray Scattering. J. Phys.: Condens. Matter 2006, 18, 11399. (32). Li, J.; Li, X.; Wu, K.; Feng, D.; Zhang, T.; Zhang, Y. Thickness and Stability of Water Film Confined inside Nanoslits and Nanocapillaries of Shale and Clay. Int. J. Coal Geol. 2017, 179, 253-268. (33). Lei, W.; McKenzie, D. R. Nanoscale Capillary Flows in Alumina: Testing the Limits of Classical Theory. J. Phys. Chem. Lett. 2016, 7, 2647-2652. (34). Berim, G. O.; Ruckenstein, E. Adsorption Isoterms and Capillary Condensation in a Nanoslit with Rough Walls: A Density Functional Theory. Langmuir 2012, 28, 11384-11390. (35). Dalla Bernardina, S.; Paineau, E.; Brubach, J.-B.; Judeinstein, P.; Rouzière, S.; Launois, P.; Roy, P. Water in Carbon Nanotubes: The Peculiar Hydrogen Bond Network Revealed by Infrared Spectroscopy. J. Am. Chem. Soc. 2016, 138, 10437-10443. (36). Li, Y.; Alibakhshi, M. A.; Zhao, Y.; Duan, C. Exploring Ultimate Water Capillary Evaporation in Nanoscale Conduits. Nano Lett. 2017, 17, 4813-4819. (37). Li, L.; Kazoe, Y.; Mawatari, K.; Sugii, Y.; Kitamori, T. Viscosity and Wetting Property of Water Confined in Extended Nanospace Simultaneously Measured from Highly-pressurized Meniscus Motion. J. Phys. Chem. Lett. 2012, 3, 2447-2452. (38). Sparreboom, W.; van den Berg, A.; Eijkel, J. C. Principles and Applications of Nanofluidic Transport. Nat. Nanotechnol. 2009, 4, 713-720.

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 20

(39). Bao, B.; Zandavi, S. H.; Li, H.; Zhong, J.; Jatukaran, A.; Mostowfi, F.; Sinton, D. Bubble Nucleation and Growth in Nanochannels. Phys. Chem. Chem. Phys. 2017, 19, 8223-8229. (40). Park, S.; Jung, Y.; Son, S. Y.; Cho, I.; Cho, Y.; Lee, H.; Kim, H.-Y.; Kim, S. J. Capillarity Ion Concentration Polarization as Spontaneous Desalting Mechanism. Nat. Commun. 2016, 7, 11223. (41). Zhong, J.; Zandavi, S. H.; Li, H.; Bao, B.; Persad, A. H.; Mostowfi, F.; Sinton, D. Condensation in One-dimensional Dead-end Nanochannels. ACS Nano 2017, 11, 304-313. (42). Li, H.; Zhong, J.; Pang, Y.; Zandavi, S. H.; Persad, A. H.; Xu, Y.; Mostowfi, F.; Sinton, D. Direct Visualization of Fluid Dynamics in Sub-10 nm Nanochannels. Nanoscale 2017, 9, 9556-9561. (43). P.J. Linstrom and W.G. Mallard, E., Nist Chemistry Webbook, Nist Standard Reference Database Number 69. National Institute of Standards and Technology, Gaithersburg MD, 20899, 2016. (44). Churaev, N.; Sobolev, V. Prediction of Contact Angles on the Basis of the FrumkinDerjaguin Approach. Adv. Colloid Interface Sci. 1995, 61, 1-16. (45). Fisher, L.; Gamble, R.; Middlehurst, J. The Kelvin Equation and the Capillary Condensation of Water. Nature 1981, 290, 575-576. (46). Awschalom, D.; Warnock, J.; Shafer, M. Liquid-film Instabilities in Confined Geometries. Phys. Rev. Lett. 1986, 57, 1607. (47). Kast, W.; Hohenthanner, C.-R. Mass Transfer within the Gas-phase of Porous Media. Int. J. Heat Mass Transfer 2000, 43, 807-823. (48). Rossi, M. P.; Ye, H.; Gogotsi, Y.; Babu, S.; Ndungu, P.; Bradley, J.-C. Environmental Scanning Electron Microscopy Study of Water in Carbon Nanopipes. Nano Lett. 2004, 4, 989993. (49). Phan, V. N.; Nguyen, N. T.; Yang, C.; Joseph, P.; Djeghlaf, L.; Bourrier, D.; Gue, A. M. Capillary Filling in Closed End Nanochannels. Langmuir 2010, 26, 13251-5. (50). Haneveld, J.; Tas, N. R.; Brunets, N.; Jansen, H. V.; Elwenspoek, M. Capillary Filling of Sub-10 nm Nanochannels. J. Appl. Phys. 2008, 104, 014309. (51). Persad, A. H.; Ward, C. A. Expressions for the Evaporation and Condensation Coefficients in the Hertz-Knudsen Relation. Chem. Rev. 2016, 116, 7727-7767. (52). Zandavi, S. H.; Ward, C. Vapour Adsorption Kinetics: Statistical Rate Theory and Zeta Adsorption Isotherm Approach. Phys. Chem. Chem. Phys. 2016, 18, 25538-25545. (53). Xie, Q.; Xiao, S.; Duan, C. Geometry-dependent Drying in Dead-end Nanochannels. Langmuir 2017, 33, 8395-8403. (54). Weijs, J. H.; Lohse, D. Why Surface Nanobubbles Live for Hours. Phys. Rev. Lett. 2013, 110, 054501. (55). Iwamatsu, M.; Horii, K. Capillary Condensation and Adhesion of Two Wetter Surfaces. J. Colloid Interface Sci. 1996, 182, 400-406. (56). Zeng, Y.; Phadungbut, P.; Do, D.; Nicholson, D. Anatomy of Adsorption in Open-end and Closed-end Slit Mesopores: Adsorption, Desorption, and Equilibrium Branches of Hysteresis Loop. J. Phys. Chem. C 2014, 118, 25496-25504. (57). Nguyen, P. T.; Do, D.; Nicholson, D. On the Irreversibility of the Adsorption Isotherm in a Closed-end Pore. Langmuir 2013, 29, 2927-2934.

ACS Paragon Plus Environment

20