# unsupervised_multilingual_alignment_using_wasserstein_barycenter__ae7e3983.pdf

Unsupervised Multilingual Alignment using Wasserstein Barycenter

Xin Lian1,3 , Kshitij Jain2 , Jakub Truszkowski2 , Pascal Poupart1,2,3 and Yaoliang Yu1,3

1University of Waterloo, Waterloo, Canada 2Borealis AI, Waterloo, Canada 3Vector Institute, Toronto, Canada {x9lian, k22jain, yaoliang.yu, ppoupart}@uwaterloo.ca, {jakub.truszkowski}@borealisai.com

We study unsupervised multilingual alignment, the problem of ﬁnding word-to-word translations between multiple languages without using any parallel data. One popular strategy is to reduce multilingual alignment to the much simpliﬁed bilingual setting, by picking one of the input languages as the pivot language that we transit through. However, it is wellknown that transiting through a poorly chosen pivot language (such as English) may severely degrade the translation quality, since the assumed transitive relations among all pairs of languages may not be enforced in the training process. Instead of going through a rather arbitrarily chosen pivot language, we propose to use the Wasserstein barycenter as a more informative mean language: it encapsulates information from all languages and minimizes all pairwise transportation costs. We evaluate our method on standard benchmarks and demonstrate state-of-the-art performances.

1 Introduction Many natural language processing tasks, such as part-ofspeech tagging, machine translation and speech recognition, rely on learning a distributed representation of words. Recent developments in computational linguistics and neural language modeling have shown that word embeddings can capture both semantic and syntactic information. This led to the development of the zero-shot learning paradigm as a way to address the manual annotation bottleneck in domains where other vector-based representations must be associated with word labels. This is a fundamental step to make natural language processing more accessible. A key input for machine translation tasks consists of embedding vectors for each word. Mikolov et al. [2013b] were the ﬁrst to release their pre-trained model and gave a distributed representation of words. After that, more software for training and using word embeddings emerged. The rise of continuous word embedding representations has revived research on the bilingual lexicon alignment problem [Rapp, 1995; Fung, 1995], where the initial goal was to learn a small dictionary of a few hundred words by leveraging statistical similarities between two languages. Mikolov et al.

[2013a] formulated bilingual word embedding alignment as a quadratic optimization problem that learns an explicit linear mapping between word embeddings, which enables us to even infer meanings of out-of-dictionary words [Zhang et al., 2016; Dinu et al., 2015; Mikolov et al., 2013a]. Xing et al. [2015] showed that restricting the linear mapping to be orthogonal further improves the result. These pioneering works required some parallel data to perform the alignment. Later on, [Smith et al., 2017; Artetxe et al., 2017; Artetxe et al., 2018a] reduced the need of supervision by exploiting common words or digits in different languages, and more recently, unsupervised methods that rely solely on monolingual data have become quite popular [Gouws et al., 2015; Zhang et al., 2017b; Zhang et al., 2017a; Lample et al., 2018; Artetxe et al., 2018b; Dou et al., 2018; Hoshen and Wolf, 2018; Grave et al., 2019]. Encouraged by the success on bilingual alignment, the more ambitious task that aims at simultaneously and unsupervisedly aligning multiple languages has drawn a lot of attention recently. A naive approach that performs all pairwise bilingual alignment separately would not work well, since it fails to exploit all language information, especially when there are low resource ones. A second approach is to align all languages to a pivot language, such as English [Smith et al., 2017], allowing us to exploit recent progresses on bilingual alignment while still using information from all languages. More recently, [Chen and Cardie, 2018; Taitelbaum et al., 2019b; Taitelbaum et al., 2019a; Alaux et al., 2019; Wada et al., 2019] proposed to map all languages into the same language space and train all language pairs simultaneously. Please refer to the related work section for more details. In this work, we ﬁrst show that the existing work on unsupervised multilingual alignment (such as [Alaux et al., 2019]) amounts to simultaneously learning an arithmetic mean language from all languages and aligning all languages to the common mean language, instead of using a rather arbitrarily pre-determined input language (such as English). Then, we argue for using the (learned) Wasserstein barycenter as the pivot language as opposed to the previous arithmetic barycenter, which, unlike the Wasserstein barycenter, fails to preserve distributional properties in word embeddings. Our approach exploits available information from all languages to enforce coherence among language spaces by enabling accurate com-

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

