We study weakly interacting Bose gases using the functional renormalization group with a hydrodynamic effective action. We use a scale-dependent parametrization of the boson fields that interpolates between a Cartesian representation at high momenta and an amplitude-phase one for low momenta. We apply this to Bose gases in two and three dimensions near the superfluid phase transition where they can be described by statistical O(2) models. We are able to give consistent physical descriptions of the infrared regime in both two and three dimensions. In particular, and in contrast to previous studies using the functional renormalization group, we find a stable superfluid phase at finite temperatures in two dimensions. We compare our results for the superfluid and boson densities with Monte-Carlo simulations, and we find they are in reasonable agreement. arXiv:1806.10373v2 [cond-mat.quant-gas] 21 Sep 2018 in the physical limit will depend on the choices made in setting up the flow equations. These include the form of the regulator, and ways of optimizing this choice have been developed [29,30].There are many versions of the FRG; in this work we use the one introduced by Wetterich [26, 27] for the average effective action, which is based on the Legendre transform of the logarithm of the partition function. The exact flow equation for this effective action has a one-loop structure. However, in practice, a truncation scheme needs to be used in order to solve for the RG flow. In this work we employ a derivative expansion; other choices can be found in the literature [31][32][33]. Truncations based on derivative expansions have been widely used in applications of the FRG and, despite the approximations introduced, they have been successful in describing phase transitions and critical exponents in a variety of systems [27,28].In particular, with fields in the Cartesian representation, the FRG has had some success in describing the weaklyinteracting Bose gas and related O(N ) models in three dimensions. In the critical regime, it yields results that are in good agreement with Monte-Carlo simulations. Away from that region, it has also been used to investigate bulk thermodynamic properties [18,27,31,[34][35][36][37][38]. As shown by Dupuis [18,39], one appealing aspect of this FRG approach is that, despite the potential IR divergences of loop diagrams as the physical limit is approached, it does respect the relevant Ward identity [17], giving a vanishing anomalous self-energy.This self-energy, and the resulting IR-divergent propagator for the longitudinal mode, arise from a momentumindependent interaction that vanishes in the physical limit. This treatment thus fails to capture the interaction between the amplitude fluctuations of the condensate that is present in Popov's effective theory. It also omits the leading interaction between the Goldstone modes, which is of second order in momenta. These limitations may be remedied by the inclusion of higher-order terms in the derivative expansion of the action. However, in current ...