Nonlinear mobility of a driven system: Temperature and disorder effects

We consider the dissipative nonlinear dynamics of a model of interacting atoms driven over a substrate potential. The substrate parameters can be suitably tuned in order to introduce disorder effects starting from two geometrically opposed ideal cases: commensurate and incommensurate interfaces. The role of temperature is also investigated through the inclusion of a stochastic force via a Langevin molecular dynamics approach. Here, we focus on the most interesting tribological case of underdamped sliding dynamics. For different values of the chain stiffness, we evaluate the static friction threshold and consider the depinning transition mechanisms as a function of the applied driving force. As experimentally observed in QCM frictional measurements of adsorbed layers, we find that disorder operates differently depending on the starting geometrical configuration. For commensurate interfaces, randomness lowers considerably the chain depinning threshold. On the contrary, for incommensurate mating contacts, disorder favors static pinning destroying the possible frictionless (superlubric) sliding states. Interestingly, thermal and disorder effects strongly influence also the occurrence of parametric resonances inside the chain, capable of converting the kinetic energy of the center-of-mass motion into internal vibrational excitations. We comment on the nature of the different dynamical states and hysteresis (due to system bi-stability) observed at different increasing and decreasing strengths of the external force.     

Optimal conditions for submonolayer nanoisland growth in ion beam assisted deposition

We study the optimal conditions for nanoisland growth in ion beam assisted deposition (IBAD). This situation occurs when adatom islands remain small enough to prevent the onset of three-dimensional growth, while at the same time preventing ion-induced surface erosion. To this end, we develop a rate equation model of IBAD, which embodies continuous deposition of adatoms and creation of vacancies, recombination of vacancies at adatom island edges, as well as recombination of adatoms at vacancy island edges. These rate equations are solved by numerical simulations based on the particle coalescence method. To determine the optimal growth condition, we find the largest mean size of the vacancy islands leading to their survival. We show that at this onset between the rough and smooth layer-by-layer growth regimes there is a simple exponential relation between the largest size of the vacancy islands and the external control parameters of the growth.     

Height and roughness distributions in thin films with Kardar-Parisi-Zhang scaling

We study height and roughness distributions of films grown with discrete Kardar-Parisi-Zhang (KPZ) models in a small time regime which is expected to parallel the typical experimental conditions. Those distributions are measured with square windows of sizes 8 ? r ? 128 gliding through a much larger surface. Results for models with weak finite-size corrections indicate that the absolute value of the skewness and the value of the kurtosis of height distributions converge to 0.2 ? ∣S∣ ? 0.3 and 0 ? Q ? 0.5, respectively. Despite the low accuracy of these results, they give additional support to a recent claim of KPZ scaling in oligomer films. However, there are significant finite-size effects in the scaled height distributions of models with large local slopes, such as ballistic deposition, which suggests that comparison of height distributions must not be used to rule out KPZ scaling. On the other hand, roughness distributions of the same models show good data collapse, with negligible dependence on time and window size. The estimates of skewness and kurtosis for roughness distributions are 1.7 ? S ? 2 and 3 ? Q ? 7. A stretched exponential tail was found, which seems to be a particular feature of KPZ systems in 2 + 1 dimensions. Moreover, the KPZ roughness distributions cannot be fitted by those of 1/f^α noise. This study suggests that the roughness distribution is the best option to test KPZ scaling in the growth regime, and provides quantitative data for future comparison with other models or experiments.     

On the flow of thermal defects from the bulk to surfaces

This paper treats flow of defects between bulk and surface sites, as a crystal passes towards equilibrium, for some practical cases. These include the realistic but quite elaborate example in which vacancy flow from the bulk is coupled to surface step edges, acting as sinks, by reaction with adatoms that are believed to dominate transport on metal surfaces. It is shown how surface processes modify the defect flow from the bulk only at short times. Lacking accurate parameters (such as concentrations) for surface defects, a crude modeling of the theoretical results is offered in order to explore likely generic behavior. The model employs a recently described approximate universality of behavior, scaled to the melting temperature, relevant mainly to fcc (1 1 1) surfaces. Under a range of conditions it is the reaction of advacancies with adatoms that provides the important channel for bulk vacancy flow. Adatom flow onto the terraces from surface step edge sinks is the bottleneck to flow above a crossover temperature (depending on step spacing) and equilibrium recombination is the bottleneck below the crossover.     

上海证券市场回报率分析——基于ARIMA模型和GARCH模型

