Generalized directed loop method for quantum Monte Carlo simulations

Monte Carlo Based Techniques for Quantum Magnets with Long-Range Interactions

Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights into quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give an overview of recent advancements to investigate quantum magnets with long-range interactions focusing on two techniques based on Monte Carlo integration. First, the method of perturbative continuous unitary transformations where classical Monte Carlo integration is applied within the embedding scheme of white graphs. This linked-cluster expansion allows extracting high-order series expansions of energies and observables in the thermodynamic limit. Second, stochastic series expansion quantum Monte Carlo integration enables calculations on large finite systems. Finite-size scaling can then be used to determine the physical properties of the infinite system. In recent years, both techniques have been applied successfully to one- and two-dimensional quantum magnets involving long-range Ising, XY, and Heisenberg interactions on various bipartite and non-bipartite lattices. Here, we summarise the obtained quantum-critical properties including critical exponents for all these systems in a coherent way. Further, we review how long-range interactions are used to study quantum phase transitions above the upper critical dimension and the scaling techniques to extract these quantum critical properties from the numerical calculations.

Numerical Interchain Mean-Field Theory for the Specific Heat of the Bimetallic Ferromagnetically Coupled Chain Compound MnNi(NO2)4(en)2 (en = Ethylenediamine)

We present a detailed study of the field-dependent specific heat of the bimetallic ferromagnetically coupled chain compound MnNi(NO2)4(en)2, en = ethylenediamine. For this material, which in zero field orders antiferromagnetically below TN=2.45 K, small fields suppress magnetic order. Instead, in such fields, a double-peak-like structure in the temperature dependence of the specific heat is observed. We attribute this behavior to the existence of an acoustic and an optical mode in the spin-wave dispersion as a result of the existence of two different spins per unit cell. We compare our experimental data to numerical results for the specific heat obtained by exact diagonalization and Quantum Monte Carlo simulations for the alternating spin-chain model, using parameters that have been derived from the high-temperature behavior of the magnetic susceptibility. The interchain coupling is included in the numerical treatment at the mean-field level. We observe remarkable agreement between experiment and theory, including the ordering transition, using previously determined parameters. Furthermore, the observed strong effect of an applied magnetic field on the ordered state of MnNi(NO2)4(en)2 promises interesting magnetocaloric properties.

Origin of the Unusual Ground-State Spin S = 9 in a Cr10 Single-Molecule Magnet

The molecular wheel [Cr₁₀(OMe)₂₀(O₂CCMe₃)₁₀], abbreviated {Cr₁₀}, with an unusual intermediate total spin S = 9 and non-negligible cluster anisotropy, D/k_B = -0.045(2) K, is a rare case among wheels based on an even number of 3d-metals, which usually present an antiferromagnetic (AF) ground state (S = 0). Herein, we unveil the origin of such a behavior. Angular magnetometry measurements performed on a single crystal confirmed the axial anisotropic behavior of {Cr₁₀}. For powder samples, the temperature dependence of the susceptibility plotted as χT(T) showed an overall ferromagnetic (FM) behavior down to 1.8 K, whereas the magnetization curve M(H) did not saturate at the expected 30 μ_B/fu for 10 FM coupled 3/2 spin Cr³⁺ ions, but to a much lower value, corresponding to S = 9. In addition, the X-ray magnetic circular dichroism (XMCD) measured at high magnetic field (170 kOe) and 7.5 K showed the polarization of the cluster moment up to 23 μ_B/fu. The magnetic results can be rationalized within a model, including the cluster anisotropy, in which the {Cr₁₀} wheel is formed by two semiwheels, each with four Cr³⁺ spins FM coupled (J_FM/k_B = 2.0 K), separated by two Cr³⁺ ions AF coupled asymmetrically (J₂₃/k_B = J₇₈/k_B = -2.0 K; J₃₄/k_B = J₈₉/k_B = -0.25 K). Inelastic neutron scattering and heat capacity allowed us to confirm this model leading to the S = 9 ground state and first excited S = 8. Single-molecule magnet behavior with an activation energy of U/k_B = 4.0(5) K in the absence of applied field was observed through ac susceptibility measurements down to 0.1 K. The intriguing magnetic behavior of {Cr₁₀} arises from the detailed asymmetry in the molecule interactions produced by small-angle distortions in the angles of the Cr-O-Cr alkoxy bridges coupling the Cr³⁺ ions, as demonstrated by ab initio and density functional theory calculations, while the cluster anisotropy can be correlated to the single-ion anisotropies calculated for each Cr³⁺ ion in the wheel.

