Source author record

Guy Cohen

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

25 published item(s)

preprint2023arXiv

Stark-Many body localization in interacting infinite dimensional systems

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization (Stark-MBL) at least up to time which scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a non-monotonic dependence of the dynamical exponent on the electric field. As the field increases towards a critical value proportional to the Hubbard interaction strength, transport slows down, becoming sub-diffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting super-diffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

preprint2022arXiv

$L^2$-Quasi-compact and hyperbounded Markov operators

A Markov operator $P$ on a probability space $(S,Σ,μ)$, with $μ$ invariant, is called {\it hyperbounded} if for some $1 \le p<q \le \infty$ it maps (continuously) $L^p$ into $L^q$. We deduce from a recent result of Glück that a hyperbounded $P$ is quasi-compact, hence uniformly ergodic, in all $L^r(S,μ)$, $1<r< \infty$. We prove, using a method similar to Foguel's, that a hyperbounded Markov operator has periodic behavior similar to that of Harris recurrent operators, and for the ergodic case obtain conditions for aperiodicity. Given a probability $ν$ on the unit circle, we prove that if the convolution operator $P_νf:=ν*f$ is hyperbounded, then $ν$ is atomless. We show that there is $ν$ absolutely continuous such that $P_ν$ is not hyperbounded, and there is $ν$ with all powers singular such that $P_ν$ is hyperbounded. As an application, we prove that if $P_ν$ is hyperbounded, then for any sequence $(n_k)$ of distinct positive integers with bounded gaps, $(n_kx)$ is uniformly distributed mod 1 for $ν$ almost every $x$ (even when $ν$ is singular).

preprint2022arXiv

Interaction expansion inchworm Monte Carlo solver for lattice and impurity models

Multi-orbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. Here, we generalize the inchworm quantum Monte Carlo method to the interaction expansion and explore its application to typical single- and multi-orbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. So far, for the systems studied here, it remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. The problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

preprint2022arXiv

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

preprint2020arXiv

Correlated nonequilibrium steady states without energy flux

Floquet engineering of closed quantum systems can lead to the formation of long-lived prethermal states that, in general, eventually thermalize to infinite temperature. Coupling these driven systems to dissipative baths can stabilize such states, establishing a true nonequilibrium steady state. We demonstrate that in a certain strongly interacting lattice model coupled to a bath and driven by an electric field, such steady states can have the remarkable property that the cycle-averaged rate of energy transfer between the lattice and the baths vanishes. Despite this, we show that these states retain a clear nonequilibrium nature.

preprint2020arXiv

Green's function methods for single molecule junctions

We present a brief pedagogical review of theoretical Green's function methods applicable to open quantum systems out of equilibrium in general, and single molecule junctions in particular. We briefly describe experimental advances in molecular electronics, then discuss different theoretical approaches. We then focus on Green's function methods. Two characteristic energy scales governing the physics are many-body interactions within the junctions, and molecule-contact coupling. We therefore discuss weak interactions and weak coupling, as two limits that can be conveniently treated within, respectively, the standard nonequilibrium Green's function (NEGF) method and its many-body flavors (pseudoparticle and Hubbard NEGF). We argue that the intermediate regime, where the two energy scales are comparable, can in many cases be efficiently treated within the recently introduced superperturbation dual fermion approach. Finally, we review approaches for going beyond these analytically accessible limits, as embodied by recent developments in numerically exact methods based on Green's functions.

preprint2020arXiv

Joint and double coboundaries of commuting contractions

Let $T$ and $S$ be commuting contractions on a Banach space $X$. The elements of $(I-T)(I-S)X$ are called {\it double coboundaries}, and the elements of $(I-T)X \cap (I-S)X$ are called {\it joint cobundaries}. For $U$ and $V$ the unitary operators induced on $L_2$ by commuting invertible measure preserving transformations which generate an aperiodic $\mathbb Z^2$-action, we show that there are joint coboundaries in $L_2$ which are not double coboundaries. We prove that if $α$,$β\in (0,1)$ are irrational, with $T_α$ and $T_β$ induced on $L_1(\mathbb T)$ by the corresponding rotations, then there are joint coboundaries in $C(\mathbb T)$ which are not measurable double cobundaries (hence not double coboundaries in $L_1(\mathbb T)$).

