Molecular simulation of gas adsorption equilibria in nanoporous materials

This work is divided into two distinct parts. The first part consists of the study of the metal organic framework UiO-66Zr, where the aim was to determine the force field that best describes the adsorption equilibrium properties of two different gases, methane and carbon dioxide. The other part of the work focuses on the study of the single wall carbon nanotube topology for ethane adsorption; the aim was to simplify as much as possible the solid-fluid force field model to increase the computational efficiency of the Monte Carlo simulations. The choice of both adsorbents relies on their potential use in adsorption processes, such as the capture and storage of carbon dioxide, natural gas storage, separation of components of biogas, and olefin/paraffin separations. The adsorption studies on the two porous materials were performed by molecular simulation using the grand canonical Monte Carlo (μ,V,T) method, over the temperature range of 298-343 K and pressure range 0.06-70 bar. The calibration curves of pressure and density as a function of chemical potential and temperature for the three adsorbates under study, were obtained Monte Carlo simulation in the canonical ensemble (N,V,T); polynomial fit and interpolation of the obtained data allowed to determine the pressure and gas density at any chemical potential. The adsorption equilibria of methane and carbon dioxide in UiO-66Zr were simulated and compared with the experimental data obtained by Jasmina H. Cavka et al. The results show that the best force field for both gases is a chargeless united-atom force field based on the TraPPE model. Using this validated force field it was possible to estimate the isosteric heats of adsorption and the Henry constants. In the Grand-Canonical Monte Carlo simulations of carbon nanotubes, we conclude that the fastest type of run is obtained with a force field that approximates the nanotube as a smooth cylinder; this approximation gives execution times that are 1.6 times faster than the typical atomistic runs.

Flow analysis-hydride generation-gas phase derivative molecular absorption spectrophotometric determination of antimony in antileishmanial drugs

In the present work, the development of a method based on the coupling of flow analysis (FA), hydride generation (HG), and derivative molecular absorption spectrophotometry (D-EAM) in gas phase (GP), is described in order to determine total antimony in antileishmanial products. Second derivative order (D²224nm) of the absorption spectrum (190 - 300 nm) is utilized as measurement criterion. Each one of the parameters involved in the development of the proposed method was examined and optimized. The utilization of the EAM in GP as detection system in a continuous mode instead of atomic absorption spectrometry represents the great potential of the analytic proposal.

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.

The molecular structure of hexamethyldigermane determined by gas-phase electron diffraction with theoretical calculations for (CH3)3M-M(CH3)3 Where M = C, Si, and Ge

Gas-phase electron diffraction (GED) data together with results from ab initio molecular orbital calculations (HF and MP2/6-311+G(d,p)) have been used to determine the structure of hexamethyldigermane ((CH3)3Ge-Ge(CH3)3). The equilibrium symmetry is D3d, but the molecule has a very low-frequency, largeamplitude, torsional mode (φCGeGeC) that lowers the thermal average symmetry. The effect of this largeamplitude mode on the interatomic distances was described by a dynamic model which consisted of a set of pseudoconformers spaced at even intervals. The amount of each pseudoconformer was obtained from the ab initio calculations (HF/6-311+G(d,p)). The results for the principal distances (ra) and angles (∠h1) obtained from the combined GED/ab initio (with estimated 1σ uncertainties) are r(Ge-Ge) ) 2.417(2) Å, r(Ge-C) ) 1.956(1) Å, r(C-H) ) 1.097(5) Å, ∠GeGeC ) 110.5(2)°, and ∠GeCH ) 108.8(6)°. Theoretical calculations were performed for the related molecules ((CH3)3Si-Si(CH3)3 and (CH3)3C-C(CH3)3).

Molecular structure of 2,5-dihydropyrrole (C4NH7), obtained by gas-phase electron diffraction and theoretical calculations

The structure of 2,5-dihydropyrrole (C4NH7) has been determined by gas-phase electron diffraction (GED), augmented by the results from ab initio calculations employing third-order Moller-Plesset (MP3) level of theory and the 6-311+G(d,p) basis set. Several theoretical calculations were performed. From theoretical calculations using MP3/6-311+G(d,p) evidence was obtained for the presence of an axial (63%) (N-H bond axial to the CNC plane) and an equatorial conformer (37%) (N-H bond equatorial to the CNC plane). The five-membered ring was found to be puckered with the CNC plane inclined at 21.8 (38)° to the plane of the four carbon atoms.

Molecular structures of trimethylchlorogermane, (CH3)(3)GeCl, and trimethylbromogermane, (CH3)(3)GeBr, obtained by gas-phase electron diffraction and theoretical calculations

The structures of trimethylchlorogermane ((CH3)(3)GeCl) and trimethylbromogermane ((CH3)(3)GeBr) have been determined by gas-phase electron diffraction (GED), augmented by the results from ab initio calculations employing second-order Moller-Plesset (MP2) level of theory and the 6-311+G(d) basis set. All the electrons were included in the correlation calculation. The results from the ab initio calculations indicated that these molecules have C-3v symmetry, and models with this symmetry were used in the electron diffraction analysis. The results for the principal distances (r(g)) and angles (angle(alpha)) from the combined GED/ab initio study of trimethylchlorogermane (with estimated 2sigma uncertainties) are: r(Ge-C) = 1.950(4) Angstrom, r(Ge-Cl) = 2.173(4) Angstrom, r(C-H) = 1.090(9) Angstrom, angleCGeC = 112.7(7)degrees, angleCGeCl = 106.0(8)degrees, angleGeCH = 107.8(12)degrees. The results for the principal distances (r(g)) and angles (angle(alpha)) from the combined GED/ab initio study of trimethylbromogermane (with estimated 2sigma uncertainties) are: r(Ge-C) = 1.952(7) Angstrom, r(Ge-Br) = 2.325(4) Angstrom, r(C-H) = 1. 140(28) Angstrom, angleCGeC = 114.2(11)degrees, angleCGeBr = 104.2(13)degrees, angleGeCH 106.9(43)degrees. Local C-3v symmetry and staggered conformation were assumed for the methyl groups.

