In many industrial applications as well as in nature, the mass transfer of CO2 at vapor-liquid interfaces in aqueous systems plays an important role. In this work, this process was studied on the atomistic level using non-equilibrium molecular dynamics simulations. In a first step, a molecular model of the system water + CO2 was developed that represents both bulk and interfacial equilibrium properties well. This system is characterized by a very large adsorption and enrichment of CO2 at the vapor-liquid interface. Then, non-equilibrium mass transfer simulations were carried out using a method that was developed recently: CO2 is inserted into the vapor phase of a simulation box which contains a liquid slab. Surprising effects are observed at the interface such as a net repulsion of CO2 particles from the interface and a complex time dependence of the amount of CO2 adsorbed at the interface.

Mass transfer through fluid interfaces is an important phenomenon in industrial applications as well as in naturally occurring processes. In this work, we investigate the mass transfer across vapor-liquid interfaces in binary mixtures using molecular dynamics simulations. We investigate the influence of interfacial properties on mass transfer by studying three binary azeotropic mixtures known to have different interfacial behaviors. Emphasis is placed on the effect of the intermolecular interactions by choosing mixtures with the same pure components but different cross-interactions such that different azeotropic behaviors are obtained. The molar flux is created by utilizing a non-stationary molecular dynamics simulation approach, where particles of one component are inserted into the vapor phase over a short period of time before the system's response to this insertion is monitored. From a direct comparison of the density profiles and the flux profiles in close proximity to the interface, we analyze the particles' tendency to accumulate in the interfacial region throughout the different stages of the simulation. We find that for mixtures with strong attractive cross-interactions, the inserted particles are efficiently transported into the liquid phase. For systems with weak attractive cross-interactions, the inserted particles show a tendency to accumulate in the interfacial region, and the flux through the system is lower. The results from this work indicate that the accumulation of particles at the interface can act as a hindrance to mass transfer, which has practical relevance in technical processes.

Fluids confined in small volumes behave differently than fluids in bulk systems. For bulk systems, a compact summary of the system's thermodynamic properties is provided by equations of state. However, there is currently a lack of successful methods to predict the thermodynamic properties of confined fluids by use of equations of state, since their thermodynamic state depends on additional parameters introduced by the enclosing surface. In this work, we present a consistent thermodynamic framework that represents an equation of state for pure, confined fluids. The total system is decomposed into a bulk phase in equilibrium with a surface phase. The equation of state is based on an existing, accurate description of the bulk fluid and uses Gibbs' framework for surface excess properties to consistently incorporate contributions from the surface. We apply the equation of state to a Lennard-Jones spline fluid confined by a spherical surface with a Weeks-Chandler-Andersen wall-potential. The pressure and internal energy predicted from the equation of state are in good agreement with the properties obtained directly from molecular dynamics simulations. We find that when the location of the dividing surface is chosen appropriately, the properties of highly curved surfaces can be predicted from those of a planar surface. The choice of the dividing surface affects the magnitude of the surface excess properties and its curvature dependence, but the properties of the total system remain unchanged. The framework can predict the properties of confined systems with a wide range of geometries, sizes, interparticle interactions, and wall-particle interactions, and it is independent of ensemble. A targeted area of use is the prediction of thermodynamic properties in porous media, for which a possible application of the framework is elaborated.

We present a new method for computing chemical potential differences of macroscopic systems by sampling fluctuations in small systems. The small system method, presented by Schnell et al. [Schnell et al., J. Phys. Chem. B, 2011, 115, 10911], is used to create small embedded systems from molecular dynamics simulations, in which fluctuations of the number of particles are sampled. The sampled fluctuations represent the Boltzmann distributed probability of the number of particles. The overlapping region of two such distributions, sampled from two different systems, is used to compute their chemical potential difference. Since the thermodynamics of small systems is known to deviate from the classical thermodynamic description, the particle distributions will deviate from the macroscopic behavior as well. We show how this can be utilized to calculate the size dependence of chemical potential differences and eventually extract the chemical potential difference in the thermodynamic limit. The macroscopic chemical potential difference is determined with a relative error of 3% in systems containing particles that interact through the truncated and shifted Lennard-Jones potential. In addition to computing chemical potential differences in the macroscopic limit directly from molecular dynamics simulation, the new method provides insights into the size dependency that is introduced to intensive properties in small systems.

Small systems have higher surface area-to-volume ratios than macroscopic systems. The thermodynamics of small systems therefore deviates from the description of classical thermodynamics. One consequence of this is that properties of small systems can be dependent on the system's ensemble. By comparing the properties in grand canonical (open) and canonical (closed) systems, we investigate how a small number of particles can induce an ensemble dependence. Emphasis is placed on the insight that can be gained by investigating ideal gases. The ensemble equivalence of small ideal gas systems is investigated by deriving the properties analytically, while the ensemble equivalence of small systems with particles interacting via the Lennard-Jones or the Weeks-Chandler-Andersen potential is investigated through Monte Carlo simulations. For all the investigated small systems, we find clear differences between the properties in open and closed systems. For systems with interacting particles, the difference between the pressure contribution to the internal energy, and the difference between the chemical potential contribution to the internal energy, are both increasing with the number density. The difference in chemical potential is, with the exception of the density dependence, qualitatively described by the analytic formula derived for an ideal gas system. The difference in pressure, however, is not captured by the ideal gas model. For the difference between the properties in the open and closed systems, the response of increasing the particles' excluded volume is similar to the response of increasing the repulsive forces on the system walls. This indicates that the magnitude of the difference between the properties in open and closed systems is related to the restricted movement of the particles in the system. The work presented in this paper gives insight into the mechanisms behind ensemble in-equivalence in small systems, and illustrates how a simple statistical mechanical model, such as the ideal gas, can be a useful tool in these investigations.