LT Scaling in Depleted Quantum Spin Ladders

Using a combination of neutron scattering, calorimetry, quantum Monte Carlo simulations, and analytic results we uncover confinement effects in depleted, partially magnetized quantum spin ladders. We show that introducing nonmagnetic impurities into magnetized spin ladders leads to the emergence of a new characteristic length L in the otherwise scale-free Tomonaga-Luttinger liquid (serving as the effective low-energy model). This results in universal LT scaling of staggered susceptibilities. Comparison of simulation results with experimental phase diagrams of prototypical spin ladder compounds bis(2,3-dimethylpyridinium)tetrabromocuprate(II) (DIMPY) and bis(piperidinium)tetrabromocuprate(II) (BPCB) yields excellent agreement.

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C9H18N2CuBr4

Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C₉H₁₈N₂CuBr₄. Single-crystal heat capacity and neutron diffraction measurements reveal that the Néel-ordered phase breaks down beyond a critical pressure of P_c ∼ 1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above P_c.

There is a class of quantum phase transitions that do not fit into the traditional Landau paradigm, but are described in terms of fractionalized degrees of freedom and emergent gauge fields. Hong et al. find evidence of such a transition in a molecular spin-1/2 antiferromagnetic ladder compound under hydrostatic pressure.

Family of Nanoclusters, Ln33 (Ln = Sm/Eu) and Gd32, Exhibiting Magnetocaloric Effects and Fluorescence Sensing for MnO4. Inorg Chem 2022;61:8861-8869. [PMID: 35653200 DOI: 10.1021/acs.inorgchem.2c00898]

A family of nanoclusters, [Ln₃₃(EDTA)₁₂(OAc)₂(CO₃)₄(μ₃-OH)₃₆(μ₅-OH)₄(H₂O)₃₈]·OAc·xH₂O (x ≈ 50, Ln = Sm for 1; x ≈ 70, Ln = Eu for 2) and [Gd₃₂(EDTA)₁₂(OAc)₂(C₂O₄)(CO₃)₂(μ₃-OH)₃₆(μ₅-OH)₄(H₂O)₃₆]·x(H₂O) (x ≈ 70 for 3; H₄EDTA = ethylene diamine tetraacetic acid), was prepared through the assembly of repeating subunits under the action of an anion template. The analysis of the structures showed that compounds 1 and 2 containing 33 Ln³⁺ ions were isostructural, which were constructed by three kinds of subunits in the presence of CO₃^2- as an anion template, while compound 3 had a slightly different structure. Compound 3 containing 32 Gd³⁺ ions was formed by three types of subunits in the presence of CO₃^2- and C₂O₄^2- as a mixed anion template. The CO₃^2- anions came from the slow fixation of CO₂ in the air. Meanwhile, one kind of high-nuclearity lanthanide clusters showed high chemical stability. The quantum Monte Carlo (QMC) calculation suggested that weak antiferromagnetic interactions were dominant between Gd³⁺ ions in 3. Magnetocaloric studies showed that compound 3 had a large entropy change of 43.0 J kg^-1 K^-1 at 2 K and 7 T. Surprisingly, compound 2 showed excellent recognition and detection effects for permanganate in aqueous solvents based on the fluorescence quenching phenomenon.

Quantum Monte Carlo study of topological phases on a spin analogue of Benalcazar-Bernevig-Hughes model

We study the higher-order topological spin phases based on a spin analogue of Benalcazar-Bernevig-Hughes model in two dimensions using large-scale quantum Monte Carlo simulations. A continuous Néel-valence bond solid quantum phase transition is revealed by tuning the ratio between dimerized spin couplings, namely, the weak and strong exchange couplings. Through the finite-size scaling analysis, we identify the phase critical points, and consequently, map out the full phase diagrams in related parameter spaces. Particularly, we find that the valence bond solid phase can be a higher-order topological spin phase, which has a gap for spin excitations in the bulk while demonstrates characteristic gapless spin modes at corners of open lattices. We further discuss the connection between the higher-order topological spin phases and the electronic correlated higher-order phases, and find both of them possess gapless spin corner modes that are protected by higher-order topology. Our result exemplifies higher-order physics in the correlated spin systems and will contribute to further understandings of the many-body higher-order topological phenomena.

