Efficient calculation of local dose distribution for response modelling in proton and ion beams

We present an algorithm for fast and accurate computation of the local dose distribution in MeV beams of protons, carbon ions or other heavy-charged particles. It uses compound Poisson-process modelling of track interaction and succesive convolutions for fast computation. It can handle mixed particle fields over a wide range of fluences. Since the local dose distribution is the essential part of several approaches to model detector efficiency or cellular response it has potential use in ion-beam dosimetry and radiotherapy.Comment: 9 pages, 3 figure

Optimal transport on wireless networks

We present a study of the application of a variant of a recently introduced heuristic algorithm for the optimization of transport routes on complex networks to the problem of finding the optimal routes of communication between nodes on wireless networks. Our algorithm iteratively balances network traffic by minimizing the maximum node betweenness on the network. The variant we consider specifically accounts for the broadcast restrictions imposed by wireless communication by using a different betweenness measure. We compare the performance of our algorithm to two other known algorithms and find that our algorithm achieves the highest transport capacity both for minimum node degree geometric networks, which are directed geometric networks that model wireless communication networks, and for configuration model networks that are uncorrelated scale-free networks.Comment: 5 pages, 4 figure

Nonequilibrium critical dynamics of the relaxational models C and D

We investigate the critical dynamics of the $n$-component relaxational models C and D which incorporate the coupling of a nonconserved and conserved order parameter S, respectively, to the conserved energy density rho, under nonequilibrium conditions by means of the dynamical renormalization group. Detailed balance violations can be implemented isotropically by allowing for different effective temperatures for the heat baths coupling to the slow modes. In the case of model D with conserved order parameter, the energy density fluctuations can be integrated out. For model C with scalar order parameter, in equilibrium governed by strong dynamic scaling (z_S = z_rho), we find no genuine nonequilibrium fixed point. The nonequilibrium critical dynamics of model C with n = 1 thus follows the behavior of other systems with nonconserved order parameter wherein detailed balance becomes effectively restored at the phase transition. For n >= 4, the energy density decouples from the order parameter. However, for n = 2 and n = 3, in the weak dynamic scaling regime (z_S <= z_rho) entire lines of genuine nonequilibrium model C fixed points emerge to one-loop order, which are characterized by continuously varying critical exponents. Similarly, the nonequilibrium model C with spatially anisotropic noise and n < 4 allows for continuously varying exponents, yet with strong dynamic scaling. Subjecting model D to anisotropic nonequilibrium perturbations leads to genuinely different critical behavior with softening only in subsectors of momentum space and correspondingly anisotropic scaling exponents. Similar to the two-temperature model B the effective theory at criticality can be cast into an equilibrium model D dynamics, albeit incorporating long-range interactions of the uniaxial dipolar type.Comment: Revtex, 23 pages, 5 eps figures included (minor additions), to appear in Phys. Rev.

Evelopment of the model of diagnosis of the risk of bankruptcy

The article presents an overview of foreign and domestic models for the diagnosis of bankruptcy risk, and gives a brief description of them. Also considered the development of our own model of bankruptcy risk diagnostics for Russian enterprise

Critical behavior of a one-dimensional monomer-dimer reaction model with lateral interactions

A monomer-dimer reaction lattice model with lateral repulsion among the same species is studied using a mean-field analysis and Monte Carlo simulations. For weak repulsions, the model exhibits a first-order irreversible phase transition between two absorbing states saturated by each different species. Increasing the repulsion, a reactive stationary state appears in addition to the saturated states. The irreversible phase transitions from the reactive phase to any of the saturated states are continuous and belong to the directed percolation universality class. However, a different critical behavior is found at the point where the directed percolation phase boundaries meet. The values of the critical exponents calculated at the bicritical point are in good agreement with the exponents corresponding to the parity-conserving universality class. Since the adsorption-reaction processes does not lead to a non-trivial local parity-conserving dynamics, this result confirms that the twofold symmetry between absorbing states plays a relevant role in determining the universality class. The value of the exponent $\delta_2$, which characterizes the fluctuations of an interface at the bicritical point, supports the Bassler-Brown's conjecture which states that this is a new exponent in the parity-conserving universality class.Comment: 19 pages, 22 figures, to be published in Phys. Rev

Effects of differential mobility on biased diffusion of two species

Using simulations and a simple mean-field theory, we investigate jamming transitions in a two-species lattice gas under non-equilibrium steady-state conditions. The two types of particles diffuse with different mobilities on a square lattice, subject to an excluded volume constraint and biased in opposite directions. Varying filling fraction, differential mobility, and drive, we map out the phase diagram, identifying first order and continuous transitions between a free-flowing disordered and a spatially inhomogeneous jammed phase. Ordered structures are observed to drift, with a characteristic velocity, in the direction of the more mobile species.Comment: 15 pages, 4 figure

