Quantum many-body physics has been studied for many decades. A lot of intriguing phenomena have been observed and theories have been developed. Yet many problems remain unsolved, which is largely due to a lack of general efficient classical computation method. One-dimensional (1D) systems have also drawn much attention over the past few decades. One reason is that many unique strongly correlated quantum phenomena only appear in low dimensions. Another reason is that many exact results can be obtained in 1D for cross benchmarking, such as Bethe Ansatz, Bosonization, Bose-Fermi mapping, etc. In addition, many numerical methods, such as Matrix-Product-State based and Monte Carlo methods, work most efficiently in 1D. In recent years, 1D quantum gases have been realized in many cold atom labs, providing experimental motivation for their studies.

One of the most mysterious assumptions about a system of identical quantum particles is that the wavefunction must be symmetric (for bosons) or anti-symmetric (for fermions). This means that no two fermions can occupy the same state, while multiple occupancy is allowed for bosons. Bosons and fermions are therefore generally very different. But in 1D this distinction could become rather subtle, provided that the bosonic multiple occupancy is suppressed, which can happen if strong repulsion exists between them. However, when considering spin degrees of freedom, things become more complicated. In this work, we develop a generalized Bose-Fermi mapping theory, under which, the charge degrees of freedom is mapped to a spinless fermions, while the spin degrees of freedom to an effective spin chain. This mapping works for arbitrary spin, arbitrary trapping potentials, arbitrary spin-dependent interaction potential, and for either spinor bosons or fermions. In the strong interaction limit, the wavefunction of the system is represented by the strong coupling ansatz wavefunction (SCAW), which takes the form of a direct product of spinless fermions wavefunction and spin wavefunction. Using this mapping technique, we study the dynamics and collective modes of the system. Furthermore, we develop a very efficient method to calculate the one-body density matrix, from which we can calculate the momentum distribution of the system.