Accurate computations of Rashba spin-orbit coupling in interacting systems:
from the Fermi gas to real materials
We describe the treatment of Rashba spin-orbit coupling (SOC) in interacting many-fermion systems within the auxiliary-field quantum Monte Carlo framework, and present a set of illustrative results. These include numerically exact calculations on the ground-state properties of the spin-balanced, attractive two-dimensional Fermi gas, as well as a study of a tight-binding Hamiltonian with repulsive interaction. These systems are formally connected via the Hubbard Hamiltonian with SOC, but cover different physics ranging from superfluidity and triplet pairing to SOC in real materials in the presence of strong interactions in localized orbitals. We carry out detailed benchmark studies of the method in the latter case when an approximation is needed to control the sign problem for repulsive Coulomb interactions. The methods presented here provide an approach for predictive computations in materials to study the interplay of SOC and strong correlation.
Spin-orbit coupling lies at the heart of a tremendous variety of physical phenomena, in contexts ranging from semiconductors and metals to 2D heterostructures and even ultracold atoms. The field has garnered renewed interest recently due to the central role that SOC plays in many exotic topological phases, including the quantum spin Hall effect, topological insulators and superconductors, and Majorana fermions reviewRashbaMaterials (); reviewRashba_nat (); reviewRashba_rep_prog_phys (). These developments underpin modern progress in spintronics Koo1515 (); reviewSpintronics (), have opened the field of spin-orbitronics, and carry important implications for quantum computation and information reviewQuantumSpintronics ().
Recent progress in the understanding and manipulation of SOC in real materials has been greatly complemented by rapid and remarkable progress in the field of cold atom physics. Two notable achievements are the ability to tune interaction strengths in ultracold gases using Feshbach resonances, and the ability to load ultracold gases into optical lattice potentials reviewBloch (); reviewLewenstein (). These techniques have ushered in a new era of quantum engineering and simulation, and provide an ideal platform to study many-body physics 2DFG_expt (); 2DFG_AFQMC (); 2DFG_DMC (); 2DFG_FINITET_PRL (). A wide variety of lattice models have already been experimentally realized, with even more likely in the near future. These clean systems, with a high degree of experimental control, offer a novel means of exploring exotic states of matter. Many of the difficulties encountered in real materials, including fixed structural and electromagnetic properties, and disorder, are surmounted in ultracold atom experiments, which can be freely tuned to probe broad parameter regimes. The third major step was the development of methods to generate “synthetic SOC” in ultracold atoms artificialgaugefieldJaksch (); PhysRevLett.95.010403 (); PhysRevLett.108.225303 (); PhysRevLett.107.255301 (); PhysRevLett.111.185302 (); GerbierArtifcialGaugeFields (), which has created new paths towards the realization of many of the exotic states long sought after in real materials, such as high- superconductivity and topological phases, including quantum spin Hall states and Majorana fermions. These techniques have enabled experimental emulations of a wide variety of fundamental models, including several with Rashba, Dresselhaus, or equal Rashba-Dresselhaus SOC SOC_BEC_2009 (); SOC_BEC_2011 (); SOC_DFG_CHEUK (); SOC_DFG_WANG (); RASHBASOC2D_EXPT ().
This expanding experimental horizon has generated tremendous opportunities for theoretical and numerical approaches to many-body systems. The delicate nature of many of the exotic states demands high-accuracy theoretical and numerical input. Careful characterization of these phases, and their behavior in the presence of strong interactions, will provide crucial benchmarks for experiments as well as essential guidelines for the design and fabrication of novel devices reviewSpinHallDevices (). These systems are typically quite challenging from a theoretical or numerical perspective. Ab initio treatment of strongly-correlated real materials remains an outstanding challenge in many-body physics. However, ultracold atom experiments, which can be measured to high accuracy, offer a unique opportunity to test and calibrate our theoretical and computational methods.
The auxiliary-field quantum Monte Carlo (AFQMC) method is a general computational approach for interacting many-fermion systems. AFQMC has been applied to a variety of electronic models and materials without SOC, including Hubbard-like models MingpuHubbardBenchmark (); Hao-symmetry (), molecular systems AFQMC_backprop_molecules (); HydrogenBenchmark (); Mo2_AFQMC (), and solids AFQMC_downfolding (). In ultracold atoms, a number of recent successes have highlighted the accuracy and capability of AFQMC, including studies of the attractive spin-balanced Fermi gas 2DFG_AFQMC (); Ettore-gaps () and fully treating interaction and Rashba SOC on the same footing 2DFG_SOC_AFQMC (); 2DFG_RASHBASOC_OPLATT_AFQMC (). The Hamiltonian for the spin-balanced system with contact interaction preserves time-reversal symmetry, thus guaranteeing that simulations are free of the sign problem nosignproblem (), so the results obtained are numerically exact. This is a particularly notable achievement given the emergence of quantum gas microscopes with single-site resolution PhysRevLett.116.235301 (); Greif953 (); Cheuk2015 (); Haller2015 (); Parsons2015 (), which provide a precise calibration of the technique via direct comparison with a clean, tunable system.
A natural next step is to treat electronic Hamiltonians with SOC, which is the topic of focus in the present paper. We present an illustrative set of results on the attractive 2D Fermi gas with Rashba SOC, as well as the first study of Rashba SOC in a real material context using AFQMC. From the ultracold atom Hamiltonians we have treated, an immediate connection to real materials can be made by considering systems with repulsive interactions. We demonstrate the capabilities of the AFQMC method by considering a repulsive Hubbard model, which captures much of the underlying physics of real materials. In this way, the two types of systems we study here, the real material and ultracold atoms, are described by a similar Hamiltonian, namely the Hubbard model, with opposite sign interactions. This Hamiltonian is of fundamental importance in condensed matter and many-body physics. It can be engineered and emulated by optical lattice experiments 2DHubbard_Hulet (), with the sign of the interaction controlled using Feshbach resonances. This would permit a direct comparison of experimental results with our calculations, and consequently a better understanding of the physics of both systems, as well as an additional calibration of the AFQMC method.
With the repulsive interaction in the Hubbard model, the sign problem reemerges in AFQMC, but as we will show, this can be systematically controlled by applying a constraint CPMC (). The accuracy achieved is comparable to what has been systematically seen in real materials without SOC. Since the AFQMC method provides access to a wide range of observables to probe the charge, spin, pairing, and transport properties, we expect this development to lead to significant applications in correlated materials with strong SOC.
We organize the remainder of the paper as follows. In section 2 we introduce the general Hamiltonian used to describe systems with Rashba SOC, and establish its connection to both real materials and the Fermi gas. Section 3 outlines the AFQMC method, highlighting several of the recent advances incorporated to treat SOC. Two different sampling approaches are discussed, for cases without and with the sign problem, respectively. Section 4 provides a demonstration of the technique, first in the context of cold atoms and then on a repulsive Hubbard model, which is closely connected to the behavior of real materials. Finally, section 5 summarizes and offers an outlook on open questions and promising directions in this burgeoning field.
We begin with a form of the SOC Hamiltonian on a lattice,
where creates a spin- ( or ) particle with momentum , and denotes the density operator in real-space on site with spin-. Although we have used the Rashba form of SOC, our discussions are general for other types. For two-dimensional lattice systems the dispersion and SOC terms are and . The parameter , set to unity throughout this work, determines the strength of nearest-neighbor hopping, the parameter controls the strength of SOC, and the parameter determines the strength of the on-site interaction. For the Fermi gas, which represents the low-density limit of the lattice model, and . In both cases, a natural description of the system in terms of helicity bands, , is obtained by diagonalizing the non-interacting Hamiltonian.
The Fermi gas Hamiltonian is defined in the dilute limit of Eq. (1), and the interaction is negative: . In order to directly compare with experiments, the parameters and are fully specified, as we further discuss in Sec. 4.1. The physics is dictated by the quantity , where the Fermi wave-vector measures the inverse of the average inter-particle spacing while is the scattering length. As decreases in 2D, the system undergoes a BCS-BEC crossover. Our study will examine the interplay between this effect and SOC, which induces triplet pairing.
For the second application, namely as a model for electronic systems, we will consider , and consider a range of values for from to about , which represents strong local interactions as is typical in models for cuprates. Most of our calculations will be at intermediate interaction , and we examine different regimes of SOC strength by varying . Our study here is mostly for testing the algorithm. In most cases we focus on small lattice sizes for which we can obtain exact results from exact diagonalization. With this description of the model for both cases, we can now outline its treatment within the AFQMC framework.
We will first provide an overview of the AFQMC method Lecture-notes (); BSS (); Koonin (), and then discuss some of the extensions necessary to treat systems with SOC. In general, ground-state QMC methods rely on imaginary-time projection to obtain the many-body ground-state of a given hamiltonian from a starting trial wavefunction . This projection proceeds according to,
provided (i.e. the trial wavefunction cannot be orthogonal to the many-body ground-state).
In order to carry out this projection numerically, we first discretize the imaginary time interval into time slices,
so that the limit in (2) can be obtained iteratively via,
To proceed, we must rewrite the many-body propagator in a single-particle form. This is accomplished by applying the Trotter-Suzuki breakup,
where contains the one-body terms and the two-body terms of the Hamiltonian in (1). This step is followed by a suitable Hubbard-Stratonovich (HS) transformation continuous-HS-transformation (). We list below four varieties of discrete HS transformation commonly used in AFQMC simulations: The charge decomposition HS_transform_discrete () is written,
with determined according to ; The spin decomposition in the - or -direction has the form,
and the spin decomposition in the direction is,
This procedure yields the following form for the propagator,
where is a set of auxiliary fields at a given time slice, with dimension equal to the size of the single-particle basis, which in the case of lattice systems is typically the number of lattice sites. Using the charge decomposition as an example, the probability density function is uniform, and,
with . The many-body propagator is now composed of single particle operators with the fermions in external auxiliary fields. The integration over auxiliary field configurations recovers the two-body interactions.
3.2 Metropolis sampling of paths in AF space
Ground-state observables are calculated according to,
The denominator in (11) is,
and we have introduced the notation,
In the above, represents an auxiliary field configuration at time slice , and the collection of auxiliary fields comprises a path in auxiliary field space.
With this shorthand, Eq. (11) can be written as a path integral over auxiliary fields,
This integral can be evaluated using standard Monte Carlo techniques, such as the Metropolis algorithm, which samples auxiliary-fields from to obtain a Monte Carlo estimate of the expectation value in (14). To accelerate the sampling procedure we employ a dynamic force bias 2DFG_AFQMC (); Lecture-notes (), which improves the acceptance ratio and consequently the efficiency of the algorithm. In addition, we remove the infinite variance problem using the bridge link method inf_var ().
3.3 Branching random walks and the constraint
In systems that preserve time-reversal symmetry, such as the spin-balanced attractive Hubbard model (even in the presence of SOC, as discussed below) or the half-filled repulsive Hubbard model, the denominator remains positive and the path-integrals in Eq. (11) can be calculated with the Monte Carlo technique outlined above, which samples the probability density function via the Metropolis algorithm. In cases where time-reversal symmetry is broken, is no longer guaranteed to be positive, and the Monte Carlo signal obtained by this straightforward sampling procedure is lost to sampling noise. This is a manifestation of the well-known sign problem Lecture-notes ().
In order to treat systems with a sign problem we recast the procedure outlined above as an open-ended random walk in Slater-determinant space, which then allows the imposition of a constraint to prevent the decay of the Monte Carlo signal CPMC (). Returning to (4), we have,
In this formulation the wave function at each step is represented by an ensemble of Slater determinants,
where denotes the weight for the -th walker at time step , and the sum runs over the entire walker population at that time step.
To implement this constraint within the random walk procedure we define an importance function,
which prevents walkers from acquiring a negative overlap with .
With the addition of importance sampling we carry out a modified version of the random walk,
with modified probability density function,
and modified wavefunction,
The true wavefunction is then,
The modified projection equation written in (20) is identical to the version in (16). In the modified version auxiliary fields are sampled from , which favors determinants with larger overlaps with , and vanishes if the overlap is zero, thus enforcing the constraint in (18). After the random walk has equilibrated, the ground state wavefunction is represented by the distribution of weighted random walkers.
It is important to note that the branching random walk approach discussed here is closely related to the Metropolis sampling discussed in Sec. 3.2. (The dynamic force bias we use in the Metropolis sampling has the same origin and the same form as that obtained from the importance sampling of in Eq. (21).) The only difference is that the Metropolis procedure, while more convenient for computing observables, encounters severe ergodicity problems when the constraint has to be imposed. In the current approach, the total energy is calculated in a manner similar to that in Sec. 3.2, except is fixed to be . For computing observables that do not commute with the Hamiltonian, we use the backpropagation technique CPMC (); AFQMC_bosons (); AFQMC_backprop_molecules ().
3.4 Generalizations to SOC
Having outlined the standard AFQMC procedures, we can now address the modifications necessary to treat SOC. In systems without SOC, there is no mixing of different spin sectors, so random walkers can be factored into - and -spin components,
For a system of particles with single particle orbitals of dimension , the walker is just a Slater determinant, , with matrix form,
The one body propagator used to project these random walkers can also be split into - and -spin components,
where contains only spin- operators, and contains only spin-. We set to be the matrix representation of , which is a square matrix.
The above applies to both the Metropolis and the branching random walk approaches. In the Metropolis approach discussed in Sec. 3.2, each path preserves the original spin sectors defined by , and both and will take the same form. In the branching random walk approach discussed in Sec. 3.3, each walker will remain in the same form.
When SOC is included the simple factorization into separate spin components is no longer possible. Instead, each random walker must be in generalized Hartree-Fock (GHF) form consisting of spin-orbitals, which has a matrix representation ,
with . In more compact notation,
where in general the are linear combinations of and , corresponding to Eq. 25. In the absence of SOC, reduces to , and . Because the one body operator couples spin- and spin-, its matrix representation, , is now ,
Note that is zero in the absence of SOC. The ensemble of random walkers, now in the form of Slater determinants, is propagated and sampled until the procedure converges to a stochastic representation of the many-body ground-state.
Observables are computed using the Green’s function between two Slater determinants, which can be measured by,
Each is a matrix, obtained according to,
As with the random walkers and one-body propagators, the application of Wick’s theorem must be modified due to the inclusion of SOC. Spin-flip terms in the Wick expansion, which make no contribution before and can be neglected, must now be included.
We take as an example the calculation of double occupancy,
Without SOC, this observable has the form,
In the presence of SOC it takes the expanded form,
We note that Hamiltonians which include pairing terms can be treated with an additional generalization, where the walker is extend to Hartree-Fock-Bogoliubov space HFB_paper ().
Having provided a brief description of the AFQMC method, we present in the following sections an illustration of its capabilities by highlighting a small set of the many observables that the method makes possible to measure to high accuracy. We discuss results for ultracold atoms and the repulsive electronic model in two separate subsections.
4.1 Fermi atomic gas with synthetic SOC
We begin with a survey of recent results on ultracold atoms, where we will focus on dilute Fermi gas systems in 2D. Further results can be found in 2DFG_SOC_AFQMC (), including results for optical lattices 2DFG_RASHBASOC_OPLATT_AFQMC (). As mentioned in Sec. 2, the Hamiltonian has here. The calculations are performed at the dilute limit, with . In order to map the results to the continuum system of Fermi atomic gases studied in experiment, we first introduce an overall energy scaling defined by the ground-state energy per particle of the corresponding non-interacting Fermi gas, , with being the number density for the 2D lattice. The interaction strength is uniquely defined WernerCastin (); 2DFG_AFQMC () by . In order to compare physically equivalent systems we introduce two dimensionless parameters:
that specify the strengths of the SOC and interaction, respectively. In the above, is the two-body binding energy at and is directly related to 2DFG_AFQMC (). The strong SOC regime, characterized by occupation of only the helicity band, occurs for . In the weak SOC regime (), both and are occupied. There is a smooth transition between the two regimes at . Small corresponds to the BCS limit of the BCS-BEC crossover, where the physics is best described in terms of weakly interacting Cooper pairs, whereas large corresponds to the BEC limit, in which strongly interacting fermions are bound into tight pairs resembling bosons.
Figures 1 and 2 plot two examples of the computed momentum distribution and pairing wavefunctions. The first figure corresponds to a weakly interacting system, with , and the second to a strongly interacting system, with . Both systems are in the modest SOC regime, with . Large lattice sizes are treated here, so that the residual discretization error and finite-size (number of particles) error are small. The results provide quantitative information for the continuum bulk 2D Fermi gas.
The momentum distribution in the weakly interacting system differs only slightly from the non-interacting case, whose regions of occupation and Fermi surfaces (one for each helicity band) are indicated by the shading and vertical dashed lines. At large interaction strength the momentum distribution reveals that there is occupation of higher momentum states that are well beyond the non-interacting region. In addition to becoming considerably broader as a function of interaction strength, the distribution also becomes smoother, showing no evidence of the cusps present at small interaction strengths.
AFQMC also provides access, at a quantitative level, to the rich pairing structure induced by the presence of SOC. The bottom row of figures 1 and 2 present two examples of the pairing wavefunction. SOC mixes the singlet and triplet channels, so the pairing wavefunction has both singlet and triplet components. At small interaction strengths the pairing wavefunction has strong peaks indicating that pairing is concentrated near the Fermi surfaces. Also evident is the fact that the singlet and triplet components of the pairing wavefunction have similar amplitude. This is no longer the case for strongly interacting systems, which favor singlet pairing over triplet pairing because of the enhanced on-site attraction, as reflected by the relative amplitude of the two components of the pairing wavefunction. As a reflection of the behavior of the momentum distribution, the well-localized peaks of the pairing wavefunction also become broad and smooth at large interaction strength, suggesting that pairing occurs across a wide range of momenta.
We can also reliably measure many other observables, including charge-charge and spin-spin correlation functions. Figures 3 and 4 illustrate the effects of SOC on the charge-charge and spin-spin correlations at two different values of interaction strength. At weak interaction (Fig. 3) and weak SOC, the charge-charge correlation exhibits a small amplitude oscillation about the square of the average density. This oscillation is suppressed as SOC strength is increased. The spin-spin correlation shows an opposite effect, with a small oscillation about zero at large SOC strength that is not present for weak SOC. With increased interaction strength (Fig. 4), both correlation functions become shorter ranged with larger amplitude, which is evident in the steeper slopes of each as they approach their asymptotic values.
4.2 Towards real materials — the repulsive Hubbard model with SOC
As a demonstration of the applicability of this approach to electronic systems, we present a study of the repulsive Hubbard model with Rashba SOC. We first examine the energetics, which highlight the capability of CP-AFQMC to achieve high-accuracy results for systems with a sign problem. To benchmark our method we focus on a supercell with electrons, for which we perform exact diagonalization calculations for comparison.
Plotted in Fig. 5 is the total energy as a function of SOC strength at an interaction strength , in comparison with exact results. In this set of simulations the free-electron wavefunction was used as the trial wavefunction (see Eq. (18)). Very good agreement is seen with exact results. The slightly larger discrepancy at small SOC strength can be improved with an improved trial wavefunction, as illustrated in Fig. 6. In this case a generalized Hartree-Fock solution was used as and the CP-AFQMC energy shows improved convergence to the exact energy. The behavior of the algorithm as a function of interaction strength is investigated in Fig. 7, which shows the total energy vs. at a fixed SOC strength of . Nearly exact results are obtained for up to intermediate interaction. Even in the strong interaction limit, the relative error is a small fraction of a percent, which is well within the accuracy of AFQMC seen in Hubbard-like models MingpuHubbardBenchmark (); Hao-symmetry () or in molecules and solids AFQMC_backprop_molecules (); HydrogenBenchmark (); AFQMC_downfolding ().
In addition to ground-state energies, high-accuracy measurements of other observables, including charge, spin and pairing properties, are possible with CP-AFQMC. As indicated in the procedure outlined above, observables can be obtained from the Green’s function, which is calculated according to (31). We show an example of the Green’s function in Fig. 8. Comparison of the CP-AFQMC result with the exact result establishes the high degree of accuracy attainable in CP-AFQMC measurements of observables.
5 Summary and Outlook
Spin-orbit coupling is an essential ingredient of many fascinating phenomena across a wide range of physical contexts. From spintronics to topological insulators to quantum information, a precise characterization of the behavior of systems subject to SOC will not only improve the fundamental understanding of these phenomena and phases, it will also inspire and guide the fabrication of new devices.
The results presented here show the applicability and the accuracy of the AFQMC method, which provides precision, and in some cases, such as the attractive 2D Fermi gas, numerically exact treatments of strongly interacting many-body systems. As we have demonstrated, AFQMC is well suited to study lattice models with either attractive or repulsive interactions. More realistic treatments of materials can be achieved by replacing the on-site Hubbard interaction with generalized (4-index) two-body matrix elements. Interactions of this generalized form can be incorporated using any single-particle basis, such as plane waves or Gaussians, but require a new decomposition and Hubbard-Stratonovich transformation, which introduces complex auxiliary fields and consequently a phase problem. This problem can be systematically controlled using a generalization of the CP approach, the phaseless approximation PRL-phaseless (), which can achieve similar accuracy to what is illustrated above, in realistic material simulations HydrogenBenchmark (); AFQMC_downfolding (). Our discussion of the treatment of SOC carries through straightforwardly, using GHF-type walkers. As we have demonstrated, this approach enables the treatment of SOC with no degradation of the exquisite accuracy achievable in electronic structure calculations with AFQMC.
As the experimental and theoretical landscape continues to expand, there will be many opportunities
and challenges in the study of SOC and strong interaction in many-body systems. These challenges,
including characterizing the effect of interaction on topological phases PhysRevB.89.195124 (); PhysRevB.90.245120 (); annurev-conmatphys-SPT (), will require complementary theoretical, numerical, and experimental progress to generate
new understanding. Given its unique ability to treat strongly-correlated many-body systems with high precision,
AFQMC will play a vital role in the combined effort.
Acknowledgments This research was supported by NSF (grant no. DMR-1409510), and the Simons Foundation. Computing was carried out at the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF grant number ACI-1053575, and the computational facilities at the College of William and Mary. The Flatiron Institute is supported by the Simons Foundation.
A. Machon, H. C. Koo, J. Nitta, S. Frolov, R. Duine,
New perspectives for Rashba
spin-orbit coupling, Nature Materials 14 (2015) 871–882.
A. Soumyanarayanan, N. Reyren, A. Fert, C. Panagopoulos,
Emergent phenomena induced by
spin-orbit coupling at surfaces and interfaces, Nature 539 (2016) 509–517.
D. Bercioux, P. Lucignano,
Quantum transport in
Rashba spinÐorbit materials: a review, Reports on Progress in Physics
78 (10) (2015) 106001.
H. C. Koo, J. H. Kwon, J. Eom, J. Chang, S. H. Han, M. Johnson,
Control of Spin
Precession in a Spin-Injected Field Effect Transistor, Science 325 (5947)
C. Chappert, A. Fert, F. N. Van Dau,
The emergence of spin electronics
in data storage, Nature Materials 6 (2007) 813–823.
D. D. Awschalom, L. C. Bassett, A. S. Dzurak, E. L. Hu, J. R. Petta,
Spintronics: Engineering and Manipulating Atom-Like Spins in
Semiconductors, Science 339 (6124) (2013) 1174–1179.
I. Bloch, J. Dalibard, W. Zwerger,
with ultracold gases, Rev. Mod. Phys. 80 (2008) 885–964.
M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen(De), U. Sen,
Ultracold atomic gases in
optical lattices: mimicking condensed matter physics and beyond, Advances
in Physics 56 (2) (2007) 243–379.
A. A. Orel, P. Dyke, M. Delehaye, C. J. Vale, H. Hu,
of a trapped two-dimensional strongly interacting Fermi gas, New Journal of
Physics 13 (11) (2011) 113032.
H. Shi, S. Chiesa, S. Zhang,
properties of strongly interacting Fermi gases in two dimensions, Phys.
Rev. A 92 (2015) 033603.
G. Bertaina, S. Giorgini,
Crossover in a Two-Dimensional Fermi Gas, Phys. Rev. Lett. 106 (2011)
E. R. Anderson, J. E. Drut,
Compressibility, and Contact of the Two-Dimensional Attractive Fermi Gas,
Phys. Rev. Lett. 115 (2015) 115301.
D. Jaksch, P. Zoller,
Creation of effective
magnetic fields in optical lattices: the Hofstadter butterfly for cold
neutral atoms, New Journal of Physics 5 (1) (2003) 56.
K. Osterloh, M. Baig, L. Santos, P. Zoller, M. Lewenstein,
Cold Atoms in
Non-Abelian Gauge Potentials: From the Hofstadter ”Moth” to Lattice Gauge
Theory, Phys. Rev. Lett. 95 (2005) 010403.
K. Jiménez-García, L. J. LeBlanc, R. A. Williams, M. C. Beeler, A. R.
Perry, I. B. Spielman,
Substitution in an Engineered Lattice Potential, Phys. Rev. Lett. 108
M. Aidelsburger, M. Atala, S. Nascimbène, S. Trotzky, Y.-A. Chen, I. Bloch,
Realization of Strong Effective Magnetic Fields in an Optical Lattice,
Phys. Rev. Lett. 107 (2011) 255301.
H. Miyake, G. A. Siviloglou, C. J. Kennedy, W. C. Burton, W. Ketterle,
Harper Hamiltonian with Laser-Assisted Tunneling in Optical Lattices, Phys.
Rev. Lett. 111 (2013) 185302.
F. Gerbier, J. Dalibard,
Gauge fields for
ultracold atoms in optical superlattices, New Journal of Physics 12 (3)
Y.-J. Lin, R. L. Compton, K. Jimenez-Garcia, J. V. Porto, I. B. Spielman,
Synthetic magnetic fields for
ultracold neutral atoms, Nature 462 (2009) 628–632.
Y.-J. Lin, K. Jimenez-Garcia, I. B. Spielman,
Bose-Einstein condensates, Nature 471 (2011) 83–86.
L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yefsah, W. S. Bakr, M. W.
Spectroscopy of a Spin-Orbit Coupled Fermi Gas, Phys. Rev. Lett. 109 (2012)
P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. Zhai, J. Zhang,
Coupled Degenerate Fermi Gases, Phys. Rev. Lett. 109 (2012) 095301.
L. Huang, Z. Meng, P. Wang, P. Peng, S.-L. Zhang, L. Chen, D. Li,
Q. Zhou, J. Zhang,
Experimental realization of a
two-dimensional synthetic spin-orbit coupling in ultracold Fermi gases,
Nature Physics 12 (2016) 540.
T. Jungwirth, J. Wunderlich, K. Olejnik,
Spin Hall effect devices, Nature
Materials 11 (2012) 382–390.
M. Qin, H. Shi, S. Zhang,
of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo
method, Phys. Rev. B 94 (2016) 085103.
H. Shi, S. Zhang,
auxiliary-field quantum Monte Carlo calculations, Phys. Rev. B 88 (2013)
- (27) M. Motta, S. Zhang, Computation of ground-state properties in molecular systems: back-propagation with auxiliary-field quantum Monte Carlo, ArXiv e-printsarXiv:1707.02684.
- (28) M. Motta, D. M. Ceperley, G. Kin-Lic Chan, J. A. Gomez, E. Gull, S. Guo, C. Jimenez-Hoyos, T. Nguyen Lan, J. Li, F. Ma, A. J. Millis, N. V. Prokof’ev, U. Ray, G. E. Scuseria, S. Sorella, E. M. Stoudenmire, Q. Sun, I. S. Tupitsyn, S. R. White, D. Zgid, S. Zhang, Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods, ArXiv e-printsarXiv:1705.01608.
W. Purwanto, S. Zhang, H. Krakauer,
Auxiliary-field quantum Monte
Carlo calculations of the molybdenum dimer, The Journal of Chemical Physics
144 (24) (2016) 244306.
F. Ma, W. Purwanto, S. Zhang, H. Krakauer,
Carlo Calculations in Solids with Downfolded Hamiltonians, Phys. Rev. Lett.
114 (2015) 226401.
- (31) E. Vitali, H. Shi, M. Qin, S. Zhang, Visualizing the BEC-BCS crossover in the two-dimensional Fermi gas: pairing gaps and dynamical response functions from ab initio computations, ArXiv e-printsarXiv:1705.07929.
H. Shi, P. Rosenberg, S. Chiesa, S. Zhang,
Spin-Orbit Coupling, Strong Interactions, and the BCS-BEC Crossover in the
Ground State of the Two-Dimensional Fermi Gas, Phys. Rev. Lett. 117 (2016)
- (33) P. Rosenberg, H. Shi, S. Zhang, Ultracold atoms in a square lattice with spin-orbit coupling: Charge order, superfluidity, and topological signatures, ArXiv e-printsarXiv:1707.02994.
C. Wu, S.-C. Zhang,
condition for absence of the sign problem in the fermionic quantum Monte
Carlo algorithm, Phys. Rev. B 71 (2005) 155115.
L. W. Cheuk, M. A. Nichols, K. R. Lawrence, M. Okan, H. Zhang, M. W. Zwierlein,
2D Fermionic Mott Insulators of with Single-Site
Resolution, Phys. Rev. Lett. 116 (2016) 235301.
D. Greif, M. F. Parsons, A. Mazurenko, C. S. Chiu, S. Blatt, F. Huber, G. Ji,
imaging of a fermionic Mott insulator, Science 351 (6276) (2016) 953–957.
- (37) L. W. Cheuk, M. A. Nichols, M. Okan, T. Gersdorf, V. V. Ramasesh, W. S. Bakr, T. Lompe, M. W. Zwierlein, Quantum-gas microscope for fermionic atoms, Physical Review Letters 114 (19) (2015) 1–5. arXiv:1503.02648, doi:10.1103/PhysRevLett.114.193001.
E. Haller, J. Hudson, A. Kelly, D. A. Cotta, B. Peaudecerf, G. D. Bruce,
S. Kuhr, Single-atom imaging of
fermions in a quantum-gas microscope, Nat Phys 11 (2015) 738–742.
M. F. Parsons, F. Huber, A. Mazurenko, C. S. Chiu, W. Setiawan,
K. Wooley-Brown, S. Blatt, M. Greiner,
imaging of fermionic in an optical lattice, Phys. Rev.
Lett. 114 (2015) 213002.
R. A. Hart, P. M. Duarte, T.-L. Yang, X. Liu, T. Paiva, E. Khatami, R. T.
Scalettar, N. Trivedi, D. A. Huse, R. G. Hulet,
antiferromagnetic correlations in the Hubbard model with ultracold atoms,
Nature 519 (2015) 211–214.
S. Zhang, J. Carlson, J. E. Gubernatis,
Monte Carlo method for fermion ground states, Phys. Rev. B 55 (1997)
- (42) S. Zhang, in: E. Pavarini, E. Koch, U. Schollwöck (Eds.), Emergent Phenomena in Correlated Matter: Modeling and Simulation, Vol. 3, Verlag des Forschungszentrum Jülich, Jülich, Germany, 2013.
- (43) F. F. Assaad, Quantum Monte Carlo Methods on Lattices: The Determinantal Approach, in: J. Grotendorst, D. Marx, A. Muramatsu (Eds.), Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms, Vol. 10, NIC, Jülich, Germany, 2002, pp. 99–156.
G. Sugiyama, S. Koonin,
field Monte-Carlo for quantum many-body ground states, Annals of Physics
168 (1) (1986) 1 – 26.
Partition Functions, Phys. Rev. Lett. 3 (1959) 77–78.
J. E. Hirsch,
Hubbard-Stratonovich transformation for fermion lattice models, Phys. Rev.
B 28 (1983) 4059–4061.
H. Shi, S. Zhang,
in fermion quantum Monte Carlo calculations, Phys. Rev. E 93 (2016) 033303.
W. Purwanto, S. Zhang,
Carlo method for the ground state of many-boson systems, Phys. Rev. E 70
H. Shi, S. Zhang,
computations by stochastic sampling in Hartree-Fock-Bogoliubov space, Phys.
Rev. B 95 (2017) 045144.
F. Werner, Y. Castin,
for quantum gases in two and three dimensions: Two-component fermions,
Phys. Rev. A 86 (2012) 013626.
S. Zhang, H. Krakauer,
Carlo Method using Phase-Free Random Walks with Slater Determinants, Phys.
Rev. Lett. 90 (2003) 136401.
C. Wang, T. Senthil,
fermionic topological insulators/superconductors in three dimensions, Phys.
Rev. B 89 (2014) 195124.
Y.-Z. You, C. Xu,
topological states of interacting fermions and bosons, Phys. Rev. B 90
Topological Phases of Quantum Matter, Annual Review of Condensed Matter
Physics 6 (1) (2015) 299–324.