positions between language mappings. We conduct extensive experiments on standard publicly available benchmark datasets and demonstrate competitive performance against current state-of-the-art alternatives.

2 Multilingual Lexicon Alignment

In this section we set up the notations and deﬁne our main problem: the multilingual lexicon alignment problem. Given m languages L1, . . . , Lm, each represented by a vocabulary Vi consisting of ni respective words. Following Mikolov et al. [2013a], we assume a monolingual word embedding Xi = [xi,1, . . . , xi,ni] Rni di for each language Li has been trained independently on its own data. We are interested in ﬁnding all pairwise mappings Ti k : Rdi Rdk that translate a word xi,ji in language Li to a corresponding word xk,jk = Ti k(xi,ji) in language Lk. In the following, for the ease of notation, we assume w.l.o.g. that ni n and di d. Note that we do not have access to any parallel data, i.e., we are in the much more challenging unsupervised learning regime. Our work is largely inspired by that of Alaux et al. [2019], which we review below ﬁrst. Along the way we point out some crucial observations that motivated our further development. Alaux et al. [2019] employ the following joint alignment approach that minimizes the total sum of mis-alignment costs between every pair of languages:

min Qi Od,Pik Pn

k=1,k =i Xi Qi Pik Xk Qk 2, (1)

where Qi Od is a d d orthogonal matrix and Pik Pn is an n n permutation matrix1. Since Qi is orthogonal, this approach ensures transitivity among word embeddings: Qi maps the i-th word embedding space Xi into a common space X, and conversely Q 1 i = Q i maps X back to Xi. Thus, Qi Q k maps Xi to Xk, and if we transit through an intermediate word embedding space Xt, we still have the desired transitive property Qi Q t Qt Q k = Qi Q k . The permutation matrix Pik serves as an inferred correspondence between words in language Li and language Lk. Naturally, we would again expect some form of transitivity in these pairwise correspondences, i.e., Pik Pkt Pit, which, however, is not enforced in (1). A simple way to ﬁx this is to decouple Pik into the product P i Pk, in the same way as how we dealt with Qi. This leads to the following variant:

argmin Qi Od,Pi Pn

k=1 Pi Xi Qi Pk Xk Qk 2 (2)

= argmin Qi Od,Pi Pn

i=1 Pi Xi Qi 1

k=1 Pk Xk Qk 2 (3)

= argmin Qi Od,Pi Pn min X Rn d

i=1 Pi Xi Qi X 2 (4)

1Alaux et al. [2019] also introduced weights αik > 0 to encode the relative importance of the language pair (i, k).

Figure 1: Comparing the Wasserstein barycenter and arithmetic mean (bottom panel) for two input distributions (top panel).

where Eq. 3 follows from the deﬁnition of variance and X in Eq. 4 admits the closed-form solution:

k=1 Pk Xk Qk. (5)

Thus, had we known the arithmetic mean language X beforehand, the joint alignment approach of Alaux et al. [2019] would reduce to a separate alignment of each language Xi to the mean language X that serves as the pivot. An efﬁcient optimization strategy would then consist of alternating between separate alignment (i.e., computing Qi and Pi) and computing the pivot language (i.e., (5)). We now point out two problems in the above formulation. First, a permutation assignment is a 1-1 correspondence that completely ignores polysemy in natural languages, that is, a word in language Li can correspond to multiple words in language Lk. To address this, we propose to relax the permutation Pi into a coupling matrix that allows splitting a word into different words. Second, the pivot language in (5), being a simple arithmetic average, may be statistically very different from any of the m given languages, see Figure 1 and below. Besides, intuitively it is perhaps more reasonable to allow the pivot language to have a larger dictionary so that it can capture all linguistic regularities in all m languages. To address this, we propose to use the Wasserstein barycenter as the pivot language. The advantage of using Wasserstein barycenter instead of the arithmetic average is that the Wasserstein metric gives a natural geometry for probability measures supported on a geometric space. In Figure 1, we demonstrate the difference between Wasserstein Barycenter and arithmetic average of two input distributions. It is intuitively clear that the Wasserstein barycenter preserves the geometry of the input distributions.

3 Our Approach We take a probabilistic approach, treating each language Li as a probability distribution over its word embeddings:

j=1 pijδxi j (6)

where pij is the probability of occurrence of the j-th word xi j in language Li (often approximated by the relative frequency of word xi j in its training documents), and δxi j is the

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

