Kinetic Monte Carlo study of electrochemical growth in a heteroepitaxial system

Structural and kinetic aspects of 2-D irreversible metal deposition under potentiostatic conditions are analyzed by means of dynamic Monte Carlo simulations employing embedded atom potentials for a model system. Three limiting models, all considering adatom diffusion, were employed to describe adatom deposition. The first model (A) considers adatom deposition on any free substrate site on the surface at the same rate. The second model (B) considers adatom deposition only on substrate sites which exhibit no neighboring sites occupied by adatoms. The third model (C) allows deposition at higher rates on sites presenting neighboring sites occupied by adatoms. Under the proper conditions, the coverage (theta) versus time (t) relationship for the three cases can be heuristically fitted to the functional form theta = 1 - exp(-betat(alpha)), where alpha and beta are parameters. We suggest that the value of the parameter alpha can be employed to distinguish experimentally between the three cases. While model A trivially delivers a = 1, models B and C are characterized by alpha 1, respectively.

The application of a Trous wave filtering and Monte Carlo analysis on secis 2001 solar eclipse observations

Eight thousand images of the solar corona were captured during the June 2001 total solar eclipse. New software for the alignment of the images and an automated technique for detecting intensity oscillations using multi-scale wavelet analysis were developed. Large areas of the images covered by the Moon and the upper corona were scanned for oscillations and the statistical properties of the atmospheric effects were determined. The a Trous wavelet transform was used for noise reduction and Monte Carlo analysis as a significance test of the detections. The effectiveness of those techniques is discussed in detail.

Is the pinning of ordinary dislocations in γ-TiAl intrinsic or extrinsic in nature? A combined atomistic and kinetic Monte Carlo approach

We address the question of the observed pinning of 1/2

Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

Nonrigid Structure-From-Motion From 2-D Images Using Markov Chain Monte Carlo

In this paper we present a new method for simultaneously determining three dimensional (3-D) shape and motion of a non-rigid object from uncalibrated two dimensional (2- D) images without assuming the distribution characteristics. A non-rigid motion can be treated as a combination of a rigid rotation and a non-rigid deformation. To seek accurate recovery of deformable structures, we estimate the probability distribution function of the corresponding features through random sampling, incorporating an established probabilistic model. The fitting between the observation and the projection of the estimated 3-D structure will be evaluated using a Markov chain Monte Carlo based expectation maximisation algorithm. Applications of the proposed method to both synthetic and real image sequences are demonstrated with promising results.

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments. (C) 2004 by Radiation Research Society.

Two-dimensional Monte Carlo simulations of H i line formation in massive young stellar object disc winds

Massive young stellar objects (YSOs) are powerful infrared Hi line emitters. It has been suggested that these lines form in an outflow from a disc surrounding the YSO. Here, new two-dimensional Monte Carlo radiative transfer calculations are described which test this hypothesis. Infrared spectra are synthesized for a YSO disc wind model based on earlier hydrodynamical calculations. The model spectra are in qualitative agreement with the observed spectra from massive YSOs, and therefore provide support for a disc wind explanation for the Hi lines. However, there are some significant differences: the models tend to overpredict the Bra/Br? ratio of equivalent widths and produce line profiles which are slightly too broad and, in contrast to typical observations, are double-peaked. The interpretation of these differences within the context of the disc wind picture and suggestions for their resolution via modifications to the assumed disc and outflow structure are discussed. © 2005 RAS.