Ann Du

Crest De
email: ducrest@math.univ-angers.fr

Paolo Rossi
email: paolo.rossi@u-bourgogne.fr

QUANTUM D 4 DRINFELD-SOKOLOV HIERARCHY AND QUANTUM SINGULARITY THEORY

In this paper we compute explicitly the double ramification hierarchy and its quantization for the D 4 Dubrovin-Saito cohomological field theory obtained applying the Givental-Teleman reconstruction theorem to the D 4 Coxeter group Frobenius manifold, or equivalently the D 4 Fan-Jarvis-Ruan-Witten cohomological field theory (with respect to the non-maximal diagonal symmetry group J = Z/3Z). We then prove its equivalence to the corresponding Dubrovin-Zhang hierarchy, which was known to coincide with the D 4 Drinfeld-Sokolov hierarchy. Our techniques provide hence an explicit quantization of the D 4 Drinfeld-Sokolov hierarchy. Moreover, since the DR hierarchy is well defined for partial CohFTs too, our approach immediately computes the DR hierarchies associated to the invariant sectors of the D 4 CohFT with respect to folding of the Dynkin diagram, the B 3 and G 2 Drinfeld-Sokolov hierarchies.

Contents

Drinfeld-Sokolov D 4 hierarchy

In [DS84], Drinfeld and Sokolov described how to associate a hierarchy of integrable PDEs to any semi-simple Lie algebra through its associated affine Kac-Moody algebra. In the most general case, the hierarchy is constructed as an infinite sequence of matrix Lax equations (zerocurvature equations). In the same article, for simple Lie algebras of types A n , B n and C n , the authors also provided a scalar Lax pair along with a bi-Hamiltonian representation using pseudo-differential operators. For the D n case however, Drinfeld and Sokolov represented one part of the hierarchy with a scalar Lax pair and bi-Hamiltonian representation ([DS84], pp. 2019-2021), this part we call the positive flows of the hierarchy, while the remaining part we call the negative flows.

In [LWZ10], the authors produced a complete picture of the hierarchy of D n type (positive and negative flows) in terms of scalar Lax pairs and bi-Hamiltonian representation. To do so, they introduced a new kind of pseudo-differential operators called of the second type (while the traditional ones are called of the first type). Roughly speaking, the operators of the second type are allowed to contain not only infinitely many nonzero terms in negative degree, but also infinitely many nonzero terms in positive degrees, along with considerations of gradation to ensure the well-definedness of the product of operators. In that same article [LWZ10], the authors also described the tau-symmetric bi-Hamiltonian structure (in the sense of [DZ05]) of the Drinfeld-Sokolov hierarchy of D n type.

In what follows, we briefly review the original Drinfeld-Sokolov reduction for the D 4 case, for it entirely defines the hierarchy [DS84]. Then we follow the approach of [LWZ10] to define the scalar Lax pairs of the positive and negative flows of the hierarchy. Finally, we describe and compute the tau-symmetric bi-Hamiltonian structure of D 4 .

Introduction

In [START_REF] Dubrovin | Differential geometry of the space of orbits of a Coxeter group[END_REF], based on the results of K. Saito [START_REF] Saito | Primitive forms for a universal unfolding of a function with an isolated critical point[END_REF][START_REF] Saito | The higher residue pairings K (k) F for a family of hypersurface singular points[END_REF][START_REF] Saito | Period mapping associated to a primitive form[END_REF] and in particular on his theory of primitive forms, Dubrovin constructed a generically semisimple Frobenius manifold structure on the space of orbits of any finite, irreducible Coxeter group or, via the relation of these with simple hypersurface singularities, on the space of miniversal deformations of any polynomial W : C n → C with an isolated critical point at the origin.

Although in general these Frobenius manifolds are not semisimple at the origin, it was proved in [START_REF] Milanov | Analyticity of the total ancestor potential in singularity theory[END_REF] that, for singularities of type ADE, using Givental-Teleman theory [START_REF]Semisimple Frobenius structures at higher genus[END_REF][START_REF] Teleman | The structure of 2D semi-simple field theories[END_REF] at a nearby semisimple point and then shifting the result back to the origin, yields a well defined conformal cohomological field theory [START_REF] Kontsevich | Gromov-Witten classes, quantum cohomology, and enumerative geometry[END_REF] which we refer to as the Dubrovin-Saito CohFT of the singularity.

In [START_REF] Fan | The Witten equation and its virtual fundamental cycle[END_REF][START_REF] Fan | The Witten equation, mirror symmetry, and quantum singularity theory[END_REF] Fan, Jarvis and Ruan, inspired by ideas of Witten, introduced another construction associating a certain moduli space of decorated curves and, through a virtual fundamental class on it, a conformal CohFT to a certain class of quasi-homogeneous polynomials with an isolated critical point at the origin, together with the choice of admissible symmetry group.

In case the symmetry group is the full automorphism group of the polynomial G max , in a typical instance of mirror symmetry, there is an isomorphism between the FJRW CohFT and the Dubrovin-Saito CohFT of the transposed singularity [START_REF] Fan | The Witten equation, mirror symmetry, and quantum singularity theory[END_REF]. When the symmetry group is smaller the mirror partner to the FJRW theory is in general not known, but for the case of ADE singularities the only admissible groups are always maximal, with the exception of singularities of type D 2n for which however the isomorphism with the Dubrovin-Saito theory still holds.

The case of D 4 is in may ways the most subtle among the ADE singularities. For instance the general method of proof for the mirror symmetry result of [START_REF] Fan | The Witten equation, mirror symmetry, and quantum singularity theory[END_REF] did not work in the D 4 case because of the specific form of the CohFT and its phase space. Indeed the proof for the mirror theorem was completed in [START_REF] Fan | Witten's D 4 Integrable Hierarchies Conjecture[END_REF]. In essence the complication (specifically the appearance of the so called broad sectors of the phase space, which complicate the description of the CohFT) originates from the peculiar symmetry of the D 4 singularity. This is apparent from the corresponding Dynkin diagram, indeed the most symmetric of the ADE diagrams.

In this paper we focus on the example of the D 4 Dubrovin-Saito CohFT, or the FJRW CohFT of the singularity D 4 : W = x 3 +xy 2 with symmetry group J 1 where J 1 (x, y) = (e 2πi 1 3 x, e 2πi 1 3 y) and approach it from the point of view of the associated integrable hierarchies.

It is well known that, given a semisimple CohFT one can associate to it the Dubrovin-Zhang hierarchy [START_REF] Dubrovin | Normal forms of hierarchies of integrable PDEs, Frobenius manifolds and Gromov-Witten invariants, a new 2005 version of[END_REF], an integrable hierarchy of Hamiltonian PDEs which controls the intersection theory of the CohFT with psi classes on the moduli space of stable curves. It was proved in [START_REF] Givental | Simple singularities and integrable hierarchies[END_REF][START_REF] Frenkel | Soliton equations, vertex operators, and simple singularities[END_REF] that the Dubrovin-Zhang hierarchy of the ADE CohFTs are the Drinfeld-Sokolov hierarchies associated to the corresponding ADE semisimple algebras.

