Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.

Adsorption of Supercritical fluids on graphitised thermal carbon black: Molecular layer structure theory versus grand canonical Monte Carlo simulation

This paper presents a detailed analysis of adsorption of supercritical fluids on nonporous graphitized thermal carbon black. Two methods are employed in the analysis. One is the molecular layer structure theory (MLST), proposed recently by our group, and the other is the grand canonical Monte Carlo (GCMC) simulation. They were applied to describe the adsorption of argon, krypton, methane, ethylene, and sulfur hexafluoride on graphitized thermal carbon black. It was found that the MLST describes all the experimental data at various temperatures well. Results from GCMC simulations describe well the data at low pressure but show some deviations at higher pressures for all the adsorbates tested. The question of negative surface excess is also discussed in this paper.

Adsorption of argon from sub- to supercritical conditions on graphitized thermal carbon black and in graphitic slit pores: A grand canonical Monte Carlo simulation study

In this paper we consider the adsorption of argon on the surface of graphitized thermal carbon black and in slit pores at temperatures ranging from subcritical to supercritical conditions by the method of grand canonical Monte Carlo simulation. Attention is paid to the variation of the adsorbed density when the temperature crosses the critical point. The behavior of the adsorbed density versus pressure (bulk density) shows interesting behavior at temperatures in the vicinity of and those above the critical point and also at extremely high pressures. Isotherms at temperatures greater than the critical temperature exhibit a clear maximum, and near the critical temperature this maximum is a very sharp spike. Under the supercritical conditions and very high pressure the excess of adsorbed density decreases towards zero value for a graphite surface, while for slit pores negative excess density is possible at extremely high pressures. For imperfect pores (defined as pores that cannot accommodate an integral number of parallel layers under moderate conditions) the pressure at which the excess pore density becomes negative is less than that for perfect pores, and this is due to the packing effect in those imperfect pores. However, at extremely high pressure molecules can be packed in parallel layers once chemical potential is great enough to overcome the repulsions among adsorbed molecules. (c) 2005 American Institute of Physics.

Refined method of potential enhancement in the equilibria characterization of activated carbon. Comparison with GCMC and DFT

In this paper, we present a technique for equilibria characterization of activated carbon having slit-shaped pores. This method was first developed by Do (Do, D. D. A new method for the characterisation of micro-mesoporous materials. Presented at the International Symposium on New Trends in Colloid and Interface Science, September 24-26, 1998 Chiba, Japan) and applied by his group and other groups for characterization of pore size distribution (PSD) as well as adsorption equilibria determination of a wide range of hydrocarbons. It is refined in this paper and compared with the grand canonical Monte Carlo (GCMG) simulation and density functional theory (DFT). The refined theory results in a good agreement between the pore filling pressure versus pore width and those obtained by GCMG and DFT. Furthermore, our local isotherms are qualitatively in good agreement with those obtained by the GCMC simulations. The main advantage of this method is that it is about 4 orders of magnitude faster than the GCMC simulations, making it suitable for optimization studies and design purposes. Finally, we apply our method and the GCMG in the derivation of the PSD of a commercial activated carbon. It was found that the PSD derived from our method is comparable to that derived from the GCMG simulations.

GCMC-surface area of carbonaceous materials with N-2 and Ar adsorption as an alternative to the classical BET method

In this paper we apply a new method for the determination of surface area of carbonaceous materials, using the local surface excess isotherms obtained from the Grand Canonical Monte Carlo simulation and a concept of area distribution in terms of energy well-depth of solid–fluid interaction. The range of this well-depth considered in our GCMC simulation is from 10 to 100 K, which is wide enough to cover all carbon surfaces that we dealt with (for comparison, the well-depth for perfect graphite surface is about 58 K). Having the set of local surface excess isotherms and the differential area distribution, the overall adsorption isotherm can be obtained in an integral form. Thus, given the experimental data of nitrogen or argon adsorption on a carbon material, the differential area distribution can be obtained from the inversion process, using the regularization method. The total surface area is then obtained as the area of this distribution. We test this approach with a number of data in the literature, and compare our GCMC-surface area with that obtained from the classical BET method. In general, we find that the difference between these two surface areas is about 10%, indicating the need to reliably determine the surface area with a very consistent method. We, therefore, suggest the approach of this paper as an alternative to the BET method because of the long-recognized unrealistic assumptions used in the BET theory. Beside the surface area obtained by this method, it also provides information about the differential area distribution versus the well-depth. This information could be used as a microscopic finger-print of the carbon surface. It is expected that samples prepared from different precursors and different activation conditions will have distinct finger-prints. We illustrate this with Cabot BP120, 280 and 460 samples, and the differential area distributions obtained from the adsorption of argon at 77 K and nitrogen also at 77 K have exactly the same patterns, suggesting the characteristics of this carbon.

Grand canonical Monte Carlo simulation study of methane adsorption at an open graphite surface and in slitlike carbon pores at 273 K

