Models of gas-liquid solubilities

A recent method for phase equilibria, the AGAPE method, has been used to predict activity coefficients and excess Gibbs energy for binary mixtures with good accuracy. The theory, based on a generalised London potential (GLP), accounts for intermolecular attractive forces. Unlike existing prediction methods, for example UNIFAC, the AGAPE method uses only information derived from accessible experimental data and molecular information for pure components. Presently, the AGAPE method has some limitations, namely that the mixtures must consist of small, non-polar compounds with no hydrogen bonding, at low moderate pressures and at conditions below the critical conditions of the components. Distinction between vapour-liquid equilibria and gas-liquid solubility is rather arbitrary and it seems reasonable to extend these ideas to solubility. The AGAPE model uses a molecular lattice-based mixing rule. By judicious use of computer programs a methodology was created to examine a body of experimental gas-liquid solubility data for gases such as carbon dioxide, propane, n-butane or sulphur hexafluoride which all have critical temperatures a little above 298 K dissolved in benzene, cyclo-hexane and methanol. Within this methodology the value of the GLP as an ab initio combining rule for such solutes in very dilute solutions in a variety of liquids has been tested. Using the GLP as a mixing rule involves the computation of rotationally averaged interactions between the constituent atoms, and new calculations have had to be made to discover the magnitude of the unlike pair interactions. These numbers have been seen as significant in their own right in the context of the behaviour of infinitely-dilute solutions. A method for extending this treatment to "permanent" gases has also been developed. The findings from the GLP method and from the more general AGAPE approach have been examined in the context of other models for gas-liquid solubility, both "classical" and contemporary, in particular those derived from equations-of-state methods and from reference solvent methods.

Visualising and controlling the flow in biomolecular systems at and between multiple scales:from atoms to hydrodynamics at different locations in time and space

A novel framework for modelling biomolecular systems at multiple scales in space and time simultaneously is described. The atomistic molecular dynamics representation is smoothly connected with a statistical continuum hydrodynamics description. The system behaves correctly at the limits of pure molecular dynamics (hydrodynamics) and at the intermediate regimes when the atoms move partly as atomistic particles, and at the same time follow the hydrodynamic flows. The corresponding contributions are controlled by a parameter, which is defined as an arbitrary function of space and time, thus, allowing an effective separation of the atomistic 'core' and continuum 'environment'. To fill the scale gap between the atomistic and the continuum representations our special purpose computer for molecular dynamics, MDGRAPE-4, as well as GPU-based computing were used for developing the framework. These hardware developments also include interactive molecular dynamics simulations that allow intervention of the modelling through force-feedback devices.

Universal duality in a Luttinger liquid coupled to a generic environment

We study a Luttinger liquid (LL) coupled to a generic environment consisting of bosonic modes with arbitrary density-density and current-current interactions. The LL can be either in the conducting phase and perturbed by a weak scatterer or in the insulating phase and perturbed by a weak link. The environment modes can also be scattered by the imperfection in the system with arbitrary transmission and reflection amplitudes. We present a general method of calculating correlation functions under the presence of the environment and prove the duality of exponents describing the scaling of the weak scatterer and of the weak link. This duality holds true for a broad class of models and is sensitive to neither interaction nor environmental modes details, thus it shows up as the universal property. It ensures that the environment cannot generate new stable fixed points of the renormalization group flow. Thus, the LL always flows toward either conducting or insulating phase. Phases are separated by a sharp boundary which is shifted by the influence of the environment. Our results are relevant, for example, for low-energy transport in (i) an interacting quantum wire or a carbon nanotube where the electrons are coupled to the acoustic phonons scattered by the lattice defect; (ii) a mixture of interacting fermionic and bosonic cold atoms where the bosonic modes are scattered due to an abrupt local change of the interaction; (iii) mesoscopic electric circuits.