Source author record

J. K. Freericks

J. K. Freericks appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el quant-ph cond-mat.quant-gas cond-mat.supr-con cond-mat.mtrl-sci cond-mat.stat-mech physics.atom-ph cond-mat.mes-hall cond-mat.other cond-mat physics.optics Computational Complexity Computational Engineering, Finance, and Science math-ph math.MP physics.chem-ph supr-con

Catalog footprint

What is connected

83works

17topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Hidden local adiabatic ramp in the modulated time evolution and the quantum approximate optimization algorithm

Adiabatic state preparation provides an analytical solution for generating the ground state of a target Hamiltonian, starting from an easily prepared ground state of the initial Hamiltonian. While effective for time-dependent Hamiltonians with an energy gap to the first coupled excited state, the process becomes exceedingly slow as the gap narrows. Rather than strictly following the adiabatic theorem, a more robust approach allows controlled diabatic excitations during the evolution and numerically optimizes the path to eliminate these excitations by the end. In this work, this is achieved via modulated time evolution, using a time dependent oscillating field lambda(t) to modulate the Hamiltonian, in conjunction with a transverse field B(t) whose optimized shape closely resembles a local adiabatic ramp. Beyond modulated time evolution, the quantum approximate optimization algorithm (QAOA), which also employs a transverse field defined as beta(t) over gamma(t), exhibits a shape similar to the local adiabatic ramp. This resemblance offers a more intuitive and physically motivated way to understand the QAOA algorithm through the lens of time evolution.

preprint2025arXiv

Out-of-time ordered correlation functions for the localized $f$ electrons in the Falicov-Kimball model

We provide an exact evaluation of the out-of-time correlation (OTOC) functions for the localized $f$-particle states in the Falicov-Kimball model within dynamical mean-field theory. Different regimes of quantum chaos and quantum scrambling are distinguished by the winding numbers of the block Toeplitz matrices used in the calculation. The similarities of these fermionic OTOCs and their logarithmic derivatives for time evolution with the OTOCs for quantum spin models with disorder are also discussed.

preprint2022arXiv

Fixed Depth Hamiltonian Simulation via Cartan Decomposition

Simulating quantum dynamics on classical computers is challenging for large systems due to the significant memory requirements. Simulation on quantum computers is a promising alternative, but fully optimizing quantum circuits to minimize limited quantum resources remains an open problem. We tackle this problem presenting a constructive algorithm, based on Cartan decomposition of the Lie algebra generated by the Hamiltonian, that generates quantum circuits with time-independent depth. We highlight our algorithm for special classes of models, including Anderson localization in one dimensional transverse field XY model, where a O(n^2)-gate circuits naturally emerge. Compared to product formulas with significantly larger gate counts, our algorithm drastically improves simulation precision. In addition to providing exact circuits for a broad set of spin and fermionic models, our algorithm provides broad analytic and numerical insight into optimal Hamiltonian simulations.

preprint2021arXiv

Decomposition of high-rank factorized unitary coupled-cluster operators using ancilla and multi-qubit controlled low-rank counterparts

The factorized form of the unitary coupled-cluster approximation is one of the most promising methodologies to prepare trial states for strongly correlated systems within the variational quantum eigensolver framework. The factorized form of the UCC ansatz can be systematically applied to a reference state to generate the desired entanglement. The difficulty associated with such an approach is the requirement of simultaneously entangling a growing number of qubits, which quickly exceeds the hardware limitations of today's quantum machines. In particular, while circuits for singles and double excitations can be performed on current hardware, higher-rank excitations require too many gate operations. In this work, we propose a set of new schemes that trade off using extra qubits for a reduced gate depth to decompose high-rank UCC excitation operators into significantly lower depth circuits. These results will remain useful even when fault-tolerant machines are available to reduce the overall state-preparation circuit depth.

preprint2021arXiv

Flexibility of the factorized form of the unitary coupled cluster ansatz

The factorized form of the unitary coupled cluster ansatz is a popular state preparation ansatz for electronic structure calculations of molecules on quantum computers. It often is viewed as an approximation (based on the Trotter product formula) for the conventional unitary coupled cluster operator. In this work, we show that the factorized form is quite flexible, allowing one to range from conventional configuration interaction, to conventional unitary coupled cluster, to efficient approximations that lie in between these two. The variational minimization of the energy often allows simpler factorized unitary coupled cluster approximations to achieve high accuracy, even if they do not accurately approximate the Trotter product formula. This is similar to how quantum approximate optimization algorithms can achieve high accuracy with a small number of levels.

preprint2021arXiv

Native multiqubit Toffoli gates on ion trap quantum computers

We examine the detailed scenario for implementing n-control-qubit Toffoli gates and select gates on ion-trap quantum computers, especially those that shuttle ions into interaction zones. We determine expected performance of these gates with realistic parameters for an ion-trap quantum computer and taking into account the time variation of the exchange integrals. This allows us to estimate the errors due to spin-phonon entanglement as well. While there are challenges with implementing these gates, because their performance always has some degree of error, they should be feasible on current hardware, but they may be too slow to be used efficiently in quantum codes on noisy intermediate scale quantum computers.

preprint2020arXiv

Algebraic derivation of Kramers-Pasternack relations based on the Schrodinger factorization method

The Kramers-Pasternack relations are used to compute the moments of r (both positive and negative) for all radial energy eigenfunctions of hydrogenic atoms. They consist of two algebraic recurrence relations, one for positive powers and one for negative. Most derivations employ the Feynman-Hellman theorem or a brute-force integration to determine the second inverse moment, which is needed to complete the recurrence relations for negative moments. In this work, we show both how to derive the recurrence relations algebraically and how to determine the second inverse moment algebraically, which removes the pedagogical confusion associated with differentiating the Hamiltonian with respect to the angular momentum quantum number l in order to find the inverse second moment.

preprint2020arXiv

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of non-interacting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final non-equilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to 5 Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. The methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

preprint2020arXiv

Many Body Thermodynamics on Quantum Computers via Partition Function Zeros

Interacting quantum systems illustrate complex phenomena including phase transitions to novel ordered phases. The universal nature of critical phenomena reduces their description to determining only the transition temperature and the critical exponents. Numerically calculating these results for systems in new universality classes is complicated due to critical slowing down, requiring increasing resources near the critical point. An alternative approach analytically continues the calculation onto the complex plane and determines the partition function via its zeros. Here we show how to robustly perform this analysis on noisy intermediate scale trapped ion quantum computers in a scalable manner, using the XXZ model as a prototype. We illustrate the transition from XY-like behavior to Ising-like behavior as a function of the anisotropy. While quantum computers cannot yet scale to the thermodynamic limit, our work provides a pathway to do so as hardware improves, allowing the determination of critical phenomena for systems that cannot be solved otherwise.

preprint2020arXiv

Quantum computation of magnon spectra

We demonstrate quantum computation of two-point correlation functions for a Heisenberg spin chain. Using the IBM Q 20 quantum machines, we find that for two sites the correlation functions produce the exact results reliably. For four sites, results from the IBM Q 20 Tokyo quantum computer are noisy due to read out errors and decoherence. Nevertheless, the correlation functions retain the correct spectral information. This is illustrated in the frequency domain by accurately extracting the magnon energies from peaks in the spectral function.

preprint2020arXiv

Test of the unitary coupled-cluster variational quantum eigensolver for a simple strongly correlated condensed-matter system

The variational quantum eigensolver has been proposed as a low-depth quantum circuit that can be employed to examine strongly correlated systems on today's noisy intermediate-scale quantum computers. We examine details associated with the factorized form of the unitary coupled-cluster variant of this algorithm. We apply it to a simple strongly correlated condensed-matter system with nontrivial behavior---the four-site Hubbard model at half filling. This work show some of the subtle issues one needs to take into account when applying this algorithm in practice, especially to condensed-matter systems.

preprint2016arXiv

Constant Matrix Element Approximation to Time-Resolved Angle-Resolved Photoemission Spectroscopy

We describe different approximations associated with employing a constant matrix element for the coupling of light to multiband electrons in the context of time-resolved angle-resolved photoemission spectroscopy (TR-ARPES). In particular, we demonstrate that the constant matrix approximation only holds for one of the coordinate representations, and changing to other bases, requires including nonconstant corrections to the matrix element. We also discuss some simplifying approximations, where a constant matrix element is employed in multiple bases, and the consequences of this further approximation (especially with respect to the TR-ARPES signal no longer being nonnegative). We also discuss issues related to gauge invariance of the final spectra.

preprint2016arXiv

Direct observation of Higgs mode oscillations in the pump-probe photoemission spectra of electron-phonon mediated superconductors

