Numerical Study of Order in a Gauge Glass Model

The XY model with quenched random phase shifts is studied by a T=0 finite size defect energy scaling method in 2d and 3d. The defect energy is defined by a change in the boundary conditions from those compatible with the true ground state configuration for a given realization of disorder. A numerical technique, which is exact in principle, is used to evaluate this energy and to estimate the stiffness exponent $\theta$. This method gives $\theta = -0.36\pm0.013$ in 2d and $\theta = +0.31\pm 0.015$ in 3d, which are considerably larger than previous estimates, strongly suggesting that the lower critical dimension is less than three. Some arguments in favor of these new estimates are given.Comment: 4 pages, 2 figures, revtex. Submitted to Phys. Rev. Let

The metastate approach to thermodynamic chaos

In realistic disordered systems, such as the Edwards-Anderson (EA) spin glass, no order parameter, such as the Parisi overlap distribution, can be both translation-invariant and non-self-averaging. The standard mean-field picture of the EA spin glass phase can therefore not be valid in any dimension and at any temperature. Further analysis shows that, in general, when systems have many competing (pure) thermodynamic states, a single state which is a mixture of many of them (as in the standard mean-field picture) contains insufficient information to reveal the full thermodynamic structure. We propose a different approach, in which an appropriate thermodynamic description of such a system is instead based on a metastate, which is an ensemble of (possibly mixed) thermodynamic states. This approach, modelled on chaotic dynamical systems, is needed when chaotic size dependence (of finite volume correlations) is present. Here replicas arise in a natural way, when a metastate is specified by its (meta)correlations. The metastate approach explains, connects, and unifies such concepts as replica symmetry breaking, chaotic size dependence and replica non-independence. Furthermore, it replaces the older idea of non-self-averaging as dependence on the bulk couplings with the concept of dependence on the state within the metastate at fixed coupling realization. We use these ideas to classify possible metastates for the EA model, and discuss two scenarios introduced by us earlier --- a nonstandard mean-field picture and a picture intermediate between that and the usual scaling/droplet picture.Comment: LaTeX file, 49 page

Quantum Spin Glasses

Ising spin glasses in a transverse field exhibit a zero temperature quantum phase transition, which is driven by quantum rather than thermal fluctuations. They constitute a universality class that is significantly different from the classical, thermal phase transitions. Most interestingly close to the transition in finite dimensions a quantum Griffiths phase leads to drastic consequences for various physical quantities: for instance diverging magnetic susceptibilities are observable over a whole range of transverse field values in the disordered phase.Comment: 10 pages LaTeX (Springer Lecture Notes style file included), 1 eps-figure; Review article for XIV Sitges Conference: Complex Behavior of Glassy System

Time reparametrization group and the long time behaviour in quantum glassy systems

We study the long time dynamics of a quantum version of the Sherrington-Kirkpatrick model. Time reparametrizations of the dynamical equations have a parallel with renormalization group transformations, and within this language the long time behaviour of this model is controlled by a reparametrization group (R$_p$G) fixed point of the classical dynamics. The irrelevance of the quantum terms in the dynamical equations in the aging regime explains the classical nature of the violation of the fluctuation-dissipation theorem.Comment: 4 page

Random quantum magnets with long-range correlated disorder: Enhancement of critical and Griffiths-McCoy singularities

We study the effect of spatial correlations in the quenched disorder on random quantum magnets at and near a quantum critical point. In the random transverse field Ising systems disorder correlations that decay algebraically with an exponent rho change the universality class of the transition for small enough rho and the off-critical Griffiths-McCoy singularities are enhanced. We present exact results for 1d utilizing a mapping to fractional Brownian motion and generalize the predictions for the critical exponents and the generalized dynamical exponent in the Griffiths phase to d>=2.Comment: 4 pages RevTeX, 1 eps-figure include

Short-Range Ising Spin Glass: Multifractal Properties

The multifractal properties of the Edwards-Anderson order parameter of the short-range Ising spin glass model on d=3 diamond hierarchical lattices is studied via an exact recursion procedure. The profiles of the local order parameter are calculated and analysed within a range of temperatures close to the critical point with four symmetric distributions of the coupling constants (Gaussian, Bimodal, Uniform and Exponential). Unlike the pure case, the multifractal analysis of these profiles reveals that a large spectrum of the $\alpha$-H\"older exponent is required to describe the singularities of the measure defined by the normalized local order parameter, at and below the critical point. Minor changes in these spectra are observed for distinct initial distributions of coupling constants, suggesting an universal spectra behavior. For temperatures slightly above T_{c}, a dramatic change in the $F(\alpha)$ function is found, signalizing the transition.Comment: 8 pages, LaTex, PostScript-figures included but also available upon request. To be published in Physical Review E (01/March 97