unit mass at xi j. We project word embeddings into a common space through the orthogonal matrix Qi Od. Taking a word xi from each language Li, we associate a cost c(Q1x1, . . . , Qmxm) R+ for bundling these words in our joint translation. To allow polysemy, we ﬁnd a joint distribution π with ﬁxed marginals πi so that the average cost Z c(Q1x1, . . . , Qmxm) dπ(x1, . . . , xm) (7)

is minimized. If we ﬁx Qi, then the above problem is known as multi-marginal optimal transport [Gangbo and Swikech, 1998]. To simplify the computation, we take the pairwise approach of Alaux et al. [2019], where we set the joint cost c as the total sum of all pairwise costs:

c(x1, . . . , xm) = P i,k xi xj 2. (8)

Interestingly, with this choice, we can signiﬁcantly simplify the numerical computation of the multi-marginal optimal transport. We recall the deﬁnition of Wasserstein barycenter ν of m given probability distributions π1, . . . , πm:

ν = arg min µ

i=1 λi W2 2(πi, µ), (9)

where λ 0 are the weights, and the (squared) Wasserstein distance W2 2 is given as:

W2 2(πi, µ)= min Πi Γ(πi,µ)

Z x y 2dΠi(x, y). (10)

The notation Γ(πi, µ) denotes all joint probability distributions (i.e. couplings) Πi with (ﬁxed) marginal distributions πi and µ. As proven by Agueh and Carlier [2011], with the pairwise distance (8), the multi-marginal problem in (7) and the barycenter problem in (9) are formally equivalent. Hence, from now on we will focus on the latter since efﬁcient computational algorithms for it exist. We use the push-forward notation (Qi)#πi to denote the distribution of Qixi when xi follows the distribution πi. Thus, we can write our approach succinctly as:

min µ min Qi Od

