Towards the Quantitative Prediction of the Phase Behavior of Polymer Solutions by Computer Simulation

被引:5
作者
Binder, Kurt [1 ]
Mognetti, Bortolo Matteo [1 ]
Macdowell, Luis Gonzalez [2 ]
Oettel, Martin [1 ]
Paul, Wolfgang [1 ]
Virnau, Peter [1 ]
Yelash, Leonid [1 ]
机构
[1] Johannes Gutenberg Univ Mainz, Inst Phys, D-55099 Mainz, Germany
[2] Univ Complutense, Dpto Quim Fis, Fac Cc Quim, E-28040 Madrid, Spain
关键词
carbon dioxide; monte carlo simulations; polar fluids; THERMODYNAMIC PERTURBATION-THEORY; ORNSTEIN-ZERNIKE EQUATION; VAPOR-LIQUID-EQUILIBRIUM; CARBON-DIOXIDE; BINARY-MIXTURES; OF-STATE; SYSTEMS; HEXADECANE; MODELS;
D O I
10.1002/masy.200950401
中图分类号
O63 [高分子化学(高聚物)];
学科分类号
070305 ; 080501 ; 081704 ;
摘要
The phase diagram of polymer solutions (cf. e.g. alkanes dissolved in supercritical carbon dioxide) is complicated, since there are four control parameters (temperature, pressure, monomer volume fraction, chain length of the polymer) and due to the interplay of liquid-vapor transitions and fluid-fluid unmixing. As a result very intricate phase diagram topologies can result. An attempt to develop coarse-grained models that can deal with this task will be described. As usual, the polymers will be modelled as off-lattice bead-spring chains, where several chemical monomers are integrated into one effective bond, torsional degrees of freedom being disregarded. But also a coarse-grained description of the solvent is needed: we show that using a simple point particle with a quadrupole moment and suitable Lennard-Jones interaction yields a very good description of pure carbon dioxide (better than fully atomistic models with potentials from ab initio quantum chemistry calculations). The strength of the quadrupole moment is taken from experiment, and the Lennard-Jones parameters are adjusted such that the experimental critical temperature and density is correctly reproduced. This procedure works for other solvents as well, such as benzene. The pure alkane phase diagram is also reproduced by similarly chosen Lennard-Jones (LJ) potentials among the monomers, and the monomer-carbon dioxide LJ parameters are chosen from Lorentz-Berthelot mixing rules. These parameters are then used as input both for Monte Carlo calculations and approximate methods to calculate the equation of state, such as density functional and perturbation theory methods. Apart from the region close to critical points, where mean-field type methods fail, the analytical calculations agree well with the simulations. The solvent model, in which quadrupolar interactions are explicitly considered, is important in order to obtain a fair agreement with mixture experimental results.
引用
收藏
页码:1 / 9
页数:9
相关论文
共 33 条
[1]   Polymer plus solvent systems:: Phase diagrams, interface free energies, and nucleation [J].
Binder, K ;
Müller, M ;
Virnau, P ;
MacDowell, LG .
ADVANCED COMPUTER SIMULATION APPROACHES FOR SOFT MATTER SCIENCES I, 2005, 173 :1-104
[2]   FINITE SIZE SCALING ANALYSIS OF ISING-MODEL BLOCK DISTRIBUTION-FUNCTIONS [J].
BINDER, K .
ZEITSCHRIFT FUR PHYSIK B-CONDENSED MATTER, 1981, 43 (02) :119-140
[3]  
Binder K., 1995, MONTE CARLO MOL DYNA
[4]   VAPOR LIQUID EQUILIBRIUM IN THE SYSTEM CARBON-DIOXIDE + N-PENTANE FROM 252-K TO 458-K AT PRESSURES TO 10-MPA [J].
CHENG, HZ ;
DEFERNANDEZ, MEP ;
ZOLLWEG, JA ;
STREETT, WB .
JOURNAL OF CHEMICAL AND ENGINEERING DATA, 1989, 34 (03) :319-323
[5]   OPTIMIZED MONTE-CARLO DATA-ANALYSIS [J].
FERRENBERG, AM ;
SWENDSEN, RH .
PHYSICAL REVIEW LETTERS, 1989, 63 (12) :1195-1198
[6]   NEW MONTE-CARLO TECHNIQUE FOR STUDYING PHASE-TRANSITIONS [J].
FERRENBERG, AM ;
SWENDSEN, RH .
PHYSICAL REVIEW LETTERS, 1988, 61 (23) :2635-2638
[7]  
Fisher M.E., 1971, CRITICAL PHENOMENA
[8]  
Flory P J., PRINCIPLES POLYM CHE
[9]   Thermodynamics of high polymer solutions [J].
Flory, PJ .
JOURNAL OF CHEMICAL PHYSICS, 1941, 9 (08) :660-661
[10]   3-PHASE LIQUID-LIQUID-VAPOR EQUILIBRIA BEHAVIOR OF CERTAIN BINARY CO2-N-PARAFFIN SYSTEMS [J].
HOTTOVY, JD ;
LUKS, KD ;
KOHN, JP .
JOURNAL OF CHEMICAL AND ENGINEERING DATA, 1981, 26 (03) :256-258