A markov chain Monte Carlo method for estimating population mixing using Y-chromosome markers: Mixing of the Han people in China

We present a new approach for estimating mixing between populations based on non-recombining markers, specifically Y-chromosome microsatellites. A Markov chain Monte Carlo (MCMC) Bayesian statistical approach is used to calculate the posterior probability

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

Bayesian inference of population expansions in domestic bovines

< p > The past population dynamics of four domestic and one wild species of bovine were estimated using Bayesian skyline plots, a coalescent Markov chain Monte Carlo method that does not require an assumed parametric model of demographic history. Four dom

Thermodynamic properties of Cu-3(CO3)(2)(OH)(2) and the spin-1/2 distorted diamond Heisenberg chain

The thermodynamic properties of the spin-1/2 diamond quantum Heisenberg chain model have been investigated by means of the transfer matrix renormalization group (TMRG) method. Considering different crystal structures, by changing the interactions among different spins and the external magnetic fields, we first investigate the magnetic susceptibility, magnetization, and specific heat of the distorted diamond chain as a model of ferrimagnetic spin systems. The susceptibility and the specific heat show different features for different ferromagnetic (F) and antiferromagnetic (AF) interactions and different magnetic fields. A 1/3 magnetization plateau is observed at low temperature in a magnetization curve. Then, we discuss the theoretical mechanism of the double-peak structure of the magnetic susceptibility and the three-peak structure of the specific heat of the compound Cu-3(CO3)(2)(OH)(2), on which an elegant measurement was performed by Kikuchi [Phys. Rev. Lett. 94, 227201 (2005)]. Our computed results are consistent with the main characteristics of the experimental data. Meanwhile, we find that the double-peak structure of susceptibility can be found in several different kinds of spin interactions in the diamond chain. Moreover, a three-peak behavior is observed in the TMRG results of magnetic susceptibility. In addition, we perform calculations relevant for some experiments and explain the characteristics of these materials. (c) 2007 American Institute of Physics.

Loschmidt echo and Berry phase of a quantum system coupled to an XY spin chain: Proximity to a quantum phase transition

We study the Loschmidt echo (LE) of a coupled system consisting of a central spin and its surrounding environment described by a general XY spin-chain model. The quantum dynamics of the LE is shown to be remarkably influenced by the quantum criticality of the spin chain. In particular, the decaying behavior of the LE is found to be controlled by the anisotropy parameter of the spin chain. Furthermore, we show that due to the coupling to the spin chain, the ground-state Berry phase for the central spin becomes nonanalytical and its derivative with respect to the magnetic parameter lambda in spin chain diverges along the critical line lambda=1, which suggests an alternative measurement of the quantum criticality of the spin chain.

城市交通中车队状态观测器的一种新算法

本文阐述了离散时间点过程理论,时变马尔科夫链及鞅差分序列在城市交通车队状态观测器中的应用。并在此基础上,改进了[5]中的估计算法。用本文提出的非线性最小方差估计算法,对提供的交通状态进行估计,所得结果比[5]中算法精度有明显提高。在大连市某交通干线计算机控制系统中初步应用,取得了令人满意的结果。

Effective elastic modulus of bone-like hierarchical materials

A shear-lag model is used to study the mechanical properties of bone-like hierarchical materials. The relationship between the overall effective modulus and the number of hierarchy level is obtained. The result is compared with that based on the tension-shear chain model and finite element simulation, respectively. It is shown that all three models can be used to describe the mechanical behavior of the hierarchical material when the number of hierarchy levels is small. By increasing the number of hierarchy level, the shear-lag result is consistent with the finite element result. However the tension-shear chain model leads to an opposite trend. The transition point position depends on the fraction of hard phase, aspect ratio and modulus ratio of hard phase to soft phase. Further discussion is performed on the flaw tolerance size and strength of hierarchical materials based on the shear-lag analysis.

Laser light scattering studies of soluble high-performance polyimides: Solution properties and molar mass distributions

