The Wigner-Lohe model for quantum synchronization and its emergent dynamics

TL;DR

This paper introduces the Wigner-Lohe model for quantum synchronization, derived from the Schrödinger-Lohe model via the Wigner formalism, and establishes conditions for complete synchronization in identical potential systems.

Contribution

It presents a new Wigner-based model for quantum synchronization and provides a theoretical framework for asymptotic complete synchronization.

Findings

L^2-distances between wave functions tend to zero asymptotically

Framework for complete synchronization in identical potential systems

Derivation of the Wigner-Lohe model from the Schrödinger-Lohe model

Abstract

We present the Wigner-Lohe model for quantum synchronization which can be derived from the Schr\"{o}dinger-Lohe model using the Wigner formalism. For identical one-body potentials, we provide a priori sufficient framework leading the complete synchronization, in which $L^2$-distances between all wave functions tend to zero asymptotically.

Equations

$$\displaystyle\begin{aligned} {\mathrm{i}}\partial_{t}\psi_{i}=-\frac{1}{2}\Delta\psi_{i}+V_{i}\psi_{i}+\frac{{\mathrm{i}}K}{2N}\sum_{k=1}^{N}\Big{(}\frac{\|\psi_{i}\|\psi_{k}}{\|\psi_{k}\|}-\frac{\langle\psi_{k},\psi_{i}\rangle\psi_{i}}{\|\psi_{i}\|\|\psi_{k}\|}\Big{)},\quad 1\leq i\leq N,\end{aligned}$$

$$w[\psi ,\varphi ](x,p)=\frac{1}{(2\pi {)}^d}{\int }_{{\mathbb{R}}^d}{e}^{iy\cdot p}\bar{\psi }\left(x+\frac{y}{2}\right)\varphi \left(x-\frac{y}{2}\right)dy.$$

$$w_j:=w[{\psi }_j],w_{jk}:=w[{\psi }_j,{\psi }_k],w_{jk}^{+}:=\mathrm{Re}{w}_{jk}\text{mbox}and{w}_{jk}^{-}:=\mathrm{Im}{w}_{jk}.$$

$$\begin{cases}\displaystyle\partial_{t}w_{j}+p\cdot\nabla_{x}w_{j}+\Theta[V](w_{j})=\frac{K}{N}\sum_{k=1}^{N}\Big{[}w^{+}_{jk}-\Big{(}\int w^{+}_{jk}dpdx\Big{)}w_{j}\Big{]},\quad x\in\mathbb{R}^{d},~{}t>0,\\ \displaystyle\partial_{t}w^{+}_{jk}+p\cdot\nabla_{x}w^{+}_{jk}+\Theta[V](w^{+}_{jk})\\ \displaystyle\hskip 28.45274pt=\frac{K}{2N}\sum_{\ell=1}^{N}\Big{[}w^{+}_{j\ell}+w^{+}_{\ell k}-\Big{(}\int(w^{+}_{j\ell}+w^{+}_{\ell k})dpdx\Big{)}w^{+}_{jk}+\Big{(}\int(w^{-}_{j\ell}+w^{-}_{\ell k})dpdx\Big{)}w^{-}_{jk}\Big{]},\\ \displaystyle\partial_{t}w^{-}_{jk}+p\cdot\nabla_{x}w^{-}_{jk}+\Theta[V](w^{-}_{jk})\\ \displaystyle\hskip 28.45274pt=\frac{K}{2N}\sum_{\ell=1}^{N}\Big{[}w^{-}_{j\ell}+w^{-}_{\ell k}-\Big{(}\int(w^{+}_{j\ell}+w^{+}_{\ell k})dpdx\Big{)}w^{-}_{jk}+\Big{(}\int(w^{-}_{j\ell}+w^{-}_{\ell k})dpdx\Big{)}w^{+}_{jk}\Big{]}.\end{cases}$$

$${\mathrm{i}}\partial_{t}\psi_{i}=-\frac{1}{2}\Delta\psi_{i}+V_{i}\psi_{i}+\frac{{\mathrm{i}}K}{2N}\sum_{k=1}^{N}\Big{(}\frac{\|\psi_{i}\|\psi_{k}}{\|\psi_{k}\|}-\frac{\langle\psi_{k},\psi_{i}\rangle\psi_{i}}{\|\psi_{i}\|\|\psi_{k}\|}\Big{)},$$