i=1 λi W2 2[(Qi)#πi, µ], (11)

where the barycenter µ serves as the pivot language in some common word embedding space. Unlike the arithmetic average in (5), the Wasserstein barycenter can have a much larger support (dictionary size) than the m given language distributions. We can again apply the alternating minimization strategy to solve (11): ﬁxing all orthogonal matrices Qi, we ﬁnd the Wasserstein barycenter using an existing algorithm of [Cuturi and Doucet, 2014] or [Claici et al., 2018]; ﬁxing the Wasserstein barycenter µ, we solve each orthogonal matrix Qi separately:

min Qi Od min Πi Γ(πi,µ)

Z Qix y 2dΠi(x, y). (12)

Algorithm 1: Barycenter Alignment

Input: Language distribution Li = (Xi, pi)m i=1, p Output: Translation for Lk and Lm for i = 1, . . . , m do

Xi Xi mean(Xi) {Ci} cosine dist(Xi,j, Xi,k) j, k Πi GW(Ci, C1, pi, p1) UΣV SVD(X i Πi X1) Qi UV Xi Xi Qi while not converged do

ν WB(π1, , πm; λ1, , λm) for i = 1, . . . , m do

Πi OT(πi, ν) UΣV SVD(X i Πi Y) Qi UV Xi Xi Qi ;

return (Π1, . . . , Πm, Q1, . . . , Qm)

For ﬁxed coupling Πi Rn s, where s is the dictionary size for the barycenter µ, the integral can be simpliﬁed as: X

jl(Πi)jl Qixi j yl 2 X i Πi Y, Qi . (13)

Thus, using the well-known theorem of Sch onemann [1966], Qi is given by the closed-form solution Ui V i , where UiΣi V i = X i Πi Y is the singular value decomposition. Our approach is presented in Algorithm 1.

4 Experiments

We evaluate our algorithm on two standard publicly available datasets: MUSE [Lample et al., 2018] and XLING [Glavas et al., 2019]. The MUSE benchmark is a high-quality dictionary containing up to 100k pairs of words and has now become a standard benchmark for cross-lingual alignment tasks [Lample et al., 2018]. On this dataset, we conducted an experiment with 6 European languages: English, French, Spanish, Italian, Portuguese, and German. The MUSE dataset contains a direct translation for any pair of languages in this set. We also conducted an experiment with the XLING dataset with a more diverse set of languages: Croatian (HR), English (EN), Finnish (FI), French (FR), German (DE), Italian (IT), Russian (RU), and Turkish (TR). In this set of languages, we have languages coming from three different Indo-European branches, as well as two non-Indo-European languages (FI from Uralic and TR from the Turkic family) [Glavas et al., 2019].

4.1 Implementation Details

To speed up the computation, we took a similar approach as Alaux et al. [2019] and initialized space alignment matrices with the Gromov-Wasserstein approach [Alvarez-Melis and Jaakkola, 2018] applied to the ﬁrst 5k vectors ( Alaux et al. [2019] used the ﬁrst 2k vectors) and with regularization parameter ϵ of 5e 5. After the initialization, we use the space alignment matrices to map all languages into the language

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

space of the ﬁrst language. Multiplying all language embedding vectors with the corresponding space alignment matrix, we realign all languages into a common language space. In the common space, we compute the Wasserstein barycenter of all projected language distributions. The support locations for the barycenter are initialized with random samples from a standard normal distribution. The next step is to compute the optimal transport plans from the barycenter distribution to all language distributions. After obtaining optimal transport plans Ti from the barycenter to every language Li, we can imply translations from Li to Lj from the coupling Ti T j . The coupling is not necessarily a permutation matrix, and indicates the probability with which a word corresponds to another. Method and code for computing accuracies of bilingual translation pairs are borrowed from Alvarez-Melis and Jaakkola [2018].

4.2 Baselines

We compare the results of our method on MUSE with the following methods: 1) Procrustes Matching with RSLS as similarity function to imply translation pairs [Lample et al., 2018]; 2) the state-of-the-art bilingual alignment method, Gromov Wasserstein alignment (GW) [Alvarez-Melis and Jaakkola, 2018]; 3) the state-of-the-art multilingual alignment method (UMH) [Alaux et al., 2019]; 4) bilingual alignment with multilingual auxiliary information (MPPA) [Taitelbaum et al., 2019b]; and 5) unsupervised multilingual word embeddings trained with multilingual adversarial training [Chen and Cardie, 2018]. We compare the results of our method on XLING dataset with Ranking-Based Optimization (RCSLS) [Joulin et al., 2018], solution to the Procrustes Problem (PROC) [Artetxe et al., 2018b; Lample et al., 2018; Glavas et al., 2019], Gromov Wasserstein alignment (GW) [Alvarez-Melis and Jaakkola, 2018], and VECMAP [Artetxe et al., 2018b]. RCLS and PROC are supervised methods, while GW and VECMAP are both unsupervised methods. The translation accuracies for Gromov-Wasserstein are computed using the source code released by [Alvarez-Melis and Jaakkola, 2018]. For the multilingual alignment method (UMH) [Alaux et al., 2019], and the two multilingual adversarial methods [Chen and Cardie, 2018], [Taitelbaum et al., 2019b], we directly compare our accuracies to previous methods as reported in [Glavas et al., 2019].

4.3 Results

Table 2 depicts precision@1 results for all bilingual tasks on the MUSE benchmark [Lample et al., 2018]. For most language pairs, our method Barycenter Alignment (BA) outperforms all current unsupervised methods. Our barycenter approach infers a potential universal language from input languages. Transiting through that universal language, we infer translation for all pairs of languages. From the experimental results in Table 2, we can see that our approach is clearly at an advantage and it beneﬁts from using the information from all languages. Our method achieves statistically signiﬁcant improvement for 22 out of 30 language pairs (p 0.05, Mc Nemar s test, one-sided).

Table 3 shows mean average precision (MAP) for 10 bilingual tasks on the XLING dataset [Glavas et al., 2019]. In Table 1, we show several German to English translations and compare the results to Gromov-Wasserstein direct bilingual alignment. Our method is capable of incorporating both the semantic and syntactic information of one word. For example, the top ten predicted English translations for the German word M unchen, are Cambridge, Oxford, Munich, London, Birmingham, Bristol, Edinburgh, Dublin, Hampshire, Baltimore . In this case, we hit the English translation Munich. What s more important in this example is that all predicted English words are the name of some city. Therefore, our method is capable of implying M unchen is a city name. Another example is the German word sollte, which means should in English. The top ﬁve words predicted for sollte are syntactically correct - would, could, will, should and might are all modal verbs. The last three examples show polysemous words, and in all these cases our method performs better than the Gromov-Wasserstein. For German word erschienen, our algorithm predicts all three words released, appeared, and published in the top ten translations as compared to Gromov Wasserstein which only predicts published .

