关于模拟法描会聚点电场的讨论

本文论证了“ＱＥ－２”型静电场描迹中模拟聚焦电场的电极在设计上有误的观点。     

Interfacial segregation in Ag-Au,Au-Pd,and Cu-Ni alloys: I. (100) surfaces

Atomistic simulations of segregation to (100) free surface in Ag–Au, Au–Pd, and Cu–Ni alloy systems have been performed for a wide range of temperatures and compositions within the solid solution region of these alloy phase diagrams. In addition to the surface segregation profiles, surface free energies, enthalpies, and entropies were determined. These simulations were performed within the framework of the free energy simulation method, in which an approximate free energy functional is minimized with respect to atomic coordinates and atomic site occupation. The effects of the relaxation with respect to either the atomic positions or the atomic concentrations are discussed. For all alloy bulk compositions (0.05 C 0.95) and temperatures (400 T(K) 1,100) examined, Ag, Au, and Cu segregates to the surface in the Ag–Au, Au–Pd, and Cu–Ni alloy systems, respectively. The present results are compared with several theories for segregation. The resultant segregation profiles in Au–Pd and Ag–Au alloys are shown to be in good agreement with an empirical segregation theory, while in Cu–Ni alloys the disagreement in Ni-rich alloys is substantial. The width of the segregation profile is limited to approximately three to four atomic planes. The surface thermodynamic properties depend sensitively on the magnitude of the surface segregation, and some of them are shown to vary linearly with the magnitude of the surface segregation.     

Environment of the P b center at the Si(111)/oxide interface

We assess several models for the environment of the P _b center (Si dangling bond center) at the interface of Si(111) with its oxide. The comparison of hyperfine constants observed with those predicted using large cluster models favors a local structure in which there is an Si-Si bond within the oxide close to the Si dangling bond. Such Si-Si bonds are also suggested by a number of other experiments and are consistent with the reactive layer model proposed to rationalize a range of oxidation studies.     

On the Landau-Teller approximation for energy exchanges with fast degrees of freedom

We revisit the Landau-Teller heuristic approach to adiabatic invariants and, following Rapp, use it to investigate the energy exchanges between the different degrees of freedom, in simple Hamiltonian systems describing the collision of fast rotating or vibrating molecules with a fixed wall. We critically compare the theoretical results with particularly accurate numerical computations (quite small energy exchanges, namely of one part over 10³⁰, are measured).     

Langevin dynamic simulation of hysteresis in a field-swept Landau potential

Numerical simulations are done of Langevin dynamics for a uniform-orderparameter, field-swept Landau model,= –|a/2|m ²+|b/4|m ⁴–mh(t) , to study hysteresis effects. The field is swept at a constant rateh(t)=h(0)+ht. The stochastic jump values of the field {h_J from an initially prepared metastable minimumm(0) are recorded, on passage to a global minimum m(). The results are: (a) The mean jump¯h _J(h) increases (hysteresis loop widens) with h, confirming a previous theoretical criterion based on rate competition between field-sweep and inverse mean first-passage time (FPT); (b) The broad jump distribution(h _J,h) is related to intrinsically large FPT fluctuations ( ²–²)/ ² O(1), and can be quantitatively understood. Possible experimental tests of the ideas are indicated.     

COMPUTER SIMULATION OF A SINGLE POLYMER CHAIN IN DIFFERENT SOLVENTS

In the present paper, the behavior of a single polymer chain under various solvent conditions was modeled by self-avoiding walks (SAW) with nearest neighbors attraction Δε on a simple cubic lattice. Determination of the θ-condition wasbased on the numerical results of the mean square radius of gyration and end-to-end distance. It was found that at the θtemperatue Δε/kT equals -0.27. The exponents a in the Mark-Houwink equation with different interaction parameters areconsistent with the results of experiments: under θ-condition, a = 0.5, and for a good solvent α= 0.74-0.84, respectively.     

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

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

Using molecular dynamics simulations to identify the key factors responsible for chiral recognition by an amino acid-based molecular micelle

Molecular dynamics (MD) simulations were used to investigate the binding of six chiral compounds to the amino acid-based molecular micelle (MM) poly-(sodium undecyl-(L)-leucine-leucine) or poly(SULL). The MM investigated is used as a chiral selector in capillary electrophoresis. The project goal was to characterize the chiral recognition mechanism in these separations and to move toward predictive models to identify the best amino acid-based MM for a given separation. Poly(SULL) was found to contain six binding sites into which chiral compounds could insert. Four sites had similar sizes, shapes, and electrostatic properties. Enantiomers of alprenolol, propranolol, 1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, 1,1′-bi-2-naphthol, chlorthalidone, or lorazepam were separately docked into each binding pocket and MD simulations with the resulting intermolecular complexes were performed. Solvent-accessible surface area calculations showed the compounds preferentially associated with binding sites where they penetrated into the MM core and shielded their non-polar atoms from solvent. Furthermore, with five of the six compounds the enantiomer with the most favorable free energy of MM association also experienced the most favorable intermolecular hydrogen bonding interactions with the MM. This result suggests that stereoselective intermolecular hydrogen bonds play an important role in chiral discrimination in separations using amino acid-based MMs.GRAPHICAL ABSTRACT     

Numerical Study of Pore Sizes Distribution in Gels

We introduce a new numerical technique for the calculation of the pore size distribution in two-dimensional disordered systems. Our method is based on a triangulation technique which allows a closer measurement of pores surface without any morphological hypothesis.In this work, we focus our calculations in simulated gels. Such materials are modeled in two different conditions: by means of the Diffusion-Limited and Reaction-Limited Cluster-cluster Aggregation algorithms, DLCA and RLCA, respectively. In both situations, when the particles concentration decreases, the average pores size increases. The more compact cluster in RLCA, compared with DLCA, is consistent with the pore size distribution we have calculated. The simulated mean pore size is quantitatively in agreement with experimental data from literature.