Grand canonical Monte Carlo (GCMC) simulation was used for the systematic investigation of the supercritical methane adsorption at 273 K on an open graphite surface and in slitlike micropores of different sizes. For both considered adsorption systems the calculated excess adsorption isotherms exhibit a maximum. The effect of the pore size on the maximum surface excess and isosteric enthalpy of adsorption for methane storage at 273 K is discussed. The microscopic detailed picture of methane densification near the homogeneous graphite wall and in slitlike pores at 273 K is presented with selected local density profiles and snapshots. Finally, the reliable pore size distributions, obtained in the range of the microporosity, for two pitch-based microporous activated carbon fibers are calculated from the local excess adsorption isotherms obtained via the GCMC simulation. The current systematic study of supercritical methane adsorption both on an open graphite surface and in slitlike micropores performed by the GCMC summarizes recent investigations performed at slightly different temperatures and usually a lower pressure range by advanced methods based on the statistical thermodynamics.

Structure of saccharose-based carbon and transport of confined fluids: hybrid reverse Monte Carlo reconstruction and simulation studies

We present results of the reconstruction of a saccharose-based activated carbon (CS1000a) using hybrid reverse Monte Carlo (HRMC) simulation, recently proposed by Opletal et al. [1]. Interaction between carbon atoms in the simulation is modeled by an environment dependent interaction potential (EDIP) [2,3]. The reconstructed structure shows predominance of sp(2) over sp bonding, while a significant proportion of sp(3) hybrid bonding is also observed. We also calculated a ring distribution and geometrical pore size distribution of the model developed. The latter is compared with that obtained from argon adsorption at 87 K using our recently proposed characterization procedure [4], the finite wall thickness (FWT) model. Further, we determine self-diffusivities of argon and nitrogen in the constructed carbon as functions of loading. It is found that while there is a maximum in the diffusivity with respect to loading, as previously observed by Pikunic et al. [5], diffusivities in the present work are 10 times larger than those obtained in the prior work, consistent with the larger pore size as well as higher porosity of the activated saccharose carbon studied here.

Gaussian quantum Monte Carlo methods for fermions and bosons

We introduce a new class of quantum Monte Carlo methods, based on a Gaussian quantum operator representation of fermionic states. The methods enable first-principles dynamical or equilibrium calculations in many-body Fermi systems, and, combined with the existing Gaussian representation for bosons, provide a unified method of simulating Bose-Fermi systems. As an application relevant to the Fermi sign problem, we calculate finite-temperature properties of the two dimensional Hubbard model and the dynamics in a simple model of coherent molecular dissociation.

Testing multiple mapping conditioning mixing for Monte Carlo probability density function simulations

Mitarai [Phys. Fluids 17, 047101 (2005)] compared turbulent combustion models against homogeneous direct numerical simulations with extinction/recognition phenomena. The recently suggested multiple mapping conditioning (MMC) was not considered and is simulated here for the same case with favorable results. Implementation issues crucial for successful MMC simulations are also discussed.

Adsorption of argon on homogeneous graphitized thermal carbon black and heterogeneous carbon surface

In this paper we investigate the effects of surface mediation on the adsorption behavior of argon at different temperatures on homogeneous graphitized thermal carbon black and on heterogeneous nongraphitized carbon black surface. The grand canonical Monte Carlo (GCMC) simulation is used to study the adsorption, and its performance is tested against a number of experimental data on graphitized thermal carbon black (which is known to be highly homogeneous) that are available in the literature. The surface-mediation effect is shown to be essential in the correct description of the adsorption isotherm because without accounting for that effect the GCMC simulation results are always greater than the experimental data in the region where the monolayer is being completed. This is due to the overestimation of the fluid–fluid interaction between particles in the first layer close to the solid surface. It is the surface mediation that reduces this fluid–fluid interaction in the adsorbed layers, and therefore the GCMC simulation results accounting for this surface mediation that are presented in this paper result in a better description of the data. This surface mediation having been determined, the surface excess of argon on heterogeneous carbon surfaces having solid–fluid interaction energies different from the graphite can be readily obtained. Since the real heterogeneous carbon surface is not the same as the homogeneous graphite surface, it can be described by an area distribution in terms of the well depth of the solid–fluid energy. Assuming a patchwise topology of the surface with patches of uniform well depth of solid–fluid interaction, the adsorption on a real carbon surface can be determined as an integral of the local surface excess of each patch with respect to the differential area. When this is matched against the experimental data of a carbon surface, we can derive the area distribution versus energy and hence the geometrical surface area. This new approach will be illustrated with the adsorption of argon on a nongraphitized carbon at 87.3 and 77 K, and it is found that the GCMC surface area is different from the BET surface area by about 7%. Furthermore, the description of the isotherm in the region of BET validity of 0.06 to 0.2 is much better with our method than with the BET equation.

Adsorption of flexible n-alkane on graphitized thermal carbon black: analysis of adsorption isotherm by means of GCMC simulation

In this paper, we investigate the suitability of the grand canonical Monte Carlo in the description of adsorption equilibria of flexible n-alkane (butane, pentane and hexane) on graphitized thermal carbon black. Potential model of n-alkane of Martin and Siepmann (J. Phys. Chem. 102 (1998) 2569) is employed in the simulation, and we consider the flexibility of molecule in the simulation. By this we study two models, one is the fully flexible molecular model in which n-alkane is subject to bending and torsion, while the other is the rigid molecular model in which all carbon atoms reside on the same plane. It is found that (i) the adsorption isotherm results of these two models are close to each other, suggesting that n-alkane model behaves mostly as rigid molecules with respect to adsorption although the isotherm for longer chain n-hexane is better described by the flexible molecular model (ii) the isotherms agree very well with the experimental data at least up to two layers on the surface.