How much carbon dioxide can be stored in the structure H clathrate hydrates?: a molecular dynamics study

The stability of structure H (sH) carbon dioxide clathrate hydrates at three temperature-pressure conditions are determined by molecular dynamics simulations on a 3x3x3 sH unit cell replica. Simulations are performed at 100 K at ambient pressure, 273 K at 100 bars and also 300 K and 5.0 kbars. The small and medium cages of the sH unit cell are occupied by a single carbon dioxide guest and large cage guest occupancies of 1-5 are considered. Radial distribution functions are given for guests in the large cages and unit cell volumes and configurational energies are studied as a function of large cage CO(2) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at three temperature/pressure conditions stated above. At the low temperature, large cage occupancy of 5 is the most stable while at the higher temperature, the occupancy of 3 is the most favored. Calculations are also performed to show that the CO(2) sH clathrate is more stable than the methane clathrate analog. Implications on CO(2) sequestration by clathrate formation are discussed.     

On quantitative structure of small Ree groups

The main aim of this article is to study quantitative structure of small Ree Groups ²G₂(q). Here, we prove that small Ree groups are uniquely determined by their orders and the set of the number of elements of the same order.     

A theoretical study of the solvent effect on the interaction of C<Subscript>20</Subscript> and N<Subscript>2</Subscript>H<Subscript>2</Subscript>

In this work, the interaction of C₂₀ and the N₂H₂ fragment is investigated at the M062X/6-311G(d,p) level of theory in both gas and solution phases. The interaction energies obtained by the standard method are corrected by the basis set superposition error (BSSE) during the geometry optimization for all molecules at the same levels of theory. The results obtained from these calculations reveal that the interaction between C₂₀ and N₂H₂ increases in the presence of more polar solvents. Values of the electrophilic charge transfer show the charge flow from C₂₀ to N₂H₂. The influence of the solvent on the hyperpolarizability indicates that β_tot values decrease on passing from vacuum to the solution phase.     

Stochastic perturbation theory: a low-scaling approach to correlated electronic energies

We discuss a stochastic implementation of M?ller-Plesset (MP) theory based upon the concept of a "graph," a set of connected Slater determinants. We show how contributions from an arbitrary level, MPn, of perturbation theory can be expressed diagrammatically in terms of graphs, and that these may be stochastically sampled to give a good estimate of the energy. We show this to be the case for Ne, Ar, N2, and H2O molecules. N-molecule chains of He atoms and H2 molecules at equilibrium and stretched geometries show an effective scaling of O[N(2.6)] and O[N(5.6)] for MP2 and MP3 theories.     

Molecular simulation of non-equilibrium methane hydrate decomposition process

We recently performed constant energy molecular dynamics simulations of the endothermic decomposition of methane hydrate in contact with water to study phenomenologically the role of mass and heat transfer in the decomposition rate [S. Alavi, J.A. Ripmeester, J. Chem. Phys. 132 (2010) 144703]. We observed that with the progress of the decomposition front temperature gradients are established between the remaining solid hydrate and the solution phases. In this work, we provide further quantitative macroscopic and molecular level analysis of the methane hydrate decomposition process with an emphasis on elucidating microscopic details and how they affect the predicted rate of methane hydrate decomposition in natural methane hydrate reservoirs. A quantitative criterion is used to characterize the decomposition of the hydrate phase at different times. Hydrate dissociation occurs in a stepwise fashion with rows of sI cages parallel to the interface decomposing simultaneously. The correlations between decomposition times of subsequent layers of the hydrate phase are discussed.     

Soybean Hulls Pretreated Using Thermo-Mechanical Extrusion—Hydrolysis Efficiency,Fermentation Inhibitors,and Ethanol Yield

Soybean hulls were subjected to thermo-mechanical extrusion pretreatment at various in-barrel moisture contents and screw speeds. Extrusion degraded the lignocellulosic structure and enhanced enzymatic hydrolysis of soybean hulls, with up to 155% increase in glucose yield as compared to untreated substrate. Greater glucose yields were observed at higher in-barrel moistures (45% and 50%) and lower screw speed (280 and 350 rpm). Maximum 74% cellulose to glucose conversion resulted from using a two-enzyme cocktail consisting of cellulase and β-glucosidase. Conversion increased to 87% when a three-enzyme cocktail having a cell wall degrading enzyme complex was used for hydrolysis. Fermentation inhibitors, such as furfural, 5-(hydroxymethyl)-2-furaldehyde (HMF), and acetic acid, were found in the extrusion pretreated soybean hulls and hydrolysate. However, their concentrations were below the known thresholds for inhibition. Fermentation of hydrolysate by Saccharomyces cerevisiae led to high yields of ethanol, with concentration ranging from 13.04 to 15.44 g/L.     

Investigation of the full configuration interaction quantum Monte Carlo method using homogeneous electron gas models

Using the homogeneous electron gas (HEG) as a model, we investigate the sources of error in the "initiator" adaptation to full configuration interaction quantum Monte Carlo (i-FCIQMC), with a view to accelerating convergence. In particular, we find that the fixed-shift phase, where the walker number is allowed to grow slowly, can be used to effectively assess stochastic and initiator error. Using this approach we provide simple explanations for the internal parameters of an i-FCIQMC simulation. We exploit the consistent basis sets and adjustable correlation strength of the HEG to analyze properties of the algorithm, and present finite basis benchmark energies for N = 14 over a range of densities 0.5 ≤ r(s) ≤ 5.0 a.u. A single-point extrapolation scheme is introduced to produce complete basis energies for 14, 38, and 54 electrons. It is empirically found that, in the weakly correlated regime, the computational cost scales linearly with the plane wave basis set size, which is justifiable on physical grounds. We expect the fixed-shift strategy to reduce the computational cost of many i-FCIQMC calculations of weakly correlated systems. In addition, we provide benchmarks for the electron gas, to be used by other quantum chemical methods in exploring periodic solid state systems.     

Firefly algorithm with chaos

A recently developed metaheuristic optimization algorithm, firefly algorithm (FA), mimics the social behavior of fireflies based on the flashing and attraction characteristics of fireflies. In the present study, we will introduce chaos into FA so as to increase its global search mobility for robust global optimization. Detailed studies are carried out on benchmark problems with different chaotic maps. Here, 12 different chaotic maps are utilized to tune the attractive movement of the fireflies in the algorithm. The results show that some chaotic FAs can clearly outperform the standard FA.