The double ramification hierarchy is newer construction, introduced in [Bur15, BR16b], associating a quantum integrable system to a CohFT (there is no need for semisimplicity in this case, and in fact the construction even works in the classical limit for partial CohFTs). This construction, inspired by ideas from symplectic field theory [EGH00, FR11, Ros10], uses intersection theory of the CohFT with the double ramification cycle, Hodge classes and psi classes, and its relation with the Dubrovin-Zhang hierarchy is the object of a conjecture (the strong DR/DZ equivalence conjecture of [START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF][START_REF] Buryak | DR/DZ equivalence conjecture and tautological relations[END_REF]) which predicts that, given a semisimple CohFT (so that the DZ side is well defined), the classical limit of the DR hierarchy and the DZ hierarchy coincide after a specific change of coordinates in the phase space of the system. This conjecture has been proved for several specific CohFTs [START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF][START_REF] Buryak | Integrable systems of double ramification type[END_REF][START_REF] Buryak | Extended r-spin theory in all genera and the discrete KdV hierarchy[END_REF] and, in particular, in [START_REF] Buryak | Towards a description of the double ramification hierarchy for Witten's r-spin class[END_REF][START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF], for the A N Dubrovin-Saito CohFTs (which are nothing but the (N + 1)-spin CohFTs) for N ≤ 5.

In this paper we show that the classical double ramification hierarchy of the D 4 CohFT coincides with the D 4 Drinfeld-Sokolov hierarchy (with no need of coordinate change in this case) and hence, via the results of [START_REF] Givental | Simple singularities and integrable hierarchies[END_REF], with the Duborinv-Zhang hierarchy. Moreover we compute the quantum double ramification hierarchy which hence provides an explicit (closed formulas) quantization of the D 4 Drinfeld-Sokolov hierarchy.

Finally we exploit the particularly symmetric nature of the D 4 CohFT, which possesses both a Z 2 and a Z 3 symmetry, to consider its invariant parts as two partial CohFTs, for which the classical DR hierarchy gives rise to the B 3 and G 2 Drinfeld-Sokolov hierarchies, and thanks to the results of [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF] we deduce that the partial CohFT potentials are tau functions of the corresponding DR hierarchy.

g = o(8) ⊗ C[λ, λ -1 ]
the loop algebra of o(8). We choose generators {e i , f i , h i | 0 ≤ i ≤ 4} of g as in [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF] (also as in [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF]). We define the principal gradation g = k∈Z g k with deg e i = -deg f i = 1 and deg h i = 0. We use the notation g >0 = k>0 g k and similarly for g <0 . We denote Λ = 4 i=0 e i ∈ g 1 the principal cyclic element of g. Then it is well known [Kos59, Kac78, DS84] that the principal Heisenberg subaglebra Ker(Λ) admits the following decomposition:

Ker(Λ) = k∈Z odd C • Λ k ⊕ k∈Z odd C • Γ k , Λ k ∈ g k , Γ k ∈ g 3k .
(See [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF] for the expression of Λ k , Γ k .) Let us denote by b (resp. n) the negative Borel (resp. negative nilpotent) subalgebra of o(8). The starting point of the Drinfeld-Sokolov hierarchy is a matrix-valued differential operator of the form

L = ∂ x + Λ + q(x), q ∈ C ∞ (R, b). (1.1)
For any function S ∈ C ∞ (R, n), the operator L = e ad S L also has the form of Equation (1.1); we say that L and L are gauge equivalent. In [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF], the authors proved the fundamental property that there exists a (non unique) function U ∈ C ∞ (R, g <0 ) such that the operator L 0 = e ad U L has the form

L 0 = ∂ x + Λ + H(x), H ∈ C ∞ (R, Ker(Λ) ∩ g <0 ).
(Note that L and L 0 are not gauge equivalent). The Drinfeld-Sokolov hierarchy of type D 4 is then defined as the two sequences of matrix Lax equations:

∂L ∂t k = -e ad U Λ k ≥0 , L , (1.2) ∂L ∂ tk = e ad U Γ k ≥0 , L , (1.3) for k ∈ Z odd +
, where ( • ) ≥0 denotes the projection onto the subspace g ≥0 . These equations eventually take the form of evolutionary PDEs on the coordinates of the matrix q. Changing the representative of L in its gauge equivalence class amounts to changing coordinates. They do not depend on the choice of the function U . The flows with respect to the variable t k (1.2) are those we call the positive flows, while those with respect to the variables tk (1.3) are the negative flows.

1.2. Salar Lax pairs. We will start by describing the scalar Lax pairs of the positive flows (1.2). We follow the formal algebraic approach to pseudo-differential operators as in [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF]. We start by setting indeterminates s 1 k , . . . , s 4 k , with k ≥ 0 and denote s α 0 = s α . Define a formal derivation

∂ x = 4 α=1 k≥0 s α k+1 ∂ ∂s α k , so that s α k+1 = ∂ x (s α k )
acting on the ring of generalized differential polynomials

A s = C[[s * ]][s * >0 ][[ε]],
where s * denotes s 1 , . . . , s 4 . In what follows we will mostly drop the index in A s , simply writing A when no confusion can arise. We define a gradation A = k≥0 A [k] by setting deg s α k = k and deg ε = -1. Changes of coordinates s α → s α will be described by transformations of the form

s α (s * ≥0 ; ε) ∈ A [0] , 0 = det ∂ s α | ε=0 ∂s β α,β=1,...,4
, called Miura transforms and should be seen as maps from A s to A s . We refer to [START_REF] Rossi | Integrability, quantization and moduli spaces of curves[END_REF] for more details on generalized differential polynomials. As usual, we define the ring of pseudodifferential operators that, here, we call of the first type (following the vocabulary of [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF]),

D -= m k=-∞ a k ∂ k x m ∈ Z, a k ∈ A [0] . (1.4)
The product of two pseudo-differential operators of the first type is defined by, for any a, b ∈ D - and any n, m ∈ Z,

a∂ n x • b∂ m x = k≥0 n k ab k ε k ∂ n+m-k x , n k = n(n -1) • • • (n -k + 1) k! , (1.5) where b k = ∂ k x (b)
, and extended by linearity. Notice that we added a factor ε k to be coherent with the introduction of generalized differential polynomials. For any operator X = k≤m a k ∂ k

x , we call the positive part (respectively the negative part) of X the operator

