用计算模拟方法研究Al_3P_4O_(16)~(3-)计量比的二维层状磷酸铝中的有机胺模板能力

将快速 Monte Carlo方法与分子动力学方法相结合 ,研究了不同种类有机分子在 Al3P4 O3- 1 6 计量比的二维层状磷酸铝形成中的模板能力 .依据主 -客体之间非键相互作用能 (包括范德华能、氢键能和库仑能 ) ,可合理地解释已知实验现象 ,并能有效地预测出适于形成某一特定无机层结构的有机胺模板剂 .通过选择理论预测的有机胺分子作为模板剂 ,成功地合成了二维层状磷酸铝化合物 Al3P4 O1 6 · 1 .5 H3NC6 H1 0 NH3.     

Computer experiments on crystalline nylons: structural analysis of nylons with large aliphatic segments

A theoretical study based on force-field calculations has been performed to investigate the structural preferences of crystalline even nylons n with large and very large aliphatic segments. Atomistic energy calculations and Monte Carlo simulations were carried out considering the conventional and forms of nylons 10, 12, 18, 24, and 32. Results indicated that the form is the most favored for nylons 10, 12, 18, and 24. However, the structure was unstable for nylon 32, a polymer in which the density of hydrogen bonds is almost negligible. In this case, the arrangement is energetically more favored than the one.     

水煤气变换反应的微观动力学分析──Cu基催化剂上H_2S中毒失活原因的 Monte Carlo 模拟研究

用ＢＯＣＭＰ方法及ＭｏｎｔｅＣａｒｌｏ模拟技术，对Ｈ２Ｓ导致Ｃｕ基催化剂失活的原因进行了计算分析。研究结果表明，当原料气中存在Ｈ２Ｓ时，ＷＧＳ反应的活化能明显高于无Ｈ２Ｓ时的活化能，随着表面Ｈ２Ｓ浓度的增大（θ＝０，０．１０，０．２５），反应的活化能也逐渐变大（其大小比为１∶１．３４∶３．３），究其原因可归结为Ｈ２Ｓ的存在使得反应物分子的吸附热减小，从而使Ｈ２Ｏ的解离吸附（ＷＧＳ反应的速控步骤）活化能增大。     

Reflection of electrons and photons from solids bombarded by 0.1 - to 100-MeV electrons

Electron and photon reflection ratios (in number and energy) for absorbers bombarded by electrons have been computed with the ITS Monte Carlo system version 3. Electrons of energies from 0.1 to 100 MeV have been assumed normally incident on an effectively semi-infinite absorber. The absorbers considered are elemental solids of atomic numbers from 4 to 92. The data on the electron reflection ratios agree rather well with the experimental data collected from literature except some discrepancies when the number-reflection ratio is small. For photons, the number-reflection ratio increases with increasing energy, but the energy-reflection ratio shows a maximum around 10 MeV. Empirical equations for the electron reflection ratios and the photon energy-reflection ratio are given (for electrons, graphs only).     

Dynamics of the formation of two-dimensional ordered structures

The growth of ordered domains in lattice gas models, which occurs after the system is quenched from infinite temperature to a state below the critical temperatureT _c, is studied by Monte Carlo simulation. For a square lattice with repulsion between nearest and next-nearest neighbors, which in equilibrium exhibits fourfold degenerate (2×1) superstructures, the time-dependent energy E(t), domain size L(t), and structure functionS(q, t) are obtained, both for Glauber dynamics (no conservation law) and the case with conserved density (Kawasaki dynamics). At late times the energy excess and halfwidth of the structure factor decrease proportional tot ^–x, whileL(t) t ^x, where the exponent x=1/2 for Glauber dynamics and x1/3 for Kawasaki dynamics. In addition, the structure factor satisfies a scaling lawS(k,t)=t ^2xS(kt^x). The smaller exponent for the conserved density case is traced back to the excess density contained in the walls between ordered domains which must be redistributed during growth. Quenches toT>T _c, T=T_c (where we estimate dynamic critical exponents) andT=0 are also considered. In the latter case, the system becomes frozen in a glasslike domain pattern far from equilibrium when using Kawasaki dynamics. The generalization of our results to other lattices and structures also is briefly discussed.     

