The Zero-Removing Property and Lagrange-Type Interpolation Series

The classical Kramer sampling theorem, which provides a method for obtaining orthogonal sampling formulas, can be formulated in a more general nonorthogonal setting. In this setting, a challenging problem is to characterize the situations when the obtained nonorthogonal sampling formulas can be expressed as Lagrange-type interpolation series. In this article a necessary and sufficient condition is given in terms of the zero removing property. Roughly speaking, this property concerns the stability of the sampled functions on removing a finite number of their zeros

Feedback-optimized parallel tempering Monte Carlo

We introduce an algorithm to systematically improve the efficiency of parallel tempering Monte Carlo simulations by optimizing the simulated temperature set. Our approach is closely related to a recently introduced adaptive algorithm that optimizes the simulated statistical ensemble in generalized broad-histogram Monte Carlo simulations. Conventionally, a temperature set is chosen in such a way that the acceptance rates for replica swaps between adjacent temperatures are independent of the temperature and large enough to ensure frequent swaps. In this paper, we show that by choosing the temperatures with a modified version of the optimized ensemble feedback method we can minimize the round-trip times between the lowest and highest temperatures which effectively increases the efficiency of the parallel tempering algorithm. In particular, the density of temperatures in the optimized temperature set increases at the "bottlenecks'' of the simulation, such as phase transitions. In turn, the acceptance rates are now temperature dependent in the optimized temperature ensemble. We illustrate the feedback-optimized parallel tempering algorithm by studying the two-dimensional Ising ferromagnet and the two-dimensional fully-frustrated Ising model, and briefly discuss possible feedback schemes for systems that require configurational averages, such as spin glasses.Comment: 12 pages, 14 figure

Lack of self-averaging of the specific heat in the three-dimensional random-field Ising model

We apply the recently developed critical minimum energy subspace scheme for the investigation of the random-field Ising model. We point out that this method is well suited for the study of this model. The density of states is obtained via the Wang-Landau and broad histogram methods in a unified implementation by employing the N-fold version of the Wang-Landau scheme. The random-fields are obtained from a bimodal distribution ($h_{i}=\pm2$), and the scaling of the specific heat maxima is studied on cubic lattices with sizes ranging from $L=4$ to $L=32$. Observing the finite-size scaling behavior of the maxima of the specific heats we examine the question of saturation of the specific heat. The lack of self-averaging of this quantity is fully illustrated and it is shown that this property may be related to the question mentioned above.Comment: 8 pages, 7 figures, extended version with two new figures, version as accepted for publication to Physical Review

Dynamics of the Wang-Landau algorithm and complexity of rare events for the three-dimensional bimodal Ising spin glass

We investigate the performance of flat-histogram methods based on a multicanonical ensemble and the Wang-Landau algorithm for the three-dimensional +/- J spin glass by measuring round-trip times in the energy range between the zero-temperature ground state and the state of highest energy. Strong sample-to-sample variations are found for fixed system size and the distribution of round-trip times follows a fat-tailed Frechet extremal value distribution. Rare events in the fat tails of these distributions corresponding to extremely slowly equilibrating spin glass realizations dominate the calculations of statistical averages. While the typical round-trip time scales exponential as expected for this NP-hard problem, we find that the average round-trip time is no longer well-defined for systems with N >= 8^3 spins. We relate the round-trip times for multicanonical sampling to intrinsic properties of the energy landscape and compare with the numerical effort needed by the genetic Cluster-Exact Approximation to calculate the exact ground state energies. For systems with N >= 8^3 spins the simulation of these rare events becomes increasingly hard. For N >= 14^3 there are samples where the Wang-Landau algorithm fails to find the true ground state within reasonable simulation times. We expect similar behavior for other algorithms based on multicanonical sampling.Comment: 9 pages, 12 figure

Location of the Multicritical Point for the Ising Spin Glass on the Triangular and Hexagonal Lattices

A conjecture is given for the exact location of the multicritical point in the phase diagram of the +/- J Ising model on the triangular lattice. The result p_c=0.8358058 agrees well with a recent numerical estimate. From this value, it is possible to derive a comparable conjecture for the exact location of the multicritical point for the hexagonal lattice, p_c=0.9327041, again in excellent agreement with a numerical study. The method is a variant of duality transformation to relate the triangular lattice directly with its dual triangular lattice without recourse to the hexagonal lattice, in conjunction with the replica method.Comment: 9 pages, 1 figure; Minor corrections in notatio

Local excitations in mean field spin glasses

We address the question of geometrical as well as energetic properties of local excitations in mean field Ising spin glasses. We study analytically the Random Energy Model and numerically a dilute mean field model, first on tree-like graphs, equivalent to a replica symmetric computation, and then directly on finite connectivity random lattices. In the first model, characterized by a discontinuous replica symmetry breaking, we found that the energy of finite volume excitation is infinite whereas in the dilute mean field model, described by a continuous replica symmetry breaking, it slowly decreases with sizes and saturates at a finite value, in contrast with what would be naively expected. The geometrical properties of these excitations are similar to those of lattice animals or branched polymers. We discuss the meaning of these results in terms of replica symmetry breaking and also possible relevance in finite dimensional systems.Comment: 7 pages, 4 figures, accepted for publicatio

"Glassy Dynamics" in Ising Spin Glasses -- Experiment and Simulation

The field-cooled magnetization (FCM) processes of Ising spin glasses under relatively small fields are investigated by experiment on Fe_{0.55}Mn_{0.45}TiO_3 and by numerical simulation on the three-dimensional Edwards-Anderson model. Both results are explained in a unified manner by means of the droplet picture. In particular, the cusp-like behavior of the FCM is interpreted as evidence, not for an equilibrium phase transition under a finite magnetic field, but for a dynamical (`blocking') transition frequently observed in glassy systems.Comment: 4 pages, 7 figure

The relationship between [OIII]5007A equivalent width and obscuration in AGN

In this paper we study the relationship between the equivalent width (EW) of the [OIII]5007A narrow emission line in AGN and the level of obscuration. To this end, we combine the results of a systematic spectral analysis, both in the optical and in the X-rays, on a statistically complete sample of ~170 X-ray selected AGN from the XMM-Newton Bright Serendipitous Source sample (XBS). We find that the observed large range of [OIII]5007A equivalent widths observed in the sample (from a few A up to 500A) is well explained as a combination of an intrinsic spread, probably due to the large range of covering factors of the Narrow Line Region, and the effect of absorption. The intrinsic spread is dominant for EW below 40-50A while absorption brings the values of EW up to ~100-150A, for moderate levels of absorption (AV~0.5-2 mag) or up to ~500A for AV>2 mag. In this picture, the absorption has a significant impact on the observed EW also in type~1 AGN. Using numerical simulations we find that this model is able to reproduce the [OIII]5007A EW distribution observed in the XBS sample and correctly predicts the shape of the EW distribution observed in the optically selected sample of QSO taken from the SDSS survey.Comment: 7 pages, 5 figures. Accepted for publication in MNRA