preprint2020arXiv

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

We present a numerically exact Inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

preprint2016arXiv

Almost mixing of all orders and CLT for some $\mathbb{Z}^d$-actions on subgroups of $\mathbb{F}\_p^{\mathbb{Z}^d}$

For N d-actions by algebraic endomorphisms on compact abelian groups, the existence of non-mixing configurations is related to "S-unit type" equations and plays a role in limit theorems for such actions. We consider a family of endomorphisms on shift-invariant subgroups of F Z d p and show that there are few solutions of the corresponding equations. This implies the validity of the Central Limit Theorem for different methods of summation.

preprint2016arXiv

The Anderson--Holstein Model in Two Flavors of the Non--Crossing Approximation

The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for non-equilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary non-crossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.

preprint2015arXiv

Absence of diffusion in an interacting system of spinless fermions on a one-dimensional disordered lattice

We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems.

preprint2015arXiv

Taming the dynamical sign problem in real-time evolution of quantum many-body problems

Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new `Inchworm Algorithm', based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.

preprint2014arXiv

CLT for random walks of commuting endomorphisms on compact abelian groups

Let $\Cal S$ be an abelian group of automorphisms of a probability space $(X, {\Cal A}, μ)$ with a finite system of generators $(A_1, ..., A_d)$. Let $A^{\el}$ denote $A_1^{\ell_1} ... A_d^{\ell_d}$, for ${\el}= (\ell_1, ..., \ell_d)$. If $(Z_k)$ is a random walk on $\Z^d$, one can study the asymptotic distribution of the sums $\sum_{k=0}^{n-1} \, f \circ A^{\,{Z_k(ω)}}$ and $\sum_{\el \in \Z^d} \PP(Z_n= \el) \, A^\el f$, for a function $f$ on $X$. In particular, given a random walk on commuting matrices in $SL(ρ, \Z)$ or in ${\Cal M}^*(ρ, \Z)$ acting on the torus $\T^ρ$, $ρ\geq 1$, what is the asymptotic distribution of the associated ergodic sums along the random walk for a smooth function on $\T^ρ$ after normalization? In this paper, we prove a central limit theorem when $X$ is a compact abelian connected group $G$ endowed with its Haar measure (e.g. a torus or a connected extension of a torus), $\Cal S$ a totally ergodic $d$-dimensional group of commuting algebraic automorphisms of $G$ and $f$ a regular function on $G$. The proof is based on the cumulant method and on preliminary results on the spectral properties of the action of $\Cal S$, on random walks and on the variance of the associated ergodic sums.

preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo

We present two methods for computing two-time correlation functions or Green's functions from real time bold-line continuous time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real time methods to analytically continued spectral functions obtained from imaginary time Monte Carlo, thus probing the limits of analytic continuation.

preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo: spectral properties of the nonequilibrium Anderson impurity model

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

preprint2014arXiv

Numerical operator method for the real time dynamics of strongly-correlated quantum impurity systems far from equilibrium

We develop a method for studying the real time dynamics of Heisenberg operators in strongly-interacting nonequilibrium quantum impurity models. Our method is applicable to a wide range of interaction strengths and to bias voltages beyond the linear response regime, works at zero temperature, and overcomes the finite-size limitations faced by other numerical methods. We compare our method with quantum Monte Carlo simulations at a strong interaction strength, at which no analytical method is applicable up to now. We find a very good coincidence of the results at high bias voltage, and in the short time period at low bias voltage. We discuss the possible reason of the deviation in the long time period at low bias voltage. We also find a good coincidence of our results with the perturbation results at weak interactions.

preprint2013arXiv

Central limit theorem for commutative semigroups of toral endomorphisms

Let $\Cal S$ be an abelian finitely generated semigroup of endomorphisms of a probability space $(Ω, {\Cal A}, μ)$, with $(T_1, ..., T_d)$ a system of generators in ${\Cal S}$. Given an increasing sequence of domains $(D_n) \subset \N^d$, a question is the convergence in distribution of the normalized sequence $|D_n|^{-\frac12} \sum_{{\k} \, \in D_n} \, f \circ T^{\,{\k}}$, for $f \in L^2_0(μ)$, where $T^{\k}= T_1^{k_1} ... T_d^{k_d}$, ${\k}= (k_1, ..., k_d) \in {\N}^d$. After a preliminary spectral study when the action of $\Cal S$ has a Lebesgue spectrum, we consider $\N^d$- or $\Z^d$-actions given by commuting toral automorphisms or endomorphisms on $\T^ρ$, $ρ\geq 1$. For a totally ergodic action by automorphisms, we show a CLT for the above normalized sequence or other summation methods like barycenters, as well as a criterion of non-degeneracy of the variance, when $f$ is regular on the torus. A CLT is also proved for some semigroups of endomorphisms. Classical results on the existence and the construction of such actions by automorphisms are recalled.

preprint2013arXiv

Generalized projected dynamics for non-system observables of non-equilibrium quantum impurity models

The reduced dynamics formalism has recently emerged as a powerful tool to study the dynamics of non-equilibrium quantum impurity models in strongly correlated regimes. Examples include the non-equilibrium Anderson impurity model near the Kondo crossover temperature and the non-equilibrium Holstein model, for which the formalism provides an accurate description of the reduced density matrix of the system for a wide range of timescales. In this work, we generalize the formalism to allow for non-system observables such as the current between the impurity and leads. We show that the equation of motion for the reduced observable of interest can be closed with the equation of motion for the reduced density matrix and demonstrate the new formalism for a generic resonant level model.

preprint2013arXiv

Multiple steady-states in nonequilibrium quantum systems with electron-phonon interactions

The existence of more than one steady-state in a many-body quantum system driven out-of-equilibrium has been a matter of debate, both in the context of simple impurity models and in the case of inelastic tunneling channels. In this Letter, we combine a reduced density matrix formalism with the multilayer multiconfiguration time-dependent Hartree method to address this problem. This allows to obtain a converged numerical solution of the nonequilibrium dynamics. Considering a generic model for quantum transport through a quantum dot with electron-phonon interaction, we prove that a unique steady-state exists regardless of the initial electronic preparation of the quantum dot consistent with the converged numerical results. However, a bistability can be observed for different initial phonon preparations. The effects of the phonon frequency and strength of the electron-phonon couplings on the relaxation to steady-state and on the emergence of bistability is discussed.

preprint2012arXiv

Numerically Exact Long Time Magnetization Dynamics at the Nonequilibrium Kondo Crossover of the Anderson Impurity Model

We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with the memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular the crossover to the Kondo domain. We provide detailed predictions for novel nonequilibrium phenomena, including non-monotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields.

preprint2012arXiv

Simulating Lattice Spin Models on Graphics Processing Units

Lattice spin models are useful for studying critical phenomena and allow the extraction of equilibrium and dynamical properties. Simulations of such systems are usually based on Monte Carlo (MC) techniques, and the main difficulty is often the large computational effort needed when approaching critical points. In this work, it is shown how such simulations can be accelerated with the use of NVIDIA graphics processing units (GPUs) using the CUDA programming architecture. We have developed two different algorithms for lattice spin models, the first useful for equilibrium properties near a second-order phase transition point and the second for dynamical slowing down near a glass transition. The algorithms are based on parallel MC techniques, and speedups from 70- to 150-fold over conventional single-threaded computer codes are obtained using consumer-grade hardware.

preprint2011arXiv

Application of a semiclassical model for the second-quantized many-electron Hamiltonian to nonequilibrium quantum transport: The resonant level model

A semiclassical (SC) approach is developed for nonequilibrium quantum transport in molecular junctions. Following the early work of Miller and White [J. Chem. Phys. 84, 5059 (1986)], the many-electron Hamiltonian in second quantization is mapped onto a classical model that preserves the fermionic character of electrons. The resulting classical electronic Hamiltonian allows for real-time molecular dynamics simulations of the many-body problem from an uncorrelated initial state to the steady state. Comparisons with exact results generated for the resonant level model reveal that a semiclassical treatment of transport provides a quantitative description of the dynamics at all relevant timescales for a wide range of bias and gate potentials, and for different temperatures. The approach opens a door to treating nontrivial quantum transport problems that remain far from the reach of fully quantum methodologies.

preprint2011arXiv

Memory Effects In Nonequilibrium Quantum Impurity Models

Memory effects play a key role in the dynamics of strongly correlated systems driven out of equilibrium. In the present study, we explore the nature of memory in the nonequilibrium Anderson impurity model. The Nakajima--Zwanzig--Mori formalism is used to derive an exact generalized quantum master equation for the reduced density matrix of the interacting quantum dot, which includes a non-Markovian memory kernel. A real-time path integral formulation is developed, in which all diagrams are stochastically sampled in order to numerically evaluate the memory kernel. We explore the effects of temperature down to the Kondo regime, as well as the role of source--drain bias voltage and band width on the memory. Typically, the memory decays on timescales significantly shorter than the dynamics of the reduced density matrix itself, yet under certain conditions it develops a smaller long tail. In addition we address the conditions required for the existence, uniqueness and stability of a steady-state.

preprint2010arXiv

Dynamic consistency for Stochastic Optimal Control problems

For a sequence of dynamic optimization problems, we aim at discussing a notion of consistency over time. This notion can be informally introduced as follows. At the very first time step $t_0$, the decision maker formulates an optimization problem that yields optimal decision rules for all the forthcoming time step $t_0, t_1, ..., T$; at the next time step $t_1$, he is able to formulate a new optimization problem starting at time $t_1$ that yields a new sequence of optimal decision rules. This process can be continued until final time $T$ is reached. A family of optimization problems formulated in this way is said to be time consistent if the optimal strategies obtained when solving the original problem remain optimal for all subsequent problems. The notion of time consistency, well-known in the field of Economics, has been recently introduced in the context of risk measures, notably by Artzner et al. (2007) and studied in the Stochastic Programming framework by Shapiro (2009) and for Markov Decision Processes (MDP) by Ruszczynski (2009). We here link this notion with the concept of "state variable" in MDP, and show that a significant class of dynamic optimization problems are dynamically consistent, provided that an adequate state variable is chosen.

preprint2010arXiv

Price decomposition in large-scale stochastic optimal control

We are interested in optimally driving a dynamical system that can be influenced by exogenous noises. This is generally called a Stochastic Optimal Control (SOC) problem and the Dynamic Programming (DP) principle is the natural way of solving it. Unfortunately, DP faces the so-called curse of dimensionality: the complexity of solving DP equations grows exponentially with the dimension of the information variable that is sufficient to take optimal decisions (the state variable). For a large class of SOC problems, which includes important practical problems, we propose an original way of obtaining strategies to drive the system. The algorithm we introduce is based on Lagrangian relaxation, of which the application to decomposition is well-known in the deterministic framework. However, its application to such closed-loop problems is not straightforward and an additional statistical approximation concerning the dual process is needed. We give a convergence proof, that derives directly from classical results concerning duality in optimization, and enlghten the error made by our approximation. Numerical results are also provided, on a large-scale SOC problem. This idea extends the original DADP algorithm that was presented by Barty, Carpentier and Girardeau (2010).