This data set contains receiver functions calculated from three-component waveform data available in IRISDMC from 2000 to 2019. The waveform data are for earthquakes greater than magnitude 6.0, depth smaller than 100 km, and epicentral distances between 30 and 95 degrees. The raw waveform data are first converted to displacement and high pass filtered with a corner frequency of 0.02 Hz. Then the waveforms are windowed 20 s before and 400 s after the P arrival for analysis. The vertical-, north- and east-components are rotated to produce the L-, Q-and T-components. Subsequently, the [-20 s, 60 s] portion around the P arrival on the L-component is cut out as a parent waveform and deconvolved from the full L- and Q-components to produce the L- and Q-receiver functions. The data set is organized by year. The receiver functions in each year are compressed into a zip file named by the year when the waveforms were recorded. The receiver functions are in SAC format (https://seiscode.iris.washington.edu/projects/sac), which can be read using the Seismic Analysis Code (https://ds.iris.edu/ds/nodes/dmc/forms/sac/) or the Obspy package (https://github.com/obspy/obspy). The filename contains information of the origin time of the earthquake, the name of seismic network, the name of the seismic station, and component. For example, in the file name ‘TA.Z59A..BHL.M.2013.059.1405.SACD.01’, ‘TA’ stands for the network name, ‘Z59A’ is the station name, ‘BHL’ means the L component, ‘2013.059.1405’ means the earthquake occurred at 14:05 on the 59th day of the year 2013.