Subtleties in data analysis related to the size of critical region

We comment on the analysis of the critical behavior of a layered driven diffusive system recently done by Achahbar and Marro. We discuss why we believe their method of taking the thermodynamic limit and determining the order-parameter exponent leads to unreliable estimates.     

Formation and morphology of latex monolayers. Computer simulation studies

The results of computer simulations of monolayers created from monodisperse latex particles are presented and discussed. Layers are characterized by the normalized coverage,NC (the actual coverage of the surface related to its maximum possible coverage with particles), and by the average number of neighbors,ANN, calculated as the number of particles being in contact with a given one and averaged over all the particles on the surface. Variable parameters used in simulations include: the rate of particles deposition, the probability of lateral movements, the probability of desorption of particles adsorbed on the surface, the probability of covalent immobilization of adsorbed particles, and the on-sphere slip parameter, OSS (characterizing the scattering of a falling particle on the particles being already attached to the surface). Morphology of monolayers is qualitatively characterized by relations betweenANN andNC. It is shown that for all monolayers formed without adhesion (and without repulsion) between the particles adsorbed on the surface the dependence ofANN vs.NC is described by a characteristic master curve (regardless of the values of probabilities of desorption and lateral movements of particles). For the monolayers created including adhesive forces between the adsorbed particles the plots ofANN vs.NC lie above the master curve, while similar plots obtained for the layers made of particles showing various types of repulsive interactions are always placed below it. Thus, the dependencies ofANN vs.NC, derived from computer simulations, can be used for the determination of the character of the interparticle interactions in the real systems.     

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.     

Estimating the complexity of a class of path-following methods for solving linear programs by curvature integrals

In this paper we study a particular class of primal-dual path-following methods which try to follow a trajectory of interior feasible solutions in primal-dual space toward an optimal solution to the primal and dual problem. The methods investigated are so-called first-order methods: each iteration consists of a long step along the tangent of the trajectory, followed by explicit recentering steps to get close to the trajectory again. It is shown that the complexity of these methods, which can be measured by the number of points close to the trajectory which have to be computed in order to achieve a desired gain in accuracy, is bounded by an integral along the trajectory. The integrand is a suitably weighted measure of the second derivative of the trajectory with respect to a distinguished path parameter, so the integral may be loosely called a curvature integral.     

手性氮磷配体与钌的配位过程的研究—原位变温 31P核磁及分子模拟计算

用原位变温 ³¹P NMR和分子模拟研究了手性氮磷配体与金属钌的配位过程。首先确定了配体及其与金属的配合物的化学位移信号分别在δ=-13.0ppm和δ=48.2ppm。配位过程在293～343K温度范围内,相继生成四配位和五配位中间体。四配位体在 ³¹P NMR上对应于30.8ppm和-15.0ppm的两个信号;五配位中间体对应于35.1ppm和-16.5ppm的两个信号;最后在343K,产物中只有六配位的配合物存在,其信号出现在48.2ppm。为了得到各种中间体的结构,用分子模拟方法进行了结构优化计算,并比较了它们的构象能。计算结果表明,四配位中间体有两种构象,能量相差9kcal·mol^-1,它们可能同时存在于平衡状态。而五配位中间体只有一种绝对优势构象,尚有一个磷原子未参与配位。六配位的钌金属配合物的结构优化结果显示,两个氯原子分别位于PNNP原子所构成平面的两侧。其构象能为162.0kcal·mol^-1,其中键角能的贡献是112.5kcal·mol^-1,而非键静电作用是-41.4kcal·mol^-1,这表明分子内的静电吸引力对于形成完全配位产物是十分有利的,但是收敛的配体分子却承受了较大的键角张力。     

Effect of enzyme concentration on the dynamic behavior of a membrane-bound enzyme system

