Wettability of Al2O3 Surface by Organic Molecules: Insights from

Sep 29, 2017 - We use molecular dynamics (MD) simulations to investigate the wettability of Al2O3 (0001) by organic molecules. Diffusion coefficients ...
1 downloads 11 Views 2MB Size
Subscriber access provided by BOSTON UNIV

Article 2

3

Wettability of AlO Surface by Organic Molecules: Insights from Molecular Dynamics Simulation Yoshitake Suganuma, Satoru Yamamoto, Tomoyuki Kinjo, Takuya Mitsuoka, and Kazuhiko Umemoto J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b07062 • Publication Date (Web): 29 Sep 2017 Downloaded from http://pubs.acs.org on September 30, 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 B 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 37

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

Wettability of Al2O3 Surface by Organic Molecules: Insights from Molecular Dynamics Simulation Yoshitake Suganuma,*,† Satoru Yamamoto,‡ Tomoyuki Kinjo, † Takuya Mitsuoka, † and Kazuhiko Umemoto† †Toyota Central R&D Labs., Inc., 41-1, Yokomichi, Nagakute, Aichi 480-1192, Japan ‡Dassault Systèmes BIOVIA, 2-1-1, Osaki, Shinagawa, Tokyo 141-6020, Japan

ABSTRACT

We use molecular dynamics (MD) simulations to investigate the wettability of Al2O3 (0001) by organic molecules. Diffusion coefficients estimated for organic molecules are clearly correlated with the contact angles observed experimentally. The results of the MD simulations suggest that molecular flexibility influences wettability. In other words, wettability owing to flexible molecules, such as an epoxy tridecamer, improves with increasing temperature because the interaction between the droplet and the surface increases due to changes in molecular conformation. Conversely, for phenylene sulfide tetramer, wettability does not change with temperature because

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

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 37

of the molecular rigidity. In addition, for epoxy monomers, we analyze the different molecular structures responsible for modifying the droplet–surface interaction. For hydrogens in aromatic rings and in methyl groups, the interaction with the surface clearly decreases with increasing temperature.

ACS Paragon Plus Environment

2

Page 3 of 37

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

Introduction The control of wettability is a challenging task for industrial processes and applications such as coating, printing, antifouling treatment, and adhesion.1–10 Numerous studies address the question of how chemistry and the topography of solid surfaces affect wettability.11–18 Surface characteristics are important parameters for controlling wettability. In addition, for liquids, several studies report a relation between molecular structure and wettability.19–21 Controlling wettability requires a thorough understanding of the interaction between the liquid and the solid surface, which is also related to molecular conformation.22,23 Recently, molecular dynamics (MD) simulations have been used to study the interactions between droplets and solid surfaces, and the results suggest that MD simulations constitute a powerful method for understanding interactions and conformational dynamics in wetting.24–30 Most such research has focused on the wettability of carbon materials such as graphite and carbon nanotubes by water. However, few studies simulate the process by which organic molecules wet metals and their oxides, which are industrially significant materials. Yet, the conformational dynamics of polymers in the wetting of metals remains poorly understood. To address this issue, we herein report the results of MD simulations for understanding the process of polymer wetting on metals. In addition, the validity of the simulations is confirmed by comparing with experimentally obtained contact angles. We study typical organic molecules such as solvents and oligomers, which are used as model polymers. The solid surface is the Al2O3 (0001) surface, which models the natural oxidation film on aluminum, a typical metal used in industry and an important alternative lightweight metal for iron. Based on the results of the

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

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 37

simulations, we discuss the interaction between the liquid and the solid surface and conformational dynamics of organic molecules in wetting of the Al2O3 (0001) surface.