$$\Vert {\psi }_i(t)\Vert =\Vert {\psi }_i^0\Vert \text{mbox}fort\ge 0,1\le i\le N.$$

$$\mathrm{i}{\partial }_t{\psi }_j=-\frac{1}{2}\Delta {\psi }_j+V_j{\psi }_j+\frac{\mathrm{i}K}{2N}\sum _{k=1}^N\left({\psi }_k-⟨{\psi }_k,{\psi }_j⟩{\psi }_j\right).$$

$$V_i(x)={\Omega }_i:\text{mbox}constant,{\psi }_i(t,x)={\psi }_i(t),(t,x)\in {\mathbb{R}}_{+}\times {\mathbb{T}}^d.$$

$${\mathrm{i}}\frac{d\psi_{i}}{dt}=\Omega_{i}\psi_{i}+\frac{{\mathrm{i}}K}{2N}\sum_{k=1}^{N}\Big{(}\frac{|\psi_{i}|\psi_{k}}{|\psi_{k}|}-\frac{\langle\psi_{i},\psi_{k}\rangle\psi_{i}}{|\psi_{i}||\psi_{k}|}\Big{)}.$$

$${\psi }_i:=e^{-\mathrm{i}{\theta }_i},1\le i\le N$$

$${\dot{\theta}}_{i}\psi_{i}=\Omega_{i}\psi_{i}+\frac{{\mathrm{i}}K}{2N}\sum_{k=1}^{N}\Big{(}\psi_{k}-e^{-{\mathrm{i}}(\theta_{i}-\theta_{k})}\psi_{i}\Big{)}.$$

$${\dot{\theta }}_i={\Omega }_i+\frac{K}{N}\sum _{k=1}^N\sin ({\theta }_k-{\theta }_i).$$

$$\exists \underset{t\to \infty }{\lim }⟨{\psi }_i,{\psi }_j⟩={\alpha }_{ij}\in \mathbb{C}.$$

$$\exists \underset{t\to \infty }{\lim }|{\theta }_i(t)-{\theta }_j(t)|={\theta }_{ij}^{\infty }.$$

$$\underset{t\to \infty }{\lim }\Vert {\psi }_i(t)-{\psi }_j(t)\Vert =0,1\le i,j\le N.$$

$$\underset{t\to \infty }{\lim }⟨{\psi }_i(t),{\psi }_j(t)⟩=1.$$

$$D(\Psi ):=\underset{i,j}{\max }||{\psi }_i-{\psi }_j||.$$

$$\frac{d}{dt}D(\Psi)\leq K(D(\Psi))\Big{(}D(\Psi)-\frac{1}{2}\Big{)},\quad t>0.$$

$$K>0,V_i=V,\Vert {\psi }_i^0{\Vert }_{L^2}=1,1\le i\le N,D({\Psi }^0)<\frac{1}{2}.$$

$$D(\Psi (t))\le \frac{D({\Psi }^0)}{D({\Psi }^0)+(1-2D({\Psi }^0))e^{Kt}},t\ge 0.$$

$$\underset{K\to \infty }{\lim }\underset{t\to \infty }{\mathrm{limsup}}\underset{i,j}{\max }||{\psi }_i-{\psi }_j||=0.$$

$$z_{jk}(t):=⟨{\psi }_j,{\psi }_k⟩(t)=r_{jk}(t)+i{s}_{jk}(t),$$

$$\sum _{k=1}^N\mathrm{Re}{z}_{jk}(0)>0,\text{for any}j=1,\dots ,N.$$

$$|1-z_{jk}(t)|\lesssim e^{-Kt},\text{as}t\to \infty .$$

$$V_j(x)=V(x),1\le j\le N.$$

$$i{\partial }_t\psi =-\frac{1}{2}\Delta \psi +V\psi ,$$

$${\partial }_{t}w+p\cdot {\nabla }_{x}w+\Theta [V]w=0,$$

$$\Theta [V](w)(x,p):=-\frac{i}{(2\pi {)}^d}\int e^{i(p-p^{\prime })\cdot y}\left(V\left(x+\frac{y}{2}\right)-V\left(x-\frac{y}{2}\right)\right)w(x,p^{\prime })d{p}^{\prime }dy.$$

$${\partial }_t{w}_j+p\cdot \nabla w_j+\Theta [V]w_j=R_j,$$

$$R_j=$$

Full text