Novel non-equilibrium critical behavior in unidirectionally coupled stochastic processes

Phase transitions from an active into an absorbing, inactive state are generically described by the critical exponents of directed percolation (DP), with upper critical dimension d_c = 4. In the framework of single-species reaction-diffusion systems, this universality class is realized by the combined processes A -> A + A, A + A -> A, and A -> \emptyset. We study a hierarchy of such DP processes for particle species A, B,..., unidirectionally coupled via the reactions A -> B, ... (with rates \mu_{AB}, ...). When the DP critical points at all levels coincide, multicritical behavior emerges, with density exponents \beta_i which are markedly reduced at each hierarchy level i >= 2. This scenario can be understood on the basis of the mean-field rate equations, which yield \beta_i = 1/2^{i-1} at the multicritical point. We then include fluctuations by using field-theoretic renormalization group techniques in d = 4-\epsilon dimensions. In the active phase, we calculate the fluctuation correction to the density exponent for the second hierarchy level, \beta_2 = 1/2 - \epsilon/8 + O(\epsilon^2). Monte Carlo simulations are then employed to determine the values for the new scaling exponents in dimensions d<= 3, including the critical initial slip exponent. Our theory is connected to certain classes of growth processes and to certain cellular automata, as well as to unidirectionally coupled pair annihilation processes. We also discuss some technical and conceptual problems of the loop expansion and their possible interpretation.Comment: 29 pages, 19 figures, revtex, 2 columns, revised Jan 1995: minor changes and additions; accepted for publication in Phys. Rev.

Order in driven vortex lattices in superconducting Nb films with nanostructured pinning potentials

Driven vortex lattices have been studied in a material with strong pinning, such as Nb films. Samples in which natural random pinning coexists with artificial ordered arrays of defects (submicrometric Ni dots) have been fabricated with different geometries (square, triangular and rectangular). Three different dynamic regimes are found: for low vortex velocities, there is a plastic regime in which random defects frustrate the effect of the ordered array; then, for vortex velocities in the range 1-100 m/s, there is a sudden increase in the interaction between the vortex lattice and the ordered dot array, independent on the geometry. This effect is associated to the onset of quasi long range order in the vortex lattice leading to an increase in the overlap between the vortex lattice and the magnetic dots array. Finally, at larger velocities the ordered array-vortex lattice interaction is suppresed again, in agreement with the behavior found in numerical simulations.Comment: 8 text pages + 4 figure

Analytical expressions for stopping-power ratios relevant for accurate dosimetry in particle therapy

In particle therapy, knowledge of the stopping-power ratios (STPRs) of the ion beam for air and water is necessary for accurate ionization chamber dosimetry. Earlier work has investigated the STPRs for pristine carbon ion beams, but here we expand the calculations to a range of ions (1 <= z <= 18) as well as spread out Bragg peaks (SOBPs) and provide a theoretical in-depth study with a special focus on the parameter regime relevant for particle therapy. The Monte Carlo transport code SHIELD-HIT is used to calculate complete particle-fluence spectra which are required for determining STPRs according to the recommendations of the International Atomic Energy Agency (IAEA). We confirm that the STPR depends primarily on the current energy of the ions rather than on their charge z or absolute position in the medium. However, STPRs for different sets of stopping-power data for water and air recommended by the International Commission on Radiation Units & Measurements (ICRU) are compared, including also the recently revised data for water, yielding deviations up to 2% in the plateau region. In comparison, the influence of the secondary particle spectra on the STPR is about two orders of magnitude smaller in the whole region up till the practical range. The gained insights enable us to propose an analytic approximation for the STPR for both pristine and SOBPs as a function of penetration depth, which parametrically depend only on the initial energy and the residual range of the ion, respectively.Comment: 21 pages, 5 figures, fixed bug with figures in v

Experiments in vortex avalanches

Avalanche dynamics is found in many phenomena spanning from earthquakes to the evolution of species. It can be also found in vortex matter when a type II superconductor is externally driven, for example, by increasing the magnetic field. Vortex avalanches associated with thermal instabilities can be an undesirable effect for applications, but "dynamically driven" avalanches emerging from the competition between intervortex interactions and quenched disorder constitute an interesting scenario to test theoretical ideas related with non-equilibrium dynamics. However, differently from the equilibrium phases of vortex matter in type II superconductors, the study of the corresponding dynamical phases - in which avalanches can play a role - is still in its infancy. In this paper we critically review relevant experiments performed in the last decade or so, emphasizing the ability of different experimental techniques to establish the nature and statistical properties of the observed avalanche behavior.Comment: To be published in Reviews of Modern Physics April 2004. 17 page