Experiment Materials Water was purified with a Milli-Q deionizing system. Diiodomethane, heptadecane, ethylene glycol, and diethylene glycol were purchased from Wako Pure Chemicals Industries Ltd. (Osaka, Japan) and used as received. 1-bromonaphtalene and octylbenzene were obtained from Tokyo Chemical Industry Co. (Tokyo, Japan) and used without further purification. The epoxy oligomers jER 828 and 1009 were obtained from Mitsubishi Chemical Co.(Tokyo, Japan) and used as bisphenol A diglycidyl ether monomer (epoxy monomer) and tridecamer (epoxy tridecamer), respectively, based on molecular weights from catalog value. Phenylene sulfide (PS) tetramer and heptamer were synthesized as model poly-PS compounds.31, 32 A highly polished single crystal (0001) wafer of α-Al2O3 with a rms roughness of about 1 Å, as determined by atomic force microscopy, was obtained from Kyocera Co. (Kyoto, Japan). The wafer was cleaned by following the procedure detailed in previous research.33 The wafer was dipped in 10 mM HNO3 for 30 min and rinsed by Milli-Q water. Next, it was heated to 350 °C in the ambient air for 30 min, washed with Milli-Q water, and finally dried in the ambient air at room temperature.

ACS Paragon Plus Environment

4

Page 5 of 37

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

Contact-Angle Measurement. Contact angles were measured by using a Dropmaster DM-500 Hi (Kyowa Interface Science Co., Saitama, Japan), which offers both probe and stage heating. Solvents were placed onto an Al2O3 (0001) surface at room temperature (298 K). For epoxy and PS oligomers, contact-angle measurements were done with the probe and stage at 423 and 498 K in addition to room temperature (298 K). Next, measurements were made to determine whether the oligomers were in the molten state, and, if so, the contact angles were measured at those temperatures. The droplet volume was about 1 μL, which is sufficiently small to neglect gravity.34 The contact angles were measured only after the droplets had reached the steady state on the substrate. The contact angles were determined by using a contour curve-fitting method implemented in FAMAS software (Kyowa Interface Science Co., Saitama, Japan).

Results Table 1 lists the liquids for which the contact angle was measured and the temperatures at which the measurements were made. The contact angles reported in Table 1 are averages of six to ten measurements using two substrates. The contact angle as a function of temperature for epoxy oligomers differs from that of PS oligomers. For epoxy oligomers, the contact angles decrease with increasing temperature. However, for PS oligomers, the contact angle changes only slightly. The different results for epoxy oligomers vis-à-vis PS oligomers were analyzed by using MD simulations in terms of changes in molecular conformation.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

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 37

Table 1. Liquids used for contact-angle measurements and temperature at which contact-angle measurements were made.

Liquid water diiodomethane 1-bromonaphthalene heptadecane octylbenzene ethylene glycol diethylene glycol epoxy monomer

epoxy tridecamer PS tetramer PS heptamer

Temperature (K) 298 298 298 298 298 298 298 298 423 493 423 493 423 493 493

Contact angle (°) 12.8 38.2 21.2 10.1 6.4 14.0 12.0 16.7 7.7 2.5 46.8 30.7 23.9 24.6 39.9

Simulations Methodology The following calculations were done by using the software package Materials Studio ver. 8.0 (Dassault Systèmes BIOVIA, San Diego, CA, 2001–2009). The model consisted of an Al2O3 (0001) surface and liquid droplets. To model the Al2O3 (0001) surface, the Al2O3 lattice was first derived from the structural database of Materials Studio. A repeat unit of Al2O3 was cleaved along the (0001) crystallographic orientation so that the surface consists of oxygen atoms based on crystal

ACS Paragon Plus Environment

6

Page 7 of 37

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

truncation rod diffraction.33 Next, a vacuum slab was created above the surface so that the lattice size along the z axis was 400 Å. The surface geometry was optimized by using the DMol3 module with the PBE functional and numerical DNP basis set. An orthorhombic lattice was prepared by redefining the vector A to (210). Next, a 24×42×1 supercell was constructed to obtain a surface with dimensions of 197.8 × 199.9 × 400 Å3. An amorphous cell was used for constructing solvent and oligomer bulks, which were optimized by MD simulations for 1000 ps in the canonical NPT ensemble. From the bulk constructions, we built droplets with a radius of 20.0 Å, which is sufficient to measure a contact angle in a simulation, as reported in previous research.29 The MD simulations were conducted by the COMPASSII35 force field, which was implemented by using the Forcite module. Finally, the solvent and oligomer droplets were placed onto the Al2O3 (0001) surface to produce the initial construction shown in Figure 1.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

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 37

