Properties of Dialkylcarbonate+ 1-Alkanol Mixtures at the Vacuum

The hunt for new antibiotics grows harder as resistance builds. Scientists keep fighting the microscopic arms race, even though it's one we might even...
0 downloads 0 Views 1005KB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Properties of Dialkylcarbonate + 1-Alkanol Mixtures at Vacuum Interface Mert Atilhan, and Santiago Aparicio J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b08952 • Publication Date (Web): 07 Dec 2016 Downloaded from http://pubs.acs.org on December 12, 2016

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 C 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 29

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

Properties of Dialkylcarbonate + 1-Alkanol Mixtures at Vacuum Interface Mert Atilhan*a and Santiago Aparicio*b a

Department of Chemical Engineering, Qatar University, P.O. Box 2713, Doha, Qatar b

Department of Chemistry, University of Burgos, 09001 Burgos, Spain

*

Corresponding authors:

Santiago Aparicio: e-mail [email protected] ; phone: +34 947 258062 Mert Atilhan: e-mail [email protected] ; phone: +974 4403 4142 ABSTRACT: The properties of solutions formed by dialkylcarbonate + 1-alkanol solutions are studied using classic molecular dynamics simulations. The effect of solutions composition, temperature and dialkylcarbonate type have been analysed. The structure at the interface was analysed and compared with the bulk liquid phases, with special attention to surface composition, molecular orientation and intermolecular forces. Evaporation rates were also calculated as a function of temperature. Reported results provide a detailed characterization of interfacial properties of the presented complex systems at the nanoscopic level for the first time.

1

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

INTRODUCTION The behaviour of liquid mixtures in general is of great relevance both for practical and basic science purposes.1,2 The role of intermolecular forces on the physicochemical properties3 for mixed fluids and their change with mixture composition4 has a great relevance in process design for industrial operations involving complex liquid mixtures, which are ubiquitous in the chemical industry. Extensive efforts have been carried out for developing models for predicting properties of mixed fluids.5-8 The characterization and understanding of the behaviour of complex liquid mixtures at the molecular level has been done both from experimental and theoretical approaches for many different types of molecules, although the development of structure-property relationships is still difficult for many fluids such as those hydrogen bonding molecules.9,10 Besides the requirement of understanding the behaviour of complex liquid mixtures in bulk liquid phases there is also a compulsory need of understanding the properties of multicomponent liquid mixtures interfacial behaviour. The interfacial behaviour of liquid mixtures at gas / vacuum interfaces is of great relevance for many industrial operations11 where the interface plays a pivotal role,12,13 such as those involving phase equilibria, gas adsorption or absorption and heat or mass transfer operations. The knowledge of interfacial properties, and their molecular – level roots, would allow improving the design an efficiency of these technologies. In spite of this relevance studies the behaviour of liquid mixtures at gas/vacuum interfaces are scarce in the literature, being almost limited to measurements of surface tension at air interface14 together with thermodynamic characterization,15 but with a reduced number of studies with the objective of obtaining a detailed picture of interfaces at the nanoscopic level.16-21 The detailed characterization of complex liquid – gas/vacuum interfaces can be carried out through several experimental approaches,22-29 although being difficult, or through molecular simulations.30-34 Classic molecular dynamics simulations have proven to be useful for the characterization of complex liquid – gas/vacuum interfaces,35,36 providing detailed information about molecular arrangements at the interfaces allowing to infer the most relevant differences between structuring at the interfaces and those at bulk liquid phases. In previous studies, our group has studied the behaviour of complex liquids at vacuum or gas interfaces for complex fluids such as ionic liquids or deep eutectic solvents.3639

Following our systematic approach on the study of the properties of complex liquid 2

ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

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

mixtures at gas / vacuum interfaces, a study on the properties of dialkylcarbonate + 1alkanol mixtures at vacuum interface using molecular dynamics simulations (MD) is reported in this work. Carbonates are considered as green solvents,40-42 due to their null toxicity and biodegradability, 43-49 which make these compound suitable for developing many different clean and sustainable applications.50-52 In previous studies, we demonstrated the behaviour of dialkylcarbonate with n-alkanes53 or 1-alkanols54 as cosolvents showing the changes in fluids’ structuring and intermolecular interactions with mixture composition. In the case of dialkylcarbonate + 1-alkanol solutions, the role developed by the evolution of hydrogen bonding upon changing mixtures molar fraction showed a pivotal role on the values of physicochemical properties.54 Therefore, in order to understand the properties of dialkylcarbonate + 1-alkanol mixed fluids at the vacuum interfaces, mixtures of dimethylcarbonate (DMC) + 1-butanol (BUT) were studied in the full composition range at isothermal conditions (303 K) using molecular dynamics simulations. The effect of temperature on the interface structuring were also considered in the 293 - 353 K range and the effect of the type of dialkylcarbonate were also analysed by doing simulations also for systems involving diethylcarbonate (DEC). The reported results give a detailed picture at the nanoscopic level of these dialkylcarbonate + 1-alkanol complex mixtures at the vacuum interface.

METHODS Molecular dynamics (MD) simulations were done with the MDynaMix v.5.2 molecular modelling package.55 Force field parameterizations for the considered molecules, Figure 1, were reported in a previous work.54 Simulation boxes for x DMC + (1-x) BUT mixtures were built in the full composition range (in 0.1 mole fraction steps) containing 1000 total molecules. Periodic boundary conditions were used. Initial simulation boxes were built using Packmol program56 for each mixture composition. Low-density boxes (∼ 0.2 g cm-3) were initially built followed by a NPT equilibration at 303 K and 0.1 MPa and then several heating and quenching steps up to 500 K were done. After heating-quenching cycles, the boxes were additionally equilibrated at 303 K and 0.1 MPa in the NPT ensemble for 10 ns, the equilibration was assured through the constancy of total potential energy with simulation time. Several low density initial boxes were built for each system and no differences 3

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

