The adsorption of water in finite carbon pores

Grand canonical Monte Carlo simulations were applied to the adsorption of SPCE model water in finite graphitic pores with different configurations of carbonyl functional groups on only one surface and several pore sizes. It was found that almost all finite pores studied exhibit capillary condensation behaviour preceded by adsorption around the functional groups. Desorption showed the reverse transitions from a filled to a near empty pore resulting in a clear hysteresis loop in all pores except for some of the configurations of the 1.0nm pore. Carbonyl configurations had a strong effect on the filling pressure of all pores except, in some cases, in 1.0nm pores. A decrease in carbonyl neighbour density would result in a higher filling pressure. The emptying pressure was negligibly affected by the configuration of functional groups. Both the filling and emptying pressures increased with increasing pore size but the effect on the emptying pressure was much less. At pressures lower than the pore filling pressure, the adsorption of water was shown to have an extremely strong dependence on the neighbour density with adsorption changing from Type IV to Type III to linear as the neighbour density decreased. The isosteric heat was also calculated for these configurations to reveal its strong dependence on the neighbour density. These results were compared with literature experimental results for water and carbon black and found to qualitatively agree.

Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

Optimum conditions for adsorptive storage

The storage of gases in porous adsorbents, such as activated carbon and carbon nanotubes, is examined here thermodynamically from a systems viewpoint, considering the entire adsorption-desorption cycle. The results provide concrete objective criteria to guide the search for the Holy Grail adsorbent, for which the adsorptive delivery is maximized. It is shown that, for ambient temperature storage of hydrogen and delivery between 30 and 1.5 bar pressure, for the optimum adsorbent the adsorption enthalpy change is 15.1 kJ/mol. For carbons, for which the average enthalpy change is typically 5.8 kJ/mol, an optimum operating temperature of about 115 K is predicted. For methane, an optimum enthalpy change of 18.8 kJ/mol is found, with the optimum temperature for carbons being 254 K. It is also demonstrated that for maximum delivery of the gas the optimum adsorbent must be homogeneous, and that introduction of heterogeneity, such as by ball milling, irradiation, and other means, can only provide small increases in physisorption-related delivery for hydrogen. For methane, heterogeneity is always detrimental, at any value of average adsorption enthalpy change. These results are confirmed with the help of experimental data from the literature, as well as extensive Monte Carlo simulations conducted here using slit pore models of activated carbons as well as atomistic models of carbon nanotubes. The simulations also demonstrate that carbon nanotubes offer little or no advantage over activated carbons in terms of enhanced delivery, when used as storage media for either hydrogen or methane.

Simulation of binary mixture adsorption of methane and CO2 at supercritical conditions in carbons

Knowledge of the adsorption behavior of coal-bed gases, mainly under supercritical high-pressure conditions, is important for optimum design of production processes to recover coal-bed methane and to sequester CO2 in coal-beds. Here, we compare the two most rigorous adsorption methods based on the statistical mechanics approach, which are Density Functional Theory (DFT) and Grand Canonical Monte Carlo (GCMC) simulation, for single and binary mixtures of methane and carbon dioxide in slit-shaped pores ranging from around 0.75 to 7.5 nm in width, for pressure up to 300 bar, and temperature range of 308-348 K, as a preliminary study for the CO2 sequestration problem. For single component adsorption, the isotherms generated by DFT, especially for CO2, do not match well with GCMC calculation, and simulation is subsequently pursued here to investigate the binary mixture adsorption. For binary adsorption, upon increase of pressure, the selectivity of carbon dioxide relative to methane in a binary mixture initially increases to a maximum value, and subsequently drops before attaining a constant value at pressures higher than 300 bar. While the selectivity increases with temperature in the initial pressure-sensitive region, the constant high-pressure value is also temperature independent. Optimum selectivity at any temperature is attained at a pressure of 90-100 bar at low bulk mole fraction of CO2, decreasing to approximately 35 bar at high bulk mole fractions. (c) 2005 American Institute of Chemical Engineers.

Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

Correct procedures for the calculation of heats of adsorption for heterogeneous adsorbents from molecular simulation