Figure 1. Initial construction for MD simulation of wetting the Al2O3 (0001) surface.

The canonical ensemble NVT was used for MD wetting simulations. The integration step was set to 1.0 fs. The Berendsen loose-coupling thermostat36 was selected for finite-temperature control of the simulated system. The van der Waals interactions and electrostatic interactions were calculated by using the group-based method, and the cutoff distance was 15.5 Å. During the simulations, the Al2O3 surface was fixed because these atoms vibrated only very slightly and so the vibration was neglected. The MD simulations were done in two steps: In the first step, the systems were evolved for 500 ps to assess the wetting process (the first MD simulation). In the second step, the systems were evolved for 1000 ps for an analysis of the steady state (the second

ACS Paragon Plus Environment

8

Page 9 of 37

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

MD simulation). These MD simulations are also conducted by the COMPASSII force field. This force field seems to be proper for calculating wettability as reported previous research.29

Results and Discussion First, to confirm the validity of the MD simulations, we tried to evaluate contact angles based on the simulation results. However, we could not obtain contact angles because the droplets were fully spread on the surface, so their contact angles were near zero for most solvents and oligomers. We attribute this result to the small scale of MD simulations compared to the experiments. Then, we calculated we calculated the diffusion coefficient D instead of contact angle, which is should be related to wettability, 29 from the first 500 ps in the second MD simulation. The diffusion coefficient is given by following the equation.37

D

lim →

1 6

〈|

0 | 〉

Where Ri(t) is the position of atom i at time t, and n is the number of atoms. This equation is used for three-dimensional systems, although molecules hardly move along the z axis because the organic molecules exist on the solid surface in the present simulations. Therefore, diffusion coefficients are underestimated; however, the degrees of underestimation seem to be almost same in all simulations, so that they are sufficient to confirm the validation of the MD simulations. We compared the diffusion coefficient to the contact angle directly rather than other measures such as the energy of adhesion, which depends on some experimental values and can cause larger

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

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 37

error bars. Figure 2 shows the experimentally determined contact angle as a function of the logarithm of the simulated diffusion coefficient. The error bars in MD simulations seems to be small because mean square displacements have good linearity with respect to the simulation time. The diffusion coefficients show a clear negative correlation with contact angle. Note that diiodomethane departs significantly from other plots because COMPASSII contains force fields only for I–, I, I2–, and I3 but not for I bonded to carbon. Excluding diiodomethane, the correlation coefficient is –0.79, which is sufficient correlation to estimate the wettability of the Al2O3 (0001) surface by MD simulation.

Figure 2. Experimentally obtained contact angles as a function of the logarithm of the diffusion coefficient obtained from MD simulation. At right are snapshots of droplets of epoxy oligomers and a PS oligomer in contact-angle measurements.

ACS Paragon Plus Environment

10

Page 11 of 37

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

To study the wettability as a function of temperature, Figure 3 shows the cohesive energy densities of epoxy oligomers and PS oligomers, which is correlated with surface tension.38,39 The cohesive energy density is given by following the equation.40 〈







Where Eintra is the intramolecular energy, and Etotal is the total energy of the system. The brackets represent an average over a NPT ensemble. The cohesive energy densities are calculated from the MD simulations to optimize the construction of oligomers, and the averages over the final 500 ps in the simulations appear in Figure 3.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

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 37

Figure 3. Cohesive energy density of (a) epoxy monomers, (b) epoxy tridecamers, and (c) PS tetramers.

ACS Paragon Plus Environment

12

Page 13 of 37

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

