Toward In Silico Prediction of CO2 Diffusion in Champagne Wines

Densities, Viscosities, Thermal Expansivities, and Isothermal Compressibilities of Carbonated Hydroalcoholic Solutions for Applications in Sparkling Beverages

Densities, viscosities, isothermal compressibilities, and thermal expansivities of carbonated hydroalcoholic solutions relevant for sparkling beverages are evaluated by molecular dynamics simulations as a function of temperature and alcoholic degree. They are compared with available experimental data, among which new measurements of densities and viscosities are performed in that respect. The OPC water model seems to yield the most accurate results, and the choice of CO₂ model has little influence on the results. Theoretical densities obtained with the OPC model typically deviate by ∼2 kg m^-3 from experimental data. At low alcoholic degrees (<9% EtOH vol), experimental viscosities lie in between theoretical values derived from the Stokes-Einstein formula and the calculation of transverse current autocorrelation functions, but at higher alcoholic degrees (≥9% EtOH vol), the Stokes-Einstein relation leads to viscosities in quantitative agreement with experiments. Isothermal compressibilities estimated with a fluctuation formula roughly extend from 0.40 to 0.49 GPa^-1 in close agreement with the experimental range of values. However, thermal expansivities are found to significantly overestimate experimental data, a behavior that is partly attributed to the low temperature of maximum density of the OPC model. Despite this discrepancy, our molecular model seems to be suitable for describing several transport and thermodynamic properties of carbonated hydroalcoholic solutions. It could therefore serve as a starting point to build more realistic models for carbonated beverages, from fizzy drinks to sparkling wines.

Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.