Several procedures for calculating the heat of adsorption from Monte Carlo simulations for a heterogeneous adsorbent are presented. Simulations have been performed to generate isotherms for nitrogen at 77 K and methane at 273.15 K in graphitic slit pores of various widths. The procedures were then applied to calculate the heat of adsorption of an activated carbon with an arbitrary pore size distribution. The consistency of the different procedures shows them to be correct in calculating interaction energy contributions to the heat of adsorption. The currently favored procedure for this type of calculation, from the literature, is shown to be incorrect and in serious error when calculating the heat of adsorption of activated carbon.

Optimization of slitlike carbon nanopores for storage of hythane fuel at ambient temperatures

Carbons with slitlike pores can serve as effective host materials for storage of hythane fuel, a bridge between the petrol combustion and hydrogen fuel cells. We have used grand canonical Monte Carlo simulation for the modeling of the hydrogen and methane mixture storage at 293 K and pressure of methane and hydrogen mixture up to 2 MPa. We have found that these pores serve as efficient vessels for the storage of hythane fuel near ambient temperatures and low pressures. We find that, for carbons having optimized slitlike pores of size H congruent to 7 angstrom ( pore width that can accommodate one adsorbed methane layer), and bulk hydrogen mole fraction >= 0.9, the volumetric stored energy exceeds the 2010 target of 5.4 MJ dm(-3) established by the U. S. FreedomCAR Partnership. At the same condition, the content of hydrogen in slitlike carbon pores is congruent to 7% by energy. Thus, we have obtained the composition corresponding to hythane fuel in carbon nanospaces with greatly enhanced volumetric energy in comparison to the traditional compression method. We proposed the simple system with added extra container filled with pure free/adsorbed methane for adjusting the composition of the desorbed mixture as needed during delivery. Our simulation results indicate that light slit pore carbon nanomaterials with optimized parameters are suitable filling vessels for storage of hythane fuel. The proposed simple system consisting of main vessel with physisorbed hythane fuel, and an extra container filled with pure free/adsorbed methane will be particularly suitable for combustion of hythane fuel in buses and passenger cars near ambient temperatures and low pressures.

When "where'' is more important than "when": Birthplace and birthdate effects on the achievement of sporting expertise

In this study, we assessed whether contextual factors related to where or when an athlete is born influence their likelihood of playing professional sport. The birthplace and birth month of all American players in the National Hockey League, National Basketball Association, Major League Baseball, and Professional Golfer's Association, and all Canadian players in the National Hockey League were collected from official websites. Monte Carlo simulations were used to verify if the birthplace of these professional athletes deviated in any systematic way from the official census population distribution, and chi-square analyses were conducted to determine whether the players' birth months were evenly distributed throughout the year. Results showed a birthplace bias towards smaller cities, with professional athletes being over-represented in cities of less than 500,000 and under-represented in cities of 500,000 and over. A birth month/relative age effect (in the form of a distinct bias towards elite athletes being relatively older than their peers) was found for hockey and baseball but not for basketball and golf. Comparative analyses suggested that contextual factors associated with place of birth contribute more influentially to the achievement of an elite level of sport performance than does relative age and that these factors are essentially independent in their influences on expertise development.

Equation of state of a superfluid Fermi gas in the BCS-BEC crossover

We present a theory for a superfluid Fermi gas near the BCS-BEC crossover, including pairing fluctuation contributions to the free energy similar to that considered by Nozieres and Schmitt-Rink for the normal phase. In the strong coupling limit, our theory is able to recover the Bogoliubov theory of a weakly interacting Bose gas with a molecular scattering length very close to the known exact result. We compare our results with recent Quantum Monte Carlo simulations both for the ground state and at finite temperature. Excellent agreement is found for all interaction strengths where simulation results are available.

Performance comparisons between greedy and lagrange algorithms in adaptive MIMO MC-CDMA systems

This paper describes two algorithms for adaptive power and bit allocations in a multiple input multiple output multiple-carrier code division multiple access (MIMO MC-CDMA) system. The first is the greedy algorithm, which has already been presented in the literature. The other one, which is proposed by the authors, is based on the use of the Lagrange multiplier method. The performances of the two algorithms are compared via Monte Carlo simulations. At present stage, the simulations are restricted to a single user MIMO MC-CDMA system, which is equivalent to a MIMO OFDM system. It is assumed that the system operates in a frequency selective fading environment. The transmitter has a partial knowledge of the channel whose properties are measured at the receiver. The use of the two algorithms results in similar system performances. The advantage of the Lagrange algorithm is that is much faster than the greedy algorithm. ©2005 IEEE