The results indicate that the cohesive energy densities become smaller with increasing temperature for all oligomers, as is well known.41 According to Young’s equation,42 contact angles are determined by the surface tension of a droplet and the interaction between the droplet and the surface, and smaller surface tension leads to a decrease of contact angles. However, for the PS oligomer, the contact angles are almost the same at 423 and 493 K, in spite of the decrease in cohesive energy density. Consequently, we also analyze the interactions of oligomers with the surface to understand the difference between epoxy oligomers and PS oligomers. To evaluate the interaction between droplet and surface, we examine the movement of hydrogen near the surface, which seems to interact with the surface by some interaction including the hydrogen bond or van der Waals. In the first snapshot of the second MD simulation, we define hydrogens less than 3 Å from the top oxygens of the Al2O3 surface as interacting hydrogens (Hint). As an example, Figure 4 shows Hint for one molecule in an epoxy monomer. You can see that most of the nearest hydrogens from the surface are selected as Hint in Figure 4. Then, the 3 Å distance seems to be sufficient to evaluate the interaction with the surface. Next, from the second MD simulation, we obtain the profile of the relative Hint concentration along the z axis over 1000 ps. Figure 5 shows these profiles for epoxy oligomers and PS oligomers.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

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 37

Figure 4. Hint (yellow atoms) of a single epoxy monomer molecule.

ACS Paragon Plus Environment

14

Page 15 of 37

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

Figure 5. Relative concentration of Hint along z axis over 1000 ps from the second MD simulations for (a) epoxy monomer, (b) epoxy tridecamer, and (c) PS tetramer.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

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 37

For the epoxy monomer, only a single peak appears at 298 K in Figure 5(a), which indicates that all Hint are bound near the surface over 1000 ps by interaction with the surface. Conversely, at 423 and 493 K, the first peak becomes smaller and additional peaks appear, which suggests that some Hint begin to move against surface interaction. In other words, the interaction strength between epoxy monomer molecules and the surface decreases with increasing temperature. According to Young’s equation,40 a weaker interaction between droplet and surface leads to greater contact angles. Thus, for epoxy monomers, the decrease in contact angles with increasing temperature is attributed not to the change in interaction strength but to the decrease in surface tension (see Figure 3). For epoxy tridecamer, the first peak at 493 K is larger than that at 423 K [see Figure 5(b)], in contrast to the case of the epoxy monomer. This indicates that more hydrogen atoms interact with the surface at 493 K than at 423 K. In addition to the decrease in surface tension shown in Figure 3, this is considered to result in smaller contact angles at 493 K. To investigate the increase in droplet–surface interaction, Figure 6 shows the molecular radius of gyration as a function of time over 500 ps from the first MD simulation. The radius of gyration of epoxy tridecamer at 423 K is less than that at 0 ps, which implies that epoxy tridecamer molecules experience intramolecular interactions at 423 K. Conversely, the radius of gyration at 493 K increases in the first ~100 ps. These results suggest that, at 493 K, thermal energy allows the epoxy tridecamer molecules to interact with the surface against intramolecular interaction. As a result, the interaction between droplet molecule and surface increases and the contact angle decreases with increasing temperature.

ACS Paragon Plus Environment

16

Page 17 of 37

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

Figure 6. Radius of gyration as a function of time for epoxy tridecamer.

Conversely, for PS tetramer, the profiles of Hint are nearly same at 423 and 493 K [see Figure 5(c)], which implies that little change occurs in the droplet–surface interaction. To assess the difference vis-à-vis the epoxy tridecamer, Figure 7 shows the radius of gyration as a function of time for PS tetramer and heptamer. The profiles of PS tetramer at 423 and 493 K are also nearly the same and remain almost constant over 500 ps. Even for the PS heptamer, which is longer than the PS tetramer, the radius of gyration remains constant in contrast to the epoxy tridecamer. These results suggest that, for the epoxy tridecamer, molecular flexibility allows it to change conformation so as to increase the interaction with the surface and decrease the contact angle with increasing temperature. However, because of its molecular rigidity, PS tetramer cannot change its conformation and interaction with the surface. Thus, increasing temperature has little effect on the contact angle for PS tetramer.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

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 37

Figure 7. Radius of gyration as a function of time for PS oligomers.

