Wang-Landau study of the 3D Ising model with bond disorder

We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagnetic interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes $L$ in the range $L=8-64$. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.Comment: 7 pages, 7 figures, to be published in Eur. Phys. J.

Evidence for softening of first-order transition in 3D by quenched disorder

We study by extensive Monte Carlo simulations the effect of random bond dilution on the phase transition of the three-dimensional 4-state Potts model which is known to exhibit a strong first-order transition in the pure case. The phase diagram in the dilution-temperature plane is determined from the peaks of the susceptibility for sufficiently large system sizes. In the strongly disordered regime, numerical evidence for softening to a second-order transition induced by randomness is given. Here a large-scale finite-size scaling analysis, made difficult due to strong crossover effects presumably caused by the percolation fixed point, is performed.Comment: LaTeX file with Revtex, 4 pages, 4 eps figure

Crossover and self-averaging in the two-dimensional site-diluted Ising model

Using the newly proposed probability-changing cluster (PCC) Monte Carlo algorithm, we simulate the two-dimensional (2D) site-diluted Ising model. Since we can tune the critical point of each random sample automatically with the PCC algorithm, we succeed in studying the sample-dependent $T_c(L)$ and the sample average of physical quantities at each $T_c(L)$ systematically. Using the finite-size scaling (FSS) analysis for $T_c(L)$, we discuss the importance of corrections to FSS both in the strong-dilution and weak-dilution regions. The critical phenomena of the 2D site-diluted Ising model are shown to be controlled by the pure fixed point. The crossover from the percolation fixed point to the pure Ising fixed point with the system size is explicitly demonstrated by the study of the Binder parameter. We also study the distribution of critical temperature $T_c(L)$. Its variance shows the power-law $L$ dependence, $L^{-n}$, and the estimate of the exponent $n$ is consistent with the prediction of Aharony and Harris [Phys. Rev. Lett. {\bf 77}, 3700 (1996)]. Calculating the relative variance of critical magnetization at the sample-dependent $T_c(L)$, we show that the 2D site-diluted Ising model exhibits weak self-averaging.Comment: 6 pages including 6 eps figures, RevTeX, to appear in Phys. Rev.

The critical amplitude ratio of the susceptibility in the random-site two-dimensional Ising model

We present a new way of probing the universality class of the site-diluted two-dimensional Ising model. We analyse Monte Carlo data for the magnetic susceptibility, introducing a new fitting procedure in the critical region applicable even for a single sample with quenched disorder. This gives us the possibility to fit simultaneously the critical exponent, the critical amplitude and the sample dependent pseudo-critical temperature. The critical amplitude ratio of the magnetic susceptibility is seen to be independent of the concentration $q$ of the empty sites for all investigated values of $q\le 0.25$. At the same time the average effective exponent $\gamma_{eff}$ is found to vary with the concentration $q$, which may be argued to be due to logarithmic corrections to the power law of the pure system. This corrections are canceled in the susceptibility amplitude ratio as predicted by theory. The central charge of the corresponding field theory was computed and compared well with the theoretical predictions.Comment: 6 pages, 4 figure

Ising model on 3D random lattices: A Monte Carlo study

We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to 128,000 approx. 503 sites which are linked together according to the Voronoi/Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and \phi^4 field theory on three-dimensional regular cubic lattices.Comment: 35 pages, LaTex, 8 tables, 8 postscript figure

First Measurement of Z/gamma* Production in Compton Scattering of Quasi-real Photons

We report the first observation of Z/gamma* production in Compton scattering of quasi-real photons. This is a subprocess of the reaction e+e- to e+e-Z/gamma*, where one of the final state electrons is undetected. Approximately 55 pb-1 of data collected in the year 1997 at an e+e- centre-of-mass energy of 183 GeV with the OPAL detector at LEP have been analysed. The Z/gamma* from Compton scattering has been detected in the hadronic decay channel. Within well defined kinematic bounds, we measure the product of cross-section and Z/gamma* branching ratio to hadrons to be (0.9+-0.3+-0.1) pb for events with a hadronic mass larger than 60 GeV, dominated by (e)eZ production. In the hadronic mass region between 5 GeV and 60 GeV, dominated by (e)egamma* production, this product is found to be (4.1+-1.6+-0.6) pb. Our results agree with the predictions of two Monte Carlo event generators, grc4f and PYTHIA.Comment: 18 pages, LaTeX, 5 eps figures included, submitted to Physics Letters