4.4 Ablation Study

In this section, we show the impact of some of our design choices and hyperparameters. One of the parameters is the number of support locations. In theory, the optimal barycenter distribution could have as many support locations as the sum of the total number of support locations for all input distributions. In Figure 2, we show the impact on translation performance when we have a different number of support locations. Let nj be the number of words we have in language Lj. We picked the three most representative cases: the average number of words = Pm j=1 nj/m, twice the average number of words = 2 Pm j=1 nj/m, and the total number of words = Pm j=1 nj. As we increase the number of support locations for the barycenter distribution, we can see in Figure 2 that the performance for language translation improves. However, when we increase the number of support locations for the barycenter, the algorithm becomes costly. Therefore, in an effort to balance accuracy and computational complexity, we decided to use 10000 support locations (twice the average number of words). We also conducted a set of experiments to determine whether the inclusion of distant languages increases bilingual translation accuracy. Excluding two non-Indo-European languages Finnish and Turkish, we calculated the barycenter of Croatian (HR), English (EN), French (FR), German (DE), and Italian (IT). Figure 3 contains results for common bilingual pairs. The red bar show the bilingual translation accuracy when translating through the barycenter for all languages including Finnish and Turkish, whereas blue bar indicate the accuracy of translations that use the barycenter of the ﬁve Indo-European languages.

5 Related Work

We brieﬂy describe related work on supervised and unsupervised techniques for bilingual and multilingual alignment.

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

German English GW Prediction BA Prediction M unchen Munich London, Dublin, Oxford, Birmingham, Wellington Cambridge, Oxford, Munich, London, Birmingham Glasgow, Edinburgh, Cambridge, Toronto, Hamilton Rristol, Edinburgh, Dublin, Hampshire, Baltimore sollte should would, could, might, will, needs, would, could, will, supposed, supposed, put, willing, wanted, meant might, meant, needs, expected, able, should erschienen released published, editions, publication, edition published, editions, volumes, publication appeared printed, volumes, compilation released, titled, printed published publications, releases, titled appeared, edition, compilation aufgenommen admitted recorded, taken, recording, selected recorded, taken, recording, admitted recorded roll, placed, performing selected, sample, included taken, included eligible, motion, assessed track, featured, mixed viel lots, lot much, lot, little, more, less much, lot, little, less, too much bit, too, plenty, than, better more, than, bit, far, lots

Table 1: German-to-English translation prediction comparing results by 1) using GW alignment to imply direct bilingual mapping and 2) using Barycenter Alignment method described in Algorithm 1. We show top-10 translations of both methods. Last three examples show the polysemous words and their corresponding translations.

Figure 2: Accuracies for language pairs using different numbers of support locations for the barycenter. In our experimental setup, we have 5000 words in each language.

5.1 Supervised Bilingual Alignment Mikolov et al. [2013a] formulated the problem of aligning word embeddings as a quadratic optimization problem to ﬁnd an explicit linear map Q between the word embeddings X1 and X2 of two languages.

min Q ||X1Q PX2||2 2 (14)

This setting is supervised since the assignment matrix P that maps words of one language to another is known. Later, [Xing et al., 2015] showed that the results can be improved by restricting the linear mapping Q to be orthogonal. This corresponds to Orthogonal Procrustes [Sch onemann, 1966].

5.2 Unsupervised Bilingual Alignment In the unsupervised setting, the assignment matrix P between words is unknown, and we resort to the joint optimization:

min Q min P ||X1Q PX2||2 2. (15)

As a result, the optimization problem becomes non-convex and therefore more challenging. The problem can be relaxed into a (convex) semideﬁnite program.This method provides high accuracy at the expense of high computation complexity.

Figure 3: This graph shows the accuracy of bilingual translation pairs. The red bar indicate translation accuracy using the barycenter of all languages (HR, EN, FI, FR, DE, RU, IT, TR), while the blue bar correspond to the barycenter of (HR, EN, FR, DE, IT, RU).

Therefore, it is not suitable for large scale problems. Another way to solve (15) is to use Block Coordinate Relaxation, where we iteratively optimize each variable with other variables ﬁxed. When Q is ﬁxed, optimizing P can be done with the Hungarian algorithm in O(n3) time (which is prohibitive since n is the number of words). Cuturi and Doucet [2014] developed an efﬁcient approximation (complexity O(n2)) achieved by adding a negative entropy regularizer.Observing that both P and Q preserve the intra-language distances, Alvarez-Melis and Jaakkola [2018] cast the unsupervised bilingual alignment problem as a Gromov-Wasserstein optimal transport problem, and give a solution with minimum hyper-parameter to tune.

