We report progress toward computing the heavy quark momentum diffusion coefficient from the correlator of two chromoelectric fields attached to a Polyakov loop in pure SU(3) gauge theory. Using a multilevel algorithm ...
详细信息
We report progress toward computing the heavy quark momentum diffusion coefficient from the correlator of two chromoelectric fields attached to a Polyakov loop in pure SU(3) gauge theory. Using a multilevel algorithm and tree-level improvement, we study the behavior of the diffusion coefficient as a function of temperature in the wide range 1.1
We study properties of Markov chain Monte Carlo simulations of classical spin models with local updates. We derive analytic expressions for the mean value of the acceptance rate of single-spin-flip algorithms for the ...
详细信息
We study properties of Markov chain Monte Carlo simulations of classical spin models with local updates. We derive analytic expressions for the mean value of the acceptance rate of single-spin-flip algorithms for the one-dimensional Ising model. We find that for the Metropolis algorithm the average acceptance rate is a linear function of energy. We further provide numerical results for the energy dependence of the average acceptance rate for the three- and four-state Potts model, and the XY model in one and two spatial dimensions. In all cases, the acceptance rate is an almost linear function of the energy in the critical region. The variance of the acceptance rate is studied as a function of the specific heat. While the specific heat develops a singularity in the vicinity of a phase transition, the variance of the acceptance rate stays finite.
The Mermin-Wagner theorem prevents the stabilization of long-range magnetic order in two-dimensional layered materials, such as the pnictide superconductors, unless the magnetism is associated with a discrete symmetry...
详细信息
The Mermin-Wagner theorem prevents the stabilization of long-range magnetic order in two-dimensional layered materials, such as the pnictide superconductors, unless the magnetism is associated with a discrete symmetry breaking. A typical known example is the discrete row and column collinear magnetic state that emerges in doped iron pnictide materials due to an order-by-disorder mechanism. In these compounds, the magnetic state competes with superconductivity and the mechanism that stabilizes magnetism remains controversial. In this work, we report the phase diagram of a doped frustrated Heisenberg model obtained through Monte Carlo simulations combined with a parallel tempering simulation technique. The emergence of long-range magnetic order is stabilized by interactions between the magnetic dopant impurities.
Square ice is a statistical mechanics model for two-dimensional ice, widely believed to have a conformally invariant scaling limit. We associate a Peano (space-filling) curve to a square ice configuration, and more ge...
详细信息
Square ice is a statistical mechanics model for two-dimensional ice, widely believed to have a conformally invariant scaling limit. We associate a Peano (space-filling) curve to a square ice configuration, and more generally to a so-called six-vertex model configuration, and argue that its scaling limit is a space-filling version of the random fractal curve SLEκ, Schramm-Loewner evolution with parameter κ, where 4<κ≤12+82. For square ice, κ=12. At the “free-fermion point” of the six-vertex model, κ=8+43. These unusual values lie outside the classical interval 2≤κ≤8.
We have examined the spin ordering of an axial next-nearest-neighbor Ising model in two dimensions (2D) near above the antiphase (〈2〉 phase). We considered an NR-replica system and calculated an overlap function qm be...
详细信息
We have examined the spin ordering of an axial next-nearest-neighbor Ising model in two dimensions (2D) near above the antiphase (〈2〉 phase). We considered an NR-replica system and calculated an overlap function qm between different replicas, having used a cluster heat bath Monte Carlo method. We determined transition temperature between the 〈2〉 phase and a floating incommensurate (IC) phase as TC2/J=0.89±0.01 with frustration ratio κ(≡−J2/J1)=0.6. We found that the spin state at T≳TC2 may be called a domain state, because the spin structure is characterized by a sequentially arranged four types of domains with different 〈2〉 structures. In the domain state, the 2D XY symmetry of the spin correlation in the IC phase weakly breaks, and the diversity of the spin arrangement increases as T→TC2. The Binder ratio gL exhibits a depression at T∼TC2 and the quasiperiodic spin structure, which is realized in the IC phase, becomes diverse at T≳TC2. We discussed that the domain state is stable against the thermal fluctuation which brings a two-stage development of the spin structure at low temperatures.
We present a solution to the problem of interpreting neural networks classifying phases of matter. We devise a procedure for reconstructing the decision function of an artificial neural network as a simple function of...
详细信息
We present a solution to the problem of interpreting neural networks classifying phases of matter. We devise a procedure for reconstructing the decision function of an artificial neural network as a simple function of the input, provided the decision function is sufficiently symmetric. In this case one can easily deduce the quantity by which the neural network classifies the input. The method is applied to the Ising model and SU(2) lattice gauge theory. In both systems we deduce the explicit expressions of the order parameters from the decision functions of the neural networks. We assume no prior knowledge about the Hamiltonian or the order parameters except Monte Carlo–sampled configurations.
The Heisenberg-Kitaev model is a paradigmatic model to describe the magnetism in honeycomb-lattice Mott insulators with strong spin-orbit coupling, such as A2IrO3 (A=Na, Li) and α−RuCl3. Here, we study in detail the ...
详细信息
The Heisenberg-Kitaev model is a paradigmatic model to describe the magnetism in honeycomb-lattice Mott insulators with strong spin-orbit coupling, such as A2IrO3 (A=Na, Li) and α−RuCl3. Here, we study in detail the physics of the Heisenberg-Kitaev model in an external magnetic field. Using a combination of Monte Carlo simulations and spin-wave theory, we map out the classical phase diagram for different directions of the magnetic field. Broken SU(2) spin symmetry renders the magnetization process rather complex, with sequences of phases and metamagnetic transitions. In particular, we find various large-unit-cell and multi-Q phases including a vortex-crystal phase for a field in the [111] direction. We also discuss quantum corrections in the high-field phase.
We have studied zero-temperature metastable minima in classical m-vector component spin glasses in the presence of m-component random fields for two models, the Sherrington-Kirkpatrick (SK) model and the Viana-Bray (V...
详细信息
We have studied zero-temperature metastable minima in classical m-vector component spin glasses in the presence of m-component random fields for two models, the Sherrington-Kirkpatrick (SK) model and the Viana-Bray (VB) model. For the SK model we have calculated analytically its complexity (the log of the number of minima) for both the annealed case where one averages the number of minima before taking the log and the quenched case where one averages the complexity itself, both for fields above and below the de Almeida-Thouless (AT) field, which is finite for m>2. We have done numerical quenches starting from a random initial state (infinite temperature state) by putting spins parallel to their local fields until there is no further decrease of the energy and found that in zero field it always produces minima that have zero overlap with each other. For the m=2 and m=3 cases in the SK model the final energy reached in the quench is very close to the energy Ec at which the overlap of the states would acquire replica symmetry-breaking features. These minima have marginal stability and will have long-range correlations between them. In the SK limit we have analytically studied the density of states ρ(λ) of the Hessian matrix in the annealed approximation. Despite the fact that in the presence of a random field there are no continuous symmetries, the spectrum extends down to zero with the usual λ form for the density of states for fields below the AT field. However, when the random field is larger than the AT field, there is a gap in the spectrum, which closes up as the AT field is approached. The VB model behaves differently and seems rather similar to studies of the three-dimensional Heisenberg spin glass in a random vector field.
A new technique is explored for the Monte Carlo sampling of complex-valued distributions. The method is based on a heat bath approach where the conditional probability is replaced by a positive representation of it on...
详细信息
A new technique is explored for the Monte Carlo sampling of complex-valued distributions. The method is based on a heat bath approach where the conditional probability is replaced by a positive representation of it on the complex plane. Efficient ways to construct such representations are also introduced. The performance of the algorithm is tested on small and large lattices with a λϕ4 theory with quadratic nearest-neighbor complex coupling. The method works for moderate complex couplings, reproducing reweighting and complex Langevin results and fulfilling various Schwinger-Dyson relations.
暂无评论