<p>In this paper we study strong coupling asymptotic expansions of <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\mathcal {N}=2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi mathvariant="script">N</mi> <mo>=</mo> <mn>2</mn> </mrow> </math></EquationSource> </InlineEquation> four-dimensional <i>SU</i>(2) gauge theory partition functions in general <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(\Omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi mathvariant="normal">Ω</mi> </math></EquationSource> </InlineEquation>-background. This is done by refining the Painlevé/gauge theory correspondence in terms of quantum Painlevé equations, obtained from <InlineEquation ID="IEq3"> <EquationSource Format="TEX">\(\mathbb {C}^2/\mathbb {Z}_2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msup> <mrow> <mi mathvariant="double-struck">C</mi> </mrow> <mn>2</mn> </msup> <mo stretchy="false">/</mo> <msub> <mi mathvariant="double-struck">Z</mi> <mn>2</mn> </msub> </mrow> </math></EquationSource> </InlineEquation> blowup relations. We present a general ansatz and a systematic analysis of the expansions of the gauge theory partition functions by solving the above equations around the strong coupling singularities, including Argyres–Douglas points. We check our results via refined holomorphic anomaly equations and compare with irregular Virasoro conformal blocks.</p>