  We employ a regularized relative trace formula to establish a second moment
estimate for twisted $L$-functions across all aspects over a number field. Our
results yield hybrid subconvex bounds for both Hecke $L$-functions and twisted
$L$-functions, comparable to the Weyl bound in suitable ranges. Moreover, we
present an application of our results to address the simultaneous nonvanishing
problem.