Using the non-equilibrium Keldysh formalism, we solve the equations of motion for electron-phonon superconductivity, including an ultrafast pump field. We present results for time-dependent photoemission spectra out of equilibrium which probes the dynamics of the superconducting gap edge. The partial melting of the order by the pump field leads to oscillations at twice the melted gap frequency, a hallmark of the Higgs or amplitude mode. Thus the Higgs mode can be directly excited through the nonlinear effects of an electromagnetic field and detected without any additional symmetry breaking.

preprint2016arXiv

Energy dissipation in the time domain governed by bosons in a correlated material

In complex materials various interactions play important roles in determining the material properties. Angle Resolved Photoelectron Spectroscopy (ARPES) has been used to study these processes by resolving the complex single particle self energy $Σ(E)$ and quantifying how quantum interactions modify bare electronic states. However, ambiguities in the measurement of the real part of the self energy and an intrinsic inability to disentangle various contributions to the imaginary part of the self energy often leave the implications of such measurements open to debate. Here we employ a combined theoretical and experimental treatment of femtosecond time-resolved ARPES (tr-ARPES) and show how measuring the population dynamics using tr-ARPES can be used to separate electron-boson interactions from electron-electron interactions. We demonstrate the analysis of a well-defined electron-boson interaction in the unoccupied spectrum of the cuprate Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8+x}$ characterized by an excited population decay time constant $τ_{QP}$ that maps directly to a discrete component of the equilibrium self energy not readily isolated by static ARPES experiments.

preprint2016arXiv

Exact time evolution of the asymmetric Hubbard dimer

We examine the time evolution of an asymmetric Hubbard dimer, which has a different on-site interaction on the two sites. The Hamiltonian has a time-dependent hopping term, which can be employed to describe an electric field (which creates a Hamiltonian with complex matrix elements), or it can describe a modulation of the lattice (which has real matrix elements). By examining the symmetries under spin and pseudospin, we show that the former case involves at most a 3 x 3 block---it can be mapped onto the time evolution of a time-independent Hamiltonian, so the dynamics can be evaluated analytically and exactly (by solving a nontrivial cubic equation). We also show that the latter case reduces to at most 2 x 2 blocks, and hence the time evolution for a single Trotter step can be determined exactly, but the time evolution generically requires a Trotter product.

preprint2016arXiv

Generalized gradient expansion for inhomogeneous dynamical mean-field theory: Application to ultracold atoms in a harmonic trap

We develop a generalized gradient expansion of the inhomogeneous dynamical mean-field theory method for determining properties of ultracold atoms in a trap. This approach goes beyond the well-known local density approximation and at higher temperatures, in the normal phase, it shows why the local density approximation works so well, since the local density and generalized gradient approximations are essentially indistinguishable from each other (and from the exact solution within full inhomogeneous dynamical mean-field theory). But because the generalized gradient expansion only involves nearest-neighbor corrections, it does not work as well at low temperatures, when the systems enter into ordered phases. This is primarily due to the problem that ordered phases often satisfy some global constraints which determine the spatial ordering pattern, and the local density and generalized gradient approximations are not able to impose those kinds of constraints; they also overestimate the tendency to order. The theory is applied to phase separation of different mass fermionic mixtures represented by the Falicov-Kimball model and to determining the entropy per particle of a fermionic system represented by the Hubbard model. The generalized gradient approximation is a useful diagnostic for the accuracy of the local density approximation---when both methods agree, they are likely accurate, when they disagree, neither is likely to be correct.

preprint2016arXiv

The role of average time dependence on the relaxation of excited electron populations in nonequilibrium many-body physics

We examine the exact equation of motion for the relaxation of populations of strongly correlated electrons after a nonequilibrium excitation by a pulsed field, and prove that the populations do not change when the Green's functions have no average time dependence. We show how the average time dependence enters into the equation of motion to lowest order and describe what governs the relaxation process of the electron populations in the long-time limit. While this result may appear, on the surface, to be required by any steady-state solution, the proof is nontrivial, and provides new critical insight into how nonequilibrium populations relax, which goes beyond the assumption that they thermalize via a simple relaxation rate determined by the imaginary part of the self-energy, or that they can be described by a quasi-equilibrium condition with a Fermi-Dirac distribution and a time-dependent temperature. We also discuss the implications of this result to approximate theories, which may not satisfy the exact relation in the equation of motion.

preprint2015arXiv

Creating thermal distributions from diabatic excitations in ion-trap-based quantum simulation

The goal of adiabatic ground-state preparation is to start a simple quantum system in its ground state and adiabatically evolve the Hamiltonian to a complex one, maintaining the ground state throughout the evolution. In ion-trap-based quantum simulations, coherence times are too short to allow for adiabatic evolution for large chains, so the system evolves diabatically, creating excitations to higher energy states. Because the probability for diabatic excitation depends exponentially on the excitation energy and because the thermal distribution also depends exponentially on the excitation energy, we investigate whether the diabatic excitation can create a thermal distribution; as this could serve as an alternative for creating thermal states of complex quantum systems without requiring contact with a heat bath. In this work, we explore this relationship and determine situations where diabatic excitation can approximately create such a thermal state.

preprint2015arXiv

Determining the ground-state probability of a quantum simulation with product-state measurements

One of the goals in quantum simulation is to adiabatically generate the ground state of a complicated Hamiltonian by starting with the ground state of a simple Hamiltonian and slowly evolving the system to the complicated one. If the evolution is adiabatic and the initial and final ground states are connected due to having the same symmetry, then the simulation will be successful. But in most experiments, adiabatic simulation is not possible because it would take too long, and the system has some level of diabatic excitation. In this work, we quantify the extent of the diabatic excitation even if we do not know {\it a priori} what the complicated ground state is. Since many quantum simulator platforms, like trapped ions, can measure the probabilities to be in a product state, we describe techniques that can employ these measurements to estimate the probability of being in the ground state of the system after the diabatic evolution. These techniques do not require one to know any properties about the Hamiltonian itself, nor to calculate its eigenstate properties. All the information is derived by analyzing the product-state measurements as functions of time.

preprint2015arXiv

Exact calculation of phonon effects on spin squeezing

Theoretical models of spins coupled to bosons provide a simple setting for studying a broad range of important phenomena in many-body physics, from virtually mediated interactions to decoherence and thermalization. In many atomic, molecular, and optical systems, such models also underlie the most successful attempts to engineer strong, long-ranged interactions for the purpose of entanglement generation. Especially when the coupling between the spins and bosons is strong---such that it cannot be treated perturbatively---the properties of such models are extremely challenging to calculate theoretically. Here, exact analytical expressions for nonequilibrium spin-spin correlation functions are derived for a specific model of spins coupled to bosons. The spatial structure of the coupling between spins and bosons is completely arbitrary, and thus the solution can be applied to systems in any number of dimensions. The explicit and nonperturbative inclusion of the bosons enables the study of entanglement generation (in the form of spin squeezing) even when the bosons are driven strongly and near-resonantly, and thus provides a quantitative view of the breakdown of adiabatic elimination that inevitably occurs as one pushes towards the fastest entanglement generation possible. The solution also helps elucidate the effect of finite temperature on spin squeezing. The model considered is relevant to a variety of atomic, molecular, and optical systems, such as atoms in cavities or trapped ions. As an explicit example, the results are used to quantify phonon effects in trapped ion quantum simulators, which are expected to become increasingly important as these experiments push towards larger numbers of ions.

preprint2015arXiv

Measuring nonequilibrium retarded spin-spin Green's functions in an ion-trap based quantum simulator

Recent work proposed a variant on Ramsey interferometry for coupled spin-$1/2$ systems that directly measures the retarded spin-spin Green's function. We expand on that work by investigating nonequilibrium retarded spin-spin Green's functions within the transverse-field Ising model. We derive the lowest four spectral moments to understand the short-time behavior and we employ a Lehmann-like representation to determine the spectral behavior. We simulate a Ramsey protocol for a nonequilibrium quantum spin system that consists of a coherent superposition of the ground state and diabatically excited higher-energy states via a temporally ramped transverse magnetic field. We then apply the Ramsey spectroscopy protocol to the final Hamiltonian, which has a constant transverse field. The short-time behavior directly relates to Lieb-Robinson bounds for the transport of many-body correlations, while the long-time behavior relates to the excitation spectra of the Hamiltonian. Compressive sensing is employed in the data analysis to efficiently extract that spectra.

preprint2015arXiv

Nonequilibrium dynamical mean-field theory for the charge-density-wave phase of the Falicov-Kimball model

Nonequilibrium dynamical mean-field theory (DMFT) is developed for the case of the charge-density-wave ordered phase. We consider the spinless Falicov-Kimball model which can be solved exactly. This strongly correlated system is then placed in an uniform external dc electric field. We present a complete derivation for nonequilibrium dynamical mean-field theory Green's functions defined on the Keldysh-Schwinger time contour. We also discuss numerical issues involved in solving the coupled equations.

