Evert Jan Meijer

Accurate Effective Pair Potentials for Polymer Solutions

Dilute or semidilute solutions of nonintersecting self-avoiding walk (SAW) polymer chains are mapped onto a fluid of “soft” particles interacting via an effective pair potential between their centers of mass. This mapping is achieved by inverting the pair distribution function of the centers of mass of the original polymer chains, using integral equation techniques from the theory of simple fluids. The resulting effective pair potential is finite at all distances, has a range of the order of the radius of gyration, and turns out to be only moderately concentration-dependent. The dependence of the effective potential on polymer length is analyzed in an effort to extract the scaling limit. The effective potential is used to derive the osmotic equation of state, which is compared to simulation data for the full SAW segment model, and to the predictions of renormalization group calculations. A similar inversion procedure is used to derive an effective wall–polymer potential from the center of mass density pro...

Chiral Induction Effects in Ruthenium(II) Amino Alcohol Catalysed Asymmetric Transfer Hydrogenation of Ketones: An Experimental and Theoretical Approach

The enantioselective outcome of transfer hydrogenation reactions that are catalysed by ruthenium(II) amino alcohol complexes was studied by means of a systematically varied series of ligands. It was found that both the substituent at the 1-position in the 2-amino-1-alcohol ligand and the substituent at the amine functionality influence the enantioselectivity of the reaction to a large extent: enantioselectivities (ee values) of up to 95% were obtained for the reduction of acetophenone. The catalytic cycle of ruthenium(II) amino alcohol catalysed transfer hydrogenation was examined at the density functional theory level. The formation of a hydrogen bond between the carbonyl functionality of the substrate and the amine proton of the ligand, as well as the formation of an intramolecular H...H bond and a planar H-Ru-N-H moiety are crucially important for the reaction mechanism. The enantioselective outcome of the reaction can be illustrated with the aid of molecular modelling by the visualisation of the steric interactions between the ketone and the ligand backbone in the ruthenium(II) catalysts.

Melting line of Yukawa system by computer simulation

We located the melting line of the Yukawa system by determining the free energy of both fluid and solid phases by computer simulations. At the high densities the fluid freezes into a body‐centered‐cubic (bcc) solid, whereas for low densities it freezes into a face‐centered‐cubic (fcc) solid. For both the bcc‐fluid and the fcc‐fluid part of the melting line two coexistence points were determined. We observed that the relative root‐mean‐square displacement in the solid (Lindeman ratio) varies considerably along the melting line. At the bcc‐fluid part of the melting line the Lindeman ratio equals 0.19, whereas at the fcc‐solid part of the melting line it is smaller than 0.16.

Polymer induced depletion potentials in polymer-colloid mixtures

The depletion interactions between two colloidal plates or between two colloidal spheres, induced by interacting polymers in a good solvent, are calculated theoretically and by computer simulations. A simple analytical theory is shown to be quantitatively accurate for the case of two plates. A related depletion potential is derived for two spheres; it also agrees very well with direct computer simulations. Theories based on ideal polymers show important deviations with increasing polymer concentration: They overestimate the range of the depletion potential between two plates or two spheres at all densities, with the largest relative change occurring in the dilute regime. They underestimate the well depth at contact for the case of two plates, but overestimate it for two spheres. Depletion potentials are also calculated using a coarse graining approach which represents the polymers as “soft colloids;” good agreement is found in the dilute regime. Finally, the effect of the polymers on colloid–colloid osmot...

Ab initio molecular dynamics study of liquid methanol

Abstract We present a density-functional theory based molecular-dynamics study of the structural, dynamical, and electronic properties of liquid methanol under ambient conditions. The calculated radial distribution functions involving the oxygen and hydroxyl hydrogen show a pronounced hydrogen bonding and compare well with recent neutron diffraction data, except for an underestimate of the oxygen–oxygen correlation. We observe that, in line with infrared spectroscopic data, the hydroxyl stretching mode is significantly red-shifted in the liquid. A substantial enhancement of the dipole moment is accompanied by significant fluctuations due to thermal motion. Our results provide valuable data for improvement of empirical potentials.

Density-functional theory-based molecular simulation study of liquid methanol

We present a density-functional theory based molecular dynamics study of the structural, dynamical, and electronic properties of liquid methanol under ambient conditions. The calculated radial distribution functions involving the oxygen and hydroxyl hydrogen show a pronounced hydrogen bonding and compare well with recent neutron diffraction data. We observe that, in line with infrared spectroscopic data, the hydroxyl-stretching mode is significantly redshifted in the liquid, whereas the hydroxyl bending mode shows a blueshift. A substantial enhancement of the molecular dipole moment is accompanied by significant fluctuations due to thermal motion. We compute a value of 32 for the relative permittivity, almost identical to the experimental value of 33. Our results provide valuable data for improvement of empirical potentials.

Location of melting point at 300 K of nitrogen by Monte Carlo simulation

We present an accurate new method to compute absolute free energies of molecular solids by computer simulations. As a first application, we computed the thermodynamic phase transition between the fluid phase and the orientational disordered solid β phase of nitrogen at 300 K, using a well tested pair potential. The computed coexistence pressure and the volume change coincides within the error margins with the experimental values. The coexistence volume differed by 2% from the experimental value. To our knowledge these results constitutes the first numerical calculation of the thermodynamic stability for a model of a realistic molecular solid.

Hydration of methanol in water. A DFT-based molecular dynamics study

Abstract We studied the hydration of a single methanol molecule in aqueous solution by first-principle DFT-based molecular dynamics simulation. The calculations show that the local structural and short-time dynamical properties of the water molecules remain almost unchanged by the presence of the methanol, confirming the observation from recent experimental structural data for dilute solutions. We also see, in accordance with this experimental work, a distinct shell of water molecules that consists of about 15 molecules. We found no evidence for a strong tangential ordering of the water molecules in the first hydration shell.

Density profiles and surface tension of polymers near colloidal surfaces

The surface tension of interacting polymers in a good solvent is calculated theoretically and by computer simulations for a planar wall geometry and for the insertion of a single colloidal hard sphere. This is achieved for the planar wall and for the larger spheres by an adsorption method, and for smaller spheres by a direct insertion technique. Results for the dilute and semidilute regimes are compared to results for ideal polymers, the Asakura–Oosawa penetrable-sphere model, and to integral equations, scaling and renormalization group theories. The largest relative changes with density are found in the dilute regime, so that theories based on noninteracting polymers rapidly break down. A recently developed “soft colloid” approach to polymer–colloid mixtures is shown to correctly describe the one-body insertion free-energy and the related surface tension.

State-of-the-art models for the phase diagram of carbon and diamond nucleation

We review recent developments in the modelling of the phase diagram and the kinetics of crystallization of carbon. In particular, we show that a particular class of bond-order potentials (the so-called LCBOP models) account well for many of the known structural and thermodynamic properties of carbon at high pressures and temperatures. We discuss the LCBOP models in some detail. In addition, we briefly review the ‘history’ of experimental and theoretical studies of the phase behaviour of carbon. Using a well-tested version of the LCBOP model (viz. LCBOPI+) we address some of the more controversial hypotheses concerning the phase behaviour of carbon, in particular: the suggestion that liquid carbon can exist in two phases separated by a first-order phase transition and the conjecture that diamonds could have formed by homogeneous nucleation in Uranus and Neptune.