5.3 Multilingual Alignment In multilingual alignment, we seek to align multiple languages together while taking advantage of inter-dependencies to ensure consistency among them. A common approach consists of mapping each language to a common space X0 by minimizing some loss function l:

min Qi Od,Pi Pn

i l(Xi Qi, Pi X0) (16)

The common space may be a pivot language such as English [Smith et al., 2017; Lample et al., 2018; Joulin et al., 2018]. Nakashole and Flauger [2017] and Alaux et al. [2019] showed that constraining coherent word alignments between triplets of nearby languages improves the quality of induced

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

it-es it-fr it-pt it-en it-de es-it es-fr es-pt es-en es-de GW 92.63 91.78 89.47 80.38 74.03 89.35 91.78 92.82 81.52 75.03 GWo - - - 75.2 - - - - 80.4 - PA 87.3 87.1 81.0 76.9 67.5 83.5 85.8 87.3 82.9 68.3 MAT+MPPA 87.5 87.7 81.2 77.7 67.1 83.7 85.9 86.8 83.5 66.5 MAT+MPSR 88.2 88.1 82.3 77.4 69.5 84.5 86.9 87.8 83.7 69.0 UMH 87.0 86.7 80.4 79.9 67.5 83.3 85.1 86.3 85.3 68.7 BA 92.32 92.54 90.14 81.84 75.65 89.38 92.19 92.85 83.5 78.25

fr-it fr-es fr-pt fr-en fr-de pt-it pt-es pt-fr pt-en pt-de GW 88.0 90.3 87.44 82.2 74.18 90.62 96.19 89.9 81.14 74.83 GWo - - - 82.1 - - - - - - PA 83.2 82.6 78.1 82.4 69.5 81.1 91.5 84.3 80.3 63.7 MAT+MPPA 83.1 83.6 78.7 82.2 69.0 82.6 92.2 84.6 80.2 63.7 MAT+MPSR 83.5 83.9 79.3 81.8 71.2 82.6 92.7 86.3 79.9 65.7 UMH 82.5 82.7 77.5 83.1 69.8 81.1 91.7 83.6 82.1 64.4 BA 88.38 90.77 88.22 83.23 76.63 91.08 96.04 91.04 82.91 76.99

en-it en-es en-fr en-pt en-de de-it de-es de-fr de-pt de-en Average GW 80.84 82.35 81.67 83.03 71.73 75.41 72.18 77.14 74.38 72.85 82.84 GWo 78.9 81.7 81.3 - 71.9 - - - - 72.8 78.04 PA 77.3 81.4 81.1 79.9 73.5 69.5 67.7 73.3 59.1 72.4 77.98 MAT+MPPA 78.5 82.2 82.7 81.3 74.5 70.1 68.0 75.2 61.1 72.9 78.47 MAT+MPSR 78.8 82.5 82.4 81.5 74.8 72.0 69.6 76.7 63.2 72.9 79.29 UMH 78.9 82.5 82.7 82.0 75.1 68.7 67.2 73.5 59.0 75.5 78.46 BA 81.45 84.26 82.94 84.65 74.08 78.09 75.93 78.93 77.18 75.85 84.24

Table 2: Pairs of languages in multilingual alignment problem results for English, German, French, Spanish, Italian, and Portuguese. All reported results are precision@1 percentage. The method achieving the highest precision for each bilingual pair is highlighted in bold. Methods we are comparing to in the table are: Procrustes Matching with CSLS metric to infer translation pairs (PA) [Lample et al., 2018]; Gromov-Wasserstein alignment (GW) [Alvarez-Melis and Jaakkola, 2018] (reproduced by us using their source code); GWo refers to the results reported by Alvarez-Melis and Jaakkola [2018] in the paper; bilingual alignment with multilingual auxiliary information (MPPA) [Taitelbaum et al., 2019b]; Multilingual pseudo-supervised reﬁnement method [Chen and Cardie, 2018]; multilingual alignment method (UMH) [Alaux et al., 2019]. Asterisks denote signiﬁcant differences between BA and GW (Mc Nemar s test, one-sided), the only methods for which predictions were available.