preprint2015arXiv

Nonequilibrium response of an electron mediated charge-density-wave-ordered material to a large dc electric field

Using the Kadanoff-Baym-Keldysh formalism, we employ nonequilibrium dynamical mean-field theory to exactly solve for the nonlinear response of an electron-mediated charge-density-wave-ordered material. We examine both the dc current and the order parameter of the conduction electrons as the ordered system is driven by the electric field. Although the formalism we develop applies to all models, for concreteness, we examine the charge-density-wave phase of the Falicov-Kimball model, which displays a number of anomalous behaviors including the appearance of subgap density of states as the temperature increases. These subgap states should have a significant impact on transport properties, particularly the nonlinear response of the system to a large dc electric field.

preprint2015arXiv

Theoretical basis for quantum simulation with a planar ionic crystal in a Penning trap using a triangular rotating wall

One of the challenges with quantum simulation in ion traps is that the effective spin-spin exchange couplings are not uniform across the lattice. This can be particularly important in Penning trap realizations where the presence of an ellipsoidal boundary at the edge of the trap leads to dislocations in the crystal. By adding an additional anharmonic potential to better control interion spacing, and a triangular shaped rotating wall potential to reduce the appearance of dislocations, one can achieve better uniformity of the ionic positions. In this work, we calculate the axial phonon frequencies and the spin-spin interactions driven by a spin-dependent optical dipole force, and discuss what effects the more uniform ion spacing has on the spin simulation properties of Penning trap quantum simulators. Indeed, we find the spin-spin interactions behave more like a power law for a wide range of parameters.

preprint2014arXiv

Creation of two-dimensional coulomb crystals of ions in oblate Paul traps for quantum simulations

We develop the theory to describe the equilibrium ion positions and phonon modes for a trapped ion quantum simulator in an oblate Paul trap that creates two-dimensional Coulomb crystals in a triangular lattice. By coupling the internal states of the ions to laser beams propagating along the symmetry axis, we study the effective Ising spin-spin interactions that are mediated via the axial phonons and are less sensitive to ion micromotion. We find that the axial mode frequencies permit the programming of Ising interactions with inverse power law spin-spin couplings that can be tuned from uniform to $r^{-3}$ with DC voltages. Such a trap could allow for interesting new geometrical configurations for quantum simulations on moderately sized systems including frustrated magnetism on triangular lattices or Aharonov-Bohm effects on ion tunneling. The trap also incorporates periodic boundary conditions around loops which could be employed to examine time crystals.

preprint2014arXiv

Diabatic ramping spectroscopy of many-body excited states for trapped-ion quantum simulators

Due to the experimental time constraints of state of the art quantum simulations with trapped ions, the direct preparation of the ground state by adiabatically ramping the field of a transverse field Ising model becomes more and more difficult as the number of particles increase. We propose a spectroscopy protocol that intentionally creates excitations through diabatic ramping of the transverse field and measures a low-noise observable as a function of time for a constant field to reveal the structure of the coherent dynamics of the resulting many-body states. To simulate the experimental data, noise from counting statistics and decoherence error are added. Compressive sensing is then applied to Fourier transform the simulated data into the frequency domain and extract the the low-lying energy excitation spectrum. By using compressive sensing, the amount of data in time needed to extract this energy spectrum is sharply reduced making such experiments feasible with current technology.

preprint2014arXiv

Effect of dynamical spectral weight redistribution on effective interactions in time-resolved spectroscopy

The redistribution of electrons in an ultrafast pump-probe experiment causes significant changes to the effective interaction between electrons and bosonic modes. We study the influence of these changes on pump-probe photoemission spectroscopy for a model electron-phonon coupled system using the nonequilibrium Keldysh formalism. We show that spectral rearrangement due to the driving field preserves an overall sum rule for the electronic self-energy, but modifies the effective electron-phonon scattering as a function of energy. Experimentally, this pump-modified scattering can be tracked by analyzing the fluence or excitation energy dependence of population decay rates and transient changes in dispersion kinks.

preprint2014arXiv

Feshbach modulation spectroscopy

In the vicinity of a Feshbach resonance, a system of ultracold atoms on an optical lattice undergoes rich physical transformations which involve molecule formation and hopping of molecules on the lattice and thus goes beyond a single-band Hubbard model description. We propose to probe the behavior of this system with a harmonic modulation of the magnetic field, and thus of the scattering length, across the Feshbach resonance, as an alternative to lattice-depth modulation spectroscopy. In the regime in which the single-band Hubbard model is still valid, we provide simulation data for this type of spectroscopy. The method may uncover a route towards the efficient creation of ultracold molecules and provides an alternate means for lattice modulation spectroscopy.

preprint2014arXiv

Gauge invariance in the theoretical description of time-resolved angle-resolved pump/probe photoemission spectroscopy

Nonequilibrium calculations in the presence of an electric field are usually performed in a gauge, and need to be transformed to reveal the gauge-invariant observables. In this work, we discuss the issue of gauge invariance in the context of time-resolved angle-resolved pump/probe photoemission. If the probe is applied while the pump is still on, one must ensure that the calculations of the observed photocurrent are gauge invariant. We also discuss the requirement of the photoemission signal to be positive and the relationship of this constraint to gauge invariance. We end by discussing some technical details related to the perturbative derivation of the photoemission spectra, which involve processes where the pump pulse photoexcites electrons due to nonequilibrium effects.

preprint2014arXiv

Intrinsic anharmonic effects on the phonon frequencies and effective spin-spin interactions in a quantum simulator made from trapped ions in a linear Paul trap

The Coulomb repulsion between ions in a linear Paul trap give rise to anharmonic terms in the potential energy when expanded about the equilibrium positions. We examine the effect of these anharmonic terms on the accuracy of a quantum simulator made from trapped ions. To be concrete, we consider a linear chain of $\text{Yb}^{171+}$ ions stabilized close to the zigzag transition. We find that for typical experimental temperatures, frequencies change by no more than a factor of $0.01\%$ due to the anharmonic couplings. Furthermore, shifts in the effective spin-spin interactions (driven by a spin-dependent optical dipole force) also tend to be small for detunings to the blue of the transverse center-of-mass frequency. However, detuning the spin interactions near other frequencies can lead to nonnegligible anharmonic contributions to the effective spin-spin interactions. We also examine an odd behavior exhibited by the harmonic spin-spin interactions for a range of intermediate detunings, where nearest neighbor spins with a larger spatial separation on the ion chain interact more strongly than nearest neighbors with a smaller spatial separation.

preprint2014arXiv

Long-lived nonequilibrium states in the Hubbard model with an electric field

We study the single-band Hubbard model in the presence of a large spatially uniform electric field out of equilibrium. Using the Keldysh nonequilibrium formalism, we solve the problem using perturbation theory in the Coulomb interaction U. We present numerical results for the charge current, the total energy of the system and the double occupancy on an infinite-dimensional hypercubic lattice with nearest-neighbor hopping. The system is isolated from an external bath and is in the paramagnetic state. We show that an electric field pulse can drive the system to a steady nonequilibrium state, which does not evolve into a thermal state. We compare results obtained within second-order perturbation theory (SOPT), self-consistent second-order perturbation theory (SCSOPT) and iterated perturbation theory (IPT). We also discuss the importance of initial conditions for a system which is not coupled to an external bath.

preprint2014arXiv

Nonequilibrium sum rules for the Holstein model

We derive frequency moment sum rules for the retarded electronic Green's function and self-energy of the Holstein model for both equilibrium and nonequilibrium cases. We also derive sum rules for the phonon propagator in equilibrium and nonequilibrium. These sum rules allow one to benchmark nonequilibrium calculations and help with interpreting the behavior of electrons driven out of equilibrium by an applied electric field. We exactly evaluate the sum rules when the system is in the atomic limit. We also discuss the application of these sum rules to pump/probe experiments like time-resolved angle-resolved photoemission spectroscopy.

preprint2014arXiv

Robust finite-temperature disordered Mott insulating phases in inhomogeneous Fermi-Fermi mixtures with density and mass imbalance

Ultracold mixtures of different atomic species have great promise for realizing novel many-body phenomena. In a binary mixture of femions with a large mass difference and repulsive interspecies interactions, a disordered Mott insulator phase can occur. This phase displays an incompressible total density, although the relative density remain compressible. We use strong-coupling and Monte Carlo calculations to show that this phase exists for a broad parameter region for ultracold gases confined in a harmonic trap on a three-dimensional optical lattice, for experimentally accessible values of the trap parameters.

preprint2014arXiv

The Nonequilibrium Many-Body Problem as a paradigm for extreme data science

