共查询到20条相似文献,搜索用时 0 毫秒
1.
Christopher G. Jesudason 《Journal of statistical physics》1996,82(3-4):1207-1211
Simple criteria for convergence of Monte Carlo algorithms not necessarily requiring detailed balance for any specified transition probability are derived and it is shown that it is possible to view the algorithm as a superimposition of a Brownian motion on configurational space coupled to the transition probabilities. As such, the error contributions due to a particular Monte Carlo algorithm and the integration limits in configuration space must be distinguished from those due to the nonuniform sampling of the Brownian motion, and criteria related to the number of steps required to distinguish these errors are provided for the simplest cases involving one dimension and symmetrical probability distributions. 相似文献
2.
In this paper a stochastic series expansion quantum Monte Carlo algorithm is used to study a frustrated spin chain with diagonal next-nearest-neighbor interactions. The detailed balance conditions are carefully analyzed to improve the efficiency of simulation process. As an application of this algothrim, the total magnetization, the static structure factor and spin-stiffness are calculated for a certain set of system parameters as a function of external field strength. 相似文献
3.
YU Hong-Yi 《理论物理通讯》2010,54(2):236-240
In this paper a stochastic series expansion quantum Monte Carlo algorithm is used to study a frustrated spin chain with diagonal next-nearest-neighbor interactions. The detailed balance conditions are carefully analyzed to improve the efficiency of simulation process. As an application of this algothrim, the total magnetization, the static structure factor and spin-stiffness arecalculated for a certain set of system parameters as a function of external field strength. 相似文献
4.
Mengjie Sun 《中国物理 B》2021,30(8):87503-087503
Recently, the family of rare-earth chalcohalides were proposed as candidate compounds to realize the Kitaev spin liquid (KSL) [Chin. Phys. Lett. 38 047502 (2021)]. In the present work, we firstly propose an effective spin Hamiltonian consistent with the symmetry group of the crystal structure. Then we apply classical Monte Carlo simulations to preliminarily study the model and establish a phase diagram. When approaching to the low temperature limit, several magnetic long range orders are observed, including the stripe, the zigzag, the antiferromagnetic (AFM), the ferromagnetic (FM), the incommensurate spiral (IS), the multi-Q, and the 120° ones. We further calculate the thermodynamic properties of the system, such as the temperature dependence of the magnetic susceptibility and the heat capacity. The ordering transition temperatures reflected in the two quantities agree with each other. For most interaction regions, the system is magnetically more susceptible in the ab-plane than in the c-direction. The stripe phase is special, where the susceptibility is fairly isotropic in the whole temperature region. These features provide useful information to understand the magnetic properties of related materials. 相似文献
5.
Frank W. Elliott Jr. Andrew J. Majda David J. Horntrop Richard M. McLaughlin 《Journal of statistical physics》1995,81(3-4):717-736
Two hierarchical Monte Carlo methods for the generation of self-similar fractal random fields are compared and contrasted. The first technique, successive random addition (SRA), is currently popular in the physics community. Despite the intuitive appeal of SRA, rigorous mathematical reasoning reveals that SRA cannot be consistent with any stationary power-law Gaussian random field for any Hurst exponent; furthermore, there is an inherent ratio of largest to smallest putative scaling constant necessarily exceeding a factor of 2 for a wide range of Hurst exponentsH, with 0.30<H<0.85. Thus, SRA is inconsistent with a stationary power-law fractal random field and would not be useful for problems that do not utilize additional spatial averaging of the velocity field. The second hierarchical method for fractal random fields has recently been introduced by two of the authors and relies on a suitable explicit multiwavelet expansion (MWE) with high-moment cancellation. This method is described briefly, including a demonstration that, unlike SRA, MWE is consistent with a stationary power-law random field over many decades of scaling and has low variance. 相似文献
6.
Weilun Jiang 《中国物理 B》2022,31(4):40504-040504
We systematically test the performance of several Monte Carlo update schemes for the (2+1)d XY phase transition of quantum rotor model. By comparing the local Metropolis (LM), LM plus over-relaxation (OR), Wolff-cluster (WC), hybrid Monte Carlo (HM), hybrid Monte Carlo with Fourier acceleration (FA) schemes, it is clear that among the five different update schemes, at the quantum critical point, the WC and FA schemes acquire the smallest autocorrelation time and cost the least amount of CPU hours in achieving the same level of relative error, and FA enjoys a further advantage of easily implementable for more complicated interactions such as the long-range ones. These results bestow one with the necessary knowledge of extending the quantum rotor model, which plays the role of ferromagnetic/antiferromagnetic critical bosons or Z2 topological order, to more realistic and yet challenging models such as Fermi surface Yukawa-coupled to quantum rotor models. 相似文献
7.
We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method,in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. Weconsider two quantum mechanical models: a symmetric one V(x) = |x|/2; and an asymmetric one V(x) = ∞, forx < 0 and V(x) = x, for x ≥ 0. The results for the spectrum, wave functions and thermodynamical observables are inagreement with the analytical or Runge-Kutta calculations. 相似文献
8.
Monte Carlo Hamiltonian:Inverse Potential 总被引:2,自引:0,他引:2
The Monte Carlo Hamiltonian method developed recently
allows to investigate the ground
state and low-lying excited states of a quantum system, using Monte
Carlo (MC) algorithm
with importance sampling. However, conventional MC algorithm
has some difficulties
when applied to inverse potentials. We propose to use effective potential and
extrapolation method to solve the problem. We present examples from the hydrogen
system. 相似文献
9.
10.
Monte Carlo simulations are performed for pure and site-diluted Ising ferro- and ferrimagnets on a simple cubic lattice with up to 403 sites and with impurity concentrationx. For the diluted ferromagnet (x=0.2) the exponent= 0.392±0.03 is definitely larger than the pure model value of=0.304±0.03. In contrast, for ferrimagnetic systems (x=0, 0.1, 0.2) the values appear to be independent ofx and within the error limits consistent with the value for the pure ferromagnet, possibly because the width of the asymptotic random critical regime (or of the crossover regime) is even smaller than in the case of ferromagnets. 相似文献
11.
Don S. LemonsB.J. Albright 《Journal of Quantitative Spectroscopy & Radiative Transfer》2002,74(6):719-729
We model radiation transport by advancing computational photons through phase space with solutions to a set of stochastic differential equations of motion. Random numbers that appear in the equations of motion are sampled with deterministically chosen Gaussian quadrature weights and abscissas. In this way, the advantages of particle Monte Carlo are realized without generating statistical noise. We demonstrate this technique by performing one- and two-dimensional test problems in which gray radiation is energetically coupled to stationary material. Scattering is accomplished with a Fokker-Planck scattering operator. Free streaming, diffusion and Marshak waves are recovered in appropriate limits. 相似文献
12.
A novel path integral Monte Carlo (PIMC) approach for correlated many‐particle systems with arbitrary pair interaction in continuous space at low temperatures is presented. It is based on a representation of the N ‐particle density operator in a basis of (anti‐)symmetrized N ‐particle states (configurations of occupation numbers). The path integral is transformed into a sum over trajectories with the same topology and, finally, the limit of M → ∞, where M is the number of high‐temperature factors, is analytically performed. This yields exact expressions for the thermodynamic quantities and allows to perform efficient simulations for fermions at low temperature and weak to moderate coupling. Our method is expected to be applicable to dense quantum plasmas in the regime of strong degeneracy where conventional PIMC fails due to the fermion sign problem (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim) 相似文献
13.
Michael Betancourt 《Annalen der Physik》2019,531(3)
From its inception in the 1950s to the modern frontiers of applied statistics, Markov chain Monte Carlo has been one of the most ubiquitous and successful methods in statistical computing. The development of the method in that time has been fueled by not only increasingly difficult problems but also novel techniques adopted from physics. Here, the history of Markov chain Monte Carlo is reviewed from its inception with the Metropolis method to the contemporary state‐of‐the‐art in Hamiltonian Monte Carlo, focusing on the evolving interplay between the statistical and physical perspectives of the method. 相似文献
14.
15.
16.
Basil T. Wong 《Journal of Quantitative Spectroscopy & Radiative Transfer》2004,84(4):437-450
Monte Carlo methods (MCMs) are the most versatile approaches in solving the integro-differential equations. They are statistical in nature and can be easily adapted for simulation of the propagation of ensembles of quantum particles within absorbing, emitting, and scattering media. In this paper, we use MCM for the solution of the Boltzmann transport equation, which is the governing equation for both radiative transfer and electron-beam processing. We briefly outline the methodology for the solution of MCMs, and discuss the similarities and differences between the two different application areas. The focus of this paper is primarily on the treatment of different scattering phase functions. 相似文献
17.
We analyse the simulation of strongly degenerate electrons at finite temperature using the recently introduced permutation blocking path integral Monte Carlo (PB‐PIMC) method [T. Dornheim et al., New J. Phys. 17 , 073017 (2015)]. As a representative example, we consider electrons in a harmonic confinement and carry out simulations for up to P = 2000 so‐called imaginary‐time propagators – an important convergence parameter within the PIMC formalism. This allows us to study the P‐dependence of different observables of the configuration space in the Monte Carlo simulations and of the fermion sign problem. We find a surprisingly persisting effect of the permutation blocking for large P, which is explained by comparing different length scales. Finally, we touch upon the uniform electron gas in the warm dense matter regime. 相似文献
18.
19.
Random walk on distant mesh points Monte Carlo methods 总被引:1,自引:0,他引:1
A new technique for obtaining Monte Carlo algorithms based on the Markov chains with a finite number of states is suggested. Instead of the classical random walk on neighboring mesh points, a general way of constructing Monte Carlo algorithms that could be called random walk on distant mesh points is considered. It is applied to solve boundary value problems. The numerical examples indicate that the new methods are less laborious and therefore more efficient.In conclusion, we mention that all Monte Carlo algorithms are parallel and could be easily realized on parallel computers. 相似文献
20.
惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升. 相似文献