Acknowledgment. The work of P. Antonelli and P. Marcati is partially supported by PRIN grant 2015YCJY3A_003 and by G.N.A.M.P.A. (I.N.d.A.M.) , The work of S.-Y. Ha and D. Kim are partially supported by the National Research Foundation of Korea (NRF2014R1A2A205002096).

1. 1 Introduction

2. 2 Preliminaries

3. 2.1 The Schrödinger-Lohe model

4. 2.2 Previous results

5. 3 From Schrödinger-Lohe to Wigner-Lohe

6. 4 Emergent dynamics of the W-L model for identical potentials

7. 5 Quantum Hydrodynamics

1. Introduction

Synchronization represents a phenomenon in which rhythms of weakly coupled oscillators are adjusted to the common frequency due to their weak interactions. It is often observed in many complex systems, e.g., the flashing of fireflies, clapping of hands in a concert hall, and heartbeat regulation by pacemaker cells, etc., [1, 6, 7, 33, 35]. However, rigorous mathematical treatment of such collective phenomena were begun only several decades ago by two scientists Winfree [37] and Kuramoto [26, 27]. For the mathematical modeling of synchronization, they adopted a continuous dynamical system approach based on their heuristic and intuitive arguments. In this paper, we are mainly interested in quantum Lohe oscillators with all-to-all couplings under one-body potential. To fix the idea, consider a classical complete network consisting of $N$ nodes, where each pair of nodes is connected with an equal capacity which is assumed to be unity. We also assume that quantum Lohe oscillators with the same unit mass are positioned on the nodes of the underlying complete network. To avoid unnecessary physical complexity, we ignore entanglement and decoherence effects inherent to the quantum many-body systems. For a better physical modeling, such genuine quantum effects need to be taken into account.

Let $\psi_{i}=\psi_{i}(t,x)$ be the wave function of the $i$-th Lohe oscillator on the spatial domain $\mathbb{R}^{d}$. Then, the dynamics of Lohe oscillators with unit mass is governed by the Schrödinger-Lohe (S-L) model: for $(t,x)\in\mathbb{R}_{+}\times\mathbb{R}^{d}$.

[TABLE]

where $\|\cdot\|:=\|\cdot\|_{L^{2}}$ and $\langle\cdot,\cdot\rangle$ are the standard $L^{2}$ norm and an inner product on $\mathbb{R}^{d}$, and $V_{i}=V_{i}(x)$ and $K$ correspond to the one-body potential and nonnegative coupling strength, respectively. The S–L model (1.1) was first introduced by Australian physicist Max Lohe [28] several years ago as an infinite state generalization of the Lohe matrix model [29]. As discussed in [28, 29], quantum synchronization has received much attention from the quantum optics community because of its possible applications in quantum computing and quantum information [14, 23, 24, 25, 30, 36, 39, 40]. The emergent dynamics of the S-L system (1.1) has been partially treated in [11, 12] for some restricted class of initial data and a large coupling strength. Recently, a new approach based on the finite-dimensional reduction has been proposed in [5, 15] which significantly improve the previous results [11, 12] by the Lyapunov functional approach. However, a complete resolution of the synchronization problem for (1.1) is still far from complete.

Our main purpose of this paper is to present a quantum kinetic analogue of the S-L model (1.1) and study its emergent dynamics. The study on the quantum kinetic model for the Schrödinger equation dates back to Wigner’s paper [38], in which Wigner considered the quantum mechanical motion of a large ensemble of electrons in a vacuum under the action of the Coulomb force generated by the charge of the electrons. For the modeling of large ensemble, he introduced a quasi one-particle distribution function, so called the Wigner function and showed that it satisfies the quantum Liouville equation [8, 9, 19, 20, 32, 41].

Before we briefly describe our main results, we first recall the Wigner transform of wave function on $\mathbb{R}^{d}.$ For more basic facts on Wigner transforms we refer the reader to [22, 21].

Definition 1.1**.**

For any two wave functions $\psi,\phi\in L^{2}$, we define the Wigner transform

[TABLE]

If we choose $\psi=\phi$, then we write $w[\psi]:=w[\psi,\psi]$.

In order to shorten the formulas, we are going to introduce the following notation: if $\psi_{j}$, $j=1,\dotsc,N$ is the solution to the S-L system (1.1), then we write

[TABLE]

Our main results of this paper are as follows. First, we show that the Wigner transforms $w_{i}$ and $w_{ij}^{\pm}$ satisfies a coupled non-local system:

[TABLE]