between the properties (density and potential energy) was inferred. Then, two vacuum slabs (200 Å long in the z-direction) were considered, one at each side of the liquid boxes (Figure S1, Supporting Information), for modelling liquid – vacuum interfaces. In this way, the interfaces are placed at xy plane with z-direction perpendicular to the interface. Simulations for the liquid – vacuum interfaces were carried out in the NVT ensemble with the temperature controlled using the Nose–Hoover thermostat and considering periodic boundary conditions. Simulations in the 293 - 353 K range were carried out. Simulations including vacuum interfaces were done for 20 ns, with the last 5 ns used for analysis. Coulombic interactions were considered with the Ewald method (15 Å cut-off radius).57 The equations of motion were solved according to the Tuckerman–Berne double time step algorithm (1 and 0.1 fs for long and short time steps)58. Lennard-Jones terms were treated according to the Lorentz-Berthelot mixing rules.

RESULTS AND DISCUSSION The analysis of interfacial properties of liquid mixtures is not a simple task, mainly raising from the difficulties for defining the exact location of the interface. Several approaches have been considered to circumvent this problem, such as the Identification of the Truly Interfacial Molecules (so-called ITIM method),20,32 although the most common approach in the literature defines the interfacial region as a slab parallel to the Gibbs dividing surface. The use of Gibbs dividing surface for defining the interface region has raised some doubts;32,36 however, it has led to a satisfactory description of the properties of many types or liquid mixtures in contact with gas/vacuum layers,16,36,59 and thus, it will be considered for the analysis along this work. In a first stage, the properties of pure solvents at vacuum interface were considered. The behavior of DMC molecules is almost equivalent for those placed in the bulk liquid phase and in DMC / vacuum interfaces, with no remarkable features at the interface, Figure 2a. On the contrary, the behavior of BUT showed remarkably more complex nature with an oscillatory behavior extending from the interface and attenuating on going to the bulk liquid phase, Figure 2b. Likewise, there is a densification at the BUT / vacuum interfaces. Up to our knowledge there are no published experimental or theoretical studies on the properties of BUT / vacuum interfaces, but there are some studies reported on methanol60 or ethanol61 / vacuum interfaces. Literature results for methanol and ethanol at vacuum interfaces did not 4

ACS Paragon Plus Environment

Page 4 of 29

Page 5 of 29

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

show enrichments at the interfaces, on the contrary density profiles for these alkanols are smooth discarding restructuring at the interface.60,61 Results for BUT / vacuum show welldefined peaks at the interface, and the increase in the length of alkyl chain in 1-alkanols leads to ordering at the vacuum interface. Likewise, results in Figure 2b show that the head group of 1-alkanols, including the hydroxyl group and the neighbor methylene group occupy the same regions at the interface close to the vacuum layer whereas the alkyl tails are placed in inner interfacial regions close to the bulk liquid phase. Remarkable changes in the structure of bulk liquid phases are obtained upon DMC and BUT mixing,54 and thus, the structure of DMC + BUT mixtures at vacuum interfaces should be remarkable different to those for pure liquids reported in Figure 2, although these results also point to prevailing effect of BUT. Atomic density profiles are reported in Figure 3a, the comparison of these results with those in Figure 2 for pure fluids shows remarkable changes upon mixing. For DMC molecules at the vacuum interface, a non-smooth oscillatory behavior at the interface is inferred. Density profiles in Figure 3a show that DMC molecules are placed in inner regions of the interface, below a layer of BUT molecules. Likewise, the enrichment for BUT molecules at the interface is in contrast with the small density peaks for DMC. Regarding the arrangement of DMC molecules, the maxima for all the considered atoms appear at almost the same distances in the vicinity of vacuum interface, thus pointing to a parallel molecular arrangement (defined according to the vector joining carbon atoms in methyl groups of DMC placed parallel to the interface). Regarding BUT molecules, they are placed forming a layer in the vicinity of vacuum, Figure 3b. Results in Figure 2b showed that tail alkyl chains were placed in inner regions of the interface for pure BUT. On the other hand results in Figure 3a show that these alkyl chains are placed close to the vacuum layer for DMC + BUT mixtures, this rearrangement seem to be forced for placing DMC molecules in the inner layer immediately below the BUT layer. Moreover, density profiles for BUT in Figure 3a show that BUT molecules are skewed at the vacuum interface, with the vector joining the BUT head oxygen atom with the tail carbon at the methyl group tending to arrange perpendicular to the surface. The results in Figure 3c show that DMC molecules form islands separated by continuous regions of BUT molecules at vacuum interface, thus confirming the prevailing role of BUT molecules for determining the interfacial properties in DMC + BUT mixtures.

5

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

Although density profiles reported in Figure 3 allow to infer the molecular orientations at the interface, this was quantified though the angle distribution as a function of mixture mole fraction reported in Figure 4. The arrangements of DMC molecules are clearly different to the ones for BUT molecules. The angle distribution for DMC molecules show a complex behavior with increasing DMC mole fraction in the mixtures, increasing mole fractions leads to a widening of the distribution showing that DMC molecules may develop different arrangements in the bulk liquid phase. At the interface regions, although results in Figure 3a showed a trend for developing parallel arrangements, distribution reported in Figure 4 show that although this arrangement is highly probable, other arrangements are also possible, confirming that the behavior of DMC molecules at the interface is similar to that in bulk regions. On the contrary, angle distributions for BUT at the vacuum interface show clearly the trend for developing skewed arrangements, as the hotspots in density distributions (for φ ∼ 150o) in the full composition range show. The temperature effect on molecular arrangements at the vacuum interface is reported in Figure S2 (Supporting Information) showing that both for DMC and BUT molecules the densification at the interface vanishes with increasing temperature, leading to interfaces with vacuum with properties equivalent to those in bulk regions. This is quantified by the comparison of mole fraction calculated at the interface region with that in bulk liquid, Figure 5. Two different values were used to define interface thickness around Gibbs dividing surface (r = 5 and 10 Å) and in both cases the concentration of DMC molecules at the interfacial region is lower than in bulk regions, in agreement with density profiles reported in Figure 3. Nevertheless, the DMC mole fraction approaches to the bulk value in a linear way with increasing temperature, showing that this heterogeneity almost vanishes for temperatures larger than 353 K. The behavior of molecules at vacuum interfaces has a relevant role on the kinetics of molecular evaporation, evaporation rates, and thus on the development of vapor - liquid equilibria, which can be analyzed at the molecular level from MD simulations reported in this work. Molecules are defined as belonging to the vapor phase when being placed in a region 10 Å above the Gibbs dividing surface, in agreement with density distributions reported in Figures 3 and S2 (Supporting Information). The results in Figures 6a and 6b show that BUT molecules are evaporated more quickly than DMC molecules, the number of BUT molecules in the vapor phase is roughly four times larger than of DMC molecules. The reason for this 6

