Hot cracks are notorious defects in casting, welding, and additive manufacturing. Grain boundaries play a key role in hot cracking formation. Still, their effects are often neglected in hot cracking susceptibility predictions, and thus a quantitative picture of hot cracking susceptibility and liquid fracture is still lacking. In this work, we construct a multi-order parameter phase-field model for binary alloy solidification to simulate attractive, neutral, and repulsive grain boundaries in the thin-interface limit. Phase-field simulations of different combinations of grain boundary and bi-crystal growth types are used to obtain solidification path data, which is then used as the input into the Rappaz, Drezet and Gremaud (RDG) model (Rappaz et al., 1999). We study the so-called Λ-shape variation of hot cracking susceptibility as a function of solute concentration, showing that the magnitude of the Λ curve peak shifts to higher values and to lower concentrations as grain boundary energy increases. Meanwhile, the peak magnitude becomes higher and shifts to a higher value and higher concentrations when convergent grain growth is applied. Furthermore, in all cases examined, the corresponding pressure drops predicted for hot cracking exceed the theoretical liquid rupture stress, consistent with experimental observations. These results demonstrate that quantitative phase-field simulations, when coupled with the RDG model, are capable of relatively accurately predicting liquid rapture states relevant to hot cracking since the large disparity between the predicted pressure drop and critical stress of liquid film is eliminated through the consideration of grain boundary effects in the solidification path calculation.