Measurement of the Michel Parameters in Leptonic Tau Decays

The Michel parameters of the leptonic tau decays are measured using the OPAL detector at LEP. The Michel parameters are extracted from the energy spectra of the charged decay leptons and from their energy-energy correlations. A new method involving a global likelihood fit of Monte Carlo generated events with complete detector simulation and background treatment has been applied to the data recorded at center-of-mass energies close to sqrt(s) = M(Z) corresponding to an integrated luminosity of 155 pb-1 during the years 1990 to 1995. If e-mu universality is assumed and inferring the tau polarization from neutral current data, the measured Michel parameters are extracted. Limits on non-standard coupling constants and on the masses of new gauge bosons are obtained. The results are in agreement with the V-A prediction of the Standard Model.Comment: 32 pages, LaTeX, 9 eps figures included, submitted to the European Physical Journal

Inclusive production of charged hadrons and Ks0K_{s}^{0} mesons in photon-photon collisions

The production of charged hadrons and K_s mesons in the collisions of quasi-real photons has been measured using the OPAL detector at LEP. The data were taken at e+e- centre-of-mass energies of 161 and 172 GeV. The differential cross-sections as a function of the transverse momentum and the pseudorapidity of the charged hadrons and K_s mesons have been compared to the leading order Monte Carlo simulations of PHOJET and PYTHIA and to perturbative next-to-leading order (NLO) QCD calculations. The distributions have been measured in the range 10-125 GeV of the hadronic invariant mass W. By comparing the transverse momentum distribution of charged hadrons measured in gamma-gamma interactions with gamma-proton and meson-proton data we find evidence for hard photon interactions in addition to the purely hadronic photon interactions.The production of charged hadrons and K_s mesons in the collisions of quasi-real photons has been measured using the OPAL detector at LEP. The data were taken at e+e- centre-of-mass energies of 161 and 172 GeV. The differential cross-sections as a function of the transverse momentum and the pseudorapidity of the charged hadrons and K_s mesons have been compared to the leading order Monte Carlo simulations of PHOJET and PYTHIA and to perturbative next-to-leading order (NLO) QCD calculations. The distributions have been measured in the range 10-125 GeV of the hadronic invariant mass W. By comparing the transverse momentum distribution of charged hadrons measured in gamma-gamma interactions with gamma-proton and meson-proton data we find evidence for hard photon interactions in addition to the purely hadronic photon interactions

A search for neutral Higgs bosons in the MSSM and models with two scalar field doublets

A search is described for the neutral Higgs bosons h^0 and A^0 predicted by models with two scalar field doublets and, in particular, the Minimal Supersymmetric Standard Model (MSSM). The search in the Z^0 h^0 and h^0 A^0 production channels is based on data corresponding to an integrated luminosity of 25 pb^{-1} from e^+e^- collisions at centre-of-mass energies between 130 and 172GeV collected with the OPAL detector at LEP. The observation of a number of candidates consistent with Standard Model background expectations is used in combination with earlier results from data collected at the Z^0 resonance to set limits on m_h and m_A in general models with two scalar field doublets and in the MSSM. For example, in the MSSM, for tan(beta) > 1, minimal and maximal scalar top quark mixing and soft SUSY-breaking masses of 1 TeV, the 95% confidence level limits m_h > 59.0 GeV and m_A > 59.5 GeV are obtained. For the first time, the MSSM parameter space is explored in a detailed scan.A search is described for the neutral Higgs bosons h^0 and A^0 predicted by models with two scalar field doublets and, in particular, the Minimal Supersymmetric Standard Model (MSSM). The search in the Z^0 h^0 and h^0 A^0 production channels is based on data corresponding to an integrated luminosity of 25 pb^{-1} from e^+e^- collisions at centre-of-mass energies between 130 and 172 GeV collected with the OPAL detector at LEP. The observation of a number of candidates consistent with Standard Model background expectations is used in combination with earlier results from data collected at the Z^0 resonance to set limits on m_h and m_A in general models with two scalar field doublets and in the MSSM. For example, in the MSSM, for tan(beta) > 1, minimal and maximal scalar top quark mixing and soft SUSY-breaking masses of 1 TeV, the 95% confidence level limits m_h > 59.0 GeV and m_A > 59.5 GeV are obtained. For the first time, the MSSM parameter space is explored in a detailed scan

Spin alignment of leading K∗(892)0K^{*}(892)^{0} mesons in hadronic Z0Z^0 decays

