Enhanced room-temperature ductility and high-temperature surface passivation are critical for the applications of body-centered cubic (BCC) refractory multicomponent alloys. Here, we devise a hierarchical workflow for the rapid initial screening of promising high-ductility single-phase compositions based on relatively simple criteria. Regression models built using data from density-functional theory calculations are used to calculate the intrinsic strength and ductility for the (1-10)[111] slip system in over 10^7 quaternary alloys, from a 13-element composition space (Ti-Zr-Hf-V-Nb-Ta-Mo-W-Re-Ru-Al-Cr-Si). Computational thermodynamics calculations are then employed to screen alloys with a single BCC phase at annealing temperatures (800°C) and surface oxidation capability at service temperatures (900°C). The corresponding manuscript has been accepted by MRS Communications (MRSC-D-22-00157R1) with the title "Automated hierarchical screening of refractory multicomponent alloys with high intrinsic ductility and surface passivation potency". Please cite our paper if you would like to use our data