ACS Paragon Plus Environment

Page 6 of 29

Page 7 of 29

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

behavior can be the prevailing presence of BUT molecules at the interface region reported in Figures 3 and 5, with DMC molecules placed in inner regions, thus hindering the evaporation of DMC. The number of molecules in the vapor phase follows an almost linear behavior with simulation time for all the studied temperatures and for both molecules in DMC + BUT mixtures, which allow to calculate the evaporation rates as a function of temperature. The temperature effect on the evaporation rate reported in Figure 6c follows a clearly non-linear behavior, especially for BUT molecules. The preferential evaporation of BUT molecules over DMC ones is clearly remarkable for large temperatures. The differences in intermolecular forces for those molecules placed in bulk liquids and those at interfaces are analyzed as a function of temperature and mixture composition in Figure 7. It should be remarked that all the forces are weaker for molecules at the interfaces, as the positive values of ΔE (values at the interface minus those at the bulk) show for any temperature and mixture composition. The bulk-surface differences for energies involving homoassociations (DMC-DMC and BUT-BUT) are lower than those for heteroassociations (DMC-BUT), which is justified considering the interface structuring involving a top layer rich in BUT and a lower layer rich in DMC. This interfacial arrangement allows the efficient development of homoassociations but hinders heteroassociations. Likewise, DMC – DMC interactions suffer less changes than BUT – BUT ones, which may be justified considering and end of aisle effect for BUT molecules (which do not have molecules for interacting above them at the interface). The evolution of ΔE with temperature reported in Figure 7 shows that although it increases linearly with temperature, the changes are almost negligible in the studied temperature range. The changes of ΔE with mixture composition are also almost negligible in the full composition range with the exception of data for DMC – BUT, which suddenly increases for xDMC > 0.5 pointing to weak DMC – BUT interactions at the interface for BUT – rich mixtures. As the regions closer to vacuum layer are richer in BUT whereas DMC molecules tend to occupy inner regions, only for DMC-rich mixtures BUT-DMC interactions at the interface start to be relevant which would justify the behavior reported in Figure 7b. The differences in intermolecular forces at the interface in comparison with bulk regions should lead also to different behavior for the dynamics of both regions, which is analyzed through the self – diffusion coefficients calculated from mean square displacements and reported in Figure 8. The calculation of self-diffusion coefficients for molecules at the interfaces was carried out considering molecules belonging to the 7

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

interface as those with the z-coordinate of their center of mass standing in the zGDS ± 10 Å region. For these molecules, the self-diffusion coefficients were calculated from the mean square displacement. The results for self-diffusion coefficients, D, for equimolar mixtures as a function of temperature are reported in Figure 8a comparing values in bulk and interface regions. The ability of the proposed forcefield parameterization for predicting self-diffusion coefficients of neat dialkylcarbonates in bulk liquid phases was previously studied,53 showing reasonable agreement between predicted data and experiments. Nevertheless, there are not experimental data for comparison for the self-diffusion coefficients in dialkylcarbonate + alkanol mixtures neither in bulk nor interfacial regions. The behavior of D for molecules in the bulk liquid phase follows a close to linear behavior with increasing temperature (R2 = 0.99 both for BUT and DMC). Likewise, D are larger for DMC than for BUT molecules in the bulk regions, as it may be expected from the larger intermolecular interaction energies for BUT-BUT than for DMC-DMC. The behavior of D at vacuum interfaces is complex. First, D are lower for BUT than for DMC molecules at the interface as in bulk regions. Second, the comparison of D values for interface and bulk regions show that at low temperatures (T < 303 K) DMC and BUT molecules diffuse faster at the interface than in the bulk liquid, which would be in agreement with the weakening in intermolecular forces reported in Figure 7. For T > 303 K, D values are lower at the interface than in the bulk regions. This can be justified considering that although intermolecular interaction energies are also weaker at the interface than in the bulk for these temperatures, with increased temperature the mobility induced by thermal motion is larger in the bulk than at the interface, in which, although mobility is also improved with increasing temperature, the presence of a vacuum layer above and the difficulties for leaving the interface region for moving to lower layers, because of the different molecular structuring, hinders molecular mobility in the direction perpendicular to the interface. Third, the effect of mixture composition on D is reported in Figure 8b showing an almost parallel evolution with mole fraction for molecules in the bulk and interface region. The molecular dynamics simulation allowed the calculation of surface tension, σ, as a function of pressure and temperature, which was obtained from the corresponding pressure tensors, Figure S3 (Supporting Information).62 The procedure for calculating surface tension from molecular dynamics simulations was validated by comparison between the simulated values and those experimental data from the literature (BUT and DMC; 19.5 and 29.9 mN m8

ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

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

1

, respectively, from MD, and 23.6363 and 27.764 and mN m-1, respectively, experimental

from literature, all at 303 K). Calculated surface tension for DMC + BUT mixtures decreases in a linear trend with increasing temperature, whereas the variation with mixture composition is clearly non-linear with negative deviations from linearity, Δσ, reaching a minimum at equimolar composition. Negative Δσ rise from the enrichment in BUT molecules at the interfacial region in comparison with bulk regions, and thus, leading to lower surface tension values as may be inferred from values for pure fluids. This is confirmed by the density profiles reported in Figure 9a as a function of mixture composition in which the enrichment in BUT molecules at the interface region is clearly inferred in the full composition range, even for highly diluted mixtures, thus decreasing surface tension in the mixtures. On the contrary, for DMC molecules the behavior of interface regions is very similar to the bulk phases, Figure 9b, and thus the properties of DMC + BUT / vacuum interfaces are heavily dependent on the trend of BUT molecules for occupying regions in the vicinity of the vacuum layer. The analysis of the mixture composition at the interface shows the enrichment in BUT molecules in the whole composition range, Figure 10, following a clearly non-linear trend with larger deviations at equimolar mixtures, which is in agreement with the behavior of surface tension reported in Figure S3 (Supporting Information). This enrichment in BUT molecules extends for a large region around the Gibbs dividing surface, even when considering a region 10 Å wide around the Gibbs dividing surface the compositions are lower than in the bulk region, although the effects are more remarkable for an interface region of ± 5 Å thickness. Regarding the charge distribution at the interface, results in Figure S4 (Supporting Information) show a region slightly positively charged above the Gibbs dividing surface and negatively charged below it, which rise from the atomic distribution of BUT molecules in that region reported in Figures 2 and 3; nevertheless, these values are very low and additional features cannot be inferred from them. The mixture composition effect on evaporation rates is reported in Figure 11 in which the preferential evaporation of BUT molecules is maintained in the full composition range. Both for BUT and DMC molecules evaporation rates decrease with increasing amount of the cosolvent in the mixtures but only for mixtures with composition larger than 0.8 DMC mole fraction rates for BUT are lower than for DMC, in agreement with results in Figures 9 and 10. In the last stage of the research, the effect of increasing alkyl chains in dialkylcarbonate were analyzed. It should be remarked that pure DMC is denser than pure 9

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