Second, we derive a sufficient condition for the complete synchronization of the coupled system (1.2). Finally, we also investigate the hydrodynamic formulation for the Schrödinger-Lohe system (1.1) and derive synchronization estimates in some special cases.

The rest of this paper is organized as follows. In Section 2, we present the Schrödinger-Lohe model for quantum synchronization and discuss previous works on the complete synchronization of the S-L model. In Section 3, we derive our augmented Wigner-Lohe model from the S-L model using the Wigner transform. In Section 4, we present a priori complete synchronization estimates for some restricted class of initial data. In Section 5, we also discuss a hydrodynamic model which can be obtained from the S-L model for two-oscillator case.

2. Preliminaries

In this section, we briefly present the Schrödinger-Lohe (S-L) model for Lohe synchronization, and review earlier results on the synchronization problem for the S-L model.

2.1. The Schrödinger-Lohe model

As a phenomenological model for the quantum synchronization generalizing classical Kuramoto synchronization, Lohe proposed a coupled Schrödinger-type model in [28]. For $(t,x)\in\mathbb{R}_{+}\times\mathbb{R}^{d}$ and $1\leq i\leq N,$

[TABLE]

where we normalized $\hbar=1$ and $m=1$.

Lemma 2.1**.**

[28]* Let $\Psi=(\psi_{1},\cdots\psi_{N})$ be a smooth solution to (2.1) with initial data $\Psi_{0}=(\psi^{0}_{1},\cdots,\psi^{0}_{N})$. Then, the $L^{2}$ norm of $\psi_{i}$ is constant along the flow (2.1):*

[TABLE]

In view of the previous lemma, from now on we will assume that $\|\psi^{0}_{i}\|=1$, $1\leq i\leq N$, so that system (2.1) becomes

[TABLE]

For the space-homogeneous case, we set the spatial domain to be a periodic domain $\mathbb{T}^{d}$ and choose a special choice of $V_{i}$:

[TABLE]

system (2.1) can be reduced to the Kuramoto model which is a prototype model for classical synchronization. In this special setting, the S-L model becomes

[TABLE]

We next simply take the ansatz for $\psi_{i}$:

[TABLE]

and substitute this ansatz into (2.3) to obtain

[TABLE]

Then, we take the inner product of the above relation with $\psi_{i}$ and compare the real part of the resulting relation to get the Kuramoto model for classical synchronization [1, 6, 13, 16, 17, 18]:

[TABLE]

Thus, the S-L model can be viewed as a quantum generalization of the Kuramoto model.

2.2. Previous results

In this subsection, we briefly review the previous results [11, 12, 15, 5] on the complete synchronization of the S-L model. For this, we first recall the definition of the complete synchronization as follows.

Definition 2.1**.**

Let $\Psi=(\psi_{1},\cdots\psi_{N})$ be a smooth solution to (2.1) with initial data $\Psi^{0}=(\psi_{1}^{0},\cdots,\psi_{N}^{0})$. Then, the S-L model exhibits an asymptotic phase-locking if the following relations holds:

[TABLE]

Remark 2.1**.**

For the classical phase models such as (2.5), asymptotic phase-locking is defined as the following condition:

[TABLE]

Via the relation (2.4), we can see that (2.6) and (2.7) are closely related. In fact, in [12] for identical potentials $V_{i}=V_{j}$, the complete synchronization is defined as

[TABLE]

Note that the condition (2.8) and normalization condition $||\psi_{i}||=1$ yield

[TABLE]

Thus, the condition (2.8) satisfies the condition (2.6). Recently, in [5, 15] the case with different one-body potentials was treated, at least for $N=2$. In this framework it is shown that, in some regimes, the limit in (2.6) is not $1$ but depends on the difference between the potentials. Hence the limit in (2.8) gives a non-zero constant. This is indeed the more general case, when the system (2.2) exhibits complete frequency synchronization but not phase synchronization. For more details we address the reader to [5, 15].

As mentioned in the Introduction, the S-L model was first considered in Lohe’s work [28] for the non-Abelian generalization of the Kuramoto model. However, the first rigorous mathematical studies of the S-L model were treated by the second author and his collaborators in [12, 15] in two different methodologies. The first methodology is to use $L^{2}$-diameters for $\{\psi_{i}\}$ as a Lyapunov functional and derive a Gronwall type differential inequality to conclude the complete synchronization with $\alpha_{ij}=1$. More precisely, we set

[TABLE]

