Abstract.  A rigorous full-vectorial modal solution approach based on the finite element method is used to find the propagation properties of terahertz (THz) waveguides, such as photonic crystal fibers, quantum cascaded lasers, plasmonic waveguides, power splitters, and narrow-band filters. Design approaches to reduce the modal loss due to the material and leakage loss in photonic crystal fibers and in metal-coated hollow-glass plasmonic waveguides have also been considered. The plasmonic confinement and gain threshold of quantum cascaded lasers used as THz sources and the chromatic dispersion in plasmonic waveguides are also presented.