Page 10 of 29

DEC,54 and thus DEC-containing mixtures are less dense than those with DMC, which is showed in the number density profiles reported in Figure S5 (Supporting Information). Likewise, the results in Figure S5 (Supporting Information) shows that on going from DMC to diethylcarbonate (DEC) + BUT mixtures the structuring of dialkylcarbonate at the interface is even more similar to the bulk regions: the small peaks appearing for DMC at the interface vanish for DEC. Regarding the behavior of BUT molecules, it is not affected by the longer alkyl chains in DEC in comparison with DMC-containing mixtures, the density peaks at the vacuum interface is maintained together with the trend of BUT molecules to occupy regions close to the vacuum layer and of dialkylcarbonate to stay in inner regions close to the bulk liquid phase. The most remarkable effect of the presence of DEC molecules is the decrease in evaporation rates when compared with mixtures with DMC, both for BUT and the especially for the dialkylcarbonate, Figure 12. This may be justified considering that increasing the number of methyl groups in the dialkylcarbonate increases the number of Lennard-Jones contacts

and

thus

increasing

the

intermolecular

interaction

energy

between

dialkylcarbonate molecules (DEC-DEC is 33 % larger than DMC-DMC for dialkylcarbonate molecules at the interface of dialkylcarbonate + BUT equimolar mixtures at 303 K) and also increasing dialkylcarbonate – BUT interaction energies (7 %), which results in large decrease of evaporation rates for DEC when compared with DMC and lower decrease for BUT on going from DMC to DEC molecules. Therefore, the increase in the dialkylcarbonate does not lead to remarkable changes in the structuring at the vacuum interfaces beyond the increase of van der Waals forces but without molecular reorganization.

CONCLUSIONS The structure of dialkylcarbonate + 1-alkanol / vacuum interfaces has been analyzed through molecular dynamics simulations as a function of mixture composition and temperature. The reported results show structuring at the interface completely different to the bulk liquid regions, with a trend of 1-alkanol molecules to develop a highly dense layer in the vicinity of the vacuum region confining dialkylcarbonate molecules to inner regions closer to the bulk fluid. This molecular arrangement at the interface is maintained in the full composition range and in the wide temperature range studied and is the molecular reason of the nonlinear behavior with mixture composition of relevant surface properties such as surface tension. Likewise, the heterogeneity at the interface determines the mechanism and kinetics 10

ACS Paragon Plus Environment

Page 11 of 29

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

of evaporation although the thermodynamics of the process was not characterized in the present work. The reported results show that large effects can be produced on liquid mixtures properties by the development of a vacuum interface, and it will be extended in future works to other families of molecules for rationalization purposes.

ACKNOWLEDGEMENT This work was funded by Ministerio de Economía y Competitividad (Spain, project CTQ201340476-R) and Junta de Castilla y León (Spain, project BU324U14). The Foundation of Supercomputing Center of Castile and León (FCSCL, Spain) for providing supercomputing facilities. The statements made herein are solely the responsibility of the authors.

SUPPORTING INFORMATION Figure S1 (model simulation box); Figure S2 (temperature effect on number density profiles); Figure S3 (surface tension); Figure S4 (charge density profiles); Figure S5 (comparison of number density profiles for DMC and DEC containing mixtures).

11

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 Captions.

Figure 1. Molecules and atom labelling used along this work.

Figure 2. Number density profiles of relevant atoms in (a) DMC and (b) BUT, in DMC / vacuum or BUT / vacuum interfaces at 303 K. z stands for the coordinate perpendicular to the vacuum interface, and zGDS for the coordinate of the Gibbs dividing surface corresponding to the top layer in contact with vacuum. Atoms labeling as in Figure 1. Figure 3. (a) Number density profiles for relevant atoms in x DMC + (1-x) BUT (x = 0.5) / vacuum interface at 303 K. Panels (b) and (c) show snapshots perpendicular and on the interface with vacuum, respectively. z stands for the coordinate perpendicular to the corresponding interfaces, and zGDS for the coordinate of the Gibbs dividing surface. Density profiles in panel a are obtained as averages in the 4 to 5 ns simulation range. Atoms labeling in panel a as in Figure 1. Color code in panels b and c: (blue) DMC and (red) BUT. Continuous lines in panels b and c show limits of periodic boundary conditions used in the simulations; vertical dashed lines in panel b show relevant positions. Figure 4. Two-dimensional density distribution of angles formed between a vector perpendicular to the interface and vector ܸ1 (vector joining C2 atoms in DMC, Figure 1) or vector ܸ2 (vector joining O3 and C4, methyl, atoms in BUT, Figure 1), for x DMC + (1-x) BUT / vacuum interface at 303 K. z stands for the coordinate perpendicular to the corresponding interfaces, and zGDS for the coordinate of the Gibbs dividing surface (for the top interface). Figure 5. Comparison between DMC mole fraction in the bulk liquid phase, xbulk, and in the interfacial region defined as zGDS ± r, xinterface, in x DMC + (1-x) BUT (xbulk = 0.5) as a function of temperature, T. z stands for the coordinate perpendicular to the corresponding interfaces, and zGDS for the coordinate of the Gibbs dividing surface. xbulk is showed as a dashed line.

