用于Monte Carlo模拟的复杂生物组织模型

Monte Carlo（MC）模拟被广泛应用于光子在生物组织中的传输研究。通常模拟时将生物组织近似为均匀的平板分层介质,当层状生物组织中含有异常物质（如肿瘤细胞等）或正常生物组织为非平板的复杂结构时,其模拟中的组织模型将会有相应的改变。通过探讨这几类生物组织的MC模拟模型,总结并分析模型建立的关键问题,对基于MC模拟的各种生物组织光学检测研究提供了指导。     

Monte Carlo Simulation of Nonionic Surfactants at the Oil-Water Interface

Abstract

A Monte Carlo simulation method has been developed for modelling amphiphiles at an oil-water interface. Properties are calculated for the mixture water, benzene and tetraoxyethylene glycol dodecyl ether.     

Monte Carlo Simulation of Polyampholyte Chains

Abstract

Polyampholyte copolymers containing both positive and negative monomers regularly dispersed along the chain were studied. The Monte Carlo method was used to simulate chains with charged monomers interacting by screened Coulomb potential. The neutral polyampholyte chains collapse due to the attractive electrostatic interactions. The nonneutral chains are in extended conformations due to the repulsive polyelectrolyte effects that dominate the attractive polyampholyte interactions. The results are in good agreement with experiment.     

A Monte Carlo Study of the Relative Stability of Protein Helical Forms

We present simulation results on a simple model to describe the hydrogen bonding in proteins with helical structures. The approximation distinguishes between ! helices, where each amino acid interacts with another one located four residues apart, 3 10 structures, where the number of amino acids in between is three, and the ? arrangement, in which that number is five. We found that the main features of the system are determined by the most stable structure (the ! helix) and that the other type of hydrogen bonds appears just below the denaturation temperature of the peptide. The probability of finding a 3 10 -type bond is greater at the beginning or at the end of the peptide chain, irrespectively of its length, while in short peptides the existence of those bonds increases appreciably the denaturation temperature, promoting stability. On the other hand, the temperature of denaturation decreases with the length of the peptide to reach a value independent of the number of amino acid residues.     

A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition

A new and simple method to determine equilibrium phase transition in adsorption systems exhibiting a hysteresis loop is presented as an alternative to methods such as multiple histogram reweighting, gauge cell method and thermodynamic integration. This method is based on the NVT-grand canonical Monte Carlo mid-density scheme to determine the coexistence chemical potential and coexistence densities of an adsorption system. We illustrate this new scheme with argon and methane adsorption in a number of model solids having slit and cylindrical pores. This method does not have a strong basis on thermodynamic ground, but it does provide a simple heuristic approach that is simpler to understand physically.     

Monte Carlo Simulation of Binary Gas Adsorption in Zeolite Cavities

Abstract

Grand canonical Monte Carlo simulations have been performed for binary adsorption of Lennard-Jones molecules with point multipole moments in zeolite cavities of type X. Fluid-solid electrostatic interactions were taken into account. Phase diagrams and total coverage were calculated for three binaries and compared with experimental measurements. MC simulations gave good agreement with experiment for two mixtures (C₂H₄-CO₂ and CO₂-CH₄) but there were discrepancies between simulation and experiment for the system i-C₄H₁₀-C₂H₄. The dependence of excess Gibbs free energy on the composition and pressure was studied. Negative deviations from ideality are due to energetic heterogeneity and size effects. Unlike liquid-vapor equilibrium, deviations from the Lorentz-Berthelot mixing rules for the adsorbates have little effect upon the phase behavior. Density distributions show that the components compete for the high energy sites inside the cavity; depending on its relative strength of adsorption, one component may be excluded from such positions (CH₄ in CO₂-CH₄), or the two species may share sites inside the cavity (C₂H₄-CO₂).     

Grand Canonical Ensemble Monte Carlo Simulation on a Transputer Array

Five methods are described for the distribution over a 3D Transputer array of the calculation of the pair interaction component of particle energy. The most efficient method, expressed in terms of the time to complete a simulation, depends on the size of the simulation and the Transputer array. This dependence is quantified, with emphasis on Grand Canonical Ensemble Monte Carlo simulation, and yields criteria for the optimum strategy for parallel implementation of GCEMC algorithms. The equations derived are generally applicable, and have implications for the programming of Molecular Dynamics simulations.     

Distance-scaled Force Biased Monte Carlo Simulation for Solutions containing a Strongly Interacting Solute