Generating big data pervades much of physics. But some problems, which we call extreme data problems, are too large to be treated within big data science. The nonequilibrium quantum many-body problem on a lattice is just such a problem, where the Hilbert space grows exponentially with system size and rapidly becomes too large to fit on any computer (and can be effectively thought of as an infinite-sized data set). Nevertheless, much progress has been made with computational methods on this problem, which serve as a paradigm for how one can approach and attack extreme data problems. In addition, viewing these physics problems from a computer-science perspective leads to new approaches that can be tried to solve them more accurately and for longer times. We review a number of these different ideas here.

preprint2014arXiv

Theoretical Description of Coherent Doublon Creation via Lattice Modulation Spectroscopy

Using a recently developed strong-coupling method, we present a comprehensive theory for doublon production processes in modulation spectroscopy of a three-dimensional system of ultracold fermionic atoms in an optical lattice with a trap. The theoretical predictions compare well to the experimental time traces of doublon production. For experimentally feasible conditions, we provide a quantitative prediction for the presence of a nonlinear "two-photon" excitation at strong modulation amplitudes.

preprint2014arXiv

Theory of pump-probe photoemission in graphene: Ultrafast tuning of Floquet bands and local pseudospin textures

The control of physical properties of solids with short laser pulses is an intriguing prospect of ultrafast materials science. Continuous-wave high-frequency laser driving with circular polarization was predicted to induce a light-matter coupled new state possessing a quasi-static band structure with an energy gap and a quantum Hall effect, coined "Floquet topological insulator". Whereas the envisioned Floquet topological insulator requires well separated Floquet bands and therefore high-frequency pumping, a natural follow-up question regards the creation of Floquet-like states in graphene with realistic pump laser pulses. Here we predict that with short low-frequency laser pulses attainable in pump-probe experiments, states with local spectral gaps at the Dirac points and novel pseudospin textures can be achieved in graphene using circular light polarization. We demonstrate that time- and angle-resolved photoemission spectroscopy can track these states by measuring sizeable energy gaps and quasi-Floquet energy bands that form on femtosecond time scales. By analyzing Floquet energy level crossings and snapshots of pseudospin textures near the Dirac points, we identify transitions to new states with optically induced nontrivial changes of sublattice mixing that can lead to Berry curvature corrections of electrical transport and magnetization.

preprint2014arXiv

Ubiquity of Linear Resistivity at Intermediate Temperature in Bad Metals

Bad metals display transport behavior that differs from what is commonly seen in ordinary metals. One of the most significant differences is a resistivity that is linear in temperature and rises to well above the Ioffe-Regel limit (where the mean-free path is equal to the lattice spacing). Using an exact Kubo formula, we show that a linear resistivity naturally occurs for many systems when they are in an incoherent intermediate-temperature state. We verify the analytic arguments with numerical calculations for a simplified version of the Hubbard model which is solved with dynamical mean-field theory. Similar features have also been seen in Hubbard models, where they can begin at even lower temperatures due to the formation of resilient quasiparticles.

preprint2013arXiv

Beyond Planck-Einstein quanta: amplitude driven quantum excitation

Linear-response quantum excitation is proportional to the amplitude of the field, with the energy of the excitation given by the driving frequency. As the amplitude is increased, there is a crossover, where the excitation energy is governed by the amplitude of the driving field, not its frequency. As the amplitude is increased even further, then complex quantum oscillations develop, whose origin is related to Wannier-Stark physics, but has no simple explanation. We illustrate this phenomena with the exact solution of the simplest model of a charge-density-wave insulator driven by a uniform time-dependent electric field.

preprint2013arXiv

Effect of defects on the phonons and the effective spin-spin interactions of an ultracold Penning trap quantum simulator

We generalize the analysis of the normal modes for a rotating ionic Coulomb crystal in a Penning trap to allow for inhomogeneities in the system. Our formal developments are completely general, but we choose to examine a crystal of Be+ ions with BeH+ defects to compare with current experimental efforts. We examine the classical phonon modes (both transverse and planar) and we determine the effective spin-spin interactions when the system is driven by an axial spin-dependent optical dipole force. We examine situations with up to approximately 15% defects. We find that most properties are not strongly influenced by the defects, indicating that the presence of a small number of defects will not significantly affect experiments.

preprint2013arXiv

Exact solution for Bloch oscillations of a simple charge-density-wave insulator

Charge-density-wave systems have a static modulation of the electronic charge at low temperatures when they enter an ordered state. While they have been studied for decades in equilibrium, it is only recently that they have been examined in nonequilibrium with time-resolved studies. Here, we present the exact solution for the nonequilibrium response of electrons (in the simplest model for a charge density wave) when the system is placed under a strong DC electric field. This allows us to examine the formation of driven Bloch oscillations and how the presence of a current modifies the nonequilibrium density of states.

preprint2013arXiv

Exact solution for high harmonic generation and the response to an ac driving field for a charge density wave insulator

We develop and exactly solve a model for electrons driven by pulsed or continuous ac fields. The theory includes both the photoexcitation process as well as the subsequent acceleration of the electrons. In the case of an ac response, we examine both the nonequilibrium density of states and the current. In the case of pulsed light for high harmonic generation, we find the radiated light assumes a nearly universal behavior, with only limited dependence on the parameters of the system, except for the amplitude of the driving field, which determines the range of high harmonics generated and a tendency toward a narrowing of the peaks in a charge density wave versus a metal. This type of high harmonic generation can potentially be used for the creation of solid-state-based ultrafast light sources.

preprint2013arXiv

Mapping of the unoccupied states and relevant bosonic modes via the time dependent momentum distribution

The unoccupied states of complex materials are difficult to measure, yet play a key role in determining their properties. We propose a technique that can measure the unoccupied states, called time-resolved Compton scattering, which measures the time-dependent momentum distribution (TDMD). Using a non-equilibrium Keldysh formalism, we study the TDMD for electrons coupled to a lattice in a pump-probe setup. We find a direct relation between temporal oscillations in the TDMD and the dispersion of the underlying unoccupied states, suggesting that both can be measured by time-resolved Compton scattering. We demonstrate the experimental feasibility by applying the method to a model of MgB$_2$ with realistic material parameters.

preprint2013arXiv

Nonequilibrium "melting" of a charge density wave insulator via an ultrafast laser pulse

We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental features in TaS2 and TbTe3 indicating a quasiuniversality for nonequilibrium "melting" with qualitative features that can be easily understood within a simple picture.

preprint2013arXiv

Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse $10^d\times 10^d$ matrix, with $d > 6$. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in non-equilibrium setups without mass transport.

preprint2013arXiv

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

We develop a strong-coupling ($t \ll U$) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to proximity effect. Our calculations for the scaled density in the vacuum to superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.

preprint2013arXiv

Theoretical description of high-order harmonic generation in solids

We consider several aspects of high-order harmonic generation in solids: the effects of elastic and inelastic scattering; varying pulse characteristics; and inclusion of material-specific parameters through a realistic band structure. We reproduce many observed characteristics of high harmonic generation experiments in solids including the formation of only odd harmonics in inversion-symmetric materials, and the nonlinear formation of high harmonics with increasing field. We find that the harmonic spectra are fairly robust against elastic and inelastic scattering. Furthermore, we find that the pulse characteristics play an important role in determining the harmonic spectra.

preprint2013arXiv

Thermalization of field driven quantum systems

There is much interest in how quantum systems thermalize after a sudden change, because unitary evolution should preclude thermalization. The eigenstate thermalization hypothesis resolves this because all observables for quantum states in a small energy window have essentially the same value; it is violated for integrable systems due to the infinite number of conserved quantities. Here, we show that when a system is driven by a DC electric field there are five generic behaviors: (i) monotonic or (ii) oscillatory approach to an infinite-temperature steady state; (iii) monotonic or (iv) oscillatory approach to a nonthermal steady state; or (v) evolution to an oscillatory state. Examining the Hubbard model (which thermalizes under a quench) and the Falicov-Kimball model (which does not), we find both exhibit scenarios (i-iv), while only Hubbard shows scenario (v). This shows richer behavior than in interaction quenches and integrability in the absence of a field plays no role.

preprint2013arXiv

Universal thermopower of bad metals

"Bad metals" have a large linear resistivity at high-T that is universally seen in oxides close to the Mott-Hubbard insulating phase. They also have an universal thermopower alpha(T): (i) at very low doping (lightly doped) alpha(T) has a pronounced low-T peak that shifts to higher-T with doping; (ii) at moderate doping (underdoped) alpha(T) has a small low-T peak that shifts to lower-T with doping and has a high-T sign change; and (iii) at the highest doping (overdoped) alpha(T) is negative and depends monotonically on T. Here we show that the simplified Hubbard model provides an easy to understand description of this phenomena due to the universal form for the chemical potential versus T for doped Mott insulators and the applicability of the Kelvin formula for the thermopower.

preprint2013arXiv