X + = k≥0 a k ∂ k x (respectively X -= k<0 a k ∂ k x ).
As already defined by Drinfeld and Sokolov [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF], the scalar Lax operator of the hierarchy of type D 4 has the form1

L = ∂ 6 x + ∂ -1 x 3 µ=1 s µ ∂ 2µ-1 x + ∂ 2µ-1 x s µ + ∂ -1 x ∂ -1 x , (1.6)
where ( )2 = s4 . As expected, L ∈ D -. The operator L satisfies the additional condition that L * ∂ x = ∂ x L, where the formal adjoint of an operator is defined by (a∂ k

x ) * = (-∂ x ) k a. In [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF], the authors give the following proposition.

Proposition 1.1. There exists a unique pseudo-differential operator P ∈ D -, called the 6-th root of L, of the form

P = ∂ x + k<0 p k ∂ k
x such that P 6 = L. The operator P satisfies [P, L] = 0 and

P * ∂ x + ∂ x P = 0. (1.7)
Moreover, Equation (1.7) is equivalent to the condition that for every k ∈ Z odd + , the free term of the operator (P k ) + vanishes, i.e. (P k ) + (1) = 0.

The above proposition implies that the following equations are well defined:

∂L ∂t k = P k + , L , k ∈ Z odd + (1.8)
These equation first appeared in [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF], they coincide with Equations (1.2). Eventually, they give evolutionary PDEs on the functions s 1 , . . . , s 4 .

Scalar Lax pairs of the negative flows. We now describe the scalar Lax pairs for the negative flows (1.3). As mentioned above, we need to introduce the pseudo-differential operators of the second type, following [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF]. First of all, we extend the additive group D -of pseudodifferential operators of the first type into the additive group

D = k∈Z a k ∂ k x a k ∈ A [0] .
We use on A [0] the gradation induced by deg

s α k = k and denote A [0] = k≥0 A [0] k . We say that an operator X ∈ D -⊂ D is homogeneous of degree k ∈ Z if X = ≤k a ∂ x , a ∈ A [0] k-.
We denote by D k the additive subgroup of D -of homogeneous operators of the first type of degree k. Then D = k∈Z D k . Similarly, for any d ∈ Z we define the additive subgroup D + (d) = k≥d D k . Now we define the additive subgroup

D + = d∈Z D + (d) .
The elements of D + are the so-called pseudo-differential operators of the second type. For these operators, we allow infinitely many nonzero terms of negative powers, as well as infinitely many nonzero terms of positive powers of ∂ x . In the end, the spaces D -and D + can be described as

D -= k≤m i≥0 a k,i ∂ k x a k,i ∈ A [0] i D + = k≥m i≥m-k a k,i ∂ k x a k,i ∈ A [0] i
Both are complete topological rings. Also, we define the subgroup

D b = D -∩ D + ,
whose elements are called bounded operators. Since for any X ∈ D k , Y ∈ D , their product satisfies X • Y ∈ D k+ , it follows that the product (1.5) in D -can be extended to D + to form a ring. The following proposition will be important for our computations [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF].

Proposition 1.2. There exists a unique operator Q ∈ D + , caleld the square root of L, of the form

Q = ∂ -1 x + k≥0 Q m • ∂ x (1.9) such that Q 2 = L. Here, Q m ∈ D b , is homogeneous of degree 2m and satisfies Q * m = Q m . Moreover, the operator Q satisfies Q * ∂ x + ∂ x Q = 0, Q * + ( ) = - m≥0 ∂ x Q m ( ) = - 1 2 ∂ x L + (1). (1.10)
Notice that in the above proposition, Q * + ( ) denotes the evaluation of the differential operator Q * + at the function . Thanks to this proposition, the following equations are well defined (we rewrite Equation (1.8) for reasons of clarity):

