We demonstrate that $$\|μ\|_{U^3([N])} \ll_{A}^{\text{ineff}} \log^{-A}(N)$$ $$\|Λ- Λ_Q\|_{U^3([N])} \ll_{A}^{\text{ineff}} \log^{-A}(N)$$ for any $A > 0$ where $Λ_Q$ is an approximant to the von Mangoldt function and will be defined below, improving upon a bound of Tao-Teräväinen (2021). As a consequence, among other things, we have the following: $$\mathbb{E}_{x, y \in [N], x + 3y \in [N]} Λ(x)Λ(x + y)Λ(x + 2y)Λ(x + 3y) = \mathfrak{S} + O_A(\log^{-A}(N))$$ where $\mathfrak{S}$ is the singular series for the configuration $(x, x + y, x + 2y, x + 3y)$. In fact, we show that $$\|μ- μ_{Siegel}\|_{U^3([N])} \ll \exp(-O(\log^{1/C}(N)))$$ $$\|Λ- Λ_{Siegel}\|_{U^3([N])} \ll \exp(-O(\log^{1/C}(N)))$$ where $μ_{Siegel}$ and $Λ_{Siegel}$ are approximants of $μ$, and $Λ$, respectively, representing the Siegel zero contribution of $μ$ and are defined in the above article. To do so, we use an improvement of the $U^3$ inverse theorem due to Sanders and we follow the approach of Green and Tao (2007), opting to use the ``old-fashioned" approach to equidistribution on two-step nilmanifolds which was also considered by Green and Tao (2017), and by Gowers and Wolf (2010). To the author's knowledge, this is the first time that quadratic Fourier analysis over $\mathbb{Z}/N\mathbb{Z}$ has achieved quasi-polynomial type bounds in applications.