One-dimensional many-body entangled open quantum systems with tensor network methods

We present a collection of methods to simulate entangled dynamics of open quantum systems governed by the Lindblad equation with tensor network methods. Tensor network methods using matrix product states have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the matrix product state design is tailored for closed one-dimensional systems governed by the Schr\"odinger equation, the next step for many-body quantum dynamics is the simulation of open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators, and locally purified tensor networks. Selected examples guide possible applications of the methods and serve moreover as a benchmark between the techniques. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose-Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, push- ing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that matrix product density operators are favorable for the exciton problem in our study. All three methods access the same library, i.e., the software package Open Source Matrix Product States, allowing us to have a meaningful comparison between the approaches based on the selected examples. For example, tensor operations are accessed from the same subroutines and with the same optimization eliminating one possible bias in a comparison of such numerical methods.Comment: 24 pages, 8 figures. Small extension of time evolution section and moving quantum simulators to introduction in comparison to v

A two-step trigonometrically fitted semi-implicit hybrid method for solving special second order oscillatory differential equation

In this paper, we derived a semi-implicit hybrid method (SIHM) which is a two-step method to solve special second order ordinary differential equations (ODEs). The SIHM which is three-stage and fourth-order is then trigonometrically fitted and denoted by TF-SIHM3(4). The method is constructed using trigonometrically fitted properties instead of using phase-lag and amplification properties. Numerical integration show that TF-SIHM3(4) is more accurate in term of accuracy compared to the existing explicit and implicit methods of the same order

Numerical Simulations of Water Waves\u27 Modulational Instability Under the Action of Wind and Dissipation

Since the work of Benjamin & Feir (1967), water waves propagating in infinite depth are known to be unstable to modulational instability. The evolution of such wave trains is well described through fully nonlinear simulations, but also by means of simplified models, such as the nonlinear Schrödinger equation. Segur et al. (2005) and Wu et al. (2006) studied theoretically and numerically the evolution of this instability, and both concluded that a long term restabilization occurs in these conditions. More recently, Kharif et al. (2010) considered wind forcing and viscous dissipation within the framework of a forced and damped nonlinear Schrödinger equation, and discussed the range of parameters for which this behavior is still valid. This work aims to demonstrate how numerical simulations are useful to analyze their theoretical predictions. Since we are dealing with long term stability, results are especially complicated to obtain experimentally. Thus, numerical simulations of the fully nonlinear equations turn out to be a very useful tool to provide a validation for the model. Here, the evolution of the modulational instability is investigated within the framework of the two-dimensional fully non linear potential equations, modified to include wind forcing and viscous dissipation. The wind model corresponds to the Miles theory. The introduction of dissipation in the equations is briefly discussed. The marginal stability curve derived from the fully nonlinear numerical simulations coincides with the curve obtained by Kharif et al. (2010) from a linear stability analysis. Furthermore, the long term evolution of the wave trains can be obtained through the numerical simulations, and it is found that the presence of wind forcing promotes the occurrence of a permanent frequency-downshifting without invoking damping due to breaking wave phenomenon

Interaction instability of localization in quasiperiodic systems

Integrable models form pillars of theoretical physics because they allow for full analytical understanding. Despite being rare, many realistic systems can be described by models that are close to integrable. Therefore, an important question is how small perturbations influence the behavior of solvable models. This is particularly true for many-body interacting quantum systems where no general theorems about their stability are known. Here, we show that no such theorem can exist by providing an explicit example of a one-dimensional many-body system in a quasiperiodic potential whose transport properties discontinuously change from localization to diffusion upon switching on interaction. This demonstrates an inherent instability of a possible many-body localization in a quasiperiodic potential at small interactions. We also show how the transport properties can be strongly modified by engineering potential at only a few lattice sites.Comment: 10 pages; (v2: additional explanations, data, and references

エネルギー関数を持つ発展方程式に対する幾何学的数値計算法

学位の種別: 課程博士審査委員会委員 : (主査)東京大学教授 松尾 宇泰, 東京大学教授 中島 研吾, 東京大学准教授 鈴木 秀幸, 東京大学准教授 長尾 大道, 東京大学准教授 齋藤 宣一University of Tokyo(東京大学

Soliton generation and control in engineered materials

Optical solitons provide unique opportunities for the control of light‐bylight. Today, the field of soliton formation in natural materials is mature, as the main properties of the possible soliton states are well understood. In particular, optical solitons have been observed experimentally in a variety of materials and physical settings, including media with cubic, quadratic, photorefractive, saturable, nonlocal and thermal nonlinearities. New opportunities for soliton generation, stability and control may become accessible in complex engineered, artificial materials, whose properties can be modified at will by, e.g., modulations of the material parameters or the application gain and absorption landscapes. In this way one may construct different types of linear and nonlinear optical lattices by transverse shallow modulations of the linear refractive index and the nonlinearity coefficient or complex amplifying structures in dissipative nonlinear media. The exploration of the existence, stability and dynamical properties of conservative and dissipative solitons in settings with spatially inhomogeneous linear refractive index, nonlinearity, gain or absorption, is the subject of this PhD Thesis. We address stable conservative fundamental and multipole solitons in complex engineered materials with an inhomogeneous linear refractive index and nonlinearity. We show that stable two‐dimensional solitons may exist in nonlinear lattices with transversally alternating domains with cubic and saturable nonlinearities. We consider multicomponent solitons in engineered materials, where one field component feels the modulation of the refractive index or nonlinearity while the other component propagates as in a uniform nonlinear medium. We study whether the cross‐phase‐modulation between two components allows the stabilization of the whole soliton state. Media with defocusing nonlinearity growing rapidly from the center to the periphery is another example of a complex engineered material. We study such systems and, in contrast to the common belief, we have found that stable bright solitons do exist when defocusing nonlinearity grows towards the periphery rapidly enough. We consider different nonlinearity landscapes and analyze the types of soliton solution available in each case. Nonlinear materials with complex spatial distributions of gain and losses also provide important opportunities for the generation of stable one‐ and multidimensional fundamental, multipole, and vortex solitons. We study onedimensional solitons in focusing and defocusing nonlinear dissipative materials with single‐ and double‐well absorption landscapes. In two‐dimensional geometries, stable vortex solitons and complexes of vortices could be observed. We not only address stationary vortex structures, but also steadily rotating vortex solitons with azimuthally modulated intensity distributions in radially symmetric gain landscapes. Finally, we study the possibility of forming stable topological light bullets in focusing nonlinear media with inhomogeneous gain landscapes and uniform twophoton absorption