Figure 6. (a,b) Number of molecules evaporated, N, for x DMC + (1-x) BUT (xbulk = 0.5) / vacuum interface as a function of temperature. Linear fits are showed for guiding purposes as dashed lines. Panel c shows the variation of the slopes of linear fits, (dN / dt), with temperature. Evaporated molecules are defined as those with z > (zGDS + 10 Å), where zGDS stands for the coordinate of the Gibbs dividing surface. Linear fits are showed for guiding purposes as dashed lines.

Figure 7. Differences between intermolecular interaction energies, ΔE, for molecules at the vacuum interface and those in the bulk liquid phases for (a) x DMC + (1-x) BUT (x = 0.5) as a function of temperature or (b) x DMC + (1-x) BUT at 303 K as a function of mole fraction. Interface region defined as zGDS ± r, with r = 10 Å, where zGDS stands for the coordinate of the Gibbs dividing surface.

12

ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29

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. Self-diffusion coefficients, D, for molecules at the vacuum interface and those in the bulk liquid phases for (a) x DMC + (1-x) BUT (x = 0.5) as a function of temperature or (b) x DMC + (1-x) BUT at 303 K as a function of mole fraction. Interface region defined as zGDS ± r, with r = 10 Å, where zGDS stands for the coordinate of the Gibbs dividing surface.

Figure 9. Mixture composition effect on number density profiles of (a) O3 atom in BUT and (b) C1 atom in DMC, in x DMC + (1-x) BUT / vacuum interface. z stands for the coordinate perpendicular to the corresponding interfaces, and zGDS for the coordinate of the Gibbs dividing surface. Density profiles are obtained as averages in the 4 to 5 ns simulation range. Atoms labeling as in Figure 1.

Figure 10. Comparison between DMC mole fraction in the bulk liquid phase, xbulk, and in the interfacial region defined as zGDS ± r, xinterface, in x DMC + (1-x) BUT at 303 K (x = xbulk). z stands for the coordinate perpendicular to the corresponding interfaces, and zGDS for the coordinate of the Gibbs dividing surface.

Figure 11. Number of molecules evaporated, N, for x DMC + (1-x) BUT / vacuum interface at 303 K. Linear fits are showed for guiding purposes as dashed lines. Panel c shows the variation of the slopes of linear fits, (dN / dt), with temperature. Evaporated molecules are defined as those with z > (zGDS + 10 Å), where zGDS stands for the coordinate of the Gibbs dividing surface. Linear fits are showed for guiding purposes as dashed lines.

Figure 12. Number of molecules evaporated, N, for x DMC + (1-x) BUT / vacuum (x = 0.5) and x DEC + (1-x) BUT / vacuum (x = 0.5) interfaces at 303 K. Linear fits are showed for guiding purposes as dashed lines. Numbers -1

inside the Figure show the variation of the slopes of linear fits, (dN / dt, ns ), with temperature. Evaporated molecules are defined as those with z > (zGDS + 10 Å), where zGDS stands for the coordinate of the Gibbs dividing surface. Linear fits are showed for guiding purposes as dashed lines.

13

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

O2

C3,H2

O

H1,C2

C2,H1

H3C

CH3

O C1 O O1

R

HO Hb,O3

O1 DMC

R = – (CH2)n – CH3

O2

O H1,C3

H3C

C4,H2

C2,H1

O C1 O O1

C2,H1 C3,H1

CH3

O1

n = 2 (1-butanol)

DEC

Figure 1.

14

ACS Paragon Plus Environment

BUT

Page 14 of 29

Page 15 of 29

0.008

(a) DMC

0.008

C2 O1 C1 O2

0.004

0.002

0

0 -50 -25 0 z - zGDS / Å

C4 C3 O3 HB

0.004

0.002

-75

(b) BUT

0.006 ρ / Å-3

0.006 ρ / Å-3

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

25

-75

-50 -25 0 z - zGDS / Å

Figure 2.

15

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

Figure 3.

16

ACS Paragon Plus Environment

Page 16 of 29

Page 17 of 29

DMC (

)

x = 0.0

x = 0.1

φ / deg

x = 0.3

x = 0.5

x = 0.7

x = 0.9

x = 1.0

180

180

180

180

180

180

150

150

150

150

150

150

120

120

120

120

120

120

90

90

90

90

90

90

60

60

60

60

60

60

30

30

30

30

30

30

0 -50 -40 -30 -20 -10 0 10

180

180

0 -50 -40 -30 -20 -10 0 10

180

0 -50 -40 -30 -20 -10 0 10

180

0 -50 -40 -30 -20 -10 0 10

180

0 -50 -40 -30 -20 -10 0 10

180

0 -50 -40 -30 -20 -10 0 10

1

high 0.9

150

150

150

150

150

150

120

120

120

120

120

120

90

90

90

90

90

90

60

60

60

60

60

60

30

30

30

30

30

30

60

0.8

)

50