Two soluble high-performance polyimides, poly(BCPOBDA/DMMDA) and poly(ODPA/DMMDA), in CHCl3 at 25 degrees C have been studied using laser light scattering. We found that the z-average radius of gyration ([R(g)]) can be scaled to the weight-average molecular weight (M(w)) as [R(g)] (nm) = 4.95 x 10(-2)M(w)(0.52) and [R(g)] (nm) = 1.25 x 10(-2)M(w)(0.66) respectively for poly(BCPOBDA/DMMDA) and poly(ODPA/DMMDA), indicating that poly(ODPA/DMMDA) in CHCl3 at 25 degrees C has a more extended chain conformation than poly(BCPOBDA/DMMDA). Using the wormlike chain model approach, we found that the Flory characteristic ratios (C*) of poly(BCPOBDA/DMMDA) and poly(ODPA/DMMDA) are similar to 20 and similar to 31, respectively, indicating that both of them have a slightly extended chain conformation in comparison with typical flexible polymer chains, such as polystyrene, whose C-infinity is similar to 10. A combination of the weight-average molar mass (M(w)) with the translational diffusion coefficient distributions (G(D)) has led to D (cm(2)/s) = 3.53 x 10(-4)M(-0.579) and D (cm(2)/s) = 4.30 x 10(-4)M(-0.613) respectively for two soluble high-performance polyimides, poly(BCPOBDA/DMMDA) and poly(ODPA/DMMTA), in CHCl3 at 25 degrees C. Using these two calibrations, we have successfully characterized the molar mass distributions of the two polyimides from their corresponding G(D)s. The exponents of these two calibrations further confirm that both of the polyimides have a slightly extended coil chain conformation in CHCl3. The chain flexibility difference between these two polyimides has also been discussed.

Properties of the superheavy nuclei (288)115 and its alpha-decay chain in a generalized liquid-drop model

The properties of the nuclei belonging to the newly observed nuclei starting from (288)115 have been studied with the generalized liquid drop model connected with WKB approximation. The calculated results have been compared with the results of the DDM3Y theory and the experimental data. The half lives of this new alpha decay chain have been well tested from the consistence of the macroscopic, microscopic and the experimental data.

MOLECULAR-GEOMETRY AND CHAIN ENTANGLEMENT - PARAMETERS FOR THE TUBE MODEL

The tube diameter in the reptation model is the distance between a given chain segment and its nearest segment in adjacent chains. This dimension is thus related to the cross-sectional area of polymer chains and the nearest approach among chains, without effects of thermal fluctuation and steric repulsion. Prior calculated tube diameters are much larger, about 5 times, than the actual chain cross-sectional areas. This is ascribed to the local freedom required for mutual rearrangement among neighboring chain segments. This tube diameter concept seems to us to infer a relationship to the corresponding entanglement spacing. Indeed, we report here that the critical molecular weight, M(c), for the onset of entanglements is found to be M(c) = 28 A/([R2]0/M), where A is the chain cross-sectional area and [R2]0 the mean-square end-to-end distance of a freely jointed chain of molecular weight M. The new, computed relationship between the critical number of backbone atoms for entanglement and the chain cross-sectional area of polymers, N(c) = A0,44, is concordant with the cross-sectional area of polymer chains being the parameter controlling the critical entanglement number of backbone atoms of flexible polymers.

Recent advance in the study of ribosomal model.

With the technological developments of cryoelectron microscope, X-ray diffraction and the growing data available on various components of ribosome, some marvelously intricate structural models of the Escherichia coli 70S ribosome have been reconstructed. The picture of the ribosomal model are detailed, including the placement of the mRNA, the arrangement of the A-site and P-site tRNAs and the peptidyltransferase within the interface gap as well as the path of nascent polypeptide chain, which results in a better understanding of the structure and function of ribosome as well as the translational process.

Quantum phase transitions in the S=1/2 distorted diamond chain

By means of the second derivative of the ground-state and first-excited energy, the quantum phase transitions (QPTs) for the distorted diamond chain (DDC) with ferromagnetic and antiferromagnetic frustrated interactions and the trimerized case are investigated, respectively. Our results show the plentiful quantum phases owing to the spin interaction competitions in the model. Meanwhile, by using the transfer-matrix renormalization-group technique, we study the two-site thermal entanglement of the DDC model in the thermodynamic limit for a further understanding of the QPTs.

Quantum critical phenomena in the XY spin chain with the Dzyaloshinski-Moriya interaction (Retracted article. See vol 80, artn 019903, 2009)

We study the effects of the Dzyaloshinski-Moriya (DM) anisotropic interaction on the ground-state properties of the Heisenberg XY spin chain by means of the fidelity susceptibility, order parameter, and entanglement entropy. Our results show that the DM interaction could influence the distribution of the regions of quantum phase transitions and cause different critical regions in the XY spin model. Meanwhile, the DM interaction has effective influence on the degree of entanglement of the system and could be used to increase the entanglement of the spin system.