en-de it-fr hr-ru en-hr de-ﬁ tr-fr ru-it ﬁ-hr tr-hr tr-ru PROC (1k) 0.458 0.615 0.269 0.225 0.264 0.215 0.360 0.187 0.148 0.168 PROC (5k) 0.544 0.669 0.372 0.336 0.359 0.338 0.474 0.294 0.259 0.290 PROC-B 0.521 0.665 0.348 0.296 0.354 0.305 0.466 0.263 0.210 0.230 RCSLS (1k)0.501 0.637 0.291 0.267 0.288 0.247 0.383 0.214 0.170 0.191 RCSLS (5k)0.580 0.682 0.404 0.375 0.395 0.375 0.491 0.321 0.285 0.324 VECMAP 0.521 0.667 0.376 0.268 0.302 0.341 0.463 0.280 0.223 0.200 GW 0.667 0.751 0.683 0.123 0.454 0.485 0.508 0.634 0.482 0.295 BA 0.683 0.799 0.667 0.646 0.508 0.513 0.512 0.601 0.481 0.355

Table 3: Mean average precision (MAP) accuracies of several current methods on XLING dataset.

bilingual lexicons. Chen and Cardie [2018] extended the work of [Lample et al., 2018] to the multilingual case using adversarial algorithms. Taitelbaum et al. extended Procrustes Matching to the multi-Pairwise case [Taitelbaum et al., 2019b], and also designed an improved representation of the source word using auxiliary languages [Taitelbaum et al., 2019a].

6 Conclusion

In this paper, we discussed previous attempts to solve the multilingual alignment problem, compared similarity between the approaches and pointed out a problem with existing formulations. Then we proposed a new method using the Wasserstein

barycenter as a pivot for the multilingual alignment problem. At the core of our algorithm lies a new inference method based on an optimal transport plan to predict the similarity between words. Our barycenter can be interpreted as a virtual universal language, capturing information from all languages. The algorithm we proposed improves the accuracy of pairwise translations compared to the current state-of-the-art method..

Acknowledgments

We thank the reviewers for their critical comments and we are grateful for funding support from NSERC and Mitacs.

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)

[Agueh and Carlier, 2011] Martial Agueh and Guillaume Carlier. Barycenters in the Wasserstein space. SIAM Journal on Mathematical Analysis, 43(2):904 924, 2011. [Alaux et al., 2019] Jean Alaux, Edouard Grave, Marco Cuturi, and Armand Joulin. Unsupervised Hyper-alignment for Multilingual Word Embeddings. In ICLR, 2019. [Alvarez-Melis and Jaakkola, 2018] David Alvarez-Melis and Tommi S Jaakkola. Gromov-Wasserstein Alignment of Word Embedding Spaces. In ACL, 2018. [Artetxe et al., 2017] Mikel Artetxe, Gorka Labaka, and Eneko Agirre. Learning bilingual word embeddings with (almost) no bilingual data. In ACL, 2017. [Artetxe et al., 2018a] Mikel Artetxe, Gorka Labaka, and Eneko Agirre. Generalizing and Improving Bilingual Word Embedding Mappings with a Multi-Step Framework of Linear Transformations. In AAAI, 2018. [Artetxe et al., 2018b] Mikel Artetxe, Gorka Labaka, and Eneko Agirre. A robust self-learning method for fully unsupervised cross-lingual mappings of word embeddings. In ACL, 2018. [Chen and Cardie, 2018] Xilun Chen and Claire Cardie. Unsupervised Multilingual Word Embeddings. EMNLP, 2018. [Claici et al., 2018] Sebastian Claici, Edward Chien, and Justin Solomon. Stochastic Wasserstein Barycenters. In ICML, 2018. [Cuturi and Doucet, 2014] Marco Cuturi and Arnaud Doucet. Fast computation of wasserstein barycenters. ICML, 2014. [Dinu et al., 2015] Georgiana Dinu, Angeliki Lazaridou, and Marco Baroni. Improving zero-shot learning by mitigating the hubness problem. In ICLR workshop, 2015. [Dou et al., 2018] Zi-Yi Dou, Zhi-Hao Zhou, and Shujian Huang. Unsupervised Bilingual Lexicon Induction via Latent Variable Models. In EMNLP, 2018. [Fung, 1995] Pascale Fung. Compiling bilingual lexicon entries from a non-parallel english-chinese corpus. In Third Workshop on Very Large Corpora, 1995.