Helicity density matrix elements for inclusive K*(892)^0 mesons from hadronic Z^0 decays have been measured over the full range of K^*0 momentum using data taken with the OPAL experiment at LEP. A preference for occupation of the helicity zero state is observed at all scaled momentum x_p values above 0.3, with the matrix element rho_00 rising to 0.66 +/- 0.11 for x_p > 0.7. The values of the real part of the off-diagonal element rho_1-1 are negative at large x_p, with a weighted average value of -0.09 +/- 0.03 for x_p > 0.3, in agreement with new theoretical predictions based on Standard Model parameters and coherent fragmentation of the qq(bar) system from the Z^0 decay. All other helicity density matrix elements measured are consistent with zero over the entire x_p range. The K^*0 fragmentation function has also been measured and the total rate determined to be 0.74 +/- 0.02 +/- 0.02 K*(892)^0 mesons per hadronic Z^0 decay.Helicity density matrix elements for inclusive K*(892)^0 mesons from hadronic Z^0 decays have been measured over the full range of K^*0 momentum using data taken with the OPAL experiment at LEP. A preference for occupation of the helicity zero state is observed at all scaled momentum x_p values above 0.3, with the matrix element rho_00 rising to 0.66 +/- 0.11 for x_p > 0.7. The values of the real part of the off-diagonal element rho_1-1 are negative at large x_p, with a weighted average value of -0.09 +/- 0.03 for x_p > 0.3, in agreement with new theoretical predictions based on Standard Model parameters and coherent fragmentation of the qq(bar) system from the Z^0 decay. All other helicity density matrix elements measured are consistent with zero over the entire x_p range. The K^*0 fragmentation function has also been measured and the total rate determined to be 0.74 +/- 0.02 +/- 0.02 K*(892)^0 mesons per hadronic Z^0 decay.Helicity density matrix elements for inclusive K*(892)^0 mesons from hadronic Z^0 decays have been measured over the full range of K^*0 momentum using data taken with the OPAL experiment at LEP. A preference for occupation of the helicity zero state is observed at all scaled momentum x_p values above 0.3, with the matrix element rho_00 rising to 0.66 +/- 0.11 for x_p > 0.7. The values of the real part of the off-diagonal element rho_1-1 are negative at large x_p, with a weighted average value of -0.09 +/- 0.03 for x_p > 0.3, in agreement with new theoretical predictions based on Standard Model parameters and coherent fragmentation of the qq(bar) system from the Z^0 decay. All other helicity density matrix elements measured are consistent with zero over the entire x_p range. The K^*0 fragmentation function has also been measured and the total rate determined to be 0.74 +/- 0.02 +/- 0.02 K*(892)^0 mesons per hadronic Z^0 decay.Helicity density matrix elements for inclusive K*(892)^0 mesons from hadronic Z^0 decays have been measured over the full range of K^*0 momentum using data taken with the OPAL experiment at LEP. A preference for occupation of the helicity zero state is observed at all scaled momentum x_p values above 0.3, with the matrix element rho_00 rising to 0.66 +/- 0.11 for x_p > 0.7. The values of the real part of the off-diagonal element rho_1-1 are negative at large x_p, with a weighted average value of -0.09 +/- 0.03 for x_p > 0.3, in agreement with new theoretical predictions based on Standard Model parameters and coherent fragmentation of the qq(bar) system from the Z^0 decay. All other helicity density matrix elements measured are consistent with zero over the entire x_p range. The K^*0 fragmentation function has also been measured and the total rate determined to be 0.74 +/- 0.02 +/- 0.02 K*(892)^0 mesons per hadronic Z^0 decay.Helicity density matrix elements for inclusive K ∗ (892) 0 mesons from hadronic Z 0 decays have been measured over the full range of K ∗ 0 momentum using data taken with the OPAL experiment at LEP. A preference for occupation of the helicity zero state is observed at all scaled momentum x p values above 0.3, with the matrix element ϱ 00 rising to 0.66 ± 0.11 for x p > 0.7. The values of the real part of the off-diagonal element ϱ 1 - 1 are negative at large x p , with a weighted average value of −0.09 ± 0.03 for x p > 0.3, in agreement with new theoretical predictions based on Standard Model parameters and coherent fragmentation of the q q system from the Z 0 decay. All other helicity density matrix elements measured are consistent with zero over the entire x p range. The K ∗ 0 fragmentation function has also been measured and the total rate determined to be 0.74 ± 0.02 ± 0.02 K ∗ (892) 0 mesons per hadronic Z 0 decay