Using an artificial electric field to create the analog of the red spot of Jupiter in light-heavy Fermi-Fermi mixtures of ultracold atoms

Time-of-flight images are a common tool in ultracold atomic experiments, employed to determine the quasimomentum distribution of the interacting particles. If one introduces a constant artificial electric field, then the quasimomentum distribution evolves in time as Bloch oscillations are generated in the system and then damped showing a complex series of patterns. Surprisingly, in different mass Fermi-Fermi mixtures, these patterns can survive for long times, and resemble the stability of the red spot of Jupiter in classical nonlinear hydrodynamics. In this work, we illustrate the rich phenomena that can be seen in these driven quantum systems.

preprint2012arXiv

Correlation tuned cross-over between thermal and nonthermal states following ultrafast transient pumping

We examine electron-electron mediated relaxation following excitation of a correlated system by an ultrafast electric field pump pulse. The results reveal a dichotomy in the temporal evolution as one tunes through a Mott metal-to-insulator transition: in the metallic regime relaxation can be characterized by evolution toward a steady-state electronic distribution well described by Fermi-Dirac statistics with an increased effective temperature; however, in the insulating regime this quasithermal paradigm breaks down with relaxation toward a nonthermal state with a more complicated electronic distribution that does not vary monotonically as a function of energy. We characterize the behavior by studying changes in the energy, photoemission response, and electronic distribution as functions of time. Qualitatively these results should be observable on short enough time scales that the electrons behave like an isolated system not in contact with additional degrees of freedom which can act as a thermal bath. Importantly, proper modeling used to analyze experimental findings should account for this behavior, especially when using strong driving fields or studying materials whose physics may manifest the effects of strong correlations.

preprint2012arXiv

Emergence and Frustration of Magnetic Order with Variable-Range Interactions in a Trapped Ion Quantum Simulator

Frustration, or the competition between interacting components of a network, is often responsible for the complexity of many body systems, from social and neural networks to protein folding and magnetism. In quantum magnetic systems, frustration arises naturally from competing spin-spin interactions given by the geometry of the spin lattice or by the presence of long-range antiferromagnetic couplings. Frustrated magnetism is a hallmark of poorly understood systems such as quantum spin liquids, spin glasses and spin ices, whose ground states are massively degenerate and can carry high degrees of quantum entanglement. The controlled study of frustrated magnetism in materials is hampered by short dynamical time scales and the presence of impurities, while numerical modeling is generally intractable when dealing with dynamics beyond N~30 particles. Alternatively, a quantum simulator can be exploited to directly engineer prescribed frustrated interactions between controlled quantum systems, and several small-scale experiments have moved in this direction. In this article, we perform a quantum simulation of a long-range antiferromagnetic quantum Ising model with a transverse field, on a crystal of up to N = 16 trapped Yb+ atoms. We directly control the amount of frustration by continuously tuning the range of interaction and directly measure spin correlation functions and their dynamics through spatially-resolved spin detection. We find a pronounced dependence of the magnetic order on the amount of frustration, and extract signatures of quantum coherence in the resulting phases.

preprint2012arXiv

Exact solution of a variety of X-ray probes in the Falicov-Kimball model with dynamical mean-field theory

We examine core-level X-ray photoemission spectroscopy (XPS), X-ray absorption near-edge spectroscopy (XANES) and X-ray emission spectroscopy (XES) in the Falicov-Kimball model by using the exact solution from dynamical mean-field theory. XPS measures the core-hole propagator, XANES measures the absorption of X-rays when the core electron is excited to an unoccupied electronic state of the solid, and not emitted, and XES measures the spectra of light emitted as electrons fill the core-hole state created via some form of X-ray excitation. These three spectra are closely related to one another and display orthogonality catastrophe behavior at T=0. We illustrate how to evaluate these spectra efficiently at finite temperature, with a primary focus on the details of XANES.

preprint2012arXiv

Intrinsic phonon effects on analog quantum simulators with ultracold trapped ions

Linear Paul traps have been used recently to simulate the transverse field Ising model with long-range spin-spin couplings. We study the intrinsic effects of phonon creation (from the initial phonon ground state) on the spin-state probability and spin entanglement for such quantum spin simulators. While it has often been assumed that phonon effects are benign because they play no role in the pure Ising model, they can play a significant role when a transverse field is added to the model. We use a many-body factorization of the quantum time-evolution operator of the system, adiabatic perturbation theory and exact numerical integration of the Schrödinger equation in a truncated spin-phonon Hilbert space followed by a tracing out of the phonon degrees of freedom to study this problem. We find that moderate phonon creation often makes the probabilities of different spin states behave differently from the static spin Hamiltonian. In circumstances in which phonon creation is minor, the spin dynamics state probabilities converge to the static spin Hamiltonian prediction at the cost of reducing the spin entanglement. We show how phonon creation can severely impede the observation of kink transitions in frustrated spin systems when the number of ions increases. Many of our results also have implications for quantum simulation in a Penning trap.

preprint2012arXiv

Many-body effects on the capacitance of multilayers made from strongly correlated materials

Recent work by Kopp and Mannhart on novel electronic systems formed at oxide interfaces has shown interesting effects on the capacitances of these devices. We employ inhomogeneous dynamical mean-field theory to calculate the capacitance of multilayered nanostructures. These multilayered nanostructures are composed of semi-infinite metallic leads coupled via a strongly correlated dielectric barrier region. The barrier region can be adjusted from a metallic regime to a Mott insulator through adjusting the interaction strength. We examine the effects of varying the barrier width, temperature, potential difference, screening length, and chemical potential. We find that the interaction strength has a relatively strong effect on the capacitance, while the potential and temperature show weaker dependence.

preprint2012arXiv

Phonon mediated quantum spin simulator employing a planar ionic crystal in a Penning trap

We derive the normal modes for a rotating Coulomb ion crystal in a Penning trap, quantize the motional degrees of freedom, and illustrate how they can by driven by a spin-dependent optical dipole force to create a quantum spin simulator on a triangular lattice with hundreds of spins. The analysis for the axial modes (oscillations perpendicular to the two-dimensional crystal plane) follow a standard normal-mode analysis, while the remaining planar modes are more complicated to analyze because they have velocity-dependent forces in the rotating frame. After quantizing the normal modes into phonons, we illustrate some of the different spin-spin interactions that can be generated by entangling the motional degrees of freedom with the spin degrees of freedom via a spin-dependent optical dipole force. In addition to the well-known power-law dependence of the spin-spin interactions when driving the axial modes blue of phonon band, we notice certain parameter regimes in which the level of frustration between the spins can be engineered by driving the axial or planar phonon modes at different energies. These systems may allow for the analog simulation of quantum spin glasses with large numbers of spins.

preprint2012arXiv

Pulsed high harmonic generation of light due to pumped Bloch oscillations in noninteracting metals

We derive a simple theory for high-order harmonic generation due to pumping a noninteracting metal with a large amplitude oscillating electric field. The model assumes that the radiated light field arises from the acceleration of electrons due to the time-varying current generated by the pump, and also assumes that the system has a constant density of photoexcited carriers, hence it ignores the dipole excitation between bands (which would create carriers in semiconductors). We examine the circumstances under which odd harmonic frequencies would be expected to dominate the spectrum of radiated light, and we also apply the model to real materials like ZnO, for which high-order harmonic generation has already been demonstrated in experiments.

preprint2012arXiv

Quasi-universal transient behavior of a nonequilibrium Mott insulator driven by an electric field

We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field, and that for a broad range of field strengths, a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments, but shows that an interesting new quasi-universal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasi-universal behavior of the current and spectral function for different values of the hopping.

preprint2012arXiv

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these moments for equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential, but require knowledge of equal time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moments to moments of spectral functions calculated from a variety of different approximations and use them to benchmark their accuracy.

preprint2012arXiv

Strongly enhanced thermal transport in a lightly doped Mott insulator at low temperature

We show how a lightly doped Mott insulator has hugely enhanced electronic thermal transport at low temperature. It displays universal behavior independent of the interaction strength when the carriers can be treated as nondegenerate fermions and a nonuniversal "crossover" region where the Lorenz number grows to large values, while still maintaining a large thermoelectric figure-of-merit. The electron dynamics are described by the Falicov-Kimball model which is solved for arbitrary large on-site correlation with a dynamical mean-field theory algorithm on a Bethe lattice. We show how these results are generic for lightly doped Mott insulators as long as the renormalized Fermi liquid scale is pushed to very low temperature and the system is not magnetically ordered.

preprint2011arXiv

Density wave patterns for fermionic dipolar molecules on a square optical lattice: Mean-field-theory analysis