本文分析中国上海证券市场回报率。分别通过APdMA模型和GARCH模型，发现若用APdMA模型分析和建立时间序列模型，一次自回归项是不够的，需要高次项，在大多数情形，若运用GARCH模型，则GARCH（1，1）就能够很好的拟合数据。     

Dynamics of surface reactions: from “hot” precursors to conventional thermal activation

Heterogeneous catalytic reactions usually contain steps, e.g., adsorption, generating energy-rich (“hot”) precursors for other steps. With increasing the rate of energy relaxation one can observe a transition from the reaction regime dominated by “hot” precursors to conventional thermal activation. To illustrate this transition in detail and to show what may happen in various situations, I present an analytical model based on the Fokker-Planck equation for energy relaxation.     

Monte Carlo simulations of oscillations, chaos and pattern formation in heterogeneous catalytic reactions

Experimental studies employing surface science methods indicate that kinetic oscillations, chaos, and pattern formation in heterogeneous catalytic reactions often result from the interplay of rapid chemical reaction steps and relatively slow complementary processes such as oxide formation or adsorbate-induced surface restructuring. In general, the latter processes should be analysed in terms of theory of phase transitions. Therefore, the conventional mean-field reaction–diffusion equations widely used to describe oscillations in homogeneous reactions are strictly speaking not applicable. Under such circumstances, application of the Monte Carlo method becomes almost inevitable. In this review, we discuss the advantages and limitations of employing this technique and show what can be achieved in this way. Attention is focused on Monte Carlo simulations of CO oxidation on (1 0 0) and (1 1 0) single-crystal Pt and polycrystal Pt, Pd and Ir surfaces and of NO reduction by CO and H₂ on Pt(1 0 0). CO oxidation on supported nanometre-sized catalyst particles and NO reduction on composite catalysts are also discussed. The results show that with current computer facilities the MC technique has become an effective tool for analysing temporal oscillations and pattern formation on the nanometre scale in catalytic reactions occurring on both single crystals and supported particles.     

Monte Carlo simulations of catalytic CO oxidation on fractal surfaces of dimension between two and three

We present a method for generating fractal surfaces of dimension between two and three. By using the method, five fractal surfaces with dimension 2.262, 2.402, 2.524, 2.631, and 2.771 are created. For each of these surfaces, the reaction of carbon monoxide and oxygen is simulated by using a Monte Carlo method based on the ZGB model [Phys. Rev. Lett. 24 (1986) 2553]. The results show that the catalytic CO oxidation proceeds more efficiently on a surface with higher fractal dimension. It is also found that as the fractal dimension of the surface becomes higher, the first-order kinetic phase transition point (y₂) is shifted to a higher partial pressure of CO. This implies that poisoning of the catalyst surface due to CO segregation sets in at a higher CO partial pressure for surfaces with more complexity.     

Roughness effects on the double-layer charge capacitance: the case of Helmholtz layer induced roughness attenuation

For electrical double layers, the presence of a Helmholtz layer could lead to electrode roughness attenuation. The latter is assumed of self-affine type which is characterized by the roughness amplitude w, the correlation length ξ, and the roughness exponent H. For sufficiently rough metal electrode surfaces (H1 and/or ratios w/ξ0.1) the diffuse/Helmholtz layer interface would not have the same roughness parameters with the metal electrode surface. If the latter is smoothened at lateral length scales smaller than a healing length Λ (ξ), the diffuse charge capacitance decreases and approaches values close to that of the Gouy–Chapman theory for flat electrodes.     

Computational study of cycloaddition reactions on the SiC(1 0 0)-c(2 × 2) surface

Cycloaddition reactions between 1,3-butadiene and the C-terminated SiC(1 0 0)-c(2 × 2) surface have been addressed using quantum-chemical methods. The c(2 × 2) structure consists of ---CC--- bridges between underlayer Si atoms which themselves form Si---Si bonds. Of various possible reaction products, the one formed by a [2 + 4] reaction with the ---CC--- bridge (giving a species resembling 1,4-cyclohexadiene) is the lowest in energy. Density functional calculations for the bare c(2 × 2) surface, using a cluster model with mechanical embedding, gave good agreement with structural parameters obtained in previous fully ab initio studies. Similar calculations for the cycloaddition product and for the transition state gave a reaction energy of −50.3 kcal/mol and an activation energy of +6.1 kcal/mol to form a planar ring structure lying normal to the surface. Detailed results for the frequency and infrared polarization behavior of adsorbate vibrational modes have also been obtained.