Novel physics arises when strongly correlated system is driven out of equilibrium by external fields. Dramatic changes in physical properties, such as conductivity, are empirically observed in strongly correlated materials under high electric field. In particular, electric-field driven metal-insulator transitions are well-known as the resistive switching effect in a variety of materials, such as VO
2, V
2O
3 and other transition metal oxides. To satisfactorily explain both the phenomenology and its underlying mechanism, it is required to model microscopically the out-of-equilibrium dissipative lattice system of interacting electrons. In this thesis, we developed a systematic method of modeling non-equilibrium steady states for dissipative lattice systems by means of Non-equilibrium Green's function and Dynamical Mean Field Theory. We firstly establish a "minimum model" to formulate the strong-field transport in non-interacting dissipative electron lattice. This model is exactly soluble and convenient for discussing energy dissipation and steady-state properties. The formalism is then combined with Dynamical Mean Field Theory to provide a systematic framework describing the nonequilibrium steady-state of correlated materials. We use the formalism to study the strong-field transport properties of correlated materials, Mott insulators and Dirac electrons in graphene. We concentrate on the microscopic description of resistive switching. Of particular interest is the filament formation during the dynamical phase transition, which has been interpreted as a result of the delicate interplay between dissipation and Mott physics. We will also examine
IV characteristics and particularly the current saturation of Dirac electrons in graphene. The arXiv version has been updated with minor modifications and corrections.
In this work, we develop a non-equilibrium steady-state non-crossing approximation (NESS-NCA) impurity solver applicable to general impurity problems. The choice of the NCA as the impurity solver enables both a more accurate description of correlation effects with larger Coulomb interaction and scalability to multi-orbital systems. Based on this development, we investigate strongly correlated non-equilibrium states of a dissipative lattice system under constant electric fields. Both the electronic Coulomb interaction and the electric field are treated non-perturbatively using dynamical mean-field theory in its non-equilibrium steady-state form (NESS-DMFT) with the NESS-NCA impurity solver. We validate our implementation using a half-filled single-band Hubbard model attached to a fictitious free Fermion reservoir, which prevents temperature divergence. As a result, we identify metallic and insulating phases as functions of the electric field and the Coulomb interaction along with a phase coexistence region amid the metal-to-insulator transition (MIT). We find that the MIT driven by the electric field is qualitatively similar to the equilibrium MIT as a function of temperature, differing from results in previous studies using the iterative perturbation theory (IPT) impurity solver. Finally, we highlight the importance of the morphology of a correlated system under the influence of an electric field.
We investigate the non-equilibrium dynamics and steady-state properties of a driven-dissipative Bose-Hubbard chain using a self-consistent Gutzwiller mean-field (GMF) approach. By employing a robust Picard iteration scheme, we solve the non-linear master equation for the non-equilibrium steady state (NESS) in the presence of strong Kerr nonlinearity. We identify two distinct dynamical regimes governed by the interplay between coherent drive, dissipation, and interaction: a regular quasilinear regime and a chaotic regime. Linear stability analysis reveals that the transition to the chaotic regime is triggered by parametric instabilities arising from the drive-induced coherence. Furthermore, we characterize the onset of quantum chaos by calculating the out-of-time-order correlator (OTOC). Our results show that in the strong coupling regime, the OTOC exhibits rapid exponential growth and saturation, providing a clear signature of information scrambling in this open quantum system. The proposed numerical framework offers an efficient pathway to explore many-body correlations in larger photonic lattices.