Random graph coloring:Statistical physics approach

The problem of vertex coloring in random graphs is studied using methods of statistical physics and probability. Our analytical results are compared to those obtained by exact enumeration and Monte Carlo simulations. We critically discuss the merits and shortcomings of the various methods, and interpret the results obtained. We present an exact analytical expression for the two-coloring problem as well as general replica symmetric approximated solutions for the thermodynamics of the graph coloring problem with p colors and K-body edges. ©2002 The American Physical Society.

Probing three-way interactions in moderated multiple regression:Development and application of a slope difference test

Researchers often use 3-way interactions in moderated multiple regression analysis to test the joint effect of 3 independent variables on a dependent variable. However, further probing of significant interaction terms varies considerably and is sometimes error prone. The authors developed a significance test for slope differences in 3-way interactions and illustrate its importance for testing psychological hypotheses. Monte Carlo simulations revealed that sample size, magnitude of the slope difference, and data reliability affected test power. Application of the test to published data yielded detection of some slope differences that were undetected by alternative probing techniques and led to changes of results and conclusions. The authors conclude by discussing the test's applicability for psychological research. Copyright 2006 by the American Psychological Association.

Discriminating depth in corrugated stereo surfaces: facilitation by a pedestal is explained by removal of uncertainty

With luminance gratings, psychophysical thresholds for detecting a small increase in the contrast of a weak ‘pedestal’ grating are 2–3 times lower than for detection of a grating when the pedestal is absent. This is the ‘dipper effect’ – a reliable improvement whose interpretation remains controversial. Analogies between luminance and depth (disparity) processing have attracted interest in the existence of a ‘disparity dipper’. Are thresholds for disparity modulation (corrugated surfaces), facilitated by the presence of a weak disparity-modulated pedestal? We used a 14-bit greyscale to render small disparities accurately, and measured 2AFC discrimination thresholds for disparity modulation (0.3 or 0.6 c/deg) of a random texture at various pedestal levels. In the first experiment, a clear dipper was found. Thresholds were about 2× lower with weak pedestals than without. But here the phase of modulation (0 or 180 deg) was varied from trial to trial. In a noisy signal-detection framework, this creates uncertainty that is reduced by the pedestal, which thus improves performance. When the uncertainty was eliminated by keeping phase constant within sessions, the dipper effect was weak or absent. Monte Carlo simulations showed that the influence of uncertainty could account well for the results of both experiments. A corollary is that the visual depth response to small disparities is probably linear, with no threshold-like nonlinearity.

Statistical mechanics of distribution networks

This thesis presents an analysis of the stability of complex distribution networks. We present a stability analysis against cascading failures. We propose a spin [binary] model, based on concepts of statistical mechanics. We test macroscopic properties of distribution networks with respect to various topological structures and distributions of microparameters. The equilibrium properties of the systems are obtained in a statistical mechanics framework by application of the replica method. We demonstrate the validity of our approach by comparing it with Monte Carlo simulations. We analyse the network properties in terms of phase diagrams and found both qualitative and quantitative dependence of the network properties on the network structure and macroparameters. The structure of the phase diagrams points at the existence of phase transition and the presence of stable and metastable states in the system. We also present an analysis of robustness against overloading in the distribution networks. We propose a model that describes a distribution process in a network. The model incorporates the currents between any connected hubs in the network, local constraints in the form of Kirchoff's law and a global optimizational criterion. The flow of currents in the system is driven by the consumption. We study two principal types of model: infinite and finite link capacity. The key properties are the distributions of currents in the system. We again use a statistical mechanics framework to describe the currents in the system in terms of macroscopic parameters. In order to obtain observable properties we apply the replica method. We are able to assess the criticality of the level of demand with respect to the available resources and the architecture of the network. Furthermore, the parts of the system, where critical currents may emerge, can be identified. This, in turn, provides us with the characteristic description of the spread of the overloading in the systems.

Reentrant condensation of proteins in solution induced by multivalent counterions

Negatively charged globular proteins in solution undergo a condensation upon adding trivalent counterions between two critical concentrations C* and C**, C*simulations (under these strong electrostatic coupling conditions) support an effective inversion of charge on surface side chains through binding of the multivalent counterions.