Diffusion on a solid surface: anomalous is normal

We present a numerical study of classical particles diffusing on a solid surface. The particles motion is modeled by an underdamped Langevin equation with ordinary thermal noise. The particle-surface interaction is described by a periodic or a random two-dimensional potential. The model leads to a rich variety of different transport regimes, some of which correspond to anomalous diffusion such as has recently been observed in experiments and Monte Carlo simulations. We show that this anomalous behavior is controlled by the friction coefficient and stress that it emerges naturally in a system described by ordinary canonical Maxwell-Boltzmann statistics.

Proton emission off nuclei induced by kaons in flight

We study the (K-, p) reaction on nuclei with a 1 GeV/c momentum kaon beam, paying special attention to the region of emitted protons having kinetic energy above 600 MeV, which was used to claim a deeply attractive kaon nucleus optical potential. Our model describes the nuclear reaction in the framework of a local density approach and the calculations are performed following two different procedures: one is based on a many-body method using the Lindhard function and the other is based on a Monte Carlo simulation. The simulation method offers flexibility to account for processes other than kaon quasielastic scattering, such as K- absorption by one and two nucleons, producing hyperons, and allows consideration of final-state interactions of the K-, the p, and all other primary and secondary particles on their way out of the nucleus, as well as the weak decay of the produced hyperons into pi N. We find a limited sensitivity of the cross section to the strength of the kaon optical potential. We also show a serious drawback in the experimental setup-the requirement for having, together with the energetic proton, at least one charged particle detected in the decay counter surrounding the target-as we find that the shape of the original cross section is appreciably distorted, to the point of invalidating the claims made in the experimental paper on the strength of the kaon nucleus optical.

Vortices in atomic Bose-Einstein condensates in the large-gas-parameter region

In this work we compare the results of the Gross-Pitaevskii and modified Gross-Pitaevskii equations with ab initio variational Monte Carlo calculations for Bose-Einstein condensates of atoms in axially symmetric traps. We examine both the ground state and excited states having a vortex line along the z axis at high values of the gas parameter and demonstrate an excellent agreement between the modified Gross-Pitaevskii and ab initio Monte Carlo methods, both for the ground and vortex states.

Does introduction of thresholds in decision aids benefit the patient?: Comparison between findings-based and threshold-based diagnostic decision aids.

PURPOSE: To assess how different diagnostic decision aids perform in terms of sensitivity, specificity, and harm. METHODS: Four diagnostic decision aids were compared, as applied to a simulated patient population: a findings-based algorithm following a linear or branched pathway, a serial threshold-based strategy, and a parallel threshold-based strategy. Headache in immune-compromised HIV patients in a developing country was used as an example. Diagnoses included cryptococcal meningitis, cerebral toxoplasmosis, tuberculous meningitis, bacterial meningitis, and malaria. Data were derived from literature and expert opinion. Diagnostic strategies' validity was assessed in terms of sensitivity, specificity, and harm related to mortality and morbidity. Sensitivity analyses and Monte Carlo simulation were performed. RESULTS: The parallel threshold-based approach led to a sensitivity of 92% and a specificity of 65%. Sensitivities of the serial threshold-based approach and the branched and linear algorithms were 47%, 47%, and 74%, respectively, and the specificities were 85%, 95%, and 96%. The parallel threshold-based approach resulted in the least harm, with the serial threshold-based approach, the branched algorithm, and the linear algorithm being associated with 1.56-, 1.44-, and 1.17-times higher harm, respectively. Findings were corroborated by sensitivity and Monte Carlo analyses. CONCLUSION: A threshold-based diagnostic approach is designed to find the optimal trade-off that minimizes expected harm, enhancing sensitivity and lowering specificity when appropriate, as in the given example of a symptom pointing to several life-threatening diseases. Findings-based algorithms, in contrast, solely consider clinical observations. A parallel workup, as opposed to a serial workup, additionally allows for all potential diseases to be reviewed, further reducing false negatives. The parallel threshold-based approach might, however, not be as good in other disease settings.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Malaria haplotype frequency estimation.

We present a Bayesian approach for estimating the relative frequencies of multi-single nucleotide polymorphism (SNP) haplotypes in populations of the malaria parasite Plasmodium falciparum by using microarray SNP data from human blood samples. Each sample comes from a malaria patient and contains one or several parasite clones that may genetically differ. Samples containing multiple parasite clones with different genetic markers pose a special challenge. The situation is comparable with a polyploid organism. The data from each blood sample indicates whether the parasites in the blood carry a mutant or a wildtype allele at various selected genomic positions. If both mutant and wildtype alleles are detected at a given position in a multiply infected sample, the data indicates the presence of both alleles, but the ratio is unknown. Thus, the data only partially reveals which specific combinations of genetic markers (i.e. haplotypes across the examined SNPs) occur in distinct parasite clones. In addition, SNP data may contain errors at non-negligible rates. We use a multinomial mixture model with partially missing observations to represent this data and a Markov chain Monte Carlo method to estimate the haplotype frequencies in a population. Our approach addresses both challenges, multiple infections and data errors.