In [12], authors derived a differential inequality for the diameter $D(\Psi)$:

[TABLE]

This leads to an exponential synchronization of the (1.1).

Theorem 2.1**.**

[12]* Suppose that the coupling strength and initial data satisfy*

[TABLE]

Then, for any solution $\Psi=(\psi_{1},\dots,\psi_{N})$ to (1.1), the diameter $D(\Psi)$ satisfies

[TABLE]

Remark 2.2**.**

For distinct one-body potentials, we do not have an asymptotic phase-locking estimate for the S-L model yet, however in [11], for some restricted class of initial data and large coupling strength, a weaker concept of synchronization, namely practical synchronization estimates have been obtained:

[TABLE]

On the other hand, at least in the two oscillator case, it is possible to improve considerably the practical synchronization result: indeed in [5, 15] a complete picture of different regimes is shown, where the system (2.2) exhibits complete synchronization or dephasing, i.e. time periodic orbits for the correlation function.

Recently, an alternative approach to prove synchronization for the S-L model was proposed both in [5] and [15], by using a finite dimensional reduction. More precisely, in both papers the authors consider the correlations between the wave functions,

[TABLE]

and they study their asymptotic behavior. Moreover, in [5] the introduction of the order parameter, defined in analogy with the classical Kuramoto model, allows to give a more general result.

Theorem 2.2**.**

[5]** Let $(\psi_{1},\dotsc,\psi_{N})\in\mathcal{C}(\mathbb{R}_{+};L^{2}(\mathbb{R}^{d}))^{N}$ be the solution to (2.2) with initial data $(\psi_{1}(0),\dotsc,\psi_{N}(0))=(\psi^{0}_{1},\dotsc,\psi^{0}_{N})\in L^{2}(\mathbb{R}^{d})^{N}$, and we assume that

[TABLE]

Then we have

[TABLE]

As we will see in the next sections, the same approach used in [5, 15] will also be exploited to infer the synchronization results for the Wigner-Lohe model (3.1) and the hydrodynamical system (5.2). More precisely, for the quantum hydrodynamic system (5.2) we are going to need also some synchronization estimates proved at the $H^{1}$ regularity level. Such estimates are proved in [5].

3. From Schrödinger-Lohe to Wigner-Lohe

In this section we present a kinetic quantum analogue “the Wigner-Lohe(W-L) model” for the quantum synchronization, which can be derived from the Schrödinger-Lohe(S-L) model [28, 29] via the Wigner transform. In this and following sections, we assume that all one-body potentials are identical

[TABLE]

Recall that for a given a solution $\psi$ to the free Schrödinger equation:

[TABLE]

then its Wigner transform $w=w[\psi]$ satisfies

[TABLE]

where the operator $\Theta[V]$ is defined by

[TABLE]

Hence, to derive the Wigner-Lohe system (1.2) we just need to see how the nonlocal coupling in (2.2) translates at the Wigner level. More precisely, let $\psi_{j}$ be a solution to (2.2), then by defining $w_{j}=w[\psi_{j}]$, we see that it satisfies

[TABLE]

where the remainder term $R_{j}$ is given by

[TABLE]

where $r_{jk}(t):=\operatorname{Re}\langle\psi_{j},\psi_{k}\rangle(t)=\int w_{jk}^{+}(t,x,p)\,dxdp$. Let us recall that this last equality comes from one of the basic properties of Wigner transforms, namely

[TABLE]

for any $f,g\in L^{2}$. Resuming, the equation for $w_{j}$ is given by

[TABLE]

We now need to derive the equation for $w_{jk}=w[\psi_{j},\psi_{k}]$. Since the linear part in the S-L model (2.2) is common for every wave functions (remember we chose identical potentials, $V_{j}\equiv V$), then the linear part in the Wigner equation for $w_{jk}$ will be exactly the same as for $w_{j}$. Consequently we also have

[TABLE]

where

[TABLE]

Let us recall that $z_{jk}$ is defined in (2.9) and we notice that $\overline{z_{jk}}=z_{kj}$. Hence we obtain

[TABLE]

and the equation for $w_{jk}$ becomes

[TABLE]

By using definitions for $w_{jk}^{\pm}$ and the linearity of operator $\Theta[V]$, we then obtain the Wigner-Lohe system

[TABLE]

4. Emergent dynamics of the W-L model for identical potentials

In this section, we focus on the Wigner-Lohe model with $N=2$. In this case, system (3.1) becomes