The force-biased extension of the Metropolis Monte Carlo method [1] improves convergence by sampling moves preferentially along the directions of force (and torque) [2]. For solvated systems it is particularly effective [3] when coupled with the preferential sampling scheme [4] that attempts to move solvents near the solute more frequently. However, in recent force-biased simulations of aqueous ionic solutions [5] some of the water molecules in the vicinity of the solute remained essentially stationary. Only significant reduction in the stepsize produced some accepted moves.     

A Monte Carlo Study of the Argon Fluid Around a Benzene Molecule

Abstract

Monte Carlo simulations have been carried out for argon fluid containing one benzene molecule at the supercritical region. The purpose in this study is to examine the effect of plate-like molecule on the structure of neighboring fluid composed of simple spherical molecules of the system. In the first neighbor shell of argon from the center of benzene molecule, the average potential energy of argon atoms is confirmed to have a large density dependence. This potential energy is relatively large in the high density region. It is found that the spatial distribution of argon fluid is significantly affected by the molecular shape of benzene and it has little direct connection with the attractive interaction between benzene and argon.     

Computer Simulation of Two-Dimensional Continuum Flows by the Direct Simulation Monte Carlo Method

Abstract

Computer simulations using particles are an attractive method to extract microscopic information of flow phenomena [1]. The molecular dynamics (MD) method, in which Newton's equations of motions are integrated, gives the temporal development of the system. In the MD simulation of fluid flows, the computational region is limited to atomistic scales [2]. On the other hand, the direct simulation Monte Carlo (DSMC) method, in which collisions of particles are made on a probabilistic basis, has a potential of treating a realistic system with a macroscopic scale length retaining the atomistic details. The DSMC method provides an efficient way to integrate the Boltzmann equation from the rarefied gas to the near-continuum region. Bird clarified the validity of the DSMC method in the near-continuum flow region [3]. However, the DSMC method has not been applied to the continuum region and compared with the continuum hydrodynamics.     

A Monte Carlo Simulation of Nematic and Discotic Ordering in a Polymeric Liquid Crystal

A cyclic polymeric liquid crystal system is simulated using the Metropolis Monte Carlo method in the NVT ensemble. The polymeric system consists of mesogenic moieties attached by alkyl chain spacers to siloxane ring polymers. In the model, the mesogenic moieties are represented individually by an anisotropic Lennard-Jones potential and the polymer ring is represented solely as a constraint on the relative motions of the attached mesogens. A transition from calamitic ordering to discotic ordering is observed as the ring-mesogen bond is varied from full flexible to rigid.     

方差组分估计方法MIVQUE和REML的模拟比较

利用ＭｏｎｔｅＣａｒｌｏ方法,对４种数据结构进行了ＭＩＶＱＵＥ和ＲＥＭＬ两种方差组分估计方法的模拟比较。方差组分估计所用的模型为奶牛育种中常用的公畜模型,它包括场年季固定效应、公牛组固定效应和公牛随机效应。４种数据结构中最大的有１２８４７个观察值,场年季效应和公牛效应水平数分别为７７８和４７,它与北京市目前可利用的奶牛头胎产奶量记录资料相当。最小的数据结构只有２００个观察值,１４８个场年季和２０头公牛。比较指标为估计值的偏差和方差（理论的或根据１０００次重复模拟所得的经验值）。结果表明,对于较大样本的数据结构,两种方法差异很小,它们间的估计值的相关接近于１,偏差小于真值的１％,方差近似相等。对于较小样本的数据结构,ＭＩＶＱＵＥ则明显优于ＲＥＭＬ。本研究还表明,对于ＲＥＭＬ来说,类似数据结构１的样本已能满足其渐近无偏性和有效性的大样本特性。     

Monte Carlo Simulation of Fluids in Curved Three-dimensional Space

Abstract

This work describes the methods required to perform computer simulations of three-dimensional fluids confined to the surface of a four-dimensional hypersphere. The use of such non-Euclidian spaces, or spherical boundary conditions, is convenient in cases where spatial inhomogeneities occur over length-scales comparable to that of the entire system. The form of the pressure equation in curved space is discussed, and the results of Monte Carlo simulations of hard spheres confined to the surface of a hypersphere are presented. Comparison of the simulation results to the Carnahan-Starling equation of state in flat space provides a basis for determining when curvature effects can be neglected.     

A Monte Carlo method for the simulation of kinetie models

