tag:blogger.com,1999:blog-5439168179960787195.post4643125417225966146..comments2024-03-28T17:13:01.117+10:00Comments on Condensed concepts: Getting the same answer from complementary numerical methodsRoss H. McKenziehttp://www.blogger.com/profile/09950455939572097456noreply@blogger.comBlogger2125tag:blogger.com,1999:blog-5439168179960787195.post-41552329464368104562012-09-11T02:05:42.003+10:002012-09-11T02:05:42.003+10:00Related to the topic of this thread and since Ross...Related to the topic of this thread and since Ross has mentioned NLC [it is probably better to call it NLCE (numerical linked cluster expansion) in the same spirit of HTE (high temperature expansion)] and FTLM as examples of numerical methods leading to the same results, I wanted to mention a recent study that involves NLCEs and quantum Monte Carlo (QMC) simulations on finite lattices. They were both used to explore thermodynamic properties and short-range correlations, of the single-band Fermi-Hubbard model on the square and honeycomb lattices. This work (arXiv:1206.0006), is a collaboration between myself, Marcos Rigol, his PhD student, Baoming Tang, and Thereza Paiva, and is part of a series of studies in which NLCEs have been used to explore finite temperature, and even ground state properties, of systems ranging from frustrated magnets to itinerant particles in optical lattices (see e.g., arXiv:1009.4700, arXiv:1104.5494, arXiv:1105.4147).<br /><br />We find that NLCE and QMC agree with each other very well for the two lattice geometries in the parameter regions that are accessible to both of them, despite the systematic errors that exist in QMC (e.g., finite size, or Trotter errors). Perhaps more interestingly, we find that the NLCE and QMC can be complementary numerical methods when it comes to Hubbard models; QMC runs into sampling problems as the ratio of the interaction strength over hopping amplitude, U/t, gets larger, while NLCE converges to lower temperatures in that regime. This can be understood as t^2/U becomes the relevant energy scale for large values of U/t, and NLCEs are very efficient up to where correlations (in this case, spin correlations) exceed the largest clusters considered. The opposite happens in the weak-coupling regime, where QMC is more efficient and NLCE's convergence temperature is typically higher. As a result of this, we have been able to, for example, map out the isentropic paths of the model at half filling in the temperature-interaction space for a wide range of interactions, up to three time the bandwidth (see Fig. 2 in arXiv:1206.0006). We also see that this agreement is not limited to half filling when we use a local density approximation to express different quantities of a systems that is subject to an external harmonic trapping potential, emulating the optical lattice experiments setup.<br /><br />Another feature of NLCE, and also an advantage over QMC-based methods, that I think is worth mentioning is that because of its exact diagonalization core, and therefore, full access to the partition function, properties such as the entropy or the specific heat can be calculated directly and without any need for numerical differentiation or integration. This has allowed us to obtain very accurate results for the specific heat of the strongly-coupled Fermi-Hubbard model in a wide temperature region (see Fig. 1 in arXiv:1204.1556).Ehsan Khatamihttps://www.blogger.com/profile/15025542501021605282noreply@blogger.comtag:blogger.com,1999:blog-5439168179960787195.post-51124055099113299982012-09-05T17:14:12.456+10:002012-09-05T17:14:12.456+10:00Rajiv Singh asked me to post his comment:
Thank y...Rajiv Singh asked me to post his comment:<br /><br />Thank you for the post. Perhaps one of Marcos' postdoc Ehsan can say something about how the method compares with QMC (even when there are no minus signs). It can be competitive or better than QMC if U is large enough.<br /><br />Another point, I would like to make is that compared with Exact Diagonalization (ED), NLC is particularly useful in 3D. Graphs are not that much harder to count as you go up in dimension, where as clusters with periodic boundary conditions needed<br />in ED will have very small linear dimension associated with them and even if doable will bring a lot of artifact through the periodic boundary condition.<br /><br />One example is the problem of Quantum Spin Ice on a pyrochlore lattice. The smallest periodic cluster for ED that has all the symmetries of the pyrochlore lattice has 16 sites, which is difficult (though not impossible) to diagonalize.<br />Most of the spin-space symmetries are absent in strong spin-orbit coupled systems. But, PBC will allow momentum as a good quantum number (for a 16 site cluster it only has 4 different values so will roughly reduce the Hilbert space by about a factor of 4). However, on this 16-site cluster, you only go about two neighbors in any direction before coming back to where you started. This is clearly very artificial.<br /><br />On the other hand, we recently carried out an NLC calculation for this problem in terms of clusters with complete tetrahedra to 4th order. The largest cluster had 13 sites. It gave reliable answers to thermodynamic properties down to a temperature of about J/2. The material Yb2Ti2O7 has largest exchange constant of about 2K, so the convergence was good down to 1K. In strong fields, NLC converges at all temperatures (down to T=0). That is because, the remainder (higher order terms) is high order in J/T and also in J/h. This covers a large range of parameters of interest to the experimental studies.<br /><br />A reference to our paper is:<br />Vindication of Yb2Ti2O7 as a Model Exchange Quantum Spin Ice<br />R. Applegate, N. R. Hayre, R. R. P. Singh, T. Lin, A. G. R. Day, and M. J. P. Gingras<br />PRL 28 August 2012 (5 pages)<br />097205Ross H. McKenziehttps://www.blogger.com/profile/09950455939572097456noreply@blogger.com