Non-Markovian 1/f noise consists a dominant source of decoherence in superconducting qubits, yet its slow nature poses a significant challenge for accurate simulation. Here we develop a hierarchical equations of motion (HEOM) framework that enables efficient and reliable modeling of qubit dynamics and gate operations under 1/f noise. By using the approach, it is first shown that perturbative quantum master equations may fail to reproduce the correct dephasing dynamics of a qubit coupled to slow baths. We then analyze dynamical decoupling sequences by including effects of finite pulse duration. It is found that different pulse sequences results in different behavior in error accumulation: all X-CPMG sequences exhibit linear scaling with parity effects, Y-CPMG follows quadratic growth, and alternating XY-type sequences can suppress the error accumulation significantly. Finally, we extend the framework to two-qubit cross-resonance (CR) gates, reconstructing the full Choi matrix and Pauli Transfer Matrix (PTM) to identify the incoherent error induced by 1/f noise. Together, these results establish HEOM as a robust methodology for simulating the environmental noise in superconducting circuits and provide new insights into error mechanisms in both single- and two-qubit gates.
Chen et al. (Tue,) studied this question.