A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

Characterization of activated carbons using liquid phase adsorption

A modification of the Dubinin-Radushkevich pore filling model by incorporation of the repulsive contribution to the pore potential, and of bulk non-ideality, is proposed in this paper for characterization of activated carbon using liquid phase adsorption. For this purpose experiments have been performed using ethyl propionate, ethyl butyrate, and ethyl isovalerate as adsorbates and the microporous-mesoporous activated carbons Filtrasorb 400, Norit ROW 0.8 and Norit ROX 0.8 as adsorbents. The repulsive contribution to the pore potential is incorporated through a Lennard-Jones intermolecular potential model, and the bulk-liquid phase non-ideality through the UNIFAC activity coefficient model. For the characterization of activated carbons, the generalized adsorption isotherm is utilized with a bimodal gamma function as the pore size distribution function. It is found that the model can represent the experimental data very well, and significantly better than when the classical energy-size relationship is used, or when bulk non-ideality is neglected. Excellent agreement between the bimodal gamma pore size distribution and DFT-cum-regularization based pore size distribution is also observed, supporting the validity of the proposed model. (C) 2001 Elsevier Science Ltd. All rights reserved.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

Boson-boson mixtures and triplet correlations

The hypernetted-chain formalism for boson-boson mixtures described by an extended Jastrow correlated wave function is derived, taking into account elementary diagrams and triplet correlations. The energy of an ideal boson 3He-4He mixture is computed for low values of the 3He concentration. The zero-3He-concentration limit provides a 3He chemical potential in good agreement with the experimental value, when a McMillan two-body correlation factor and the Lennard-Jones potential are adopted. If the Euler equations for the two-body correlation factors are solved in presence of triplet correlations, the agreement is again improved. At the experimental 4He equilibrium density, the 3He chemical potential turns out to be -2.58 K, to be compared with the experimental value, -2.79 K.

Model-independent bounds for the potential and kinetic energy of liquid 4He at zero temperature

Using the experimental values of the chemical potentials of liquid 4He and of a 3He impurity in liquid 4He, we derive a model-independent lower (upper) bound to the kinetic (potential) energy per particle at zero temperature. The values of the bounds at the experimental saturation density are 13.42 K for the kinetic energy and -20.59 K for the potential energy. All the theoretical calculations based on the Lennard-Jones potential violate the upper-bound condition for the potential energy.

Collective dynamic properties in simple crystals: Influence of the structural disorder

Collective dynamic properties in Lennard-Jones crystals are investigated by molecular dynamics simulation. The study is focused on properties such as the dynamic structure factors, the longitudinal and transverse currents and the density of states. The influence on these properties of the structural disorder is analyzed by comparing the results for one-component crystals with those for liquids and supercooled liquids at analogous conditions. The effects of species-disorder on the collective properties of binary crystals are also discussed.

Termodinâmica estatística de líquidos com o método de Monte Carlo. I. Metodologia

Statistical mechanics Monte Carlo simulation is reviewed as a formalism to study thermodynamic properties of liquids. Considering the importance of free energy changes in chemical processes, the thermodynamic perturbation theory implemented in the Monte Carlo method is discussed. The representation of molecular interaction by the Lennard-Jones and Coulomb potential functions is also discussed. Charges derived from quantum molecular electrostatic potential are also discussed as an useful methodology to generate an adequate set of partial charges to be used in liquid simulation.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Um modelo para a superfície líquida no estudo da dinâmica do espalhamento de Xe e Ne pelo esqualano

In this work we present a theoretical model to investigate the scattering of Xe and Ne by a liquid squalane surface. The liquid surface is modeled as a grid of harmonic oscillators with frequencies adjusted to experimental vibration as frequencies of the liquid squalane and the atom-surface interaction potential is modeled by a Lennard-Jones function. The three dimensional description of the dynamics of the process which occurs at the gas-liquid interface is obtained by the classical trajectory method. The general characteristics of the dynamics of the scattering process are in good agreement with experimental data.

An application of cell-cluster theory to a rare gas crystal /|n[by] K. Westera. - 260 St. Catharines [Ont.] : Dept. of Physics, Brock University,

The Lennard-Jones Devonshire 1 (LJD) single particle theory for liquids is extended and applied to the anharmonic solid in a high temperature limit. The exact free energy for the crystal is expressed as a convergent series of terms involving larger and larger sets of contiguous particles called cell-clusters. The motions of all the particles within cell-clusters are correlated to each other and lead to non-trivial integrals of orders 3, 6, 9, ... 3N. For the first time the six dimensional integral has been calculated to high accuracy using a Lennard-Jones (6-12) pair interaction between nearest neighbours only for the f.c.c. lattice. The thermodynamic properties predicted by this model agree well with experimental results for solid Xenon.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

A treatment of atomic mean square displacement in higher order perturbation theory

A general derivation of the anharmonic coefficients for a periodic lattice invoking the special case of the central force interaction is presented. All of the contributions to mean square displacement (MSD) to order 14 perturbation theory are enumerated. A direct correspondance is found between the high temperature limit MSD and high temperature limit free energy contributions up to and including 0(14). This correspondance follows from the detailed derivation of some of the contributions to MSD. Numerical results are obtained for all the MSD contributions to 0(14) using the Lennard-Jones potential for the lattice constants and temperatures for which the Monte Carlo results were calculated by Heiser, Shukla and Cowley. The Peierls approximation is also employed in order to simplify the numerical evaluation of the MSD contributions. The numerical results indicate the convergence of the perturbation expansion up to 75% of the melting temperature of the solid (TM) for the exact calculation; however, a better agreement with the Monte Carlo results is not obtained when the total of all 14 contributions is added to the 12 perturbation theory results. Using Peierls approximation the expansion converges up to 45% of TM• The MSD contributions arising in the Green's function method of Shukla and Hubschle are derived and enumerated up to and including 0(18). The total MSD from these selected contributions is in excellent agreement with their results at all temperatures. Theoretical values of the recoilless fraction for krypton are calculated from the MSD contributions for both the Lennard-Jones and Aziz potentials. The agreement with experimental values is quite good.

Molecular dynamics calculation of mean square displacement in alkali metals and rare gas solids and comparison with lattice dynamics

Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.