Indeed, based on the results for the epoxy tridecamer and PS tetramer, the flexibility of molecules seems to effect the droplet–surface interaction and lead to the different changes in contact angle with temperature for the epoxy tridecamer and PS tetramer. Subsequently, the effect of molecular flexibility on the interaction is also examined for epoxy monomer. Figure 8 shows the radius of gyration as a function of time for the epoxy monomer at 298, 423, and 493 K. All profiles are almost the same and remain constant because the epoxy monomer molecule is so short that it is at maximum extension, as shown in Figure 4. Accordingly, for the epoxy monomer, molecular flexibility does not affect the droplet–surface interaction strength.

ACS Paragon Plus Environment

18

Page 19 of 37

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

Figure 8. Radius of gyration as a function of time for epoxy monomer.

We now analyze the decrease in interaction strength between epoxy monomer molecules and the surface by evaluating the changes in molecular conformation. For this purpose, all hydrogens (Hall) and Hint in Figure 4 are classified as three types: hydrogens bonded to a carbon next to an oxygen (HnearO), hydrogens in aromatic rings (Hring), and hydrogen in methyl groups (Hmethyl). Next, as for the profiles shown in Figure 5, Figure 9 shows profiles of relative concentration along the z axis over 1000 ps from the second MD simulation for the three types of hydrogens. The profiles of HnearO in Hall and HnearO in Hint are almost the same [Figures 9(a) and 9(b)], which indicates that most of the HnearO exists near the surface and interacts with it. Moreover, temperature has little effect on the profiles of HnearO, which suggests that most of the HnearO are bound close to the surface regardless of temperature; that is, HnearO seems to interact strongly with the surface.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

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 37

Figure 9. Relative concentration profiles of Hall and Hint, classified by their bonding position, along z axis over 1000 ps from the second MD simulation.

ACS Paragon Plus Environment

20

Page 21 of 37

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

Two peaks appear for Hring in Hall in Figure 9(c), which suggests that aromatic rings in the epoxy monomer molecule are not parallel to the surface but tilted, as shown in Figure 4. In addition, the profile of Hring in Hint has only one peak at 298 K [Figure 9(d)], which indicates that Hring is bound to the surface due to interactions with the surface. In other words, the aromatic rings hardly move at 298 K. However, the profiles at 423 and 498 K change significantly compared with that at 298 K: the small-distance peaks shrink and, simultaneously, the new peaks grow [Figure 9(d)]. This suggests that the hydrogens in aromatic rings begin to move over the temperature transition, which seems to happen because they cannot resist thermal motion and so interact more weakly with the surface. As a result, the aromatic rings begin to tilt to the other side at 423 and 493 K. Finally, two split peaks appear in the profiles of Hmethyl in Hall [Figure 9(e)] that seem to correspond to two methyl groups in the molecule and their rotational movement. The profiles of Hmethyl in Hint show only a single split peak [Figure 9(f)], in contrast with the profiles in Figure 9(e). This means that only one methyl group in the molecule is adsorbed onto the surface. Furthermore, the first peak shrinks and the second peak increases in Figure 9(f), which indicates that the methyl groups rotate more freely at 423 and 493 K than at 298 K. On the whole, for the epoxy monomer, the interaction of Hring and Hmethyl with the surface decreases significantly with increasing temperature, so that aromatic rings begin to tilt to the other side or the methyl groups begin to rotate freely at the higher temperatures. Thus, the wettability of Al2O3 (0001) by organic molecules can be simulated by MD simulation. In addition, the change in conformation of organic molecules influences the liquid-surface interaction and contact angle.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

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 37

Conclusion The process by which organic molecules wet an Al2O3 (0001) surface is analyzed by using MD simulations. The diffusion coefficients calculated by these simulations provide satisfactory estimates of the wettability of the Al2O3 (0001) surface. Moreover, the effect of changes in molecular conformation is evaluated for several oligomers. The results indicate that the molecular flexibility affects the molecule–surface interaction. In addition, for the epoxy monomer, the interaction with the surface of hydrogens in aromatic rings and in methyl groups decreases dramatically with increasing temperature. Thus, this method can aid in understanding the consequences of molecular structure for the wettability of Al2O3 (0001) surfaces and thereby help to control the wettability of this surface.