We model a system of ultracold fermionic dipolar molecules on a two-dimensional square lattice. Assuming that the molecules are in their nondegenerate hyperfine ground state, and that the dipole moment is polarized perpendicular to the plane (as in the recent experiments on 40K-87Rb molecules), we approximate these molecules as spinless fermions with long range repulsive dipolar interactions. We use mean field theory to obtain the restricted phase diagram as a function of the filling, the strength of interaction and the temperature. We find a number of ordered density wave phases in the system, as well as phase separation between these phases. A Monte Carlo analysis shows that the higher-period phases are usually suppressed in the exact solution.

preprint2011arXiv

Efficiency for preforming molecules from mixtures of light Fermi and heavy Bose atoms in optical lattices: the strong-coupling-expansion method

We discuss the application of a strong-coupling expansion (perturbation theory in the hopping) for studying light-Fermi-heavy-Bose (like $^{40}$K-$^{87}$Rb) mixtures in optical lattices. We use the strong-coupling method to evaluate the efficiency for pre-forming molecules, the entropy per particle and the thermal fluctuations. We show that within the strong interaction regime (and at high temperature), the strong-coupling expansion is an economical way to study this problem. In some cases, it remains valid even down to low temperatures. Because the computational effort is minimal, the strong-coupling approach allows us to work with much larger system sizes, where boundary effects can be eliminated, which is particularly important at higher temperatures. Since the strong-coupling approach is so efficient and accurate, it allows one to rapidly scan through parameter space in order to optimize the pre-forming of molecules on a lattice (by choosing the lattice depth and interspecies attraction). Based on the strong-coupling calculations, we test the thermometry scheme based on the fluctuation-dissipation theorem and find the scheme gives accurate temperature estimation even at very low temperature. We believe this approach and the calculation results will be useful in the design of the next generation of experiments, and will hopefully lead to the ability to form dipolar matter in the quantum degenerate regime.

preprint2011arXiv

Momentum distribution and ordering in mixtures of ultracold light and heavy fermionic atoms

The momentum distribution is one of the most important quantities which provides information about interactions in many-body systems. At the same time it is a quantity that can easily be accessed in experiments on ultracold atoms. In this paper, we consider mixtures of light and heavy fermionic atoms in an optical lattice described effectively by the Falicov-Kimball model. Using a Monte Carlo method, we study how different ordered density-wave phases can be detected by measurement of the momentum distribution of the light atoms. We also demonstrate that ordered phases can be seen in Bragg scattering experiments. Our results indicate that the main factor that determines the momentum distribution of the light atoms is the trap confinement. On the other hand, the pattern formed by the heavy atoms seen in the Bragg scattering experiments is very sensitive to the temperature and possibly can be used in low-temperature thermometry.

preprint2011arXiv

Onset of a Quantum Phase Transition with a Trapped Ion Quantum Simulator

A quantum simulator is a well controlled quantum system that can simulate the behavior of another quantum system which may require exponentially large classical computing resources to understand otherwise. In the 1980s, Feynman proposed the use of quantum logic gates on a standard controllable quantum system to efficiently simulate the behavior of a model Hamiltonian. Recent experiments using trapped ions and neutral atoms have realized quantum simulation of Ising model in presence of external magnetic fields, and showed almost arbitrary control in generating non-trivial Ising coupling patterns. Here we use laser-cooled trapped 171-Yb+ ions to simulate the emergence of magnetism in a system of interacting spins by implementing a fully-connected non-uniform ferromagnetic Ising model in a transverse magnetic field. To link this quantum simulation to condensed matter physics, we measure scalable correlation functions and order parameters appropriate for the description of larger systems, such as various moments of the magnetization. By increasing the Ising coupling strengths compared with the external field, the crossover from paramagnetism to ferromagnetic order sharpens as the system is scaled up from N = 2 to 9 trapped ion spins. This points toward the onset of a quantum phase transition that should become infinitely sharp as the system approaches the macroscopic scale. We compare the measured ground state order to theory, which may become intractable for non-uniform Ising couplings as the number of spins grows beyond 20- 30 and even NP complete for a fully-connected frustrated Ising model, making this experiment an important benchmark for large-scale quantum simulation.

preprint2011arXiv

Resonant inelastic X-ray scattering in a Mott insulator

We calculate the resonant inelastic X-ray scattering (RIXS) response in a Mott insulator which is described by the Falicov-Kimball model. The model can be solved exactly within the single site dynamical mean-field theory (DMFT) approximation and the calculated RIXS response is accurate up to a local background correction. We find that on resonance the RIXS response is greatly enhanced over various other non-resonant background effects and the response systematically evolves from a single peak structure, arising due to relaxation processes within the lower Hubbard band, to a two peak structure, arising due to relaxation processes within the upper Hubbard band as well as across the Mott gap into the lower Hubbard band, as we vary the incident photon frequency to access states from the bottom of the lower Hubbard band to the top of the upper Hubbard band. The charge transfer excitations are found to disperse monotonically outwards (as a function of transfered energy) as we go from the center of the Brillouin zone towards the zone corner. These correlation induced features have been observed by Hasan {\it et. al.} (Science {\bf 288}, 1811 (2000)) and many other experimentalists in RIXS measurements over various transition metal oxide compounds and are found to be robust and survive even for large Auger lifetime broadening effects. As a comparison, we also calculate the dynamic structure factor for this model, which is proportional to the nonresonant part of the response, and does not show these specific signatures.

preprint2010arXiv

Effect of vertex corrections on the longitudinal transport through multilayered nanostructures: Exact dynamical mean-field theory approach applied to the inhomogeneous Falicov-Kimball model

Inhomogeneous dynamical mean-field theory is employed to calculate the vertex-corrected electronic charge transport for multilayered devices composed of semi-infinite metallic lead layers coupled through a strongly correlated material barrier region. The barrier region can be tuned from a metal to a Mott insulator through adjusting the interaction strength and the particle filling. We use the Falicov-Kimball model to describe the barrier region because an exact expression for the vertex corrections is known, allowing us to determine their effect on transport. The dc conductivity is calculated and we find the effects of the vertex corrections are relatively small, manifesting themselves in a small reduction in the resistance-area product. This reduction saturates in absolute magnitude as the barrier layer becomes thick, as expected due to the vanishing nature of the vertex corrections in bulk. The vertex corrections have a larger relative effect on the resistance-area product for more metallic and thinner devices.

preprint2010arXiv

Quantum Simulation and Phase Diagram of the Transverse Field Ising Model with Three Atomic Spins

We perform a quantum simulation of the Ising model with a transverse field using a collection of three trapped atomic ion spins. By adiabatically manipulating the Hamiltonian, we directly probe the ground state for a wide range of fields and form of the Ising couplings, leading to a phase diagram of magnetic order in this microscopic system. The technique is scalable to much larger numbers of trapped ion spins, where phase transitions approaching the thermodynamic limit can be studied in cases where theory becomes intractable.

preprint2010arXiv

Resonant Raman scattering effects in a nesting driven charge-density-wave insulator: exact analysis of the spinless Falicov-Kimball model with dynamical mean-field theory

We calculate the total electronic Raman scattering spectrum for a system with a charge density wave on an infinite-dimensional hypercubic lattice. The problem is solved exactly for the spinless Falicov-Kimball model with dynamical mean-field theory. We include the nonresonant, mixed, and resonant contributions in three common experimental polarizations, and analyze the response functions for representative values of the energy of the incident photons. The complicated scattering response can be understood from the significant temperature dependence of the many-body density of states, and includes a huge enhancement for photon frequencies near the charge-density-wave gap energy.

preprint2010arXiv

Temporal response of nonequilibrium correlated electrons

In this work we examine the time-resolved, instantaneous current response for the spinless Falicov-Kimball model at half-filling, on both sides of the Mott-Hubbard metal-insulator transition, driven by a strong electric field pump pulse. The results are obtained using an exact, nonequilibrium, many-body impurity solution specifically designed to treat the out-of-equilibrium evolution of electrons in time-dependent fields. We provide a brief introduction to the method and its computational details. We find that the current develops Bloch oscillations, similar to the case of DC driving fields, with an additional amplitude modulation, characterized by beats and induced by correlation effects. Correlations primarily manifest themselves through an overall reduction in magnitude and shift in the onset time of the current response with increasing interaction strength.

preprint2010arXiv

Time-resolved photoemission of correlated electrons driven out of equilibrium

We describe the temporal evolution of the time-resolved photoemission response of the spinless Falicov-Kimball model driven out of equilibrium by strong applied fields. The model is one of the few possessing a metal-insulator transition and admitting an exact solution in the time domain. The nonequilibrium dynamics, evaluated using an extension of dynamical mean-field theory, show how the driven system differs from two common viewpoints - a quasiequilibrium system at an elevated effective temperature (the "hot" electron model) or a rapid interaction quench ("melting" of the Mott gap) - due to the rearrangement of electronic states and redistribution of spectral weight. The results demonstrate the inherent trade-off between energy and time resolution accompanying the finite width probe pulses, characteristic of those employed in pump-probe time-domain experiments, which can be used to focus attention on different aspects of the dynamics near the transition.