[TABLE]

where we have

[TABLE]

Let us remark that the system (4.1), complemented with the definitions (4.2) above, can be considered independently on the S-L system (2.2). For such a system we will prescribe initial data $w_{1}^{0},w_{2}^{0},w_{12}^{0}$ such that $w_{1}^{0}$ and $w_{2}^{0}$ are real valued, $\int w_{1}^{0}\,dxdp=\int w_{2}^{0}\,dxdp=1$, $w_{12}^{0}$ is complex valued, $|\int w_{12}^{0}\,dxdp|\leq 1$. Let us also notice that the last equation is complex valued, so that we don’t split it into two coupled equations for $w_{12}^{+}$ and $w_{12}^{-}$ as in (3.1).

Let us now prove the synchronization for (4.1). First of all we remark that, by integrating the last equation over the whole phase space, we find the following ODE

[TABLE]

for which it is straightforward to give its asymptotic behavior.

Lemma 4.1**.**

Let $z_{12}(0)\in\mathbb{C}$ be such that $|z_{12}(0)|\leq 1$ and $z_{12}(0)\neq-1$, then the solution $z_{12}(t)$ to (4.3) satisfies

[TABLE]

Proof.

By integrating (4.3) we obtain

[TABLE]

∎

By using the Lemma above it is then possible to show the complete synchronization for the W-L model (4.1).

Theorem 4.1**.**

Let $(w_{1},w_{2},w_{12})$ be a solution to (4.1) with initial data $(w_{1}(0),w_{2}(0),w_{12}(0))=(w_{1}^{0},w_{2}^{0},w_{12}^{0})$ such that

[TABLE]

and

[TABLE]

Then we have

[TABLE]

Proof.

It follows from (4.1) that it is possible to write down the equation for the difference $w_{1}-w_{2}$,

[TABLE]

By multiplying it by $2(w_{1}-w_{2})$ and by integrating over the whole phase space we obtain

[TABLE]

By Lemma 4.1 we know that $|1-r_{12}(t)|\lesssim e^{-Kt}$, hence by Gronwall’s inequality we obtain the synchronization result. ∎

5. Quantum Hydrodynamics

In this Section, we derive the hydrodynamic equations associated to the Schrödinger-Lohe model (2.2). Here we follow the approach developed in [2, 3, 4] where a polar factorisation method is exploited in order to define the hydrodynamical quantities also in the vacuum region. In order to simplify the exposition we mainly focus on the case of two identical oscillators. In this case the Schrödinger-Lohe model reads

[TABLE]

The case with $N$ non-identical oscillators can be treated similarly with obvious modifications, but the study of this special case will simplify substantially the exposition.

In order to derive the hydrodynamics associated to system (5.1), we first need to ensure that it is globally well-posed in $H^{1}({\mathbf{R}}^{d})$. This is indeed a strightforward application of the standard for nonlinear Schrödinger equations [10], see for example Proposition 2.1 in [5]. Furthermore, let us notice that by defining $z_{12}(t)=\langle\psi_{1},\psi_{2}\rangle(t)$, then this function satisfies the ODE (4.3). This is not surprising because the W-L model was indeed derived from (5.1) and because of the property $\int w[\psi_{1},\psi_{2}]\,dxdp=\langle\psi_{1},\psi_{2}\rangle$. This implies that, under the same assumptions of Lemma 4.1, in this case we also have

[TABLE]

However, this synchronization result is too weak to be exploited for quantum hydrodynamic system derived from (5.1). Indeed, as we already remarked above, the natural space for the hydrodynamics is the finite energy space, namely $H^{1}$ for the wave functions. Hence we need to improve the result in the space of energy. Here we will make use of Theorem 4.5 in [5], where we address the reader for more general results in this direction.

Let us now consider the solution $(\psi_{1},\psi_{2})\in\mathcal{C}(\mathbb{R}_{+};H^{1})$ to system (5.1), given by Proposition 2.1 in [5]. To derive the hydrodynamic system associated with (5.1), we first define the mass densities, namely $\rho_{1}=|\psi_{1}|^{2}$ and $\rho_{2}=|\psi_{2}|^{2}$. By differentiating those quantities with respect to time and by using the equations above, we obtain

[TABLE]

where the associated current densities are respectively given by

[TABLE]