Critical Behavior and Griffiths-McCoy Singularities in the Two-Dimensional Random Quantum Ising Ferromagnet

We study the quantum phase transition in the two-dimensional random Ising model in a transverse field by Monte Carlo simulations. We find results similar to those known analytically in one-dimension. At the critical point, the dynamical exponent is infinite and the typical correlation function decays with a stretched exponential dependence on distance. Away from the critical point there are Griffiths-McCoy singularities, characterized by a single, continuously varying exponent, z', which diverges at the critical point, as in one-dimension. Consequently, the zero temperature susceptibility diverges for a RANGE of parameters about the transition.Comment: 4 pages RevTeX, 3 eps-figures include

Non-Fermi liquid behavior and Griffiths phase in {\it f}-electron compounds

We study the interplay among disorder, RKKY and Kondo interactions in {\it f}-electron alloys. We argue that the non-Fermi liquid behavior observed in these systems is due to the existence of a Griffiths phase close to a quantum critical point. The existence of this phase provides a unified picture of a large class of materials. We also propose new experiments that can test these ideas.Comment: 4 pages, 1 Figure. NEW version of the original manuscript. A single framework for NFL behavior in different kinds of alloys is presented. Final version finally allowed to appear on the glorious Physical Review Letter

Domain Wall Renormalization Group Study of XY Model with Quenched Random Phase Shifts

The XY model with quenched random disorder is studied by a zero temperature domain wall renormalization group method in 2D and 3D. Instead of the usual phase representation we use the charge (vortex) representation to compute the domain wall, or defect, energy. For the gauge glass corresponding to the maximum disorder we reconfirm earlier predictions that there is no ordered phase in 2D but an ordered phase can exist in 3D at low temperature. However, our simulations yield spin stiffness exponents $\theta_{s} \approx -0.36$ in 2D and $\theta_{s} \approx +0.31$ in 3D, which are considerably larger than previous estimates and strongly suggest that the lower critical dimension is less than three. For the $\pm J$ XY spin glass in 3D, we obtain a spin stiffness exponent $\theta_{s} \approx +0.10$ which supports the existence of spin glass order at finite temperature in contrast with previous estimates which obtain $\theta_{s}< 0$. Our method also allows us to study renormalization group flows of both the coupling constant and the disorder strength with length scale $L$. Our results are consistent with recent analytic and numerical studies suggesting the absence of a re-entrant transition in 2D at low temperature. Some possible consequences and connections with real vortex systems are discussed.Comment: 14 pages, 9 figures, revtex

Quality Assurance Assessment of Intra-Acquisition Diffusion-Weighted and T2-Weighted Magnetic Resonance Imaging Registration and Contour Propagation for Head and Neck Cancer Radiotherapy

BACKGROUND/PURPOSE: Adequate image registration of anatomical and functional magnetic resonance imaging (MRI) scans is necessary for MR-guided head and neck cancer (HNC) adaptive radiotherapy planning. Despite the quantitative capabilities of diffusion-weighted imaging (DWI) MRI for treatment plan adaptation, geometric distortion remains a considerable limitation. Therefore, we systematically investigated various deformable image registration (DIR) methods to co-register DWI and T2-weighted (T2W) images. MATERIALS/METHODS: We compared three commercial (ADMIRE, Velocity, Raystation) and three open-source (Elastix with default settings [Elastix Default], Elastix with parameter set 23 [Elastix 23], Demons) post-acquisition DIR methods applied to T2W and DWI MRI images acquired during the same imaging session in twenty immobilized HNC patients. In addition, we used the non-registered images (None) as a control comparator. Ground-truth segmentations of radiotherapy structures (tumour and organs at risk) were generated by a physician expert on both image sequences. For each registration approach, structures were propagated from T2W to DWI images. These propagated structures were then compared with ground-truth DWI structures using the Dice similarity coefficient and mean surface distance. RESULTS: 19 left submandibular glands, 18 right submandibular glands, 20 left parotid glands, 20 right parotid glands, 20 spinal cords, and 12 tumours were delineated. Most DIR methods tookcase, with the exception of Elastix 23 which took ∼458 s to execute per case. ADMIRE and Elastix 23 demonstrated improved performance over None for all metrics and structures (Bonferroni-corrected p \u3c 0.05), while the other methods did not. Moreover, ADMIRE and Elastix 23 significantly improved performance in individual and pooled analysis compared to all other methods. CONCLUSIONS: The ADMIRE DIR method offers improved geometric performance with reasonable execution time so should be favoured for registering T2W and DWI images acquired during the same scan session in HNC patients. These results are important to ensure the appropriate selection of registration strategies for MR-guided radiotherapy