Self-organized criticality and synchronization in a coupled map lattice model of integrate-and-fire oscillators

We introduce two coupled map lattice models with nonconservative interactions and a continuous nonlinear driving. Depending on both the degree of conservation and the convexity of the driving we find different behaviors, ranging from self-organized criticality, in the sense that the distribution of events (avalanches) obeys a power law, to a macroscopic synchronization of the population of oscillators, with avalanches of the size of the system.

Potts fully frustrated model: Thermodynamics, percolation and dynamics in two dimensions

We consider a Potts model diluted by fully frustrated Ising spins. The model corresponds to a fully frustrated Potts model with variables having an integer absolute value and a sign. This model presents precursor phenomena of a glass transition in the high-temperature region. We show that the onset of these phenomena can be related to a thermodynamic transition. Furthermore, this transition can be mapped onto a percolation transition. We numerically study the phase diagram in two dimensions (2D) for this model with frustration and without disorder and we compare it to the phase diagram of (i) the model with frustration and disorder and (ii) the ferromagnetic model. Introducing a parameter that connects the three models, we generalize the exact expression of the ferromagnetic Potts transition temperature in 2D to the other cases. Finally, we estimate the dynamic critical exponents related to the Potts order parameter and to the energy.

Disordered backgammon model

In this paper we consider an exactly solvable model that displays glassy behavior at zero temperature due to entropic barriers. The new ingredient of the model is the existence of different energy scales or modes associated with different relaxational time scales. Low-temperature relaxation takes place by partial equilibration of successive lower-energy modes. An adiabatic scaling solution, defined in terms of a threshold energy scale e*, is proposed. For such a solution, modes with energy ee* are equilibrated at the bath temperature, modes with ee* remain out of equilibrium, and relaxation occurs in the neighborhood of the threshold e~e*. The model is presented as a toy example to investigate the conditions related to the existence of an effective temperature in glassy systems and its possible dependence on the energy sector is probed by the corresponding observable.

Continuous phase transition in a spin-glass model without time-reversal symmetry

We investigate the phase transition in a strongly disordered short-range three-spin interaction model characterized by the absence of time-reversal symmetry in the Hamiltonian. In the mean-field limit the model is well described by the Adam-Gibbs-DiMarzio scenario for the glass transition; however, in the short-range case this picture turns out to be modified. The model presents a finite temperature continuous phase transition characterized by a divergent spin-glass susceptibility and a negative specific-heat exponent. We expect the nature of the transition in this three-spin model to be the same as the transition in the Edwards-Anderson model in a magnetic field, with the advantage that the strong crossover effects present in the latter case are absent.

Liquid-liquid phase transition for soft-core attractive potential

Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potential¿s parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.

Stationary states and phase diagram for a model of the Gunn effect under realistic boundary conditions

A general formulation of boundary conditions for semiconductor-metal contacts follows from a phenomenological procedure sketched here. The resulting boundary conditions, which incorporate only physically well-defined parameters, are used to study the classical unipolar drift-diffusion model for the Gunn effect. The analysis of its stationary solutions reveals the presence of bistability and hysteresis for a certain range of contact parameters. Several types of Gunn effect are predicted to occur in the model, when no stable stationary solution exists, depending on the value of the parameters of the injecting contact appearing in the boundary condition. In this way, the critical role played by contacts in the Gunn effect is clearly established.

Dynamical properties of the Zhang model of self-organized criticality

Critical exponents of the infinitely slowly driven Zhang model of self-organized criticality are computed for d=2 and 3, with particular emphasis devoted to the various roughening exponents. Besides confirming recent estimates of some exponents, new quantities are monitored, and their critical exponents computed. Among other results, it is shown that the three-dimensional exponents do not coincide with the Bak-Tang-Wiesenfeld [Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 364 (1988)] (Abelian) model, and that the dynamical exponent as computed from the correlation length and from the roughness of the energy profile do not necessarily coincide, as is usually implicitly assumed. An explanation for this is provided. The possibility of comparing these results with those obtained from renormalization group arguments is also briefly addressed.

Asymmetric Little spin-Glass model

We study the static properties of the Little model with asymmetric couplings. We show that the thermodynamics of this model coincides with that of the Sherrington-Kirkpatrick model, and we compute the main finite-size corrections to the difference of the free energy between these two models and to some clarifying order parameters. Our results agree with numerical simulations. Numerical results are presented for the symmetric Little model, which show that the same conclusions are also valid in this case.

Expanding scroll rings and negative tension turbulence in a model of excitable media

Scroll waves in excitable media, described by the Barkley model, are studied. In the parameter region of weak excitability, negative tension of wave filaments is found. It leads to expansion of scroll rings and instability of wave filaments. A circular filament tends to stretch, bend, loop, and produce an expanding tangle that fills up the volume. The filament does not undergo fragmentation before it touches the boundaries. Statistical properties of such Winfree turbulence of scroll waves are numerically investigated.

Surface structure and stability of PdZn and PtZn alloys: Density-functional slab model studies

Geometric parameters of binary (1:1) PdZn and PtZn alloys with CuAu-L10 structure were calculated with a density functional method. Based on the total energies, the alloys are predicted to feature equal formation energies. Calculated surface energies of PdZn and PtZn alloys show that (111) and (100) surfaces exposing stoichiometric layers are more stable than (001) and (110) surfaces comprising alternating Pd (Pt) and Zn layers. The surface energy values of alloys lie between the surface energies of the individual components, but they differ from their composition weighted averages. Compared with the pure metals, the valence d-band widths and the Pd or Pt partial densities of states at the Fermi level are dramatically reduced in PdZn and PtZn alloys. The local valence d-band density of states of Pd and Pt in the alloys resemble that of metallic Cu, suggesting that a similar catalytic performance of these systems can be related to this similarity in the local electronic structures.

Probe-specific mixed-model approach to detect copy number differences using multiplex ligation-dependent probe amplification (MLPA)

Background: MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions. In this paper, we propose a method for the normalization procedure based on a non-linear mixed-model, as well as a new approach for determining the statistical significance of altered probes based on linear mixed-model. This method establishes a threshold by using different tolerance intervals that accommodates the specific random error variability observed in each test sample.Results: Through simulation studies we have shown that our proposed method outperforms two existing methods that are based on simple threshold rules or iterative regression. We have illustrated the method using a controlled MLPA assay in which targeted regions are variable in copy number in individuals suffering from different disorders such as Prader-Willi, DiGeorge or Autism showing the best performace.Conclusion: Using the proposed mixed-model, we are able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.