ACS Paragon Plus Environment

22

Page 23 of 37

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

AUTHOR INFORMATION Corresponding Author [email protected]

REFERENCES (1) Son, Y.; Kim, C.; Yang, D. H.; Ahn, D. J. Spreading of an Inkjet Droplet on a Solid Surface with a Controlled Contact Angle at Low Weber and Reynolds Numbers. Langmuir 2008, 24, 2900– 2907. (2) Perelaer, J.; Hendriks, C. E.; de Laat, A. W. M.; Schubert, U. S. One-step Inkjet Printing of Conductive Silver Tracks on Polymer Substrates. Nanotechnology 2009, 20, 165303. (3) Sakai, M.; Yanagisawa, T.; Nakajima, A.; Kameshima, Y.; Okada, K. Effect of Surface Structure on the Sustainability of an Air Layer on Superhydrophobic Coatings in a Water–Ethanol Mixture. Langmuir 2009, 25, 13–16. (4) Gudipati, C. S.; Finlay, J. A.; Callow, J. A.; Callow, M. E.; Woolley, K. L. The Antifouling and Fouling-Release Perfomance of Hyperbranched Fluoropolymer (HBFP)–Poly(ethylene glycol) (PEG) Composite Coatings Evaluated by Adsorption of Biomacromolecules and the Green Fouling Alga Ulva. Langmuir 2005, 21, 3044–3053. (5) Fan, X.; Lin, L.; Dalsin, J. L.; Messersmith, P. B. Biomimetic Anchor for Surface-Initiated Polymerization from Metal Substrates. J. Am. Chem. Soc. 2005, 127, 15843–15847. (6) Ravichandran, V.; Obendorf, S. K. Wettability and Adhesion Studies of Grafted Poly(pphenylene terephthalamide) Fiber Surfaces. J. Adhes. Sci. Technol. 1992, 6, 1303–1323. (7) Coelho, M. A. N.; Vieira, E. P.; Motschmann, H.; Möhwald, H.; Thünemann, A. F. Human Serum Albumin on Fluorinated Surfaces. Langmuir 2003, 19, 7544–7550.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

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 24 of 37

(8) Zeng, H.; Huang, J.; Tian, Y.; Li, L.; Tirrell, M. V.; Israelachvili, J. N. Adhesion and Detachment

Mechanisms

between

Polymer

and

Solid

Substrate

Surfaces:

Using

Polystyrene−Mica as a Model System. Macromolecules 2016, 49, 5223−5231. (9) Yerushalmi-Rozen, R.; Klein, J.; Fetters, L. J. Suppression of Rupture in Thin, Nonwetting Liquid Films. Science 1994, 263, 793–795. (10) Tsukruk, V. V.; Bliznyuk, V. N. Adhesive and Friction Forces between Chemically Modified Silicon and Silicon Nitride Surfaces. Langmuir 1998, 14, 446–455 (11) Luzinov, I.; Minko, S.; Senkovsky, V.; Voronov, A.; Hild, S.; Marti, O.; Wilke, W. Synthesis and Behavior of the Polymer Covering on a Solid Surface. 3. Morphology and Mechanism of Formation of Grafted Polystyrene Layers. Macromolecules 1998, 31, 3945–3952. (12) Öner, D.; McCarthy, T. J. Ultrahydrophobic Surfaces. Effects of Topography Length Scales on Wettability. Langmuir 2000, 16, 7777–7782. (13) Ming, W.; Wu, D.; van Benthem, R.; de With, G. Superhydrophobic Films from Raspberrylike Particles. Nano Lett. 2005, 5, 2298–2301. (14) Grundke, K.; Zschoche, S.; Pöschel, K.; Gietzelt, T.; Michel, S.; Friedel, P.; Jehnichen, D.; Neumann, A. W. Wettability of Maleimide Copolymer Films: Effect of the Chain Length of nAlkyl Side Groups on the Solid Surface Tension. Macromolecules 2001, 34, 6768–6775. (15) Synytska, A.; Appelhans, D.; Wang, Z. G.; Simon, F.; Lehmann, F.; Stamm, M.; Grundke, K. Perfluoroalkyl End-Functionalized Oligoesters: Correlation between Wettability and EndGroup Segregation. Macromolecules 2007, 40, 297–305. (16) Lee, J. A.; McCarthy, T. J. Polymer Surface Modification: Topography Effects Leading to Extreme Wettability Behavior. Macromolecules 2007, 40, 3965–3969.