Furthermore, in the equation for the mass density we also find the interaction term $\rho_{12}=\operatorname{Re}(\bar{\psi}_{1}\psi_{2})$, so that $r_{12}=\operatorname{Re}\langle\psi_{1},\psi_{2}\rangle=\int\rho_{12}\,dx$.

Let us notice that $\rho_{12}$ is not a mass density, since in general it can also be negative. By using those definitions we can derive the evolution equations for the current densities $J_{1}$ and $J_{2}$. For instance, by differentiating $J_{1}$ with respect to time we find that

[TABLE]

where the new interaction term here is given by

[TABLE]

Next, we use the polar factorisation Lemma in [2, 3] to infer that, for $\psi_{1}\in H^{1}(\mathbb{R}^{d})$, we have

[TABLE]

where $\sqrt{\rho_{1}}=|\psi_{1}|,\Lambda_{1}=\operatorname{Im}(\bar{\phi}_{1}\nabla\psi_{1})$, $\phi_{1}$ is the polar factor for the wave function $\psi_{1}$ and we have $\sqrt{\rho_{1}}\Lambda_{1}=J_{1}$, see [2, 3, 4] for more details on the polar factorisation. In this way we can write down the following equation for the current density $J_{1}$:

[TABLE]

By using the equation for $\psi_{2}$ we obtain an analogous equation for $J_{2}$:

[TABLE]

Resuming, by defining the hydrodynamical quantities $\rho_{1},J_{1},\rho_{2},J_{2}$ associated to $\psi_{1},\psi_{2}$, respectively, we can derive the following system:

[TABLE]

Note that the above hydrodynamical system is not closed, as we need to derive also the evolution equations for the quantities $\rho_{12},J_{12}$. However, it is quite troublesome to derive a hydrodynamical equation for the quantity $J_{12}$. For this reason we consider the following auxiliary variables

[TABLE]

By using those variables, it is straightforward to derive their dynamical equations,

[TABLE]

where we denoted $\sigma_{12}=\operatorname{Im}(\bar{\psi}_{1}\psi_{2})$, so that

[TABLE]

Define $\psi_{d}:=\psi_{1}-\psi_{2}$, then by following some similar calculations as before we find out

[TABLE]

After some simple algebra, we obtain that

[TABLE]

where $G_{12}:=\frac{1}{2}\operatorname{Re}(\bar{\psi}_{2}\nabla\psi_{1}-\bar{\psi}_{1}\nabla\psi_{2})$. Hence the equation for $J_{d}$ is given by

[TABLE]

Once again, to close the hydrodynamic equations, we still need to determine the evolution for $\sigma_{12},G_{12}$. As before, the equation derived for $G_{12}$ would be too involved, for this reason we alternatively define $\psi_{a}=\psi_{1}-{\mathrm{i}}\psi_{2}$ and its hydrodynamical quantities $\rho_{a}=\frac{1}{2}|\psi_{a}|^{2}$, $J_{a}=\frac{1}{2}\operatorname{Im}(\bar{\psi}_{a}\nabla\psi_{a})$. If we write down the equation for $\psi_{a}$,

[TABLE]

we can then infer the equations for $\rho_{a}$ and $J_{a}$. By proceeding as before with some straightforward but long calculations, we find out

[TABLE]

We can now resume and write down the whole set of hydrodynamic equations associated to the Schrödinger-Lohe system (5.1):

[TABLE]

where

[TABLE]

By considering the system (5.2) above we can now prove the synchronization property. In view of the previosu synchronization results we expect that

[TABLE]

and furthermore

[TABLE]

To show the above properties we are going to use a synchronization result in $H^{1}$ for system (5.1) given in [5], which will be stated in the following Theorem. The result below actually holds in a more general case, see [5] for more details, however here we will state the synchronization property we are going to use for our system (5.2).

Theorem 5.1**.**

[5]** Let $(\psi^{0}_{1},\psi^{0}_{2})\in H^{1}$ be such that

[TABLE]

Then, for the solution $(\psi_{1},\psi_{2})\in\mathcal{C}(\mathbb{R}_{+};H^{1})$ emanated from such initial data, we have

[TABLE]

We apply now this result for the synchronization of system (5.2), First of all, from (5.3) we then infer that

[TABLE]

This and the polar factorisation Lemma 3 in [3] then readily implies that

[TABLE]

Furthermore, from the fact that $\psi_{d}\to 0$ in $H^{1}$ as $t\to\infty$ and the polar factorisation Lemma, again, we can also show that

[TABLE]