preprint2009arXiv

Dynamical mean-field theory for light fermion--heavy boson mixtures on optical lattices

We theoretically analyze Fermi-Bose mixtures consisting of light fermions and heavy bosons that are loaded into optical lattices (ignoring the trapping potential). To describe such mixtures, we consider the Fermi-Bose version of the Falicov-Kimball model on a periodic lattice. This model can be exactly mapped onto the spinless Fermi-Fermi Falicov-Kimball model at zero temperature for all parameter space as long as the mixture is thermodynamically stable. We employ dynamical mean-field theory to investigate the evolution of the Fermi-Bose Falicov-Kimball model at higher temperatures. We calculate spectral moment sum rules for the retarded Green's function and self-energy, and use them to benchmark the accuracy of our numerical calculations, as well as to reduce the computational cost by exactly including the tails of infinite summations or products. We show how the occupancy of the bosons, single-particle many-body density of states for the fermions, momentum distribution, and the average kinetic energy evolve with temperature. We end by briefly discussing how to experimentally realize the Fermi-Bose Falicov-Kimball model in ultracold atomic systems.

preprint2009arXiv

Improving the efficiency of ultracold dipolar molecule formation by first loading onto an optical lattice

Ultracold ground state dipolar 40K-87Rb molecules have recently been produced in a loose harmonic trap by employing a magnetic field sweep across a Feshbach resonance followed by stimulated Raman adiabatic passage [K.-K. Ni et al., Science 322, 231 (2008)]. The overall experimental efficiency for molecule formation was around 20%. We show that the efficiency can be increased to nearly 100% if one first loads the atomic gases into an optical lattice of the appropriate depth and tunes the interspecies attraction to have exactly one atom of each species at an occupied lattice site. Our proposed scheme provides a large enhancement to the dipolar molecule density even at relatively high temperatures, and avoids three-body recombination loss by preventing lattice sites from being triply occupied.

preprint2009arXiv

Impurity problems for steady-state nonequilibrium dynamical mean-field theory

The mapping of steady-state nonequilibrium dynamical mean-field theory from the lattice to the impurity is described in detail. Our focus is on the case with current flow under a constant dc electric field of arbitrary magnitude. In addition to formulating the problem via path integrals and functional derivatives, we also describe the distribution function dependence of the retarded and advanced Green's functions. Our formal developments are exact for the Falicov-Kimball model. We also show how these formal developments are modified for more complicated models (like the Hubbard model).

preprint2009arXiv

Inhomogeneous spectral moment sum rules for the retarded Green function and self-energy of strongly correlated electrons or ultracold fermionic atoms in optical lattices

Spectral moment sum rules are presented for the inhomogeneous many-body problem described by the fermionic Falicov-Kimball or Hubbard models. These local sum rules allow for arbitrary hoppings, site energies, and interactions. They can be employed to quantify the accuracy of numerical solutions to the inhomogeneous many-body problem like strongly correlated multilayered devices, ultracold atoms in an optical lattice with a trap potential, strongly correlated systems that are disordered, or systems with nontrivial spatial ordering like a charge density wave or a spin density wave. We also show how the spectral moment sum rules determine the asymptotic behavior of the Green function, self-energy, and dynamical mean field, when applied to the dynamical mean-field theory solution of the many body problem. In particular, we illustrate in detail how one can dramatically reduce the number of Matsubara frequencies needed to solve the Falicov-Kimball model, while still retaining high precision, and we sketch how one can incorporate these results into Hirsch-Fye quantum Monte Carlo solvers for the Hubbard (or more complicated) models. Since the solution of inhomogeneous problems is significantly more time consuming than periodic systems, efficient use of these sum rules can provide a dramatic speed up in the computational time required to solve the many-body problem. We also discuss how these sum rules behave in nonequilibrium situations as well, where the Hamiltonian has explicit time dependence due to a driving field or due to the time-dependent change of a parameter like the interaction strength or the origin of the trap potential.

preprint2009arXiv

Momentum distribution of the insulating phases of the extended Bose-Hubbard model

We develop two methods to calculate the momentum distribution of the insulating (Mott and charge-density-wave) phases of the extended Bose-Hubbard model with on-site and nearest-neighbor boson-boson repulsions on $d$-dimensional hypercubic lattices. First we construct the random phase approximation result, which corresponds to the exact solution for the infinite-dimensional limit. Then we perform a power-series expansion in the hopping $t$ via strong-coupling perturbation theory, to evaluate the momentum distribution in two and three dimensions; we also use the strong-coupling theory to verify the random phase approximation solution in infinite dimensions. Finally, we briefly discuss possible implications of our results in the context of ultracold dipolar Bose gases with dipole-dipole interactions loaded into optical lattices.

preprint2009arXiv

Strong-coupling perturbation theory for the extended Bose-Hubbard model

We develop a strong-coupling perturbation theory for the extended Bose-Hubbard model with on-site and nearest-neighbor boson-boson repulsions on ($d > 1$)-dimensional hypercubic lattices. Analytical expressions for the ground-state phase boundaries between the incompressible (Mott or charge-density-wave insulators) and the compressible (superfluid or supersolid) phases are derived up to third order in the hopping $t$. We also briefly discuss possible implications of our results in the context of ultracold dipolar Bose gases with dipole-dipole interactions loaded into optical lattices.

preprint2008arXiv

F-electron spectral function of the Falicov-Kimball model and the Wiener-Hopf sum equation approach

We derive an alternative representation for the $f$-electron spectral function of the Falicov-Kimball model from the original one proposed by Brandt and Urbanek. In the new representation, all calculations are restricted to the real time axis, allowing us to go to arbitrarily low temperatures. The general formula for the retarded Green's function involves two determinants of continuous matrix operators that have the Toeplitz form. By employing the Wiener-Hopf sum equation approach and Szegö's theorem, we can derive exact analytic formulas for the large-time limits of the Green's function; we illustrate this for cases when the logarithm of characteristic function (which defines the continuous Toeplitz matrix) does and does not wind around the origin. We show how accurate these asymptotic formulas are to the exact solutions found from extrapolating matrix calculations to the zero discretization size limit.

preprint2008arXiv

Perfect dc conductance of a finite width Mott insulator sandwiched between metallic leads at zero temperature: a quantum emergent phenomenon in strongly correlated multilayers

Using inhomogeneous dynamical mean-field theory, we argue that the normal-metal proximity effect forces any finite number of "barrier" planes that are described by the (paramagnetic) Hubbard model and sandwiched between semi-infinite metallic leads to always be a Fermi liquid at T=0. This then implies that the inhomogeneous system restores lattice periodicity at zero frequency, has a well-defined Fermi surface, and should display perfect (ballistic) conductivity or "transparency". These results are, however, fragile with respect to finite frequency, V, T, disorder, or magnetism, all of which restore the expected quantum tunneling regime through a finite-width Mott insulator. Our formal results are complemented by numerical renormalization group studies on small thickness barriers that illustrate under what circumstances this behavior might be seen in real experimental systems.

preprint2006arXiv

Nonequilibrium dynamical mean-field theory

The many-body formalism for dynamical mean-field theory is extended to treat nonequilibrium problems. We illustrate how the formalism works by examining the transient decay of the oscillating current that is driven by a large electric field turned on at time t=0. We show how the Bloch oscillations are quenched by the electron-electron interactions, and how their character changes dramatically for a Mott insulator.

preprint2001arXiv

Reduction of Josephson critical current in short ballistic SNS weak links

We present fully self-consistent calculations of the thermodynamic properties of three-dimensional clean SNS Josephson junctions, where S is an s-wave short-coherence-length superconductor and N is a clean normal metal. The junction is modeled on an infinite cubic lattice such that the transverse width of the S is the same as that of the N, and its thickness is tuned from the short to long limit. Both the reduced order parameter near the SN boundary and the short coherence length depress the critical Josephson current $I_c$, even in short junctions. This is contrasted with recent measurements on SNS junctions finding much smaller $I_cR_N$ products than expected from the standard (non-self consistent and quasiclassical) predictions. We also find unusual current-phase relations, a ``phase anti-dipole'' spatial distribution of the self-consistently determined contribution to the macroscopic phase, and an ``unexpected'' minigap in the local density of states within the N region.

preprint1998arXiv

Reevaluating electron-phonon coupling strengths: Indium as a test case for ab initio and many-body-theory methods