ACS Paragon Plus Environment

24

Page 25 of 37

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

(17) Takei, Y. G.; Aoki, T.; Sanui, K.; Ogata, N.; Sakurai, Y.; Okano, T. Dynamic Contact Angle Measurement of Temperature-Responsive Surface Properties for Poly(N-isopropylacrylamide) Grafted Surfaces. Macromolecules 1994, 27, 6163–6166. (18) Reiter, G.; Khanna, R. Kinetics of Autophobic Dewetting of Polymer Films. Langmuir 2000, 16, 6351–6357. (19) Reiter, G.; Auroy, P.; Auvray, L. Instabilities of Thin Polymer Films on Layers of Chemically Identical Grafted Molecules. Macromolecules 1996, 29, 2150–2157. (20) van Oss, C. J.; Good, R. J.; Chaudhury, M. K. Additive and Nonadditive Surface Tension Components and the Interpretation of Contact Angles. Langmuir 1988, 4, 884–891. (21) Kwok, D. Y.; Neumann, A. W. Contact Angle Measurement and Contact Angle Interpretation. Adv. Colloid Interface Sci. 1999, 81, 167–249. (22) Tao, Y. T. Structural Comparison of Self-Assembled Monolayers of n-Alkanoic Acids on the Surfaces of Silver, Copper, and Aluminum. J. Am. Chem. Soc. 1993, 115, 4350–4358. (23) Steiner, T.; Saenger, W. Geometry of C-H-O Hydrogen Bonds in Carbohydrate Crystal Structures. Analysis of Neutron Diffraction Data. J. Am. Chem. Soc. 1992, 114, 10146–10154. (24) Cao, Z.; Stevens, M. J.; Dobrynin, A. V. Elastocapillarity: Adhesion and Wetting in Soft Polymeric Systems. Macromolecules 2014, 46, 6515–6521. (25) Nieto, D. R.; Santese, F.; Toth, R.; Posocco, P.; Pricl, S.; Fermeglia, M. Simple, Fast, and Accurate In silico Estimations of Contact Angle, Surface Tension, and Work of Adhesion of Water and Oil Nanodroplets on Amorphous Polypropylene Surfaces. Appl. Mater. Interfaces 2012, 4, 2855−2859.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

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 26 of 37

(26) Werder, T.; Walther, J. H.; Jaffe, R. L.; Halicioglu, T.; Noca, F.; Koumoutsakos, P. Molecular Dynamics Simulation of Contact Angles of Water Droplets in Carbon Nanotubes. Nano Lett. 2001, 1, 697–702. (27) Giovambattista, N.; Debenedetti, P. G.; Rossky, P. J. Effect of Surface Polarity on Water Contact Angle and Interfacial Hydration Structure. J. Phys. Chem. B 2007, 111, 9581–9587. (28) de Ruijter, M. J.; Blake, T. D.; De Coninck, J. Dynamic Wetting Studied by Molecular Modeling Simulations of Droplet Spreading. Langmuir 1999, 15, 7836–7847. (29) Zhang, J.; Zhong, J.; Li, W.; Wang, M.; Liu, B.; Li, Z.; Yan, Y. Molecular Insight into the Dynamical Adsorption Behavior of Nanoscale Water Droplets on a Heterogeneous Surface. RSC Adv. 2015, 5, 52322–52329. (30) Maheshwari, S.; van der Hoef, M.; Lohse, D. Line Tension and Wettability of Nanodrops on Curved Surfaces. Langmuir 2016, 32, 316−321. (31) Vicente, J.; Abad, J. A.; López-Nicolás, R. M. A New Approach to the Synthesis of Oligomers. Application to the Synthesis of P-phenylene Thioether Wires. Tetrahedron Lett. 2005, 46, 5839–5840. (32) Reddy, V. P.; Swapna, K.; Kumar, A. V.; Rao, K. R. Indium-Catalyzed C-S Cross-Coupling of Aryl Halides with Thiols. J. Org. Chem. 2009, 74, 3189–3191. (33) Eng, P. J.; Trainor, T. P.; Brown Jr., G. E.; Waychunas, G. A.; Newville, M.; Sutton, S. R.; Rivers, M. L. Structure of the Hydrated α-Al2O3 (0001) Surface. Science 2000, 288, 1029–1033. (34) Bonn, D.; Eggers, J.; Indekeu, J.; Meunier, J.; Rolley, E. Wetting and Spreading. Rev. Mod. Phys. 2009, 81, 739−805.