BUT (

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

0.7 40

0.6 0.5

30 0.4 0.3

20

0.2 10 0.1

0 -50 -40 -30 -20 -10 0 10

0 -50 -40 -30 -20 -10 0 10

0 -50 -40 -30 -20 -10 0 10

0 -50 -40 -30 -20 -10 0 10

0 -50 -40 -30 -20 -10 0 10

z-zGDS / Å

Figure 4.

17

ACS Paragon Plus Environment

0 -50 -40 -30 -20 -10 0 10

low

0 0

0.5

1

The Journal of Physical Chemistry

0.6

r=5 Å r=10 Å

0.5 xinterface

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 29

0.4

0.3 280

300

320 T/K

340

Figure 5.

18

ACS Paragon Plus Environment

360

Page 19 of 29

250

(a) BUT

353 K 343 K 333 K 323 K 313 K 303 K 293 K

200 N

150 100 50 0 50

(b) DMC

40 N

30 20 10 0 0

1

2

3

t / ns

80 (dN / dt) / ns-1

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

BUT DMC

60 40 20

(c)

0 280

300

320 T/K

340

Figure 6.

19

ACS Paragon Plus Environment

360

The Journal of Physical Chemistry

5

16

DMC-DMC DMC-BUT BUT-BUT

(a)

(b)

12 ∆E / kJ mol-1

4 ∆E / kJ mol-1

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

3 2

8 4

1

0

0

-4 280

300

Page 20 of 29

320 T/K

340

360

0

0.2

0.4

0.6 x

Figure 7.

20

ACS Paragon Plus Environment

0.8

1

Page 21 of 29

3

1.2 DMC @ bulk DMC @ interface BUT @ bulk BUT @ interface

(a)

(b)

1 109 D / m2 s-1

2.5 109 D / m2 s-1

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

2 1.5

0.8

0.6

1 0.5

0.4 280

300

320 T/K

340

360

0

0.2

0.4

0.6 x

Figure 8.

21

ACS Paragon Plus Environment

0.8

1

The Journal of Physical Chemistry

0.008

(a) O3 (BUT)

ρ / Å-3

0.006 0.004

x=0.0 x=0.1 x=0.3 x=0.5 x=0.7 x=0.9 x=1.0

0.002 0 0.008

(b) C1 (DMC)

0.006 ρ / Å-3

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 29

0.004 0.002 0 -20

-10 0 z - zGDS / Å

Figure 9.

22

ACS Paragon Plus Environment

10

Page 23 of 29

1

r=5 Å r=10 Å

0.8

xinterface

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

0.6 0.4 0.2 0 0

0.2

0.4

0.6

0.8

xbulk

Figure 10.

23

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

40

(a) BUT

N

30 20 10 0 40

(b) DMC

x=0.0 x=0.1 x=0.3 x=0.5 x=0.7 x=0.9 x=1.0

N

30 20 10 0 0

1

2

3

t / ns

16 (dN / dt) / ns-1

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 29

BUT DMC

12 8 4

(c)

0 0

0.2

0.4

0.6

0.8

x Figure 11.

24

ACS Paragon Plus Environment

1

Page 25 of 29

20

DMC in DMC + BUT DEC in DEC + BUT BUT in DMC + BUT BUT in DEC + BUT

15 N

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

4.83 4.04

10

2.22

5 0.48

0 0

1

2 t / ns

Figure 12.

25

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

REFERENCES (1) Rowlinson, J. S.; Swinton, F. L. Liquids and Liquid Mixtures. Ed. Butterworth Scientific. London, 1982. (2) Gmehling, J.; Kolbe, B.; Kleiber, M.; Rarey, J. Chemical Thermodynamics for Processs Simulation. Ed. WileyVCH, Weinheim, 2012. (3) Nabi, F.; Malik, M. A.; Jesudason, C. G.; Al-Thabaiti, S. A. A Review of Molecular Interactions in Organic Binary Mixtures. Korean J. Chem. Eng. 2014, 31, 1505-1517. (4) Stone, A. The Theory of Intermolecular Forces. Ed. Oxford University Press, Oxford, 2013. (5) Kontogeorgis, G. M.; Folas, G. K. Thermodynamic Models for Industrial Applications. From Classical and Advanced Mixing Rules to Association Theories. Ed.John Wiley and Sons, New York, 2010. (6) Bell, I. H.; Wronski, J.; Quoilin, S., Lemort, V. Pure and Pseudo-Pure Fluid Thermophysical Property Evaluation and the Open-Source Thermophysical Property Library CoolProp. Ind. Eng. Chem. Res. 2014, 53, 2498-2508. (7) Xiong, R.; Sandler, S. I.; Burnett, R. I. An Improvement to COSMO-SAC for Predicting Thermodynamic Properties. Ind. Eng. Chem. Res. 2014, 53, 8265-8278. (8) Dufal, S.; Papaioammou, V.; Sadeqzadeh, M.; Pogiatzis, T., Chremos, A.; Adjiman, C. S., Jackson, G., Galindo, A. Precition of Thermodynamic Properties and Phase Behavior of Fluids and Mixtures with the SAFT-γ Mie Group-Contribution Equaiton of State. J. Chem. Eng. Data 2014, 59, 3272-3288. (9) De Villiers, A. J.; Schwarz, C. E.; Burger, A. J.; Kontogeorgis, G. M. Evalution of the PC-SAFT, SAFT and CPA Equations of State in Predicting Derivative Properties of Selected Non-polar and Hydrogen-bonding Compounds. Fluid Phase Equilibr. 2013, 338, 1-15. (10) Fouad, W. A.; Wang, L.; Haghm,oradi, A.; Asthagari, D.; Chapman, W. G. Understanding the Thermodynamics of Hydrogen Bonding in Alcohol-Containing Mixtures: Cross-Association. J. Phys. Chem. B 2016, 120, 3388-3402. (11) Als-Nielsen, J.; Jacquemain, D.; Kjaer, K.; Leveiller, F.; Lahav, M.; Leiserowitz, L. Principles and Applications of Grazing Incidence X-Ray and Neutron Scattering from Ordered Molecular Monolayers at the Air-Water Interface. Phys. Rep. 1994, 246, 251-313. (12) Martinez, I. S., Santos, C.; Baldelli, S. Structural Study at the Gas-Liquid Interface of 1-Alkyl-3methylimidazolium Alkylsulfates Using Surface Potential Measurements. ChemPhysChem 2012, 13, 18181824. (13) Oliveira, M. B.; Dominguez-Perez, M.; Freire, M. G.; Llovell, F.; Cabeza, O.; Lopes-da-Silva, J. A.; Vega, L. F.; Coutinho, J. A. P. Surface Tension of Binary Mixtures Of 1-Alkyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide Ionic Liquids: Experimental Measurements and Soft-SAFT Modeling. J. Phys. Chem. B 2012, 116, 12133-12143. (14) Abbasi, A.; Gitifar, V., Setoodeh, P. QSPR Strategy to Model and Analyze Surface Tension of Binary – Liquid Mixtures: A Large – Data – Set Case. Chemometr. Intell. Lab. 2016, 155, 36-45. (15) Soledade, M.; Santos, C. S., Reis, J. C. R. Activity Coefficients in the Surface Phase of Liquid Mixtures. ChemPhysChem 2015, 16, 470-475. (16) Ishihara, T.; Ishiyama, T.; Morita, A. Surface Structure of Methanol / Water Solutions vis Sum Frequency Orientational Analysis and Molecular Dynamics Simulation. J. Phys. Chem. C 2015, 119, 9879-9889. (17) Dominguez, H. The Non-Ideal Behaviour of the Interfacial Tension of the n-Heptane + Perfluoro-n-hexane Mixture: A Computational Study. Chem. Phys. Lett. 2015, 627, 77-81. (18) Bermudez-Salguero, C.; Garcia-Fadrique, J. Gibbs Excess and the Calculation of the Absolute Surface Composition of Liquid Binary Mixtures. J. Phys. Chem. B 2015, 119, 5598-5608. (19) Chilukoti, H. K., Kikugawa, G.; Ohara, T. Self-diffusion Coefficient and Structure of Binary n-Alkane Mixtures at the Liquid – Vapor Interfaces. J. Phys. Chem. B 2015, 119, 13177-13184. (20) Fábián, B., Jójárt, B.; Horvai, G.; Jedlovszky, P. Properties of the Liquid-Vapor Interface of Acetone-Water Mixtures. A Computer Simulation and ITIM Analysis Study. J. Phys. Chem. C 2015, 119, 12473-12487. (21) Goufi, A.; Malfreyt, P.; Tildesley, D. J. Computer Modelling of the Surface Tension of the Gas-Liquid and Liquid-Liquid Interface. Chem. Soc. Rev. 2016, 45, 1387-1409. (22) Eisenthal, K. B. Liquid Interfaces Probed by Second-Harmonic and Sum-Frequency Spectroscopy. Chem. Rev. 1996, 96, 1343-1360. (23) Miranda, P. B.; Shen, Y. R. Liquid Interfaces:  A Study by Sum-Frequency Vibrational Spectroscopy. J. Phys. Chem. B 1999, 103, 3292-3307. (24) Chang, T. M.; Dang, L. X. Liquid-Vapor Interface of Methanol-Water Mixtures: A Molecular Dynamics Study. J. Phys. Chem. B 2005, 109, 5759-5765. 26

ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

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

(25) Jubb, A. M.; Hua, W.; Allen, H. C. Environmental Chemistry at Vapor / Water Interfaces: Insights from Vibrational Sum Frequency Generation Spectroscopy. Ann. Rev. Phys. Chem. 2012, 63, 107-130. (26) Geissler, P. L. Water Interfaces, Solvation, and Spectroscopy. Ann. Rev. Phys. Chem. 2013, 64, 317-337. (27) Nagata, Y.; Hsieh, C. S.; Hasegawa, T., Voll, J.; Backus, E. H. G.; Bonn, M. Water Bending Mode at the Water-Vapor Interface Probed by Sum-Frequency Generation Spectroscopy: A Combines Molecular Simulation and Experimental Study. J. Phys. Chem. Lett. 2013, 4, 1872-1877. (28) Tesa-Serrate, M. A.; Marshall, B. C.; Smoll, E. J.; Purcell, S. M., Costen, M. L.; Slattery, J. M., Minton, T. K.; McKendrick, K. G. Ionic Liquid – Vacuum Interfaces probed by Reactive Atom Scattering: Influence of Alkyl Chain Length and Anion Volume. J. Phys. Chem. C, 2015, 119, 5491–5505. (29) Ni, Y., Skinner, J. L. IR and SFG Vibrational Spectroscopy of the Water Bend in the Bulk Liquid and at the Liquid-Vapor Interface, Respectively. J. Chem. Phys.2015, 143, 014502. (30) Kessler, J.; Elgabarty, H., Spura, T.; Karham, K.; Partovi-Azar, P.; Hassanali, A. A.; Kühne, T. D. Structure and Dynamics of the Instantaneous Water / Vapor Interface Revisited by Path – Integral and Ab-Initio Molecular Dynamics Simulations. J. Phys. Chem. B, 2015, 119, 10079–10086 (31) Delaunay, F.; Rodríguez, A. S.; Couvert, A.; Amrane, A.; Biard, P. F.; Szymczyk, A.; Malfreyt, P.; Ghoufi, A. Interfacial Structure of Toluene at an Ionic Liquid / Vapor Interface: A Molecular Dynamics Simulation Investigation. J. Phys. Chem. C 2015, 119, 9966-9972. (32) Idrissi, A.; Hantal, G.; Properties of the Liquid–Vapor Interface of Acetone–Methanol Mixtures, as Seen From Computer Simulation and ITIM Surface Analysis. Phys. Chem. Chem. Phys. 2015, 17, 8913-8926. (33) Li, S.; Zhao, W.; Feng, G.; Cummings, P. T. A Computational Study of Dicationic Ionic Liquids/CO2 Interfaces Langmuir 2015, 31, 2447– 2454. (34) Konieczny, J.; Szefzyk, B. Structure of Alkylimidazolium-Based Ionic Liquids at the Interface with Vacuum and Water – A Molecular Dynamics Study. J. Phys. Chem. B 2015, 119, 3795–3807. (35) Paredes, X.; Fernández, J.; Padua, A. H.; Malfrey, P.; Malberg, F.; Kirchner, B.; Pensado, A. S. Bulk and Liquid-Vapor Interface of Pyrrolidinium-Bases Ionic Liquids: A Molecular Simulation Study. J. Phys. Chem. B 2014, 118, 731–742. (36) García, G., Atilhan, M.; Aparicio, S. Interfacial Properties of Double Salt Ionic Liquids: A Molecular Dynamics Study. J. Phys. Chem. C 2015, 119, 28405–28416. (37) Aparicio, S.; Atilhan, M. On the Properties of CO2 and Flue Gas at the Piperazinium – Based Ionic Liquids Interface: A Molecular Dynamics Study. J. Phys. Chem. C 2013, 117, 15061–15074. (38) García, G.; Atilhan, M.; Aparicio, S. Interfacial Properties of Deep Eutectic Solvents Regarding to CO2 Capture. J. Phys. Chem. C, 2015, 119, 21413–21425. (39) Herrera, C.; García, G.; Alcalde, R.; Atilhan, M.; Aparicio, S. Interfacial Properties of 1-Ethyl-3methylimidazolium Glycinate Ionic Liquid Regarding CO2, SO2 and Water from Molecular Dynamics. J. Mol. Liq. 2016, 220, 910-917. (40) Tundo, P.; Arico, F.; Rosamilia, A. E.; Rossi, G. L. Dimethyl Carbonate: Green Solvent and Ambident Reagent. In: Tundo, P.; Esposito, V. (Eds.) Green Chemical Reactions, Ed. Springer, 2008; pp. 213-232. (41) Parker, H. L.; Sherwood, J.; Hunt, A.J.; Clark, J. H. Cyclic Carbonates as Green Alternative Solvents for the Heck Reaction. ACS Sus. Chem. Eng. 2014, 2, 1739-1742. (42) Bellardita, M.; Loddo, V.; Mele, A.; Panzeri, W., Parrino, F.; Pibiri, I.; Palmisano, L. Photocatalysis in dimethyl carbonate green solvent: degradation and partial oxidation of phenanthrene on supported TiO2. RSC Adv. 2014, 4, 40859-40864. (43) Brown, D.; Gaunt, I. F.; Kiss, I. S.; Butterworth, K. R. Long-term Toxicity of Diethyl Carbonate in Mice. Toxicol. 1978, 10, 291-295. (44) Tundo, P.; New Developments in Dimethyl Carbonate Chemistry. Pure Appl. Chem. 2001, 73, 1117-1124. (45) Bilde, M.; Mogelberg, T. E.; Sehested, J.; Nielsen O. J.; Wallington, T. J.; Hurley, M. D.; Japar, S. M.; Dill, M.; Orkin, V. L.; Buckley, T. J. et al. Atmospheric Chemistry of Dimethyl Carbonate:  Reaction with OH Radicals, UV Spectra of CH3OC(O)OCH2 and CH3OC(O)OCH2O2 Radicals, Reactions of CH3OC(O)OCH2O2 with NO and NO2, and Fate of CH3OC(O)OCH2O Radicals. J. Phys. Chem. A 1997, 101, 3514-3525. (46) Tundo, P.; Selva, M. The Chemistry of Dimethyl Carbonate. Acc. Chem. Res. 2002, 35, 706-716. (47) Righi, G.; Bovicelli, P.; Barontini, P.; Tirotta, I. Dimethyl Carbonate in the Regio- and Stereocontrolled Opening of Three-membered Heterocyclic Rings. Green Chem. 2012, 14, 495-502. (48) Arico, F.; Tundo, P. Dimethyl Carbonate: a Modern Green Reagent and Solvent. Russ. Chem. Rev. 2010, 79, 479-489. (49) Tundo, P.; Selva, M. The Chemistry of Dimethyl Carbonate. Acc. Chem. Res. 2002, 35, 706-716.

27

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

(50) Arteconi, A.; Mazzarini, A.; Di Nicola, G. Emissions from Ethers and Organic Carbonate Fuel Additives: A Review. Water Air Soil Pollut. 2011, 221, 405-423. (51) Schaffner, B.; Scha ̈ ffner, F.; Verevkin, S. P.; Borner, A. Organic Carbonates as Solvents in Synthesis and Catalysis. Chem. Rev. 2010, 110, 4554-4581. (52) Santos, B. A. V.; Silva, V. M. T.; Loureiro, J. M.; Rodrigues, A. E. Review of the Direct Synthesis of Dimethyl Carbonate. ChemBioEng Rev. 2014, 1, 214-229. (53) García, G.; Trenzado, J. L.; Alcalde, R., Rodríguez-Delgado, A.; Atilhan, M.; Aparicio, S. Structure of Alkylcarbonate + n-Alkane Mixed Fluids. J. Phys. Chem. B 2014, 118, 11310-11322. (54) Alcalde, R.; Atilhan, M.; Trenzado, J. L.; Aparicio, S. Physicochemical Insights on Alkylcarbonate – Alkanol Solutions. J. Phys. Chem. B 2016, 120, 5015-5028. (55) Lyubartsev, A. P.; Laaksonen, A. MDynaMix - A Scalable Portable Parallel MD Simulation Package for Arbitrary Molecular Mixtures. Comput. Phys. Commun. 2000, 128, 565-589. (56) Martínez, L.; Andrade, R.; Birgin, E. G.; Martínez, J. M. Packmol: A Package for Building Initial Configurations for Molecular Dynamics Simulations. Journal of Computational Chemistry, 2009, 30, 21572164. (57) Essmann, U. L.; Perera, M. L.; Berkowitz, T.; Darden, H.; Lee,H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577-8593. (58) Tuckerman, M.; Berne, B. J.; Martyna, G. J. Reversible Multiple Time Scale Molecular Dynamics. J. Chem. Phys. 1992, 97, 1990-2001. (59) Schöttl, S.; Touraud, D.; Kunz, W.; Zemb, T.; Horinek, D. Consistent Definitions of “the Interface” in Surfactant-Free Micellar Aggregates. Colloid. Surface A 2015, 480, 222-227. (60) Dang, L. X.; Chang, T. M. Many-body Interactions in Liquid Methanol and Its Liquid / Vapor Interface: A Molecular Dynamics Study. J. Chem. Phys. 2003, 119, 9851-9857. (61) Taylor, R.; Shields, R. L. Molecular-dynamics Simulations of the Ethanol Liquid – Vapor Interface. J. Chem. Phys. 2003, 119, 12569 - 12576. (62) Alejandre, J.; Tildesley, D. J.; Chapela, G. A. Molecular Dynamics Simulation of the Orthobaric Densities and Surface Tension of Water. J. Chem. Phys. 1995, 102, 4574−4583. (63) Jiménez, E.; Cabanas, M.; Segade, L.; García-Garabal, S.; Casas, H. Excess Volume, Changes of Refractive Index and Surface Tension of Binary 1,2-Ethanediol + 1-Propanol or 1-Butanol Mixtures at Several Temperatures. Fluid Phase Equilibr. 2001, 180, 151-164. (64) Wang, F.; Wu, J.; Liu, Z. Surface Tension of Dimethyl Carbonate (C3H6O3). Fluid Phase Equilibr. 2004, 220, 121-124.

28

ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29

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

Graphic Table of Contents

29

ACS Paragon Plus Environment