Single-cluster dynamics for the random-cluster model

We formulate a single-cluster Monte Carlo algorithm for the simulation of the random-cluster model. This algorithm is a generalization of the Wolff single-cluster method for the $q$-state Potts model to non-integer values $q>1$. Its results for static quantities are in a satisfactory agreement with those of the existing Swendsen-Wang-Chayes-Machta (SWCM) algorithm, which involves a full cluster decomposition of random-cluster configurations. We explore the critical dynamics of this algorithm for several two-dimensional Potts and random-cluster models. For integer $q$, the single-cluster algorithm can be reduced to the Wolff algorithm, for which case we find that the autocorrelation functions decay almost purely exponentially, with dynamic exponents $z_{\rm exp} =0.07 (1), 0.521 (7)$, and $1.007 (9)$ for $q=2, 3$, and 4 respectively. For non-integer $q$, the dynamical behavior of the single-cluster algorithm appears to be very dissimilar to that of the SWCM algorithm. For large critical systems, the autocorrelation function displays a range of power-law behavior as a function of time. The dynamic exponents are relatively large. We provide an explanation for this peculiar dynamic behavior.Comment: 7 figures, 4 table

Nebular Spectra of SN 1998bw Revisited: Detailed Study by One and Two Dimensional Models

Refined one- and two-dimensional models for the nebular spectra of the hyper-energetic Type Ic supernova (SN) 1998bw, associated with the gamma-ray burst GRB980425, from 125 to 376 days after B-band maximum are presented. One dimensional, spherically symmetric spectrum synthesis calculations show that reproducing features in the observed spectra, i.e., the sharply peaked [OI] 6300\AA doublet and MgI] 4570\AA emission, and the broad [FeII] blend around 5200\AA, requires the existence of a high-density O-rich core expanding at low velocities (\lsim 8,000 km s$^{-1}$) and of Fe-rich material moving faster than the O-rich material. Synthetic spectra at late phases from aspherical (bipolar) explosion models are also computed with a two-dimensional spectrum synthesis code. The above features are naturally explained by the aspherical model if the explosion is viewed from a direction close to the axis of symmetry ($\sim 30^{\rm o}$), since the aspherical model yields a high-density O-rich region confined along the equatorial axis. By examining a large parameter space (in energy and mass), our best model gives following physical quantities: the kinetic energy $E_{51} \equiv E_{\rm K}/10^{51}$ ergs \gsim 8 - 12 and the main-sequence mass of the progenitor star M_{\rm ms} \gsim 30 - 35 \Msun. The temporal spectral evolution of SN 1998bw also indicates mixing among Fe-, O-, and C-rich regions, and highly clumpy structure.Comment: 38 pages, 22 figures. ApJ, 640 (01 April 2006 issue), in pres

Effective g-factor in Majorana Wires

We use the effective g-factor of subgap states, g*, in hybrid InAs nanowires with an epitaxial Al shell to investigate how the superconducting density of states is distributed between the semiconductor core and the metallic shell. We find a step-like reduction of g* and improved hard gap with reduced carrier density in the nanowire, controlled by gate voltage. These observations are relevant for Majorana devices, which require tunable carrier density and g* exceeding the g-factor of the proximitizing superconductor. Additionally, we observe the closing and reopening of a gap in the subgap spectrum coincident with the appearance of a zero-bias conductance peak

Crossing bonds in the random-cluster model

We derive the scaling dimension associated with crossing bonds in the random-cluster representation of the two-dimensional Potts model, by means of a mapping on the Coulomb gas. The scaling field associated with crossing bonds appears to be irrelevant, on the critical as well as on the tricritical branch. The latter result stands in a remarkable contrast with the existing result for the tricritical O(n) model that crossing bonds are relevant. In order to obtain independent confirmation of the Coulomb gas result for the crossing-bond exponent, we perform a finite-size-scaling analysis based on numerical transfer-matrix calculations.Comment: 2 figure

Nonlocality of Majorana modes in hybrid nanowires