ACS Paragon Plus Environment

26

Page 27 of 37

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

(35) Sun, H.; Jin, Z.; Yang, C.; Akkermans, R. L. C.; Robertson, S. H.; Spenley, N. A.; Miller, S.; Todd, S. M. COMPASS II: Extended Coverage for Polymer and Drug-like Molecule Databases. J. Mol. Model. 2016, 22, 47–57. (36) Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690. (37) Salles, F.; Jobic, H.; Devic, T.; Llewellyn, P. L.; Serre, C.; Férey, G.; Maurin, G. Self and Transport Diffusivity of CO2 in the Metal–Organic Framework MIL-47(V) Explored by Quasielastic Neutron Scattering Experiments and Molecular Dynamics Simulations. ACS Nano 2010, 4, 143–152. (38) Schonhorn, H. Theoretical Relationship between Surface Tension and Cohesive Energy Density. J. Chem. Phys. 1965, 43, 2041–2043. (39) Becher, P. The Calculation of Cohesive Energy Density from the Surface Tension of Liquids. J. Colloid Interface Sci. 1972, 38, 291–293. (40) Hildebrand, J. H., Prausnitz, J. M., Scott, R. L. Regular and Related Solutions; Van Nostrand Reinhold Company: New York, 1970. (41) Sauer, B. B.; Dee, G. T. Su rface Tension and Melt Cohesive Energy Density of Polymer Melts Including High Melting and High Glass Transition Polymers. Macromolecules, 2002, 35, 7024–7030. (42) Young, T. An Essay on the Cohesion of Fluids. Phil. Trans. R. Soc. Lond. 1805, 95, 65–87.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

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 28 of 37

TOC Graphic

ACS Paragon Plus Environment

28

Page 29 of 37

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

Figure 1. Initial construction for MD simulation of wetting the Al2O3 (0001) surface. 68x103mm (150 x 150 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 2. Experimentally obtained contact angles as a function of the logarithm of the diffusion coefficient obtained from MD simulation. At right are snapshots of droplets of epoxy oligomers and a PS oligomer in contact-angle measurements. 177x80mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 30 of 37

Page 31 of 37

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

Figure 3. Cohesive energy density of (a) epoxy monomers, (b) epoxy tridecamers, and (c) PS tetramers. 79x147mm (150 x 150 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 4. Hint (yellow atoms) of a single epoxy monomer molecule. 80x47mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37

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

Figure 5. Relative concentration of Hint along z axis over 1000 ps from the second MD simulations for (a) epoxy monomer, (b) epoxy tridecamer, and (c) PS tetramer. 76x170mm (150 x 150 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 6. Radius of gyration as a function of time for epoxy tridecamer. 81x55mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37

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

Figure 7. Radius of gyration as a function of time for PS oligomers. 81x57mm (150 x 150 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 8. Radius of gyration as a function of time for epoxy monomer. 83x57mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37

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

Figure 9. Relative concentration profiles of Hall and Hint, classified by their bonding position, along z axis over 1000 ps from the second MD simulation. 175x171mm (150 x 150 DPI)

ACS Paragon Plus Environment