Static and dynamic aspects of the fission process of 226 Th are analyzed in a self-consistent framework based on relativistic energy density functionals. Constrained relativistic mean-field (RMF) calculations in the collective space of axially symmetric quadrupole and octupole deformations, based on the energy density functional PC-PK1 and a δ-force pairing, are performed to determine the potential energy surface of the fissioning nucleus, the scission line, the single-nucleon wave functions, energies and occupation probabilities, as functions of deformation parameters. Induced fission dynamics is described using the time-dependent generator coordinate method in the Gaussian overlap approximation. A collective Schrödinger equation, determined entirely by the microscopic single-nucleon degrees of freedom, propagates adiabatically in time the initial wave packet built by boosting the ground-state solution of the collective Hamiltonian for 226 Th. The position of the scission line and the microscopic input for the collective Hamiltonian are analyzed as functions of the strength of the pairing interaction. The effect of static pairing correlations on the pre-neutron emission charge yields and total kinetic energy of fission fragments is examined in comparison with available data, and the distribution of fission fragments is analyzed for different values of the initial excitation energy. * Electronic address: zpliphy@swu.edu.cn A microscopic description of fission presents one of the most complex problems in lowenergy theoretical nuclear physics [1, 2]. For a comprehensive recent review and an exhaustive list of references, we refer the reader to Ref. [1]. The spontaneous or induced fission process in which a heavy nucleus splits into fragments is out of reach for ab initio methods and, therefore, modern microscopic approaches are based on the framework of nuclear energy density functionals (NEDFs). Nuclear density functional theory (DFT) and its timedependent (TD) generalization have enabled a self-consistent treatment of both static and dynamic aspects of fission [3-11]. The slow large-amplitude collective motion of the compound system that eventually leads to the formation of the final fragments can be described, in a first approximation, as an adiabatic process in which the intrinsic nucleonic degrees of freedom are decoupled from macroscopic collective degrees of freedom such as multipole moments (deformations) of the mass distribution and pairing fields [1].Numerous studies of spontaneous fission, based on NEDFs, have analyzed the effects of the choice of collective coordinates (shape degrees of freedom), approximations used to calculate the collective inertia, and coupling between shape and pairing degrees of freedom on fission half-lives [12][13][14][15][16][17][18]. A quantitative description of induced fission is, in this framework, conceptually and computationally more challenging and this process has been explored less systematically. In particular, several recent studies have used the time-dependent gener...