We present a new spectral problem, a generalization of the D-Kaup-Newell spectral problem, associated with the Lie algebra sl(2,R). Zero curvature equations furnish the soliton hierarchy. The trace identity produces the Hamiltonian structure for the hierarchy. Lastly, a reduction of the spectral problem is shown to have a different soliton hierarchy with a bi-Hamiltonian structure. The first major motivation of this dissertation is to present spectral problems that generate two soliton hierarchies with infinitely many commuting conservation laws and high-order symmetries, i.e., they are Liouville integrable. We use the soliton hierarchies and a non-seimisimple matrix loop Lie algebra in order to construct integrable couplings. An enlarged spectral problem is presented starting from a generalization of the D-Kaup-Newell spectral problem. Then the enlarged zero curvature equations are solved from a series of Lax pairs producing the desired integrable couplings. A reduction is made of the original enlarged spectral problem generating a second integrable coupling system. Next, we discuss how to compute bilinear forms that are symmetric, ad-invariant, and non-degenerate on the given non-semisimple matrix Lie algebra to employ the variational identity. The variational identity is applied to the original integrable couplings of a generalized D-Kaup-Newell soliton hierarchy to furnish its Hamiltonian structures. Then we apply the variational identity to the reduced integrable couplings. The reduced coupling system has a bi-Hamiltonian structure. Both integrable coupling systems retain the properties of infinitely many commuting high-order symmetries and conserved densities of their original subsystems and, again, are Liouville integrable. In order to find solutions to a generalized D-Kaup-Newell integrable coupling system, a theory of Darboux transformations on integrable couplings is formulated. The theory pertains to a spectral problem where the spectral matrix is a polynomial in lambda of any order. An application to a generalized D-Kaup-Newell integrable couplings system is worked out, along with an explicit formula for the associated Bäcklund transformation. Precise one-soliton-like solutions are given for the m-th order generalized D-Kaup-Newell integrable coupling system.

An investigation into structures of bi-integrable and tri-integrable couplings is undertaken. Our study is based on semi-direct sums of matrix Lie algebras. By introducing new classes of matrix loop Lie algebras, we form new Lax pairs and generate several new bi-integrable and tri-integrable couplings of soliton hierarchies through zero curvature equations. Moreover, we discuss properties of the resulting bi-integrable couplings, including infinitely many commuting symmetries and conserved densities. Their Hamiltonian structures are furnished by applying the variational identities associated with the presented matrix loop Lie algebras. The goal of this dissertation is to demonstrate the efficiency of our approach and discover rich structures of bi-integrable and tri-integrable couplings by manipulating matrix Lie algebras.

In chapter 2, we study two Kaup-Newell-type matrix spectral problems, derive their soliton hierarchies within the zero curvature formulation, and furnish their bi-Hamiltonian structures by the trace identity to show that they are integrable in the Liouville sense. In chapter 5, we obtain the Riemann theta function representation of solutions for the first hierarchy of generalized Kaup-Newell systems. In chapter 3, using Hirota bilinear forms, we discuss positive quadratic polynomial solutions to generalized bilinear equations, which generate lump or lump-type solutions to nonlinear evolution equations, and propose an algorithm for computing higher-order lump or lump-type solutions. In chapter 4, we study mixed exponential and trigonometric wave solutions (called complexitons) to general bilinear equations, and propose two methods to find complexitons to generalized bilinear equations. We also succeed in proving that by choosing suitable complex coefficients in soliton solutions, multi-complexitons are actually real wave solutions from complex soliton solutions and establish the linear superposition principle for complexion solutions. In each chapter, we present computational examples.

We derive two hierarchies of 1+1 dimensional soliton-type integrable systems from two spectral problems associated with the Lie algebra of the special orthogonal Lie group SO(3,R). By using the trace identity, we formulate Hamiltonian structures for the resulting equations. Further, we show that each of these equations can be written in Hamiltonian form in two distinct ways, leading to the integrability of the equations in the sense of Liouville. We also present finite-dimensional Hamiltonian systems by means of symmetry constraints and discuss their integrability based on the existence of sufficiently many integrals of motion.