The purpose of this note is to illustrate the feasibility of simulating kinetic systems, such as commonly encountered in photosynthesis research, using the Monte Carlo (MC) method. In this approach, chemical events are considered at the molecular level where they occur randomly and the macroscopic kinetic evolution results from averaging a large number of such events. Their repeated simulation is easily accomplished using digital computing. It is shown that the MC approach is well suited to the capabilities and resources of modern microcomputers. A software package is briefly described and discussed, allowing a simple programming of any kinetic model system and its resolution. The execution is reasonably fast and accurate; it is not subject to such instabilities as found with the conventional analytical approach.Abbreviations MC Monte Carlo - RN random number - PSU photosynthetic unit Dedicated to Prof. L.N.M. Duysens on the occasion of his retirement.     

基于时域的时间分辨荧光光谱的蒙特卡洛模拟

利用蒙特卡洛方法构建了生物组织基于时域的时间分辨荧光光谱的仿真模型,并将应用该模型获得的模拟结果与生物组织的实验光谱进行了比较。结果表明：吸收系数与散射系数分别影响光谱的不同区域;低浓度情况下,模拟结果与实验光谱符合得较好（x^2〈1．2）;高浓度下,实验光谱强度会被浓度效应削弱,以致影响其与模拟光谱的吻合程度。该方法为研究生物组织荧光光谱提供了一种新思路。     

用蒙特卡罗方法模拟光在多层组织中的吸收特性

在讨论目前新颖的组织功能成像打骂能性（例如光声成像）时,光子在组织中的吸收和散射特性是一个很重要的问题,鉴于这一点,本文利用一个多层模型研究了光子在皮肤,脂肪和肌肉组织中的吸收和散射特性,得到了在组织中某一深度处光子在一个平面上的吸收分布,以及在不同吸收系数和散射系数的情况下,光子的反射,吸收和透射几率,结果表明在经过多次散射后,大部分的光子被吸收,在本文的模型中只有7.3%的光子从表面反射（包括镜面反射和漫反射）,还讨论了不同光学参灵敏对参流分布的影响。     

Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures

Abstract

We have developed a new technique, based on the standard Monte Carlo simulation method with Markov chain sampling, in which a set of three dimensional particle configurations are generated that are consistent with the experimentally measured structure factor. A(Q), and radial distribution function, g(r), of a liquid or other disordered system. Consistency is determined by a standard χ² test using the experimental errors. No input potential is required, we present initial results for liquid argon. Since the technique can work directly from the structure factor it promises to be useful for modelling the structures of glasses or amorphous materials. It also has other advantages in multicomponent systems and as a tool for experimental data analysis.     

GATE:一种用于辐射层析成像的蒙特卡罗模拟技术

描述了GATE的主要特点以及实现这些特点的GATE所具有的层级结构。通过测试由GATE软件包自带的GATE PET Benchmark的模拟结果,验证了GATE模拟引擎的正确性以及不同初始随机数下多次模拟实验的稳定性与可靠性。     

Monte Carlo Simulation Study of Adsorption Characteristics in Slit-Like Micropores Under Supercritical Conditions

Abstract

Adsorption characteristics of a solute diluted in supercritical fluids has been investigated by using the Monte Carlo simulation techniques. The Lennard-Jones potential function is used for describing interactions for a model system of CO₂ + benzene in slit-like micropores with infinite graphitic carbon walls. A modified μVT ensemble method with particle exchange proposed by Cracknell, Nicholson and Quirke (1993) is found to be much superior to the conventional μVT ensemble method especially for dense mixtures in a pore. Adsorption isotherms of CO₂ and benzene, in equilibrium with a dilute benzene mixture in CO₂ (mole fraction of benzene = 0.001), are computed by varying pressure, temperature, the benzene–surface interaction potential, and the slitwidth. Adsorption isotherm curve of CO₂ increases with an increase in pressure while that of benzene shows a maximum at a pressure far below the critical pressure of CO₂ and then it decreases with increasing pressure. The decrease in benzene adsorption with increasing pressure is attributable to both the enhanced solubility in supercritical CO₂ and the competitive adsorption of CO₂. The isotherm curves of each component at two temperatures, 313.2 K and 323.2 K, show to cross at a pressure near the critical pressure due to the “density effect” on the chemical potentials of a solute at supercritical fluid conditions. When the interaction between a solute and a surface increases, the adsorption isotherm increases. Narrowing the slitwidth results in the increase in the adsorption of solute since the external potential from two walls becomes deeper.