Spatial separation of Majorana zero modes distinguishes trivial from topological midgap states and is key to topological protection in quantum computing applications. Although signatures of Majorana zero modes in tunneling spectroscopy have been reported in numerous studies, a quantitative measure of the degree of separation, or nonlocality, of the emergent zero modes has not been reported. Here, we present results of an experimental study of nonlocality of emergent zero modes in superconductor-semiconductor hybrid nanowire devices. The approach takes advantage of recent theory showing that nonlocality can be measured from splitting due to hybridization of the zero mode in resonance with a quantum dot state at one end of the nanowire. From these splittings as well as anticrossing of the dot states, measured for even and odd occupied quantum dot states, we extract both the degree of nonlocality of the emergent zero mode, as well as the spin canting angles of the nonlocal zero mode. Depending on the device measured, we obtain either a moderate degree of nonlocality, suggesting a partially separated Andreev subgap state, or a highly nonlocal state consistent with a well-developed Majorana modeThis research was supported by Microsoft, the Danish National Research Foundation, the European Commission, and the Spanish Ministry of Economy and Competitiveness through Grants No. FIS2015-65706-P, No. FIS2015-64654-P, and No. FIS2016-80434-P (AEI/FEDER, EU), the Ramón y Cajal programme Grant No. RYC-2011-09345, and the María de Maeztu Programme for Units of Excellence in R&D (Grant No. MDM-2014-0377). C.M.M. acknowledges support from the Villum Foundation. M.-T.D. acknowledges support from State Key Laboratory of High Performance Computing, Chin

Surface and bulk transitions in three-dimensional O(n) models

Using Monte Carlo methods and finite-size scaling, we investigate surface criticality in the O$(n)$ models on the simple-cubic lattice with $n=1$, 2, and 3, i.e. the Ising, XY, and Heisenberg models. For the critical couplings we find $K_{\rm c}(n=2)=0.454 1655 (10)$ and $K_{\rm c}(n=3)= 0.693 002 (2)$. We simulate the three models with open surfaces and determine the surface magnetic exponents at the ordinary transition to be $y_{h1}^{\rm (o)}=0.7374 (15)$, $0.781 (2)$, and $0.813 (2)$ for $n=1$, 2, and 3, respectively. Then we vary the surface coupling $K_1$ and locate the so-called special transition at $\kappa_{\rm c} (n=1)=0.50214 (8)$ and $\kappa_{\rm c} (n=2)=0.6222 (3)$, where $\kappa=K_1/K-1$. The corresponding surface thermal and magnetic exponents are $y_{t1}^{\rm (s)} =0.715 (1)$ and $y_{h1}^{\rm (s)} =1.636 (1)$ for the Ising model, and $y_{t1}^{\rm (s)} =0.608 (4)$ and$y_{h1}^{\rm (s)} =1.675 (1)$ for the XY model. Finite-size corrections with an exponent close to -1/2 occur for both models. Also for the Heisenberg model we find substantial evidence for the existence of a special surface transition.Comment: TeX paper and 10 eps figure

Shear and Bulk Viscosities of a Gluon Plasma in Perturbative QCD: Comparison of Different Treatments for the gg<->ggg Process

The leading order contribution to the shear and bulk viscosities, \eta and \zeta, of a gluon plasma in perturbative QCD includes the gg -> gg (22) process, gg ggg (23) process and multiple scattering processes known as the Landau-Pomeranchuk-Migdal (LPM) effect. Complete leading order computations for \eta and \zeta were obtained by Arnold, Moore and Yaffe (AMY) and Arnold, Dogan and Moore (ADM), respectively, with the inelastic processes computed by an effective g gg gluon splitting. We study how complementary calculations with 22 and 23 processes and a simple treatment to model the LPM effect compare with the results of AMY and ADM. We find that our results agree with theirs within errors. By studying the contribution of the 23 process to \eta, we find that the minimum angle \theta among the final state gluons in the fluid local rest frame has a distribution that is peaked at \theta \sim \sqrt{\alpha_{s}}, analogous to the near collinear splitting asserted by AMY and ADM. However, the average of \theta is much bigger than its peak value, as its distribution is skewed with a long tail. The same \theta behavior is also seen if the 23 matrix element is taken to the soft gluon bremsstrahlung limit in the center-of-mass (CM) frame. This suggests that the soft gluon bremsstrahlung in the CM frame still has some near collinear behavior in the fluid local rest frame. We also generalize our result to a general SU(N_c) pure gauge theory and summarize the current viscosity computations in QCD.Comment: ReVTex 4, 18 pages, 7 figures, accepted version in Phys. Rev.