[Gangbo and Swikech, 1998] Wilfrid Gangbo and Andrzej Swikech. Optimal maps for the multidimensional Monge Kantorovich problem. Communications on Pure and Applied Mathematics, 51(1):23 45, 1998. [Glavas et al., 2019] Goran Glavas, Robert Litschko, Sebastian Ruder, and Ivan Vulic. How to (Properly) Evaluate Cross-Lingual Word Embeddings: On Strong Baselines, Comparative Analyses, and Some Misconception. In ACL, 2019. [Gouws et al., 2015] Stephan Gouws, Yoshua Bengio, and Greg Corrado. Bil BOWA: Fast Bilingual Distributed Representations without Word Alignments. In ICML, 2015. [Grave et al., 2019] Edouard Grave, Armand Joulin, and Quentin Berthet. Unsupervised Alignment of Embeddings with Wasserstein Procrustes. In AISTATS, 2019.

[Hoshen and Wolf, 2018] Yedid Hoshen and Lior Wolf. Non Adversarial Unsupervised Word Translation. In EMNLP, 2018. [Joulin et al., 2018] Armand Joulin, Piotr Bojanowski, Tomas Mikolov, Herv e J egou, and Edouard Grave. Loss in Translation: Learning Bilingual Word Mapping with a Retrieval Criterion. In EMNLP, 2018. [Lample et al., 2018] Guillaume Lample, Alexis Conneau, Marc Aurelio Ranzato, Ludovic Denoyer, and Herv e J egou. Word translation without parallel data. In ICLR, 2018. [Mikolov et al., 2013a] Tomas Mikolov, Quoc V Le, and Ilya Sutskever. Exploiting Similarities among Languages for Machine Translation. Co RR, abs/1309.4168, 2013. [Mikolov et al., 2013b] Tomas Mikolov, Ilya Sutskever, Kai Chen, Greg S Corrado, and Jeff Dean. Distributed Representations of Words and Phrases and their Compositionality. In NIPS, 2013. [Nakashole and Flauger, 2017] Ndapandula Nakashole and Raphael Flauger. Knowledge Distillation for Bilingual Dictionary Induction. In EMNLP, 2017. [Rapp, 1995] Reinhard Rapp. Identifying Word Translations in Non-parallel Texts. In ACL, 1995. [Sch onemann, 1966] Peter H Sch onemann. A generalized solution of the orthogonal procrustes problem. Psychometrika, 31(1):1 10, Mar 1966. [Smith et al., 2017] Samuel L Smith, David HP Turban, Steven Hamblin, and Nils Y Hammerla. Ofﬂine bilingual word vectors, orthogonal transformations and the inverted softmax. In ICLR, 2017. [Taitelbaum et al., 2019a] Hagai Taitelbaum, Gal Chechik, and Jacob Goldberger. Multilingual word translation using auxiliary languages. In EMNLP, 2019. [Taitelbaum et al., 2019b] Hagai Taitelbaum, Gal Chechik, and Jacob Goldberger. A multi-pairwise extension of procrustes analysis for multilingual word translation. In EMNLP, 2019. [Wada et al., 2019] Takashi Wada, Tomoharu Iwata, and Yuji Matsumoto. Unsupervised Multilingual Word Embedding with Limited Resources using Neural Language Models. In ACL, 2019. [Xing et al., 2015] Chao Xing, Dong Wang, Chao Liu, and Yiye Lin. Normalized word embedding and orthogonal transform for bilingual word translation. In NAACL, 2015. [Zhang et al., 2016] Yuan Zhang, David Gaddy, Regina Barzilay, and Tommi Jaakkola. Ten pairs to tag-multilingual pos tagging via coarse mapping between embeddings. In ACL, 2016. [Zhang et al., 2017a] Meng Zhang, Yang Liu, Huanbo Luan, and Maosong Sun. Earth Mover s Distance Minimization for Unsupervised Bilingual Lexicon Induction. In EMNLP, 2017. [Zhang et al., 2017b] Meng Zhang, Yang Liu, Huanbo Luan, and Maosong Sun. Adversarial Training for Unsupervised Bilingual Lexicon Induction. In ACL, 2017.

Proceedings of the Twenty-Ninth International Joint Conference on Artiﬁcial Intelligence (IJCAI-20)