A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

Polydifferentials and the deformation functor of curves with automorphisms

We give a relation between the dimension of the tangent space of the deformation functor of curves with automorphisms and the Galois module structure of the space of 2-holomorphic differentials. We prove a homological version of the local-global principle similar to the one of J. Bertin and A. Mézard. Let G be a cyclic subgroup of the group of automorphisms of a curve X, so that the order of G is equal to the characteristic. By using the results of S. Nakajima on the Galois module structure of the space of 2-holomorphic differentials, we compute the dimension of the tangent space of the deformation functor.     

Reactions of dimethyldivinylsilane,dimethyldivinyltin and allyltrimethyltin with diethylene (tertiary) phosphine)platinum complexes

The compounds [Pt(C₂H₄)₂(PR₃)] [PR₃ = P-tBu₂Me, P(C₆H₁₁)₃, PPh₃] react dimethyldivinylsilane or dimethyldivinyltin to give chelate complexes [Pt{(CH₂CH)₂MMe₂} (PR₃)] (M = Si or Sn). allyltrimethyltin reacts with various diethylene (tertiary phosphine)platinum compounds with cleavage of the allyl group to afford complexes [Pt(SnMe₃)(η³-C₃H₅)(PR₂)]. The NMR spectra (¹³C, ¹H and ³¹P) of the new compounds have been recorded, and the data are discussed in terms of the structures proposed.     

The dipole moments of the excited states of FeC

With the purpose of comparing expectation dipole moment values mu with finite-field obtained dipole moments mu(FF), we recalculated by the finite-field method previously reported mu values of 38 excited states of FeC. In most of the cases mu(FF) is significantly larger than mu.     

The electronic structure of vanadium carbide, VC

Within an energy range of 2.4 eV, we have explored 29 of the 36 states of the diatomic molecule VC that arise from the atoms in their ground state, V(4s23d3;4F)+C(2s2 2p2;3P). We use multireference methods with large atomic natural orbital basis sets. The ground state is of 2Delta symmetry with the first two excited states, 4Delta and 2Sigma+, located 4.2 and 7.0 kcal/mol above the X state. All the states examined in this work are relatively strongly bound and show significant charge transfer from V to C. The binding energy of the X 2Delta state is estimated to be 95.3 kcal/mol in good agreement with the experimental value.     

A multireference coupled-cluster potential energy surface of diazomethane, CH(2)N(2)

The intrinsically multireference dissociation of the C-N bond in ground-state diazomethane (CH(2)N(2)) at different angles has been studied with the multireference Brillouin-Wigner coupled-cluster singles and doubles (MRBWCCSD) method. The morphology of the calculated potential energy surface (PES) in C(s)() symmetry is similar to a multireference perturbational (CASPT3) PES. The MRBWCCSD/cc-pVTZ H(2)C-N(2) dissociation energy with respect to the asymptotic CH(2)(?(1)A(1)) + N(2)(X(1)Sigma(g)(+)) products is D(e) = 35.9 kcal/mol, or a zero-point corrected D(0) = 21.4 kcal/mol with respect to the ground-state CH(2)(X(3)B(1)) + N(2)(X(1)Sigma(g)(+)) fragments.     

Partition coefficients of organic molecules in squalane and water/ethanol mixtures by molecular dynamics simulations

Accurate partition coefficient data of migrants between a polymer and a solvent are of paramount importance for estimating the migration of the migrant over time, including the concentration of the migrant at infinite time in the two solvents. In this article it is shown how this partition coefficient can be estimated for both a small hydrophilic and a hydrophobic organic molecules between squalane (used here to mimic low density poly ethylene) and water/ethanol solutes using thermodynamic integration to calculate the free energy of solvation. Molecular dynamics simulations are performed, using the GROMACS software, by slowly decoupling of firstly the electrostatic and then the Lennard-Jones interactions between molecules in the simulation box. These calculations depend very much on the choice of force field. Two force fields have been tested in this work, the TraPPE-UA (united-atom) and the OPLS-AA (all-atom). The computational cheaper TraPPE-UA force field showed to be more accurate over the whole range of systems compared to the OPLS-AA force field. Moreover, some of the calculations were done with five different water models to investigate the influence of the specific water model on the calculations. It was found that the combination of the TraPPE-UA force field and the TIP4p water model gave the best results. Based on the methodology proposed in this article, it is possible to obtain good partition coefficients only knowing the chemical structure of the molecules in the system.     

Molecular Mixing and Flowfield Measurements in a Recirculating Shear Flow. Part I: Subsonic Flow

The mixing and flowfield of a complex geometry, similar to a rearward-facing step flow but with injection, is studied. A subsonic top-stream is expanded over a perforated ramp at an angle of 30°, through which a secondary stream is injected. The mass flux of the second stream is chosen to be insufficient to provide the entrainment requirements of the shear layer, which, as a consequence, attaches to the lower guidewall. Part of the flow is directed upstream forming a re-entrant jet within the recirculation zone that enhances mixing and flameholding. A control-volume model of the flow is found to be in good agreement with the variation of the overall pressure coefficient of the device with variable mass injection. The flowfield response to changing levels of heat release is also quantified. While increased heat release acts somewhat analogously to increased mass injection, fundamental differences in the flow behaviour are observed. The hypergolic hydrogen-fluorine chemical reaction employed allows the level of molecular mixing in the flow to be inferred. The amount of mixing is found to be higher in the expansion-ramp geometry than in classical free-shear layers. As in free-shear layers, the level of mixing is found to decrease with increasing top-stream velocity. Results for a similar configuration with supersonic flow in the top stream are reported in Part II of this two-part series.     

Use of monomer fraction data in the parametrization of association theories

Association theories such as the CPA (cubic-plus-association), NRHB (non-random hydrogen bonding) equations of state and the various variants of SAFT (statistical associating fluid theory) have been extensively applied to phase equilibrium calculations. Such models can also be used for estimating the monomer fraction of hydrogen bonding compounds and their mixtures. Monomer fraction data are obtained from spectroscopic measurements and they are available for a few compounds such as pure water and alcohols as well as for some alcohol–alkane and similar mixtures. These data are useful for an understanding of the capabilities and limitations of association models. The purpose of this work is two-fold: (i) to compare the performance of three models, CPA, NRHB and sPC-SAFT, in predicting the monomer fraction of water, alcohols and mixtures of alcohol-inert compounds and (ii) to investigate whether “improved” model parameters can be obtained if monomer fraction data are included in the parameter estimation together with vapor pressures and liquid densities. The expression “improved” implies parameters which can represent several pure compound properties as well as monomer fraction data for pure compounds and mixtures. The accuracy of experimental monomer fraction data is discussed, as well as the role of monomer fraction data in clarifying which association scheme should be used in these equations of state. The results reveal that the investigated association models (CPA, sPC-SAFT and NRHB) can predict, at least qualitatively correct, monomer fractions of associating compounds and mixtures. Only, small differences are observed between the models. In addition, it has been shown that, using a suitable association scheme, a single set of parameters can describe satisfactorily vapor pressures, liquid densities and monomer fractions of water and alcohols. The 4C scheme is the best choice for water, while for methanol there is small difference between the 2B and 3B association schemes.     

What is your “birthday elliptic curve”?

In this article, Ramanujan–Weber class invariants and its analogue are used to derive birthday elliptic curves.