Monolayer CrCl_{3} as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX_{3} (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. Here, we show that monolayer CrCl_{3} has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl_{3} for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl_{3}, as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl_{3} to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

Quantum phase transitions of interacting bosons on hyperbolic lattices

The effect of many-body interaction in curved space is studied based on the extended Bose-Hubbard model on hyperbolic lattices. Using the mean-field approximation and quantum Monte Carlo simulation, the phase diagram is explicitly mapped out, which contains the superfluid, supersolid and insulator phases at various fillings. Particularly, it is revealed that the sizes of the Mott lobes shrink and the supersolid is stabilized at smaller nearest-neighbor interaction asqin the Schläfli symbol increases. The underlying physical mechanism is attributed to the increase of the coordination number, and hence the kinetic energy and the nearest-neighbor interaction. The results suggest that the hyperbolic lattices may be a unique platform to study the effect of the coordination number on quantum phase transitions, which may be relevant to the experiments of ultracold atoms in optical lattices.

Amplitude Mode in Quantum Magnets via Dimensional Crossover

We investigate the amplitude (Higgs) mode associated with longitudinal fluctuations of the order parameter at the continuous spontaneous symmetry breaking phase transition. In quantum magnets, due to the fast decay of the amplitude mode into low-energy Goldstone excitations, direct observation of this mode represents a challenging task. By focusing on a quasi-one-dimensional geometry, we circumvent the difficulty and investigate the amplitude mode in a system of weakly coupled spin chains with the help of quantum Monte Carlo simulations, stochastic analytic continuation, and a chain-mean field approach combined with a mapping to the field-theoretic sine-Gordon model. The amplitude mode is observed to emerge in the longitudinal spin susceptibility in the presence of a weak symmetry-breaking staggered field. A conventional measure of the amplitude mode in higher dimensions, the singlet bond mode, is found to appear at a lower than the amplitude mode frequency. We identify these two excitations with the second (first) breather of the sine-Gordon theory, correspondingly. In contrast to higher-dimensional systems, the amplitude and bond order fluctuations are found to carry significant spectral weight in the quasi-1D limit.

Particle-Hole Symmetry Breaking in a Spin-Dimer System TlCuCl_{3} Observed at 100 T

The entire magnetization process of TlCuCl_{3} has been experimentally investigated up to 100 T employing the single-turn technique. The upper critical field H_{c2} is observed to be 86.1 T at 2 K. A convex slope of the M-H curve between the lower and upper critical fields (H_{c1} and H_{c2}) is clearly observed, which indicates that a particle-hole symmetry is broken in TlCuCl_{3}. By quantum Monte Carlo simulation and the bond-operator theory method, we find that the particle-hole symmetry breaking results from strong interdimer interactions.

Higgs Mode of Planar Coupled Spin Ladders and its Observation in C_{9}H_{18}N_{2}CuBr_{4}

Polarized inelastic neutron scattering experiments recently identified the amplitude (Higgs) mode in C_{9}H_{18}N_{2}CuBr_{4}, a two-dimensional near-quantum-critical spin-1/2 two-leg ladder compound, which exhibits a weak easy-axis exchange anisotropy. Here, we theoretically examine the dynamic spin structure factor of such planar coupled spin-ladder systems using large-scale quantum Monte Carlo simulations. This allows us to provide a quantitative account of the experimental neutron scattering data within a consistent quantum spin model. Moreover, we trace the details of the continuous evolution of the amplitude mode from a two-particle bound state of coupled ladders in the classical Ising limit all the way to the quantum spin-1/2 Heisenberg limit with fully restored SU(2) symmetry, where it gets overdamped by the two-magnon continuum in neutron scattering.

Thermal Critical Points and Quantum Critical End Point in the Frustrated Bilayer Heisenberg Antiferromagnet

We consider the finite-temperature phase diagram of the S=1/2 frustrated Heisenberg bilayer. Although this two-dimensional system may show magnetic order only at zero temperature, we demonstrate the presence of a line of finite-temperature critical points related to the line of first-order transitions between the dimer-singlet and -triplet regimes. We show by high-precision quantum Monte Carlo simulations, which are sign-free in the fully frustrated limit, that this critical point is in the Ising universality class. At zero temperature, the continuous transition between the ordered bilayer and the dimer-singlet state terminates on the first-order line, giving a quantum critical end point, and we use tensor-network calculations to follow the first-order discontinuities in its vicinity.

Giant Pressure Dependence and Dimensionality Switching in a Metal-Organic Quantum Antiferromagnet

We report an extraordinary pressure dependence of the magnetic interactions in the metal-organic system [CuF_{2}(H_{2}O)_{2}]_{2}pyrazine. At zero pressure, this material realizes a quasi-two-dimensional spin-1/2 square-lattice Heisenberg antiferromagnet. By high-pressure, high-field susceptibility measurements we show that the dominant exchange parameter is reduced continuously by a factor of 2 on compression. Above 18 kbar, a phase transition occurs, inducing an orbital re-ordering that switches the dimensionality, transforming the quasi-two-dimensional lattice into weakly coupled chains. We explain the microscopic mechanisms for both phenomena by combining detailed x-ray and neutron diffraction studies with quantitative modeling using spin-polarized density functional theory.

Diagnosing Fractionalization from the Spin Dynamics of Z_{2} Spin Liquids on the Kagome Lattice by Quantum Monte Carlo Simulations

Based on large-scale quantum Monte Carlo simulations, we examine the dynamical spin structure factor of the Balents-Fisher-Girvin kagome lattice spin-1/2 model, which is known to harbor an extended Z_{2} quantum spin liquid phase. We use a correlation-matrix sampling scheme combined with a stochastic analytic continuation method to resolve the spectral functions of this anisotropic quantum spin model with a three-site unit cell. Based on this approach, we monitor the spin dynamics throughout the phase diagram of this model, from the XY-ferromagnetic region to the Z_{2} quantum spin liquid regime. In the latter phase, we identify a gapped two-spinon continuum in the transverse scattering channel, which is faithfully modeled by an effective spinon tight-binding model. Within the longitudinal channel, we identify gapped vison excitations and exhibit indications for the translational symmetry fractionalization of the visons via an enhanced spectral periodicity.

Benchmarking the Variational Reduced Density Matrix Theory in the Doubly Occupied Configuration Interaction Space with Integrable Pairing Models

The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called [Formula: see text] N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.

Directed-Loop Quantum Monte Carlo Method for Retarded Interactions

The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out the bosons. This yields an exact algorithm that combines the highly efficient loop updates available in the stochastic series expansion representation with the advantages of avoiding a direct sampling of the bosons. The application to electron-phonon models reveals that the method overcomes the previously detrimental issues of long autocorrelation times and exponentially decreasing acceptance rates. For example, the resulting dramatic speedup allows us to investigate the Peierls quantum phase transition on chains of up to 1282 sites.

Spin-singlet Quantum Ground State in Zigzag Spin Ladder Cu(CF3 COO)2

The copper salt of trifluoroacetic acid, Cu(CF₃ COO)₂ , offers a new platform to investigate the quantum ground states of low-dimensional magnets. In practice, it realizes the ideal case of a solid hosting essentially isolated magnetic monolayers. These entities are constituted by well-separated two-leg half-integer spin ladders organized in a zigzag fashion. The ladders are comprised of dimeric units of edge-sharing tetragonal pyramids coupled through carbon ions. The spin-gap state in this compound was revealed by static and dynamic magnetic measurements. No indications of long range magnetic ordering down to liquid helium temperature were obtained in specific heat measurements. First principles calculations allow estimation of the main exchange interaction parameters, J_⊥ =176 K and J_∥ =12 K, consistent with the weakly interacting dimers model.

Excitation-Gap Scaling near Quantum Critical Three-Dimensional Antiferromagnets

By means of large-scale quantum Monte Carlo simulations, we examine the quantum critical scaling of the magnetic excitation gap (the triplon gap) in a three-dimensional dimerized quantum antiferromagnet, the bicubic lattice, and identify characteristic multiplicative logarithmic scaling corrections atop the leading mean-field behavior. These findings are in accord with field-theoretical predictions that are based on an effective description of the quantum critical system in terms of an asymptotically free field theory, which exhibits a logarithmic decay of the renormalized interaction strength upon approaching the quantum critical point. Furthermore, using bond-based singlet spectroscopy, we identify the amplitude (Higgs) mode resonance within the antiferromagnetic region. We find a Higgs mass scaling in accord with field-theoretical predictions that relate it by a factor of sqrt[2] to the corresponding triplon gap in the quantum disordered regime. In contrast to the situation in lower-dimensional systems, we observe in this three-dimensional coupled-dimer system a distinct signal from the amplitude mode also in the dynamical spin structure factor. Its width is observed to vanish proportional to the Higgs mass in the accessible proximity to the quantum critical point.

Numerical stabilization of entanglement computation in auxiliary-field quantum Monte Carlo simulations of interacting many-fermion systems

In the absence of a fermion sign problem, auxiliary-field (or determinantal) quantum Monte Carlo (DQMC) approaches have long been the numerical method of choice for unbiased, large-scale simulations of interacting many-fermion systems. More recently, the conceptual scope of this approach has been expanded by introducing ingenious schemes to compute entanglement entropies within its framework. On a practical level, these approaches, however, suffer from a variety of numerical instabilities that have largely impeded their applicability. Here we report on a number of algorithmic advances to overcome many of these numerical instabilities and significantly improve the calculation of entanglement measures in the zero-temperature projective DQMC approach, ultimately allowing us to reach similar system sizes as for the computation of conventional observables. We demonstrate the applicability of this improved DQMC approach by providing an entanglement perspective on the quantum phase transition from a magnetically ordered Mott insulator to a band insulator in the bilayer square lattice Hubbard model at half filling.

Sign-Problem-Free Monte Carlo Simulation of Certain Frustrated Quantum Magnets

We introduce a quantum Monte Carlo (QMC) method for efficient sign-problem-free simulations of a broad class of frustrated S=1/2 antiferromagnets using the basis of spin eigenstates of clusters to avoid the severe sign problem faced by other QMC methods. We demonstrate the utility of the method in several cases with competing exchange interactions and flag important limitations as well as possible extensions of the method.

Solution to sign problems in models of interacting fermions and quantum spins

We show that solutions to fermion sign problems that are found in the formulation where the path integral is expanded in powers of the interaction in continuous time can be extended to systems involving fermions interacting with dynamical quantum spins. While these sign problems seem unsolvable in the auxiliary field approach, solutions emerge in the world-line representation of quantum spins. Combining this idea with meron-cluster methods, we are able to further extend the class of models that are solvable. We demonstrate these solutions to sign problems by considering several examples of strongly correlated systems that contain the physics of semimetals, insulators, superfluidity, and antiferromagnetism.

First-principles investigation of exchange interactions in quasi-one-dimensional antiferromagnet CaV2O4

The effects of orbital degrees of freedom on the exchange interactions in a quasi-one-dimensional spin-1 antiferromagnet CaV2O4 are systematically studied. For this purpose a realistic low-energy electron model with the parameters derived from the first-principles calculations is constructed in the Wannier basis for the t2g bands. The exchange interactions are calculated using both the theory of infinitesimal spin rotations near the mean-field ground state and the superexchange model, which provide a consistent description. The obtained behaviour of exchange interactions differs substantially from the previously proposed phenomenological picture based on magnetic measurements and structural considerations, namely: (i) despite the quasi-one-dimensional character of the crystal structure, consisting of the zigzag chains of the edge-sharing VO6 octahedra, the electronic structure is essentially three-dimensional, that leads to finite interactions between the chains; (ii) the exchange interactions along the legs of the chains appear to dominate; and (iii) there is a substantial difference in exchange interactions in two crystallographically inequivalent chains. The combination of these three factors successfully reproduces the behaviour of experimental magnetic susceptibility.

Evidence of a field-induced Berezinskii-Kosterlitz-Thouless scenario in a two-dimensional spin-dimer system

Two-dimensional (2D) systems with continuous symmetry lack conventional long-range order because of thermal fluctuations. Instead, as pointed out by Berezinskii, Kosterlitz and Thouless (BKT), 2D systems may exhibit so-called topological order driven by the binding of vortex-antivortex pairs. Signatures of the BKT mechanism have been observed in thin films, specially designed heterostructures, layered magnets and trapped atomic gases. Here we report on an alternative approach for studying BKT physics by using a chemically constructed multilayer magnet. The novelty of this approach is to use molecular-based pairs of spin S=½ ions, which, by the application of a magnetic field, provide a gas of magnetic excitations. On the basis of measurements of the magnetic susceptibility and specific heat on a so-designed material, combined with density functional theory and quantum Monte Carlo calculations, we conclude that these excitations have a distinct 2D character, consistent with a BKT scenario, implying the emergence of vortices and antivortices.

Temperature dependent magnon-phonon coupling in bcc Fe from theory and experiment

An ab initio based framework for quantitatively assessing the phonon contribution due to magnon-phonon interactions and lattice expansion is developed. The theoretical results for bcc Fe are in very good agreement with high-quality phonon frequency measurements. For some phonon branches, the magnon-phonon interaction is an order of magnitude larger than the phonon shift due to lattice expansion, demonstrating the strong impact of magnetic short-range order even significantly above the Curie temperature. The framework closes the previous simulation gap between the ferro- and paramagnetic limits.

Reliability evaluation of thermophysical properties from first-principles calculations

Thermophysical properties, such as heat capacity, bulk modulus and thermal expansion, are of great importance for many technological applications and are traditionally determined experimentally. With the rapid development of computational methods, however, first-principles computed temperature-dependent data are nowadays accessible. We evaluate various computational realizations of such data in comparison to the experimental scatter. The work is focussed on the impact of different first-principles codes (QUANTUM ESPRESSO and VASP), pseudopotentials (ultrasoft and projector augmented wave) as well as phonon determination methods (linear response and direct force constant method) on these properties. Based on the analysis of data for two pure elements, Cr and Ni, consequences for the reliability of temperature-dependent first-principles data in computational thermodynamics are discussed.

Adiabatic loading of one-dimensional SU(N) alkaline-earth-atom fermions in optical lattices

Ultracold fermionic alkaline earth atoms confined in optical lattices realize Hubbard models with internal SU(N) symmetries, where N can be as large as ten. Such systems are expected to harbor exotic magnetic physics at temperatures below the superexchange energy scale. Employing quantum Monte Carlo simulations to access the low-temperature regime of one-dimensional chains, we show that after adiabatically loading a weakly interacting gas into the strongly interacting regime of an optical lattice, the final temperature decreases with increasing N. Furthermore, we estimate the temperature scale required to probe correlations associated with low-temperature SU(N) magnetism. Our findings are encouraging for the exploration of exotic large-N magnetic states in ongoing experiments.

Solving condensed-matter ground-state problems by semidefinite relaxations

We present a generic approach to the condensed-matter ground-state problem which is complementary to variational techniques and works directly in the thermodynamic limit. Relaxing the ground-state problem, we obtain semidefinite programs (SDP). These can be solved efficiently, yielding strict lower bounds to the ground-state energy and approximations to the few-particle Green's functions. As the method is applicable for all particle statistics, it represents, in particular, a novel route for the study of strongly correlated fermionic and frustrated spin systems in D>1 spatial dimensions. It is demonstrated for the XXZ model and the Hubbard model of spinless fermions. The results are compared against exact solutions, quantum Monte Carlo calculations, and Anderson bounds, showing the competitiveness of the SDP method.

Pair superfluidity of three-body constrained bosons in two dimensions

We examine the equilibrium properties of lattice bosons with attractive on-site interactions in the presence of a three-body hard-core constraint that stabilizes the system against collapse and gives rise to a dimer superfluid phase. Employing quantum Monte Carlo simulations, the ground state phase diagram of this system on the square lattice is analyzed. In particular, we study the quantum phase transition between the atomic and dimer superfluid regime and analyze the nature of the superfluid-insulator transitions. Evidence is provided for the existence of a tricritical point along the saturation transition line, where the transition changes from being first order to a continuous transition of the dilute Bose gas of holes. The Berzinskii-Kosterlitz-Thouless transition from the dimer superfluid to the normal fluid is found to be consistent with an anomalous stiffness jump, as expected from the unbinding of half-vortices.

Distribution of NMR relaxations in a random Heisenberg chain

NMR measurements of the (29)Si spin-lattice relaxation time T(1) were used to probe the spin-1/2 random Heisenberg chain compound BaCu(2)(Si(1-x)Ge(x))(2)O(7). Remarkable differences between the pure (x=0) and the fully random (x=0.5) cases are observed, indicating that randomness generates a distribution of local magnetic relaxations. This distribution, which is reflected in a stretched exponential NMR relaxation, exhibits a progressive broadening with decreasing temperature, caused by a growing inequivalence of magnetic sites. Compelling independent evidence for the influence of randomness is also obtained from magnetization data and Monte Carlo calculations. These results suggest the formation of random-singlet states in this class of materials, as previously predicted by theory.

Measuring entanglement in condensed matter systems

We show how entanglement may be quantified in spin and cold atom many-body systems using standard experimental techniques only. The scheme requires no assumptions on the state in the laboratory, and a lower bound to the entanglement can be read off directly from the scattering cross section of neutrons deflected from solid state samples or the time-of-flight distribution of cold atoms in optical lattices, respectively. This removes a major obstacle which so far has prevented the direct and quantitative experimental study of genuine quantum correlations in many-body systems: The need for a full characterization of the state to quantify the entanglement contained in it. Instead, the scheme presented here relies solely on global measurements that are routinely performed and is versatile enough to accommodate systems and measurements different from the ones we exemplify in this work.

Markov chain Monte Carlo method without detailed balance

We present a specific algorithm that generally satisfies the balance condition without imposing the detailed balance in the Markov chain Monte Carlo. In our algorithm, the average rejection rate is minimized, and even reduced to zero in many relevant cases. The absence of the detailed balance also introduces a net stochastic flow in a configuration space, which further boosts up the convergence. We demonstrate that the autocorrelation time of the Potts model becomes more than 6 times shorter than that by the conventional Metropolis algorithm. Based on the same concept, a bounce-free worm algorithm for generic quantum spin models is formulated as well.

Quantum Monte Carlo simulations of fidelity at magnetic quantum phase transitions

When a system undergoes a quantum phase transition, the ground-state wave function shows a change of nature, which can be monitored using the fidelity concept. We introduce two quantum Monte Carlo schemes that allow the computation of fidelity and its susceptibility for large interacting many-body systems. These methods are illustrated on a two-dimensional Heisenberg model, where fidelity estimators show marked behavior at two successive quantum phase transitions. We also develop a scaling theory which relates the divergence of the fidelity susceptibility to the critical exponent of the correlation length. A good agreement is found with the numerical results.

Solids and supersolids of three-body interacting polar molecules on an optical lattice

We study the physics of cold polar molecules loaded into an optical lattice in the regime of strong three-body interactions, as put forward recently by Büchler et al. [Nature Phys. 3, 726 (2007)]. To this end, quantum Monte Carlo simulations, exact diagonalization, and a semiclassical approach are used to explore hard-core bosons on the 2D square lattice which interact solely by long-ranged three-body terms. The resulting phase diagram shows a sequence of solid and supersolid phases. Our findings are directly relevant for future experimental implementations and open a new route towards the discovery of a lattice supersolid phase in experiment.

Evidence for an unconventional universality class from a two-dimensional dimerized quantum heisenberg model

The two-dimensional J-J' dimerized quantum Heisenberg model is studied on the square lattice by means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio alpha=J'/J. The critical point of the order-disorder quantum phase transition in the J-J' model is determined as alpha_c=2.5196(2) by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimerized models we show, contrary to the current belief, that the critical exponents of the J-J' model are not in agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion of nontrivial critical excitations at the quantum critical point.

Supersolid phase induced by correlated hopping in spin-1/2 frustrated quantum magnets

We show that correlated hopping of triplets, which is often the dominant source of kinetic energy in dimer-based frustrated quantum magnets, produces a remarkably strong tendency to form supersolid phases in a magnetic field. These phases are characterized by simultaneous modulation and ordering of the longitudinal and transverse magnetization, respectively. Using quantum Monte Carlo and a semiclassical approach for an effective hard-core boson model with nearest-neighbor repulsion on a square lattice, we prove, in particular, that a supersolid phase can exist even if the repulsion is not strong enough to stabilize an insulating phase at half-filling. Experimental implications for frustrated quantum antiferromagnets in a magnetic field at zero and finite temperature are discussed.

Short-loop algorithm for quantum Monte Carlo simulations

We present an algorithmic framework for a variant of the quantum Monte Carlo operator-loop algorithm, where nonlocal cluster updates are constructed in a way that makes each individual loop smaller. The algorithm is designed to increase simulation efficiency in cases where conventional loops become very large, do not close altogether, or otherwise behave poorly. We demonstrate and characterize some aspects of the short loop on a square lattice spin-1/2 XXZ model where, remarkably, a significant increase in simulation efficiency is observed in some parameter regimes. The simplicity of the model provides a prototype for the use of short loops on more complicated quantum systems.

Neél order in square and triangular lattice Heisenberg models

We show that the density matrix renormalization group can be used to study magnetic ordering in two-dimensional spin models. Local quantities should be extrapolated with the truncation error, not with its square root. We introduce sequences of clusters, using cylindrical boundary conditions with pinning fields, which provide for rapidly converging finite-size scaling. We determine the magnetization for both the square and triangular Heisenberg lattices with errors comparable to the best alternative approaches.

Universal temperature dependence of the magnetization of gapped spin chains

A Haldane chain under applied field is analyzed numerically, and a clear minimum of magnetization is observed as a function of temperature. We elucidate its origin using the effective theory near the critical field and propose a simple method to estimate the gap from the magnetization at finite temperatures. We also demonstrate that there exists a relation between the temperature dependence of the magnetization and the field dependence of the spin-wave velocity. Our arguments are universal for general axially symmetric one-dimensional spin systems.

Quench dynamics and nonequilibrium phase diagram of the bose-hubbard model

We investigate the time evolution of correlations in the Bose-Hubbard model following a quench from the superfluid to the Mott insulator. For large values of the final interaction strength the system approaches a distinctly nonequilibrium steady state that bears strong memory of the initial conditions. In contrast, when the final interaction strength is comparable to the hopping, the correlations are rather well approximated by those at thermal equilibrium. The existence of two distinct nonequilibrium regimes is surprising given the nonintegrability of the Bose-Hubbard model. We relate this phenomenon to the role of quasiparticle interactions in the Mott insulator.

Bose-glass phases in disordered quantum magnets

In disordered spin systems with antiferromagnetic Heisenberg exchange, transitions into and out of a magnetic-field-induced ordered phase pass through unique regimes. Using quantum Monte Carlo simulations to study the zero-temperature behavior, these intermediate regions are determined to be Bose-glass phases. The localization of field-induced triplons causes a finite compressibility and, hence, glassiness in the disordered phase.

Ultracold atoms in optical lattices with random on-site interactions

We consider the physics of lattice bosons affected by disordered on-site interparticle interactions. Characteristic qualitative changes in the zero-temperature phase diagram are observed when compared to the case of randomness in the chemical potential. The Mott-insulating regions shrink and eventually vanish for any finite disorder strength beyond a sufficiently large filling factor. Furthermore, at low values of the chemical potential both the superfluid and Mott insulator are stable towards formation of a Bose glass leading to a possibly nontrivial tricritical point. We discuss feasible experimental realizations of our scenario in the context of ultracold atoms on optical lattices.

Supersolid hard-core bosons on the triangular lattice

We determine the phase diagram of hard-core bosons on a triangular lattice with nearest-neighbor repulsion, paying special attention to the stability of the supersolid phase. Similar to the same model on a square lattice we find that for densities rho<1/3 or rho>2/3 a supersolid phase is unstable and the transition between a commensurate solid and the superfluid is of first order. At intermediate fillings 1/3<rho<2/3 we find an extended supersolid phase even at half filling rho=1/2. The emergence of the supersolid on the triangular lattice reflects a novel and interesting way for a quantum system to avoid classical frustration, similar to an order-by-disorder mechanism. It also offers an exciting possibility of realizing such phenomena in ultracold atoms on optical lattices.

Stochastic series expansion algorithm for the S = 1/2 XY model with four-site ring exchange

We describe a stochastic series expansion quantum Monte Carlo method for a two-dimensional S = 1/2 XY model (or, equivalently, hard-core bosons at half filling) which in addition to the standard pair interaction J includes a four-site term K that flips spins on a square plaquette. The model has three ordered ground state phases; for K/J approximately < or = 8 it has long-range xy spin order (superfluid bosons), for K/J approximately > or = 15 it has staggered spin order in the z direction (charge-density wave), and between these phases it is in a state with columnar order in the bond and plaquette energy densities. We discuss an implementation of directed-loop updates for the SSE simulations of this model and also introduce a "multibranch" cluster update which significantly reduces the autocorrelation times for large K/J. In addition to the pure J-K model, which in the z basis has only off-diagonal terms, we also discuss modifications of the algorithm needed when various diagonal interactions are included.