Critical Behavior and Lack of Self Averaging in the Dynamics of the Random Potts Model in Two Dimensions

We study the dynamics of the q-state random bond Potts ferromagnet on the square lattice at its critical point by Monte Carlo simulations with single spin-flip dynamics. We concentrate on q=3 and q=24 and find, in both cases, conventional, rather than activated, dynamics. We also look at the distribution of relaxation times among different samples, finding different results for the two q values. For q=3 the relative variance of the relaxation time tau at the critical point is finite. However, for q=24 this appears to diverge in the thermodynamic limit and it is ln(tau) which has a finite relative variance. We speculate that this difference occurs because the transition of the corresponding pure system is second order for q=3 but first order for q=24.Comment: 9 pages, 13 figures, final published versio

Short-time dynamics and magnetic critical behavior of two-dimensional random-bond Potts model

The critical behavior in the short-time dynamics for the random-bond Potts ferromagnet in two-dimensions is investigated by short-time dynamic Monte Carlo simulations. The numerical calculations show that this dynamic approach can be applied efficiently to study the scaling characteristic, which is used to estimate the critical exponents theta, beta/nu and z for the quenched disorered systems from the power-law behavior of the kth moments of magnetizations.Comment: 10 pages, 4 figures Soft Condensed Matte

Learning patient-specific parameters for a diffuse interface glioblastoma model from neuroimaging data

Parameters in mathematical models for glioblastoma multiforme (GBM) tumour growth are highly patient specific. Here we aim to estimate parameters in a Cahn-Hilliard type diffuse interface model in an optimised way using model order reduction (MOR) based on proper orthogonal decomposition (POD). Based on snapshots derived from finite element simulations for the full order model (FOM) we use POD for dimension reduction and solve the parameter estimation for the reduced order model (ROM). Neuroimaging data are used to define the highly inhomogeneous diffusion tensors as well as to define a target functional in a patient specific manner. The reduced order model heavily relies on the discrete empirical interpolation method (DEIM) which has to be appropriately adapted in order to deal with the highly nonlinear and degenerate parabolic PDEs. A feature of the approach is that we iterate between full order solves with new parameters to compute a POD basis function and sensitivity based parameter estimation for the ROM problems. The algorithm is applied using neuroimaging data for two clinical test cases and we can demonstrate that the reduced order approach drastically decreases the computational effort

Watersheds are Schramm-Loewner Evolution curves

We show that in the continuum limit watersheds dividing drainage basins are Schramm-Loewner Evolution (SLE) curves, being described by one single parameter $\kappa$. Several numerical evaluations are applied to ascertain this. All calculations are consistent with SLE$_\kappa$, with $\kappa=1.734\pm0.005$, being the only known physical example of an SLE with $\kappa<2$. This lies outside the well-known duality conjecture, bringing up new questions regarding the existence and reversibility of dual models. Furthermore it constitutes a strong indication for conformal invariance in random landscapes and suggests that watersheds likely correspond to a logarithmic Conformal Field Theory (CFT) with central charge $c\approx-7/2$.Comment: 5 pages and 4 figure

Large-q asymptotics of the random bond Potts model

We numerically examine the large-q asymptotics of the q-state random bond Potts model. Special attention is paid to the parametrisation of the critical line, which is determined by combining the loop representation of the transfer matrix with Zamolodchikov's c-theorem. Asymptotically the central charge seems to behave like c(q) = 1/2 log_2(q) + O(1). Very accurate values of the bulk magnetic exponent x_1 are then extracted by performing Monte Carlo simulations directly at the critical point. As q -> infinity, these seem to tend to a non-trivial limit, x_1 -> 0.192 +- 0.002.Comment: 9 pages, no figure

Ignition and chemical kinetics of acrolein-oxygen-argon mixtures behind reflected shock waves

In order to address increasing greenhouse gas emissions, the future fossil fuel shortage and increasingly stringent pollutant emission regulations, a variety of biofuels are being progressively incorporated into conventional transportation fuels. Despite the beneficial impact of biofuels on most regulated pollutants, their combustion induces the increase of a variety of aldehydes that are being considered for specific regulations due to their high toxicity. One of the most hazardous aldehyde compounds is acrolein, C_2H_3CHO. Despite its high toxicity and increased formation during bioalcohol and biodiesel combustion, no experimental data are available for acrolein combustion. In the present study, we have investigated the ignition of acrolein–oxygen–argon mixtures behind reflected shock wave using three simultaneous emission diagnostics monitoring OH∗, CH∗ and CO_2∗. Experiments were performed over a range of conditions: Φ = 0.5–2; T_5 = 1178–1602 K; and P_5 = 173–416 kPa. A tentative detailed reaction model, which includes sub-mechanisms for the three measured excited species, was developed to describe the high-temperature chemical kinetics of acrolein oxidation. Reasonable agreement was found between the model prediction and experimental data

The Random-bond Potts model in the large-q limit

We study the critical behavior of the q-state Potts model with random ferromagnetic couplings. Working with the cluster representation the partition sum of the model in the large-q limit is dominated by a single graph, the fractal properties of which are related to the critical singularities of the random Potts model. The optimization problem of finding the dominant graph, is studied on the square lattice by simulated annealing and by a combinatorial algorithm. Critical exponents of the magnetization and the correlation length are estimated and conformal predictions are compared with numerical results.Comment: 7 pages, 6 figure

Nonequilibrium critical dynamics of the two-dimensional Ising model quenched from a correlated initial state

The universality class, even the order of the transition, of the two-dimensional Ising model depends on the range and the symmetry of the interactions (Onsager model, Baxter-Wu model, Turban model, etc.), but the critical temperature is generally the same due to self-duality. Here we consider a sudden change in the form of the interaction and study the nonequilibrium critical dynamical properties of the nearest-neighbor model. The relaxation of the magnetization and the decay of the autocorrelation function are found to display a power law behavior with characteristic exponents that depend on the universality class of the initial state.Comment: 6 pages, 5 figures, submitted to Phys. Rev.

Symmetry relation for multifractal spectra at random critical points

Random critical points are generically characterized by multifractal properties. In the field of Anderson localization, Mirlin, Fyodorov, Mildenberger and Evers [Phys. Rev. Lett 97, 046803 (2006)] have proposed that the singularity spectrum $f(\alpha)$ of eigenfunctions satisfies the exact symmetry $f(2d-\alpha)=f(\alpha)+d-\alpha$ at any Anderson transition. In the present paper, we analyse the physical origin of this symmetry in relation with the Gallavotti-Cohen fluctuation relations of large deviation functions that are well-known in the field of non-equilibrium dynamics: the multifractal spectrum of the disordered model corresponds to the large deviation function of the rescaling exponent $\gamma=(\alpha-d)$ along a renormalization trajectory in the effective time $t=\ln L$. We conclude that the symmetry discovered on the specific example of Anderson transitions should actually be satisfied at many other random critical points after an appropriate translation. For many-body random phase transitions, where the critical properties are usually analyzed in terms of the multifractal spectrum $H(a)$ and of the moments exponents X(N) of two-point correlation function [A. Ludwig, Nucl. Phys. B330, 639 (1990)], the symmetry becomes $H(2X(1) -a)= H(a) + a-X(1)$, or equivalently $\Delta(N)=\Delta(1-N)$ for the anomalous parts $\Delta(N) \equiv X(N)-NX(1)$. We present numerical tests in favor of this symmetry for the 2D random $Q-$state Potts model with various $Q$.Comment: 15 pages, 3 figures, v2=final versio