Slow Dynamics of Ising Models with Energy Barriers

Using Monte Carlo simulations we study the dynamics of three-dimensional Ising models with nearest-, next-nearest-, and four-spin (plaquette) interactions. During coarsening, such models develop growing energy barriers, which leads to very slow dynamics at low temperature. As already reported, the model with only the plaquette interaction exhibits some of the features characteristic of ordinary glasses: strong metastability of the supercooled liquid, a weak increase of the characteristic length under cooling, stretched-exponential relaxation, and aging. The addition of two-spin interactions, in general, destroys such behavior: the liquid phase loses metastability and the slow-dynamics regime terminates well below the melting transition, which is presumably related with a certain corner-rounding transition. However, for a particular choice of interaction constants, when the ground state is strongly degenerate, our simulations suggest that the slow-dynamics regime extends up to the melting transition. The analysis of these models leads us to the conjecture that in the four-spin Ising model domain walls lose their tension at the glassy transition and that they are basically tensionless in the glassy phase.

Diffusion of passive scalars under stochastic convection

The diffusion of passive scalars convected by turbulent flows is addressed here. A practical procedure to obtain stochastic velocity fields with well¿defined energy spectrum functions is also presented. Analytical results are derived, based on the use of stochastic differential equations, where the basic hypothesis involved refers to a rapidly decaying turbulence. These predictions are favorable compared with direct computer simulations of stochastic differential equations containing multiplicative space¿time correlated noise.

Elastic scattering of fast electrons and positrons by atoms

Elastic scattering of relativistic electrons and positrons by atoms is considered in the framework of the static field approximation. The scattering field is expressed as a sum of Yukawa terms to allow the use of various approximations. Accurate phase shifts have been computed by combining Bühring¿s power-series method with the WKB and Born approximations. This combined procedure allows the evaluation of differential cross sections for kinetic energies up to several tens of MeV. Numerical results are used to analyze the validity of several approximate methods, namely the first- and second-order Born approximations and the screened Mott formula, which are frequently adopted as the basis of multiple scattering theories and Monte Carlo simulations of electron and positron transport.

Density matrices and momentum distributions in 3He-4He mixtures

We report variational calculations, in the hypernetted-chain (HNC)-Fermi-HNC scheme, of one-body density matrices and one-particle momentum distributions for 3He-4He mixtures described by a Jastrow correlated wave function. The 4He condensate fractions and the 3He strength poles are examined and compared with the Monte Carlo available results. The agreement has been found to be very satisfactory. Their density dependence is also studied.

Lattice-gas model of orientable molecules: Application to liquid crystals

We have analyzed a two-dimensional lattice-gas model of cylindrical molecules which can exhibit four possible orientations. The Hamiltonian of the model contains positional and orientational energy interaction terms. The ground state of the model has been investigated on the basis of Karl¿s theorem. Monte Carlo simulation results have confirmed the predicted ground state. The model is able to reproduce, with appropriate values of the Hamiltonian parameters, both, a smectic-nematic-like transition and a nematic-isotropic-like transition. We have also analyzed the phase diagram of the system by mean-field techniques and Monte Carlo simulations. Mean-field calculations agree well qualitatively with Monte Carlo results but overestimate transition temperatures.

n=1/4 domain-growth universality class: Crossover to the n=1/2 class

The kinetic domain-growth exponent is studied by Monte Carlo simulation as a function of temperature for a nonconserved order-parameter model. In the limit of zero temperature, the model belongs to the n=(1/4 slow-growth unversality class. This is indicative of a temporal pinning in the domain-boundary network of mixed-, zero-, and finite-curvature boundaries. At finite temperature the growth kinetics is found to cross over to the Allen-Cahn exponent n=(1/2. We obtain that the pinning time of the zero-curvature boundary decreases rapidly with increasing temperature.

Cut-off importance sampling of bole volume

Summary

Fragile-glass behavior of a short-range p-spin model

We propose a short-range generalization of the p-spin interaction spin-glass model. The model is well suited to test the idea that an entropy collapse is at the bottom line of the dynamical singularity encountered in structural glasses. The model is studied in three dimensions through Monte Carlo simulations, which put in evidence fragile glass behavior with stretched exponential relaxation and super-Arrhenius behavior of the relaxation time. Our data are in favor of a Vogel-Fulcher behavior of the relaxation time, related to an entropy collapse at the Kauzmann temperature. We, however, encounter difficulties analogous to those found in experimental systems when extrapolating thermodynamical data at low temperatures. We study the spin-glass susceptibility, investigating the behavior of the correlation length in the system. We find that the increase of the relaxation time is accompanied by a very slow growth of the correlation length. We discuss the scaling properties of off-equilibrium dynamics in the glassy regime, finding qualitative agreement with the mean-field theory.