Molecular structures of benzoic acid and 2-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of benzoic acid (C6H5COOH) and 2-hydroxybenzoic acid (C6H4OHCOOH) have been determined in the gas phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d, p)) predict two conformers for benzoic acid, one which is 25.0 kJ mol(-1) (MP2) lower in energy than the other. In the low-energy form, the carboxyl group is coplanar with the phenyl ring and the O-H group eclipses the C=O bond. Theoretical calculations (HF and MP2/6-311+ G(d, p)) carried out for 2-hydroxybenzoic acid gave evidence for seven stable conformers but one low-energy form (11.7 kJ mol-1 lower in energy (MP2)) which again has the carboxyl group coplanar with the phenyl ring, the O-H of the carboxyl group eclipsing the C=O bond and the C=O of the carboxyl group oriented toward the O-H group of the phenyl ring. The effects of internal hydrogen bonding in 2-hydroxybenzoic acid can be clearly observed by comparison of pertinent structural parameters between the two compounds. These differences for 2-hydroxybenzoic acid include a shorter exocyclic C-C bond, a lengthening of the ring C-C bond between the substituents, and a shortening of the carboxylic single C-O bond.

Molecular structures of 3-hydroxybenzoic acid and 4-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of 3-hydroxybenzoic acid and 4-hydroxybenzoic acid have been determined by gas-phase electron diffraction using results from quantum chemical calculations to inform the choice of restraints applied to some of the structural parameters. The results from the study presented here demonstrate that resonance hybrids are not as helpful in rationalizing the structures of 2-, 3-, and 4-hydroxybenzoic acids as are models based upon electrostatic effects.

Gas-phase rate coefficients for the reactions of O(3P), S(3P), Se(3P), and Te(3P) with alkenes: application of perturbation frontier molecular orbital theory, correlations, and structure-activity relations (SARs)

The kinetics of the reactions of the atoms O(P-3), S(P-3), Se(P-3), and Te((3)p) with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. These correlations may be employed to predict rate coefficients from the calculated HOMO energy of any other alkene of interest. The rate coefficients obtained from the correlations were used to formulate structure-activity relations (SARs) for reactions of O((3)p), S(P-3), Se (P-3), and Te((3)p) with alkenes. A comparison of the values predicted by both the correlations and the SARs with experimental data where they exist allowed us to assess the reliability of our method. We demonstrate the applicability of perturbation frontier molecular orbital theory to gas-phase reactions of these atoms with alkenes. The correlations are apparently not applicable to reactions of C(P-3), Si(P-3), N(S-4), and Al(P-2) atoms with alkenes, a conclusion that could be explained in terms of a different mechanism for reaction of these atoms.

The molecular structure of propylene sulphide (methylthiirane) by gas-phase electron diffraction and theoretical calculations: a molecule used in MOCVD

Gas-phase electron-diffraction (GED) data together with results from ab initio molecular orbital calculations have been used to determine the structure of propylene sulphide. Values found for the main structural parameters for the molecule are consistent with those obtained from microwave studies and are compared here with those found for similar sulphur containing rings of general formula S(CH2)n (n = 2–5). A high ring strain enthalpy was calculated for propylene sulphide which is consistent with the small C–S–C angle (48.2(6)degrees) and the relatively long C–S bond lengths (ra = 1.831(2) Å). This is thought to account for the ease of ring opening in propylene sulphide observed in MOCVD reactions and the ready polymerisation of the molecule.

Applicability of molecular simulations for modelling the adsorption of the greenhouse gas CF4 on carbons

Tetrafluoromethane, CF4, is powerful greenhouse gas, and the possibility of storing it in microporous carbon has been widely studied. In this paper we show, for the first time, that the results of molecular simulations can be very helpful in the study of CF4 adsorption. Moreover, experimental data fit to the results collected from simulations. We explain the meaning of the empirical parameters of the supercritical Dubinin–Astakhov model proposed by Ozawa and finally the meaning of the parameter k of the empirical relation proposed by Amankwah and Schwarz.

Hydrogen bonding in the gas-phase: the molecular structures of 2‑hydroxybenzamide (C7H7NO2) and 2‑methoxybenzamide (C8H9NO2), obtained by gas-phase electron diffraction and theoretical calculations

The structures of 2-hydroxybenzamide(C7H7NO2) and 2-methoxybenzamide (C8H9NO2) have been determined in the gas-phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d,p)) predict four stable conformers for both 2-hydroxybenzamide and 2-methoxybenzamide. For both compounds, evidence for intramolecular hydrogen bonding is presented. In 2-hydroxybenzamide, the observed hydrogen bonded fragment is between the hydroxyl and carbonyl groups, while in 2-methoxybenzamide, the hydrogen bonded fragment is between one of the hydrogen atoms of the amide group and the methoxy oxygen atom.