(1.11) ∂L ∂t k = P k + , L , ∂L ∂ tk = Q k + , L , k ∈ Z odd + .
Finally, we can state the following theorem (

L = ∂ x •L (denoted L in [LWZ10]
), which we write down using coordinates v α , namely,

L = ∂ 7 x + 3 µ=1 s µ ∂ 2µ-1 x + ∂ 2µ-1 x s µ + ∂ -1 x = ∂ 7 x + 3 µ=1 v µ ∂ 2µ-1 x + ṽµ ∂ 2µ-2 x + ∂ -1
x , (1.12) and we set v 4 = s 4 = ( ) 2 . The coordinates ṽα are related to the coordinates v α via the condition L * + L = 0 (see Equation (1.21)). In turns, the coordinates v α are related to the coordinates s α via a Miura transform by identifying the two expressions of L.

We call local functionals, the elements of the quotient space

Λ = A/ (Im(∂ x ) ⊕ C[[ε]]) .
Given a differential polynomial f ∈ A, we denote by f = f dx its class in Λ. For a local functional f ∈ Λ in formal variables v * * , we can define its variational derivatives by

δf δv α = k≥0 (-1) k ∂ k x ∂f ∂v α k ∈ A. It is well known that if f ∈ A is such that f (0) = 0 (no constant term), then f ∈ Im(∂ x )
if and only if δf δv α = 0 for any α ∈ {1, 2, 3, 4} (see e.g. [GKMZ70], Lemma 2). In particular, the variational derivatives are well defined on Λ. Next we define the variational differential (or variational derivative w.r.t. L) by

δf δ L = δf δv 4 + 1 2 3 µ=1 δf δv µ ∂ -2µ x + ∂ -2µ x δf δv µ ∈ D -.
In [START_REF] Drinfeld | Lie algebras and equations of Kortewegde Vries type[END_REF], Drinfeld and Sokolov gave the following Poisson brackets: Let two local functionals f , g ∈ Λ and their variational differentials

X = δf δ L , Y = δg δ L , then {f , g} 1 = ε -1 res X (∂ x Y + L) --( LY + ∂ x ) --(∂ x Y -L) + + ( LY -∂ x ) + dx, {f , g} 2 = ε -1 res X ( LY ) + L -L(Y L) + dx.
These brackets are compatible in the sense that for any λ, µ ∈ C, the map λ{•, •} 1 + µ{•, •} 2 still satisfies Jacobi's identity. In [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF], the authors proved the following theorem.

Theorem 1.4. The hierarchy (1.11) admits the following bi-Hamiltonian representation: for any local functional f ∈ Λ,

∂f ∂t k = {f , H k+6 } 1 = {f , H k } 2 , ∂f ∂ tk = {f , Ĥk+2 } 1 = {f , Ĥk } 2 , (1.13)
where the Hamiltonian functionals are given by

H k = 6 k res P k dx, Ĥk = 2 k res Q k dx
1.4. Tau structure. Finally, we present the tau structure (in the sense of [START_REF] Dubrovin | Normal forms of hierarchies of integrable PDEs, Frobenius manifolds and Gromov-Witten invariants, a new 2005 version of[END_REF]) of the Drinfeld-Sokolov hierarchy of type D 4 . The first step is to define the so-called topological variables: for µ ∈ {1, 2, 3} and p ≥ 0,

t µ p = 6Γ p + 1 + 2µ-1 6 Γ 2µ-1 6 t 6p+2µ-1 , t 4 p = 2Γ p + 3 2 Γ( 1 2 ) t2p+1 .
In particular, t 1 0 = t 1 = x. In the same fashion, we define the following Hamiltonian densities:

h µ,p-1 = Γ 2µ-1 6 6Γ p + 1 + 2µ-1 6 res P 6p+2µ-1 , h 4,p-1 = Γ( 1 2 ) 2Γ p + 3 2 res Q 2p+1 , (1.14) = 6 p p j=0 (2µ -1 + 6j) -1 res P 6p+2µ-1 , = 2 p p j=0 (1 + 2j) -1 res Q 2p+1 .
We denote by h α,p = h α,p dx the associated functionals. Then the Hamiltonian equations (1.13) read

∂f ∂t α p = {f , h α,p } 1 = p + 1 2 + µ α -1 {f , h α,p-1 } 2 ,
where µ α 's are the spectrum of the underlying Frobenius manifold [START_REF] Dubrovin | Geometry of 2D topological field theories[END_REF][START_REF] Dubrovin | Frobenius manifolds and central invariants for the Drinfeld-Sokolov bihamiltonian structures[END_REF]; they read µ ν = 2ν-4 6 , for ν ∈ {1, 2, 3}, and µ 4 = 0. These Hamiltonian densities satisfy the so-called tau symmetry [START_REF] Dubrovin | Normal forms of hierarchies of integrable PDEs, Frobenius manifolds and Gromov-Witten invariants, a new 2005 version of[END_REF], [START_REF] Liu | On the Drinfeld-Sokolov hierarchies of D type[END_REF]: for any α, β ∈ {1, 2, 3, 4} and p, q ≥ 1,

∂h α,p-1 ∂t β q = ∂h β,q-1 ∂t α p . (1.15)
If we define differential polynomials Ω α,p;β,q via ∂ x Ω α,p;β,q = ∂ t β q h α,p-1 , then we can find integration constants such that they satisfy Ω α,p;β,q = Ω β,q;α,p . This is a sufficient condition to define the tau function τ of the hierarchy by setting [DZ05]

∂ 2 log τ ∂t α p ∂t β q
= Ω α,p;β,q .

1.5. Explicit tau-symmetric Hamiltonian densities. In this section, we give the explicit fomulae for the Hamiltonian densities h α,-1 , h α,0 , for 1 ≤ α ≤ 4, and h 1,1 . As proven in [START_REF] Buryak | Towards a description of the double ramification hierarchy for Witten's r-spin class[END_REF], even h 1,1 (and the Poisson structure) suffices to confirm the equivalence with the double ramification hierarchy (see Section 2). Our computations will be expressed in a special set of coordinates called normal and defined by

u α = η αβ h β,-1
(summation over 1 ≤ β ≤ 4 is implicit) where the matrix (η αβ ) is given by (see e.g. [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF])

(η αβ ) =     0 0 6 0 0 6 0 0 6 0 0 0 0 0 0 2    
By Equation (1.14), the normal coordinates read, for µ ∈ {1, 2, 3},

u µ = 6 7 -2µ res P 7-2µ , u 4 = 2 res Q.
Because of the form of the operator Q (1.9), its residue is straightforwardly res Q = , which means that u 4 = 2 = 2 √ s 4 . For the remaining normal coordinates, we compute the residues in the coordinates s α and then invert the system. We find

             s 1 = 1 2 u 1 + 1 12 u 2 u 3 + 1 216 ( u 3 ) 3 + -1 8 ( u 3 1 ) 2 -1 6 u 2 2 -1 9 u 3 u 3 2 ε 2 + 23 90 u 3 4 ε 4 , s 2 = 1 2 u 2 + 1 8 ( u 3 ) 2 -1 2 u 3 2 ε 2 , s 3 = 1 2 u 3 , s 4 = 1 4 ( u 4 ) 2 .
(1.16)

(We have performed the substitution ∂ k x (f ) → ε/ √ 2 k ∂ k x (f ).
) This agrees with the expressions found in [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF], p. 751. Then the Hamiltonian densities h α,0 are given by, for µ ∈ {1, 2, 3},

h µ,0 = 6 (2µ -1)(2µ + 5) res P 2µ+5 , h 4,0 = 2 3 res Q 3
For the computation res Q 3 we write

res Q 3 = res QL = res ∂ -1 x L + res m≥0 Q m ∂ x L , = res ∂ -1 x L + res m≥0 Q m ∂ -1 x ,
where in the last equation we used the fact that

∂ x L = L = L+ + ∂ -1 x .
To compute the rightmost term, we write

Q m = k≤2m q m,k ∂ k x , then res m≥0 Q m ∂ -1 x = res m≥0 k≤2m q m,k ∂ k x ∂ -1 x = res m≥0 k≤2m q m,k k =0 k ε √ 2 ∂ k-1- x = m≥0 k≤2m q m,k ε √ 2 k k = m≥0 Q m ( ).
In the above equations, all matters of convergence are resolved by the grading of A and the fact that Q ∈ D + . Now thanks to Equation (1.10), it follows that

m≥0 Q m ( ) = 1 2 L + (1) = 3 µ=1 s µ 2µ-2 ε 2µ-2 2 µ-1 .
Using u 4 = 2 , we finally find that

h 4,0 = 1 3 res ∂ -1 x u 4 L + 1 3 u 4 3 µ=1 s µ 2µ-2 ε 2µ-2 2 µ-1 .
1.6. Explicit first Poisson structure. Finally, we compute the components of the first Poisson structure given by, for two local functionals f , g ∈ Λ and their variational differentials

X = δf δ L , Y = δg δ L , {f , g} 1 = 1 ε res X (∂ x Y + L) --( LY + ∂ x ) --(∂ x Y -L) + + ( LY -∂ x ) + dx. (1.18)
Working with normal coordinate u * * , we aim to express this bracket through an Hamiltonian operator

K αβ u = k≥0 K αβ k ∂ k x , where K αβ k ∈ A [-k+1] u , via {f , g} 1 = δf δ u α K αβ u δg δ u β .
Under a Miura transformation u * * = u * * (u * * , ε), the Hamiltonian operator transforms accoring to

K αβ u = (L * ) α µ • K µν u • L β ν , where (L * ) α µ = s≥0 ∂ u α ∂u µ s ∂ s x , L β ν = s≥0 (-∂ x ) s • ∂ u β ∂u ν s
. When clear from the context which formal variables are being used, as usual, we will drop the corresponding index in K αβ u or K αβ u and simply write K αβ . We use Dirac's notation of the components K αβ of the bracket and rewrite the above expression in the equivalent form

{f , g} 1 = δf (x) δ u α (x) { u α (x), u β (y)} 1 δg(y) δ u β (y) dxdy, (1.19) where { u α (x), u β (y)} 1 := k≥0 K αβ k (x)δ (k) (x -y). (1.20)
In Equation (1.20), the expression "{ u α (x), u β (y)} 1 " really is only a notation and does not mean that we evaluate the bracket on the pair of functions ( u α (x), u β (y)). However, one can evaluate the bracket on the pair of functions ( u α (z)δ(x -z), u β (z)δ(y -z)) (sometimes called the coordinate functionals in physics literature) and find that, almost tautologically,

{ u α (z)δ(x -z), u β δ(y -z)} 1 = δ u α (z)δ(x -z) δ u µ (z) k≥0 K µν k δ (k) (z -z ) δ u α (z )δ(x -z ) δ u µ (z ) dzdz = k≥0 K αβ k (x)δ (k) (x -y).
Hence,

{ u α (x), u β (y)} 1 = { u α (z)δ(x -z), u β δ(y -z)} 1 ,
where in the above equation, the left-hand side is purely notational, while the right-hand side can actually be computed. Now the Poisson bracket (1.18) is given in terms of the variational differentials, so that we first need to compute those of the coordinate functionals. To do so, we need to transform the normal coordinates u α into the coordinates v α of Equation (1.12). The coordinates ṽα are related to the v α 's by the condition L * + L = 0; they read

       ṽ1 = 1 2 √ 2 v 1 1 ε -1 8 √ 2 v 2 3 ε 3 + 1 8 √ 2 v 3 5 ε 5 , ṽ2 = 3 2 √ 2 v 2 1 ε -5 4 √ 2 v 3 3 ε 3 , ṽ3 = 5 2 √ 2 v 3 1 ε.
(1.21)

Next we will show that this D 4 DR hierarchy coincides with the D 4 Drinfeld-Sokolov hierarchy described in the previous section. Indeed, in [START_REF] Givental | Simple singularities and integrable hierarchies[END_REF][START_REF] Frenkel | Soliton equations, vertex operators, and simple singularities[END_REF] (see also [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF]) the authors proved that the D 4 Drinfeld-Sokolov hierarchy is the Dubrovin-Zhang hierarchy of the aforementioned D 4 CohFT and, as proven in [START_REF] Buryak | Towards a description of the double ramification hierarchy for Witten's r-spin class[END_REF], for a hierarchy of DR type to coincide with a DZ hierarchy with the same dispersionless (i.e. genus 0) limit, it is enough that their Hamiltonian operators and Hamiltonians h 1,1 coincide, which is what we will prove.

2.1. Double ramification hierarchy. We will not spell here the full definition of the classical double ramification hierarchy and its quantization, referring the reader to the papers [START_REF] Buryak | Double ramification cycles and integrable hierarchies[END_REF][START_REF] Buryak | Recursion relations for double ramification hierarchies[END_REF][START_REF] Buryak | Double ramification cycles and quantum integrable systems[END_REF][START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF], or the review [START_REF] Rossi | Integrability, quantization and moduli spaces of curves[END_REF], instead.

The construction has, as input, a cohomological field theory c g,n : V ⊗n → H * (M g,n , C) on the moduli space of stable curves M g,n , with metric η on the vector space V and unit e 1 ∈ V (see [START_REF] Kontsevich | Gromov-Witten classes, quantum cohomology, and enumerative geometry[END_REF] for the definition) and, from its intersection theory, produces an integrable system of tau-symmetric Hamiltonian PDEs with Hamiltonian densities g α,d ∈ A u , 1 ≤ α ≤ N , d ≥ -1 where only even powers of ε appear, and Hamiltonian operator

K αβ u = η αβ ∂ x .
The densities g α,d satisfy the property ∂g α,d+1 ∂u 1 = g α,d with, in particular, g α,-1 = η αµ u µ . As remarked in [START_REF] Buryak | Integrable systems of double ramification type[END_REF], an integrable system with this property and the above Poisson structure automatically possesses a tau-structure given by h α,d = δg α,d+1 δu 1 , with normal coordinates given by u α = η αµ δg µ,0 δu 1 , 1 ≤ α ≤ N .

In case the vector space V is graded, with deg e 1 = 0, the cohomological field theory is homogeneous if there exists δ ∈ C such that the maps c g,n : V ⊗n → H * (M g,n , C) and η :

V ⊗2 → C have degree deg c g,n = δ(g -1), deg η = -δ, δ ∈ C, 1 ≤ α ≤ dim V.
In this case the formal variables u α k and ε acquire a grading too and the Hamiltonian densities of the resulting DR hierarchy are homogeneous:

|u α k | = 1 -deg e α , |ε| = 1 -δ 2 , |g α,d | = d + 3 -δ -|u α |.
This grading | • | is not related to the differential grading deg(•) on A and Λ introduced in section 1.2.

The DR hierarchy coincides with the Dubrovin-Zhang hierarchy [DZ05, BPS12a, BPS12b] in genus 0 (i.e. its dispersionless limit ε → 0 is the principal hierarchy of the corresponding Frobenius manifold, see [START_REF] Dubrovin | Normal forms of hierarchies of integrable PDEs, Frobenius manifolds and Gromov-Witten invariants, a new 2005 version of[END_REF]) and the two hierarchies are conjectured to be equivalent up to a Miura transformation preserving the tau-structure. This is the strong DR/DZ equivalence conjecture of [START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF], which has been proven for several CohFTs but is open in general. In this paper we will prove the conjecture for the D 4 Dubrovin-Saito CohFT. N consider the operator

D h : A[[z]] → A[[z]] defined by D h = ∂ x • (D -1) -z{•, h}, D := k≥0 u α k ∂ ∂u α k + ε ∂ ∂ε . Suppose there exist N = dim V solutions g α (z) ∈ A [0] [[z]], α = 1, . . . , N , to D h g α (z) = 0 with initial conditions g α (z = 0) = η αµ u µ .
Then a new vector of solutions in the same class can be found by the following transformation (2.1)

g α (z) → a µ α (z)g µ (z) + b α (z), where a µ α (z) = δ µ α + i>0 a µ α,i z i ∈ C[[z]] and b α (z) = i>0 b α,i z i ∈ C[[z]].
The following result from [START_REF] Buryak | Integrable systems of double ramification type[END_REF] will constitute the main technical tool in our computation the DR hierarchy of the D 4 CohFT.

Theorem 2.1 ( [START_REF] Buryak | Integrable systems of double ramification type[END_REF]). Assume that h ∈ Λ [0] has the following properties:

(a) there exist N independent solutions g α (z) = p≥0 g α,p-1 z p ∈ A [0] [[z]], α = 1, . . . , N , to the equation (2.2) D h g α (z) = 0 with the initial conditions g α (z = 0) = η αµ u µ , (b) δh δu 1 = 1 2 η µν u µ u ν + ∂ 2 x r, r ∈ A [-2] .
Then, up to a transformation of type (2.1), we have

(i) g 1,0 = 1 2 η µν u µ u ν + ∂ 2 x (D -1) -1 r, (ii) g 1,1 = h, (iii) {g α,p , g β,q } = 0, α, β = 1, . . . , N, p, q ≥ -1, (iv) {g α,p , g β,0 } = ∂ x ∂g α,p+1 ∂u , β = 1, . . . , N, p ≥ -1, (v) ∂g α,p ∂u 1 = g α,p-1 , α = 1, . . . , N, p ≥ -1,
hence in particular h is part of an integrable tau-symmetric hierarchy.

We call a system of densities originating from an Hamiltonian h = g 1,1 as in the theorem above an integrable system of double ramification type. As proved in [START_REF] Buryak | Integrable systems of double ramification type[END_REF], the DR hierarchy of any given cohomological field theory, is always and integrable system of DR type. This fact implies in particular that the entire hierarchy of DR Hamiltonian densities can be reconstructed from g 1,1 ∈ Λ [0] alone by means of the DR recursion equation [START_REF] Buryak | Recursion relations for double ramification hierarchies[END_REF] (2.3)

∂ x (D -1)g α,d+1 = {g α,d , g 1,1 }, g α,-1 = η αµ u µ , 1 ≤ α ≤ N, d ≥ -1.
The above result can sometimes be used to effectively compute the Hamiltonian g 1,1 itself, starting from a limited amount of information on the CohFT, as we will show in the next section.

2.3. D 4 Dubrovin-Saito CohFT. In [START_REF] Dubrovin | Differential geometry of the space of orbits of a Coxeter group[END_REF], based on the work of K. Saito [Sai81, Sai83a, Sai83b], Dubrovin constructs a structure of Frobenius manifold [START_REF] Dubrovin | Geometry of 2D topological field theories[END_REF] on the space of orbits of any finite irreducible Coxeter group. This space is biholomorphic to the space of miniversal unfoldings of the corresponding simple hypersurface singularity and the Frobenius structure is given by the Milnor ring multiplicative structure and the residue pairing at the corresponding deformation. The resulting Frobenius manifold is generically semisimple and conformal, so Givental-Teleman [START_REF]Semisimple Frobenius structures at higher genus[END_REF][START_REF] Teleman | The structure of 2D semi-simple field theories[END_REF] (see also [START_REF] Pandharipande | Relations on M g,n via 3-spin structures[END_REF]) theory can be applied to produce a uniquely defined homogeneous cohomological field theory.

We are interested in this construction for the case of the D 4 simple singularity, W = x 3 +xy 2 . The resulting CohFT c D 4 g,n has phase space V = e 1 , . which, in particular, gives the metric

(2.5) η =     0 0 1 6 0 0 1 6 0 0 1 6 0 0 0 0 0 0 1 2    
Notice that this CohFT was proved to be isomorphic to the quantum singularity theory of Fan-Jarvis-Ruan-Witten [START_REF] Fan | The Witten equation and its virtual fundamental cycle[END_REF][START_REF] Fan | The Witten equation, mirror symmetry, and quantum singularity theory[END_REF] for the simple singularity W = x 3 + xy 2 , with respect to the non-maximal diagonal symmetry group J = Z/3Z [START_REF] Fan | Witten's D 4 Integrable Hierarchies Conjecture[END_REF].

2.4. D 4 double ramification hierarchy. In [START_REF] Givental | Simple singularities and integrable hierarchies[END_REF][START_REF] Frenkel | Soliton equations, vertex operators, and simple singularities[END_REF] it was proved that the Dubrovin-Zhang hierarchy for the Dubrovin-Saito CohFT associated to a Coxeter group coincides with the Drinfeld-Sokolov hierarchy of the corresponding semisimple Lie algebra. As explained above, the ε → 0 limit of both the DZ and DR hierarchies of any (semisimple) CohFT coincides with the principal hierarchy of the Frobenius manifold which is, consequently, completely determined, thanks to the recursion equation (2.3), by the Hamiltonian g

[0] 1,1 = g 1,1 | ε=0 .
To compute the latter in the D 4 case we can use (see [START_REF] Dubrovin | Normal forms of hierarchies of integrable PDEs, Frobenius manifolds and Gromov-Witten invariants, a new 2005 version of[END_REF]) the dilaton equation g

[0] 1,1 = (D -2)(F (t * )| t * =u * )dx, to obtain g [0] 1,1 = u 1 (u 2 ) 2 12 + (u 1 ) 2 u 3 12 - (u 2 ) 3 u 3 108 + (u 2 ) 2 (u 3 ) 3 432 + (u 3 ) 7 326592 + u 1 (u 4 ) 2 4 + u 2 u 3 (u 4 ) 2 12 + (u 3 ) 3 (u 4 ) 2 144 dx.
(2.6) Moreover, thanks to homogeneity of the CohFT, we know that, for any k ≥ 0,

(2.7)

|u 1 k | = 1, |u 2 k | = |u 4 k | = 2 3 , |u 3 k | = 1 3 , |ε| = 1 6 , |g 1,1 | = 7 3 .
Notice here how all the variables have positive degree, guaranteeing the polynomiality of g 1,1 and, in fact, of all the Hamiltonian densities g α,d , 1 ≤ α ≤ N , d ≥ -1.

Theorem 2.2. The double ramification hierarchy for the D 4 Dubrovin-Saito cohomological field theory is the integrable hierarchy with Hamiltonian operator

K DR u =     0 0 6∂ x 0 0 6∂ x 0 0 6∂ x 0 0 0 0 0 0 2∂ x     e 1 → e 1 e 2 → e 2 e 3 → e 3 e 4 → -e 4 , Z 3 :        e 1 → e 1 e 2 → -1 2 e 2 -3 2 e 4 e 3 → e 3 e 4 → 1 2 e 2 -1 2 e 4 .
In fact one can directly verify that the corresponding actions on the coordinates (t 1 , . . . , t 4 ) leave the Frobenius potential (2.4) unchanged and, by unicity of the Givental-Teleman reconstruction, this symmetry is inherited by the Dubrovin-Saito CohFT at all genera.

It was proved in [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF] that the restriction of a CohFT with a finite symmetry to the invariant subspace of V is a partial CohFT (i.e. a system of linear maps satisfying all the axioms of a CohFT with the exception of the gluing axiom at non-separating nodes). The authors then remark that, although the Dubrovin-Zhang hierarchy cannot be directly constructed for partial CohFTs, the two restrictions of the D 4 Dubrovin-Zhang hierarchy to the invariant sectors with respect to the two above symmetries give well defined integrable systems, isomorphic respectively to the B 3 and G 2 Drinfeld-Sokolov hierarchies.

In contrast, in [START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF] it was shown how that the double ramification hierarchy construction works for partial CohFTs too. All this then proves the following theorem.

Theorem 2.4. The double ramification hierarchy associated to the restriction of the D 4 Dubrovin-Saito CohFT to the invariant subspace with respect to the Z 2 action described in (2.9) is the integrable hierarchy with Hamiltonian operator The double ramification hierarchy associated to the restriction of the D 4 Dubrovin-Saito CohFT to the invariant subspace with respect to the Z 3 action described in (2.9) is the integrable hierarchy with Hamiltonian operator In case the CohFT is homogeneous of degree deg c g,n = δ(g -1), then we have the following grading for the quantum DR hierarchy:

K DR u =   0 0 6∂ x 0 6∂ x 0 6∂ x 0 0   and
K DR u = 0 6∂ x 6∂ x 0 and
|u α k | = 1 -deg e α , |ε| = 1 -δ 2 , | | = 2 -δ, |G α,d | = d + 3 -δ -|u α |
In light of the results of section 2, the DR quantization procedure can be applied to obtain a quantization of the D 4 Drinfeld-Sokolov hierarchy, similarly to what was achieved in [START_REF] Buryak | Towards a description of the double ramification hierarchy for Witten's r-spin class[END_REF] for the A n Gelfand-Dickey hierarchies, for n ≤ 4.

3.2. Quantum systems of DR type and the quantum D 4 hierarchy. In [START_REF] Buryak | Integrable systems of double ramification type[END_REF] the quantum version of Theorem 2.1 was also proved. We recall it here as, again, it will be applied to compute the D 4 quantum DR hierarchy. Let us consider the quantum Hamiltonian system defined by a Hamiltonian H ∈ ( Λ ) [≤0] with respect to the standard quantum commutator introduced above. Consider the operator D H : A

[[z]] → A [[z]] defined by D H f (z) = ∂ x • (D -1) - z [•, H], D := k≥0 u α k + ε ∂ ∂ε + 2 ∂ ∂ .
We warn the reader that we have upgraded the definition of the operator D. Suppose there exist N solutions G α (z) ∈ ( A ) Again, we call a system of densities like the one described in the theorem above a quantum integrable system of double ramification type and we recall from [START_REF] Buryak | Integrable systems of double ramification type[END_REF] that the quantum

1. 1 .

 1 Matrix Lax equations of type D 4 . Let us briefly explain what these positive and negative flows are in terms of Lie algebras. Type D 4 is the Dynkin diagram of the simple Lie algebra o(8). We denote by

  [LWZ10], Theorem 4.11), which gives the complete picture of the scalar Lax pairs representation of the Drinfeld-Sokolov hierarchy of type D 4 . Theorem 1.3. The flows (1.2), (1.3) of the Drinfeld-Sokolov hierarchy of type D 4 coincide with the flows of Equation (1.11) 1.3. Bi-Hamiltonian structure. Here we give the two compatible Poisson brackets and Hamiltonian densities for the first Hamiltonian structure of the Drinfeld-Sokolov hierarchy of type D 4 . They were given in [DS84], Proposition 8.3. To do so, we need to introduce the operator