The dynamic behavior of the reaction-diffusion system, composed of glucose oxidase (EC 1.1.3.4) immobilized at a uniform concentration in a membrane, used as a glucose electrode is represented by a diffusion equation with a nonlinear reaction-term in one-dimensional space. The mathematical model is analyzed by computer simulation, that is, numerical integration of the equation under various initial and boundary conditions, to examine the effect of enzyme concentration on the response characteristics (responsiveness and linearity in response) of the electrode. The analysis of the responses of the system to stepwise changes in the boundary value (glucose concentration in simple solution) infers that the enzyme concentration governs the patterns of the spatial distributions of the substrates (glucose and dissolved oxygen) in steady states and transient responses. It is also revealed that the response characteristics of the electrode are optimized with concentration of immobilized enzyme and that the system establishes the steady states at the same spatial distributions of the substrates, regardless of the boundary value. The diffusion of the substrates and the oxygen concentration also have significant effects on the response characteristics of the electrode.     

Computational studies of carbon nanotube-hydrocarbon bond strengths at nanotube ends: effect of link heteroatom and hydrocarbon structure

Semiempirical and density functional electronic structure theory methods were used to study SWNT-X--R bond strengths, where the single-walled carbon nanotube (SWNT) had an armchair or zigzag structure, the link heteroatom X was O, N(H), or S and the hydrocarbon chain R was CH(2)CH(3), CH(OH)CH(3), CHCH(2), or CH(CF(3))CH(3). In all systems the hydrocarbon was bonded to the end of the nanotube. The SWNT-X--R bond (that is, the bond joining the link atom to the hydrocarbon) is more than 0.4 eV stronger for armchair than for zigzag nanotubes with the same diameters, irrespective of whether O, N, or S are used as link atoms or whether OH, C==C, or CF(3) groups are present in the hydrocarbon chain. This raises the possibility for selective manipulation of armchair/zigzag nanotubes using a variety of link atoms and hydrocarbon structures. The SWNT-O--CH(CF(3))CH(3) bond is weaker than the SWNT-O--CH(2)CH(3) bond (for both armchair and zigzag nanotubes), while inclusion of a double bond in the ethyl chain increases the bond strengths. Also, SWNT-S--CH(2)CH(3) and SWNT-N(H)--CH(2)CH(3) bonds are stronger than SWNT-O--CH(2)CH(3) bonds.     

Clusters,metastability, and nucleation: Kinetics of first-order phase transitions

We describe and interpret computer simulations of the time evolution of a binary alloy on a cubic lattice, with nearest neighbor interactions favoring like pairs of atoms. Initially the atoms are arranged at random; the time evolution proceeds by random interchanges of nearest neighbor pairs, using probabilities compatible with the equilibrium Gibbs distribution at temperatureT. For temperatures 0.59T_c, 0.81 T_c, and 0.89T _c, with density of A atoms equal to that in the B-rich phase at coexistence, the density C₁ of clusters ofl A atoms approximately satisfies the following empirical formulas: C₁ w(1 –)³ andC ₁, (1 –)⁴Q₁w¹ (2 l 10). Herew is a parameter and we defineQ _l = _K e ^–E(K) , where the sum goes over all translationally nonequivalentl-particle clusters andE(K) is the energy of formation of the clusterK. Forl > 10,Q ₁ is not known exactly; so we use an extrapolation formulaQ _l Aw _s ^–l l ^– exp(–bl ), wherew _s is the value ofw at coexistence. The same formula (withw > w _s) also fits the observed values of C, (for small values ofl) at densities greater than the coexistence density (forT=0.59T_c): When the supersaturation is small, the simulations show apparently metastable states, a theoretical estimate of whose lifetime is compatible with the observations. For higher supersaturation the system is observed to undergo a slow process of segregation into two coexisting phases (andw therefore changes slowly with time). These results may be interpreted as a more quantitative formulation (and confirmation) of ideas used in standard nucleation theory. No evidence for a spinodal transition is found.Supported by AFOSR Grant No. 73-2430D and by ERDA Contract No. EY-76-C-02-3077*000.