Using indium as a test case, we investigate the accuracy of the electron-phonon coupling calculated with state-of-the-art ab initio and many-body theory methods. The ab initio calculations -- where electrons are treated in the local-density approximation, and phonons and the electron-phonon interaction are treated within linear response -- predict an electron-phonon spectral function alpha^2 F(omega) which translates into a relative tunneling conductance that agrees with experiment to within one part in 1000. The many-body theory calculations -- where alpha^2 F(omega) is extracted from tunneling data by means of the McMillan-Rowell tunneling inversion method -- provide spectral functions that depend strongly on details of the inversion process. For the the most important moment of alpha^2 F(omega), the mass-renormalization parameter lambda, we report 0.9 +/- 0.1, in contrast to the value 0.805 quoted for nearly three decades in the literature. The ab initio calculations also provide the transport electron-phonon spectral function alpha_{tr}^2 F(omega), from which we calculate the resistivity as a function of temperature in good agreement with experiment.

preprint1996arXiv

Vertex-corrected tunneling inversion in superconductors: Pb

The McMillan-Rowell tunneling inversion program, which extracts the electron-phonon spectral function $α^2F(Ω)$ and the Coulomb pseudopotential $μ^*$ from experimental tunneling data, is generalized to include the lowest-order vertex correction. We neglect the momentum dependence of the electron-phonon matrix elements, which is equivalent to using a local approximation. The perturbation theory is performed on the imaginary axis and then an exact analytic continuation is employed to produce the density of states on the real axis. Comparison is made with the experimental data for Pb.

preprint1995arXiv

Anomalous Normal-State Properties of High-T$_c$ Superconductors -- Intrinsic Properties of Strongly Correlated Electron Systems?

A systematic study of optical and transport properties of the Hubbard model, based on Metzner and Vollhardt's dynamical mean-field approximation, is reviewed. This model shows interesting anomalous properties that are, in our opinion, ubiquitous to single-band strongly correlated systems (for all spatial dimensions greater than one), and also compare qualitatively with many anomalous transport features of the high-T$_c$ cuprates. This anomalous behavior of the normal-state properties is traced to a ``collective single-band Kondo effect,'' in which a quasiparticle resonance forms at the Fermi level as the temperature is lowered, ultimately yielding a strongly renormalized Fermi liquid at zero temperature.

preprint1995arXiv

Competition between electron-phonon attraction and weak Coulomb repulsion

The Holstein-Hubbard model is examined in the limit of infinite dimensions. Conventional folklore states that charge-density-wave (CDW) order is more strongly affected by Coulomb repulsion than superconducting order because of the pseudopotential effect. We find that both incommensurate CDW and superconducting phases are stabilized by the Coulomb repulsion, but, surprisingly, the commensurate CDW transition temperature is more robust than the superconducting transition temperature. This puzzling feature is resolved by a detailed analysis of perturbation theory.

J. K. Freericks

What is connected

Connect this record

See the researcher in context

Building this map preview

83 published item(s)

Hidden local adiabatic ramp in the modulated time evolution and the quantum approximate optimization algorithm

Out-of-time ordered correlation functions for the localized $f$ electrons in the Falicov-Kimball model

Fixed Depth Hamiltonian Simulation via Cartan Decomposition

Decomposition of high-rank factorized unitary coupled-cluster operators using ancilla and multi-qubit controlled low-rank counterparts

Flexibility of the factorized form of the unitary coupled cluster ansatz

Native multiqubit Toffoli gates on ion trap quantum computers

Algebraic derivation of Kramers-Pasternack relations based on the Schrodinger factorization method

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

Many Body Thermodynamics on Quantum Computers via Partition Function Zeros

Quantum computation of magnon spectra

Test of the unitary coupled-cluster variational quantum eigensolver for a simple strongly correlated condensed-matter system

Constant Matrix Element Approximation to Time-Resolved Angle-Resolved Photoemission Spectroscopy

Direct observation of Higgs mode oscillations in the pump-probe photoemission spectra of electron-phonon mediated superconductors

Energy dissipation in the time domain governed by bosons in a correlated material

Exact time evolution of the asymmetric Hubbard dimer

Generalized gradient expansion for inhomogeneous dynamical mean-field theory: Application to ultracold atoms in a harmonic trap

The role of average time dependence on the relaxation of excited electron populations in nonequilibrium many-body physics

Creating thermal distributions from diabatic excitations in ion-trap-based quantum simulation

Determining the ground-state probability of a quantum simulation with product-state measurements

Exact calculation of phonon effects on spin squeezing

Measuring nonequilibrium retarded spin-spin Green's functions in an ion-trap based quantum simulator

Nonequilibrium dynamical mean-field theory for the charge-density-wave phase of the Falicov-Kimball model

Nonequilibrium response of an electron mediated charge-density-wave-ordered material to a large dc electric field

Theoretical basis for quantum simulation with a planar ionic crystal in a Penning trap using a triangular rotating wall

Creation of two-dimensional coulomb crystals of ions in oblate Paul traps for quantum simulations

Diabatic ramping spectroscopy of many-body excited states for trapped-ion quantum simulators

Effect of dynamical spectral weight redistribution on effective interactions in time-resolved spectroscopy

Feshbach modulation spectroscopy

Gauge invariance in the theoretical description of time-resolved angle-resolved pump/probe photoemission spectroscopy

Intrinsic anharmonic effects on the phonon frequencies and effective spin-spin interactions in a quantum simulator made from trapped ions in a linear Paul trap

Long-lived nonequilibrium states in the Hubbard model with an electric field

Nonequilibrium sum rules for the Holstein model

Robust finite-temperature disordered Mott insulating phases in inhomogeneous Fermi-Fermi mixtures with density and mass imbalance

The Nonequilibrium Many-Body Problem as a paradigm for extreme data science

Theoretical Description of Coherent Doublon Creation via Lattice Modulation Spectroscopy

Theory of pump-probe photoemission in graphene: Ultrafast tuning of Floquet bands and local pseudospin textures

Ubiquity of Linear Resistivity at Intermediate Temperature in Bad Metals

Beyond Planck-Einstein quanta: amplitude driven quantum excitation

Effect of defects on the phonons and the effective spin-spin interactions of an ultracold Penning trap quantum simulator

Exact solution for Bloch oscillations of a simple charge-density-wave insulator

Exact solution for high harmonic generation and the response to an ac driving field for a charge density wave insulator

Mapping of the unoccupied states and relevant bosonic modes via the time dependent momentum distribution

Nonequilibrium "melting" of a charge density wave insulator via an ultrafast laser pulse

Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

Theoretical description of high-order harmonic generation in solids

Thermalization of field driven quantum systems

Universal thermopower of bad metals

Using an artificial electric field to create the analog of the red spot of Jupiter in light-heavy Fermi-Fermi mixtures of ultracold atoms

Correlation tuned cross-over between thermal and nonthermal states following ultrafast transient pumping

Emergence and Frustration of Magnetic Order with Variable-Range Interactions in a Trapped Ion Quantum Simulator

Exact solution of a variety of X-ray probes in the Falicov-Kimball model with dynamical mean-field theory

Intrinsic phonon effects on analog quantum simulators with ultracold trapped ions

Many-body effects on the capacitance of multilayers made from strongly correlated materials

Phonon mediated quantum spin simulator employing a planar ionic crystal in a Penning trap

Pulsed high harmonic generation of light due to pumped Bloch oscillations in noninteracting metals

Quasi-universal transient behavior of a nonequilibrium Mott insulator driven by an electric field

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

Strongly enhanced thermal transport in a lightly doped Mott insulator at low temperature

Density wave patterns for fermionic dipolar molecules on a square optical lattice: Mean-field-theory analysis

Efficiency for preforming molecules from mixtures of light Fermi and heavy Bose atoms in optical lattices: the strong-coupling-expansion method

Momentum distribution and ordering in mixtures of ultracold light and heavy fermionic atoms

Onset of a Quantum Phase Transition with a Trapped Ion Quantum Simulator

Resonant inelastic X-ray scattering in a Mott insulator

Effect of vertex corrections on the longitudinal transport through multilayered nanostructures: Exact dynamical mean-field theory approach applied to the inhomogeneous Falicov-Kimball model

Quantum Simulation and Phase Diagram of the Transverse Field Ising Model with Three Atomic Spins

Resonant Raman scattering effects in a nesting driven charge-density-wave insulator: exact analysis of the spinless Falicov-Kimball model with dynamical mean-field theory

Temporal response of nonequilibrium correlated electrons

Time-resolved photoemission of correlated electrons driven out of equilibrium

Dynamical mean-field theory for light fermion--heavy boson mixtures on optical lattices

Improving the efficiency of ultracold dipolar molecule formation by first loading onto an optical lattice

Impurity problems for steady-state nonequilibrium dynamical mean-field theory

Inhomogeneous spectral moment sum rules for the retarded Green function and self-energy of strongly correlated electrons or ultracold fermionic atoms in optical lattices

Momentum distribution of the insulating phases of the extended Bose-Hubbard model