No Arabic abstract
We present a product formula to approximate the exponential of a skew-Hermitian operator that is a sum of generators of a Lie algebra. The number of terms in the product depends on the structure factors. When the generators have large norm with respect to the dimension of the Lie algebra, or when the norm of the effective operator resulting from nested commutators is less than the product of the norms, the number of terms in the product is significantly less than that obtained from well-known results. We apply our results to construct product formulas useful for the quantum simulation of some continuous-variable and bosonic physical systems, including systems whose potential is not quadratic. For many of these systems, we show that the number of terms in the product can be sublinear or subpolynomial in the dimension of the relevant local Hilbert spaces, where such a dimension is usually determined by an energy scale of the problem. Our results emphasize the power of quantum computers for the simulation of various quantum systems.
The Trotter-Suzuki decomposition is one of the main approaches for realization of quantum simulations on digital quantum computers. Variance-based global sensitivity analysis (the Sobol method) is a wide used method which allows to decompose output variance of mathematical model into fractions allocated to different sources of uncertainty in inputs or sets of inputs of the model. Here we developed a method for application of the global sensitivity analysis to the optimization of Trotter-Suzuki decomposition. We show with a proof-of-concept example that this approach allows to reduce the number of exponentiations in the decomposition and provides a quantitative method for finding and truncation unimportant terms in the system Hamiltonian.
The Trotter-Suzuki decomposition is an important tool for the simulation and control of physical systems. We provide evidence for the stability of the Trotter-Suzuki decomposition. We model the error in the decomposition and determine sufficiency conditions that guarantee the stability of this decomposition under this model. We relate these sufficiency conditions to precision limitations of computing and control in both classical and quantum cases. Furthermore we show that bounded-error Trotter-Suzuki decomposition can be achieved by a suitable choice of machine precision.
The Trotter-Suzuki approximation leads to an efficient algorithm for solving the time-dependent Schrodinger equation. Using existing highly optimized CPU and GPU kernels, we developed a distributed version of the algorithm that runs efficiently on a cluster. Our implementation also improves single node performance, and is able to use multiple GPUs within a node. The scaling is close to linear using the CPU kernels, whereas the efficiency of GPU kernels improve with larger matrices. We also introduce a hybrid kernel that simultaneously uses multicore CPUs and GPUs in a distributed system. This kernel is shown to be efficient when the matrix size would not fit in the GPU memory. Larger quantum systems scale especially well with a high number nodes. The code is available under an open source license.
Given the Hamiltonian, the evaluation of unitary operators has been at the heart of many quantum algorithms. Motivated by existing deterministic and random methods, we present a hybrid approach, where Hamiltonians with large amplitude are evaluated at each time step, while the remaining terms are evaluated at random. The bound for the mean square error is obtained, together with a concentration bound. The mean square error consists of a variance term and a bias term, arising respectively from the random sampling of the Hamiltonian terms and the operator splitting error. Leveraging on the bias/variance trade-off, the error can be minimized by balancing the two. The concentration bound provides an estimate on the number of gates. The estimates are verified by using numerical experiments on classical computers.
Symmetry lies at the heart of todays theoretical study of particle physics. Our manuscript is a tutorial introducing foundational mathematics for understanding physical symmetries. We start from basic group theory and representation theory. We then introduce Lie Groups and Lie Algebra and their properties. We next discuss with detail two important Lie Groups in physics Special Unitary and Lorentz Group, with an emphasis on their applications to particle physics. Finally, we introduce field theory and its version of the Noether Theorem. We believe that the materials cover here will prepare undergraduates for future studies in mathematical physics.