Succession of stacking disorder in hard-sphere crystal under gravity by Monte Carlo simulation

We have performed Monte Carlo simulations of hard spheres under gravity. Vertical boundaries are hard walls, which are well separated with each other. On the other hand, the periodic boundary condition is imposed in the horizontal direction. While we previously reported enhancement of crystallinity as well as crystallization due to gravity, we present here the results that demonstrate the succession of a defect. In case that the crystal formed at the bottom of the system includes kinds of stacking disorders for the (0 0 1) growth, twin band structure develops as mediated by a stacking disorder succeeded in the crystal formed in the fluid region which lies on the bottom crystal. In case that the stacking structure along horizontal direction changes from the (1 1 1) stacking to the (0 0 1) stacking, twin band structure in the (0 0 1) stacking region develops as succeeded in the crystal transformed. The twin band structure also becomes large with its upward growth.     

Quasi‐random integration in quantum chemistry: Efficiency,stability, and application to the study of small atoms and molecules constrained in spherical boxes

This project consists of two parts. In the first part, a series of test calculations is performed to verify that the integrals involved in the determination of atomic and molecular properties by standard self‐consistent field (SCF) methods can be obtained through Halton, Korobov, or Hammersley quasi‐random integration procedures. Through these calculations, we confirm that all three methods lead to results that meet the levels of precision required for their use in the calculation of properties of small atoms or molecules at least at a Hartree–Fock level. Moreover, we have ensured that the efficiency of quasi‐random integration methods that we have tested is Halton=Korobov>Hammersley?pseudo‐random. We also find that these results are comparable to those yielded by ordinary Monte Carlo (pseudo‐random) integration, with a calculation effort of two orders of smaller magnitude. The second part, which would not have been possible without the integration method previously analyzed, contains a first study of atoms constrained in spherical boxes through SCF calculations with basis functions adapted to the features of the problem: Slater‐type orbitals (STOs) trimmed by multiplying them by a function that yields 1 for 0 < r < (R‐δ), polynomial values for (R‐δ) < r < R and null for r > R, R being the radius of the box and δ a variationally determined interval. As a result, we obtain a equation of state for electrons of small systems, valid just in the limit of low temperatures, but fairly simple. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The index of ideality of correlation: improvement of models for toxicity to algae

Toxicity to algae is important characteristic of substances from ecologic point of view. The CORAL software (http://www.insilico.eu/coral) gives possibility to build up model of toxicity to algae using data on the molecular architecture and experimental toxicity, without additional data on physicochemical and/or biochemical parameters. Considerable improvement of the model is observed in the case of using the index of ideality of correlation (IIC) in the role of additional criterion of predictive potential. The IIC is calculated with using of the correlation coefficient between experimental and calculated values of endpoint for the calibration set, with taking into account the positive and negative dispersions between experimental and calculated values. The best model calculated with use the IIC is characterized (the validation set) by n?=?50, r²?=?0.947, RMSE?=?0.401 whereas, model calculated without use the IIC is characterized by n?=?50, r²?=?0.805, and RMSE?=?0.539. The suggested models are built up in accordance to five OECD principles.

     

Effect of chemical potential on the computer simulation of hydrogen storage in single walled carbon nanotubes

Hydrogen is a kind of clean, sustainable and renewable energy carrier. Of the problems to be solved for the utilization of hydrogen energy, how to store and transport hydrogen has been given high priority on the research agenda. Recently, carbon nanotubes (CNTs) were reported to be very promising candidates for hydrogen uptake[1], which may have possibility to satisfy the benchmark set by the US Department of Energy (DOE) Hydrogen Plan for fuel cell powered vehicles: a gravimetric density …