2. 2 .

 2 Integrable systems of DR type. Let {•, •} denote the Poisson structure associated to the Hamiltonian operator η µν ∂ x . For a local functional h ∈ Λ [0]

(2

  Hamiltonian densities defined by equation (2.3), where g 1,1 is obtained from the D 4 Hamiltonian (2.8) by imposing u 2 * = u 4 * = 0. It is equivalent to the G 2 Drinfeld-Sokolov hierarchy. 3. Quantum double ramification hierarchy for the D 4 Dubrovin-Saito CohFT 3.1. Quantum double ramification hierarchy. In [BR16b] Buryak's original definition of the integrable hierarchy associated to a CohFT was upgraded to include a quantization of the classical double ramification hierarchy. It consists of a system of quantum differential polynomials G α,d ∈ ( A ) [≤0] , 1 ≤ α ≤ N , d ≥ -1 , where A = A[[ ]] and deg = -2, such that G α,d | =0 = g α,d . Similarly we define Λ = A /(Im∂ x ⊕ C[[ε, ]]) and the images G α,d = G α,d dx of G α,d in this quotient commute with respect to a star product that canonically quantizes the Hamiltonian operatorK DR = η∂ x , i.e. [G α,p , G β,q ] = 0, 1 ≤ α, β ≤ N, p, q ≥ -1where, for any f ∈ A andg ∈ Λ , [f , g] = [f, g]dx and [f, g] = n≥1 r 1 ,...,rn≥0 s 1 ,...,sn≥0 (n-1+ a i -j) C a 1 ,...,an j , if j = n -1 + n i=1 a i (mod2), 0, otherwise. (3.2) and k i=1 Li -d i (z) = k-1+ d i j=1 C d 1 ,...,d k j Li -j (z), Li -d (z) := k≥0 k d z k . (3.3)

  [≤0] [[z]], α = 1, . . . , N , to D H G α (z) = 0 with the initial conditions G α (z = 0) = η αµ u µ . Then a new vector of solutions can be found by the following transformation(3.4) G α (z) → A µ α (z)G µ (z) + B α (z), where A µ α (z) = δ µ α + i>0 A µ α,i z i ∈ C[[z]] and B α (z) = i>0 B α,i (ε, )z i ∈ C[[ε, , z]]. Theorem 3.1 ([[START_REF] Buryak | Integrable systems of double ramification type[END_REF]). Assume that H ∈ ( Λ )[≤0] has the following properties:(a) there exist N independent solutions G α (z) = p≥0 G α,p-1 z p ∈ ( A ) [≤0] [[z]], α = 1, . . . , N , to the equation (3.5) D H G α (z) = 0 with the initial conditions G α (z = 0) = η αµ u µ , u µ u ν + ∂ x R + c(ε, ), R ∈ ( A ) [≤-1] , c(ε, ) ∈ C[[ε, ]], (c) G 1,1 = H.Then, up to a transformation of type (3.4), we have(i) G 1,0 = 1 2 η µν u µ u ν dx, (ii) [G α,p , G β,q ] = 0, α, β = 1, . . . , N, p, q ≥ -1, (iii) 1 [G α,p , G β,0 ] = ∂ x ∂G α,p+1 ∂u β , β=1, . . . , N, p ≥ -1, (iv) ∂G α,p ∂u 1 = G α,p-1 , α = 1, . . . , N, p ≥ -1,

  . . , e 4 with deg e 1 = 0, deg e 2 = deg e 4 = 1 3 , deg e 3 = 2 3 , deg c D 4 g,n = δ = 2 3 . In genus 0 the corresponding Frobenius structure can described by the Frobenius potential (see for instance [LRZ15])

	(2.4)
	F =	t 1 t 2 2 12	+	t 1 2 t 3 12	-	t 2 3 t 3 216	+	t 2 2 t 3 3 1296	+	t 3 7 1632960	+	t 1 t 4 2 4	+	t 2 t 3 t 4 2 24	+	t 3 3 t 4 2 432	,

  Hamiltonian densities defined by equation (2.3), where g 1,1 is obtained from the D 4 Hamiltonian (2.8) by imposing u 4 k = 0. It is equivalent to the B 3 Drinfeld-Sokolov hierarchy.

Since the coordinate s 4 is treated differently then the others, we use letters µ, ν for coordinates s 1 , s

, s

and letters α, β for coordinates s 1 , s 2 , s 3 and s

.

Acknowledgments. We would like to express our gratitude to A. Buryak as well as to M. Cafasso and V. Roubtsov for useful discussions. A. du Crest de Villeneuve also thanks the University of Angers, France, the University of Bourgogne, France and the European Marie Skodowska-Curie action 'IPaDEGAN' for offering them the best working conditions while working on the present paper.

We find:

+ 5 2592 ( u 3 1 ) 4 + 13 720 ( u 2 2 ) 2 + 55 3888 u 3 ( u 3 1 ) 2 u 3 2 + 29 1080 u 1 2 u 3 2 + 67 15552 ( u 3 ) 2 ( u 3 2 ) 2 + 13 240 ( u 4 2 )

2

Finally, the density h 1,1 is given by h 1,1 = 36 91 res P 13 .

We find:

(1.17)

Then by identifying the two operators in Equation (1.18) and inverting Equation (1.16), we find the following Miura transformation from the normal coordinates u α to the coordinates v α :

Now recall that the variational differentials are given by

It follows that the variational differentials of the coordinates functionals are given by, for µ ∈ {1, 2, 3},

As we can see in Equations (1.22) and (1.23), for any µ ∈ {1, 2, 3}, we have Y µ (z) + = 0, and Y 4 (z) -= 0. It follows that { u 4 (x), u 4 (y)} 1 = 2δ (x -y) (Note that this computation does not depend on the integer n of D n , meaning that similarly, in the D n case, { u n (x), u n (y)} 1 = 2δ (x -y).) It is easy to see that for any µ ∈ {1, 2, 3}, { u µ (x), u 4 (y)} 1 = 0. The remaining components are computed in a similar fashion, this time using

for µ, ν ∈ {1, 2, 3}. Now we equivalently write ∂ k x instead of δ (k) (x -y), for their action is identical. We write down the components { u α (x), u β (y)} 1 in the following matrix:

(1.24)

2. Classical double ramification hierarchy for the D 4 Dubrovin-Saito CohFT

In this section we compute the double ramification hierarchy for the D 4 Dubrovin-Saito or Fan-Jarvis-Ruan-Witten cohomological field theory (the latter with respect to the non-maximal diagonal symmetry group J = Z/3Z) in the framework of homogeneous integrable systems of double ramification type. This means that we will find the unique dispersive deformation of double ramification type, according to [START_REF] Buryak | Integrable systems of double ramification type[END_REF], of the principal hierarchy associated to the Frobenius manifold for the D 4 Coxeter group [START_REF] Dubrovin | Differential geometry of the space of orbits of a Coxeter group[END_REF], compatible with the natural grading of this homogeneous CohFT.

and Hamiltonian densities defined by equation (2.3), where

Proof. As explained above we know that g 1,1 |ε = 0 is given by equation (2.6) and that, by definition, only even powers of ε appear in the Hamiltonian densities of any DR hierarchy. Now remark from (2.7) that all the variables have positive degree, guaranteeing the polynomiality of g 1,1 and, in fact, of all the Hamiltonian densities g α,d , 1 ≤ α ≤ N , d ≥ -1. In particular, since |ε| = 1 6 and in each monomials there are as many x-derivatives as powers of ε, we see that any term where the power of ε is bigger than 10 is either trivial or a total x derivative.

One can then verify, by direct computation, that up to rescaling of the variable ε, there exist a unique local functional g 1,1 ∈ Λ [0] such that: g 1,1 | ε=0 is given by (2.6), contains only even powers of ε, |g 1,1 | = 7 3 and g 1,1 is of double ramification type. In practice one writes down the most general polynomial deformation of (2.6) with the given degree and notices that imposing conditions (a) and (b) of theorem 2.1 determines, up to rescalings of ε, all the coefficients. Finally, to determine the correct normalization of ε, it is sufficient to compute the coefficient of a single monomial of g 1,1 containing ε. We can do this by recalling (see for instance [START_REF] Buryak | Tau-structure for the Double Ramification Hierarchies[END_REF]) that, for any CohFT,

Recall that the DR hierarchy possesses a natural tau-structure given by the Hamiltonian densities h α,d = δg α,d+1 δu 1 . The normal coordinates for this tau structures are u α = η αµ δg µ,0 δu 1 . Explicitly one obtains

Applying this change of coordinates one obtains the Hamiltonian operator

which coincides with (1.24), and the Hamiltonian

which, up to a ∂ x -exact term, agrees with (1.17). 2.5. B 3 and G 2 double ramification hierarchies. In [START_REF] Fan | Witten's D 4 Integrable Hierarchies Conjecture[END_REF] it is proved that the D 4 Dubrovin-Saito CohFT and the FJRW theory of the singularity D 4 : W = x 3 + xy 2 with symmetry group J 1 where J 1 (x, y) = (e 2πi 1 3 x, e 2πi 1 3 y) and of the singularity D T 4 : W = x 3 y + y 2 with symmetry group G max = J 2 where J 2 (x, y) = (e 2πi 1 6 x, e 2πi 1 2 ) are all isomorphic.

In [START_REF] Liu | BCFG Drinfeld-Sokolov hierarchies and FJRW-Theory[END_REF] it was shown that the FJRW theory for (D 4 , J 1 ) and (D T 4 , G max ) carry are invariant with respect to the action of further symmetry groups Z 3 and Z 2 respectively.

Via the isomorphism and using our presentation from section 2.3 for the Dubrovin-Saito CohFT, we can express the action of the generators of these symmetries on the phase space V in the following way: 3 the quantum correction to the classical DR Hamiltonian g 1,1 is a polynomial and imposing the conditions of theorem 3.1 determine it uniquely up to a normalization constant for . The latter is fixed by recalling that, by definition (see [START_REF] Buryak | Double ramification cycles and quantum integrable systems[END_REF]), for any CohFT we have Coef u 1 i G 1,1 = -1 24 dim V .