ARTS ET CULTUREBUSINESSLes ExpertsOlivier EzrattyTECH

Comprendre l’informatique quantique – qubits

Par Olivier Ezratty, expert FrenchWeb

On peut comprendre le fonctionnement d’un ordinateur quantique sans trop se plonger dans la mécanique quantique au-delà de la compréhension de ses mécanismes de base, vus dans dans l’épisode précédent de cette série d’articles, et surtout celui de l’intrication. Par contre, il faut se plonger un peu dans quelques éléments de mathématiques, d’algèbre linéaire et de trigonométrie, ce que nous allons commencer à faire ici.

Le premier élément de base d’un ordinateur quantique est l’inévitable qubit. Vous avez certainement déjà entendu parler de cet objet mystérieux capable d’être simultanément dans la valeur 0 et 1. Les explications courantes s’arrêtent le plus souvent là et vous tombez immédiatement dans l’expectative, vous demandant comment cela peut ensuite bien fonctionner.

Nous allons donc commencer ici par expliquer le fonctionnement logique, mathématique et matériel de ces qubits. Dans la partie suivante de cette série d’articles, nous irons plus loin en décrivant tour à tour les registres, les portes, l’organisation et l’architecture complète d’un ordinateur quantique. A chaque fois, lorsque nécessaire, nous ferons le parallèle avec les ordinateurs traditionnels.

Pour pouvoir suivre cette partie, il faut connaître quelques basiques mathématiques : la trigonométrie, les vecteurs et matrices et avoir déjà entendu parler des nombres complexes.

(ci-dessus, illustration de la conférence “Le quantique, c’est fantastique” du 14 juin 2018 à Nantes au Web2day par Hélène Pouille, c.f. son site et la vidéo de la conférence)

Le principe des qubits

Les qubits sont les éléments de manipulation de base de l’information dans les ordinateurs quantiques. Ils s’opposent aux bits de l’informatique traditionnelle. Avec eux, on passe d’un monde déterministe à un monde probabiliste.

Dans l’informatique traditionnelle, les bits correspondent à des charges électriques circulantes qui traduisent le passage d’un courant électrique ou son absence. Un bit est de valeur 1 si le courant passe soit de 0 si le courant ne passe pas. La lecture d’un bit donne 1 ou 0. Elle est déterministe, à savoir que si l’on répète l’opération de lecture plusieurs fois, ou l’opération de lecture après une réédition du calcul, on obtiendra normalement le même résultat.

C’est vrai aussi bien pour le stockage de l’information que pour son transport et pour les traitements dans des processeurs. Ceci est valable modulo les erreurs qui peuvent intervenir dans le parcours. Celles-ci interviennent le plus souvent au niveau de la mémoire et sont corrigées via des systèmes… de correction d’erreurs utilisant de la redondance.

Dans un qubits, rien à voir ! Si les qubits sont généralement initialisés à 0, les opérations portant dessus vont généralement les amener à avoir un état de superposition entre 0 et 1. Ces états correspondent à l’état de base ‘'(“ground state”) et à l’état excité (“excited state”) d’un système quantique à deux états possibles. Ces qubits peuvent donc être à la fois à la valeur 0 et 1, et dans une proportion qui est variable et qui correspond à la notion de superposition d’états évoquée dans la partie précédente sur les fondements de la mécanique quantique. A la fin des calculs, lorsque l’on lit la valeur d’un qubit, on retrouve 0 ou 1. La richesse des valeurs du qubits se manifeste donc uniquement pendant les calculs et non pas à leur initialisation où lors de leur lecture à la fin des calculs. C’est un concept que vous ne comprendrez complètement que lorsque nous aurons décrit quelques algorithmes quantiques lors d’une partie suivante de cette série.

Ici, nous allons d’abord creuser le modèle mathématique de représentation des qubits et comprendre comment on peut se le représenter physiquement et mentalement.

Nous ferons alors un tour des différents types de qubits physiques. Les modèles mathématiques de représentation des qubits ne dépendent pas de leur type physique. Seules les caractéristiques de l’ordinateur sont affectées comme le taux d’erreur et la nature des portes quantiques physiques de base dites “universelles” agissant sur les qubits sachant que toutes les portes quantiques sont exécutables sur les ordinateurs quantiques.

La sphère de Bloch de représentation de l’état des qubits

Dans un modèle probabiliste classique, un pbit ou bit probabiliste aurait une probabilité p d’avoir la valeur 0 et 1-p d’avoir la valeur 1. Ce serait un modèle probabiliste linéaire.

Dans un qubit, c’est bien différent ! Le modèle de représentation mathématique de l’état d’un qubit s’appuie sur la fameuse sphère de Bloch, qui m’a donné bien du fil à retordre en termes de compréhension. Je vais partager avec vous ce que j’ai pu en comprendre après des mois de recherche dans des dizaines de cours, articles scientifiques et livres sur le calcul quantique.

Ce modèle est lié à la représentation de l’état d’un qubit ou de tout quantum à deux états par un vecteur à deux dimensions dont la longueur dite “norme” est toujours de 1. Ce vecteur a la particularité de comporter deux éléments : un nombre réel (alpha) et un nombre complexe (beta).

Dans la sphère de Bloch, l’état |0> d’un quantum à deux états est figuré par la position d’un vecteur de longueur 1 allant du centre de la sphère vers le pôle Nord de la sphère et l’état |1> est un vecteur allant du centre de la sphère à son pôle Sud.

Les états intermédiaires sont représentés par des vecteurs partant du centre de la sphère qui sont toujours de longueur 1 avec un angle theta par rapport à la verticale z et un angle psi par rapport à l’axe x situé allant du centre de la sphère à son équateur et autour de l’axe z. Histoire de simplifier les choses, les états |0> et |1> qui sont opposés dans la sphère de Bloch sont dits “orthogonaux” alors qu’ils sont opposés dans la sphère. Nous allons rapidement comprendre pourquoi.

Règle de Max Born et probabilités

Les équations décrivant l’état d’un qubit indiquent que celui-ci est la superposition de l’état |0> et de l’état |1>. Dans les équations, alpha est un nombre réel qui décrit la probabilité d’obtenir l’état |0> et beta est un nombre complexe qui décrit celle d’avoir l’état |1>. La somme des probabilités des deux états doit donner 1. Mais ce n’est pas alpha + beta. C’est alpha au carré plus beta au carré qui donnent 1. Pourquoi donc ? C’est une question qui m’a taraudé des mois pendant la préparation de ce texte.

Ce modèle probabiliste a été élaboré par Max Born en 1926. Il donne au carré du module de la fonction d’onde d’un quantum la signification d’une densité de probabilité de présence d’une particule élémentaire. C’est lié au fait que l’état |0> et l’état |1> correspondent non pas à une position précise d’une particule mais sont représentés par la fonction d’onde de Schrödinger qui décrit la distribution probabiliste de l’état du quantum dans le temps et dans l’espace. Elle est ici appliquée dans l’espace.

L’état d’un qubit est représenté par un vecteur à deux dimensions, ce qui est une information bien plus riche qu’un 0 ou un 1 d’un simple bit ou même qu’une probabilité linéaire entre 0 et 1. Ce vecteur à deux dimensions comprend les deux composantes alpha et beta que nous venons de définir.

Le passage de la fonction d’onde de Schrödinger à la représentation alpha/beta du qubit avec un nombre réel et un nombre complexe m’échappe encore. Je tente une explication qui fera peut-être s’étrangler les physiciens. Si on prend par exemple deux orbites possibles d’un électron autour du noyau d’un atome dans un modèle d’hydrogène simple (il en existe plein…), la fonction d’onde de l’orbite basse dans l’espace sera une sorte de courbe de Gauss décalée de celle de l’orbite haute. Intuitivement, je me dis que la probabilité d’avoir l’électron dans l’orbite haute ou basse est proportionnelle à l’intégrale des fonctions d’onde correspondantes à ces orbites. D’où le carré qui serait une approximation de la valeur du pic de la gaussienne, à supposer que ce pic corresponde à la taille des vecteurs alpha et beta de représentation des états des qubits. M’enfin, je ne suis pas sûr !

Pour en savoir plus et avec une meilleure exactitude scientifique, vous pouvez consulter la fiche Wikipedia de la fonction d’onde ainsi que celle de la probabilité d’amplitude. On trouve d’autres explications dans l’exemple des niveaux d’orbites d’électrons dans l’atome d’hydrogène dans Quantum Mechanics and the hydrogen atom (19 slides). L’interprétation physique de la règle statistique de Max Born reste en tout cas ouverte si l’on en juge par ce papier de juin 2018 d’Arkady Bolotin, Quantum probabilities and the Born rule in the intuitionistic interpretation of quantum mechanics (14 pages).

Trigonométrie dans la sphère de Bloch

Second mystère à résoudre, pourquoi donc l’angle theta est-il divisé par deux dans les équations décrivant un état quantique dans la sphère de Bloch dans les calculs de sinus et de cosinus des formules donnant alpha et beta ? Cela vient de ce que l’état |1> est placé en bas de la sphère pour que l’espace des états du qubit occupe toute la sphère et pas seulement son hémisphère nord.

C’est décrypté dans The Bloch Sphere de Ian Glendinning, 2005 (33 slides) qui explique cela par l’orthogonalité mathématique des deux états |0> et |1> qui sont pourtant opposés dans la sphère de Bloch. C’est encore mieux expliqué dans Why is theta/2 used for a Bloch sphere instead of theta ? qui a définitivement éclaircit ce mystère pour moi. Je m’en suis inspiré pour créer le schéma ci-dessous. C’est un peu tarabiscoté et vous n’êtes pas obligés d’y passer du temps.

Comme l’état |1> est mathématiquement orthogonal par construction à l’état |0>, leur angle “mathématique” calculé doit être de 90° (pi/2). Or, l’angle theta est du double de cet angle dans la sphère car il est de pi (pour mémoire pi*2 = circonférence de 360°). Donc, on divise theta par deux pour relier la représentation géométrique dans la sphère avec la représentation mathématique de l’état du qubit. Et surtout, pour permettre un étalement de tous les états d’un qubits sur l’ensemble de la sphère.

J’illustre ceci dans le schéma composite ci-dessus qui sert à comprendre la chose pour les lecteurs patients et comprenant les bases de la trigonométrie et des nombres complexes. Dans l’histoire, alpha est toujours un nombre réel car c’est un simple cosinus. Seul beta peut-être un nombre complexe. Il l’est dès lors que le qubit n’est pas dans le plan croisant l’axe x (theta = 0) et l’axe z (phi = 0) de la sphère de Bloch. Ce nombre complexe associe une partie réelle pour la direction x et une partie complexe pour la dimension y qui est orthogonale à x.

Un point important à noter est que la représentation sur la sphère de Bloch est un modèle mathématique probabiliste. Elle ne correspond pas à un modèle physique, comme l’angle de polarisation d’un photon ou le spin d’un électron, malgré les similitudes. Sachant néanmoins que ce modèle a été créé à l’origine pour représenter le spin d’électrons.

Lorsque le vecteur d’état du qubit est horizontal dans la sphère, c’est-à-dire qu’il va jusqu’à son équateur (schéma ci-dessous), nous sommes dans un état superposant l’état 0 et l’état 1 à égalité, mais avec une phase variable qui est liée à l’angle horizontal du vecteur phi par rapport à l’axe z comme dans le schéma ci-dessous. Pourquoi une phase ? Parce que ce vecteur est lié à la fonction d’onde de Schrödinger et que l’état d’un quantum est une fonction d’onde additionnant les fonctions d’onde de ses états de base !

Cette information riche d’un qubit est modifiée ensuite par des portes quantiques. Une porte quantique unitaire, à savoir, qui s’applique à un seul qubit, applique une rotation à l’état du qubit dans sa sphère de Bloch. Cette rotation est appliquée via une matrice de nombres complexes 2×2, dite matrice orthogonale de déterminant 1. De déterminant 1 car elle ne va pas modifier la longueur du vecteur après son application. Cette longueur restera toujours de 1. Nous examinerons la diversité de ces portes quantiques dans la partie suivante de cette série.

J’ai l’impression qu’en règle générale, les portes quantiques ne génèrent pas toutes les positions de vecteurs dans la sphère de Bloch. ce sont souvent des quarts de tours. Les points de la sphère les plus souvent utilisés sont les points cardinaux : le 0, le 1, puis les quatre points correspondant à la superposition 0 et 1 qui sont sur l’équateur de la sphère.

On doit finalement cette sphère de Bloch à trois personnes : Erwin Schrödinger pour sa fonction d’onde de 1926, Max Born pour son modèle probabiliste associé, créé la même année et à Felix Bloch (1903-1983, Suisse) qui a représenté l’état d’un quantum à deux niveaux sur la sphère en 1946. En optique et pour décrire la polarisation d’un photon, la sphère de Bloch s’appelle la sphère de Poincaré, du nom du mathématicien français Henri Poincaré, mort en 1912 et cousin germain du président Raymond Poincaré.

Voici quelques sources d’information associées à cette partie : Lectures on Quantum Computing de Dan C. Marinescu et Gabriela M. Marinescu, 2003 (274 pages), The Bloch Sphere de Ian Glendinning, 2005 (33 slides), The statistical interpretation of quantum mechanics, discours d’acceptation du prix Nobel de physique de Max Born en 1954 (12 pages) ainsi que l’excellent livre The mathematics of quantum mechanics de Martin Laforest, 2015 (111 pages), qui décrit les basiques mathématiques de l’informatique quantique avec les nombres complexes, les vecteurs, les matrices et tout le toutim.

Cycle de vie d’un qubit

Au bilan, le qubit est un objet mathématique un peu particulier :

  • On l’initialise toujours à 0, correspondant à l’état de base du qubit.
  • On le modifie ensuite de manière programmatique avec des portes quantiques pour lui faire prendre des valeurs qui sont des vecteurs dans la sphère de Bloch. La porte de Hadamard est l’une des plus courantes et elle créé un état de superposition entre un 0 et un 1. Les manipulations mathématiques de ce vecteur consistent ensuite à le faire tourner dans la sphère de Bloch avec des portes quantiques unitaires que nous verrons dans la partie suivante consacrée à la description du fonctionnement d’un ordinateur quantique. Ces manipulations reviennent à multiplier le vecteur représentant le qubit [alpha, beta] par une matrice de deux lignes et deux colonnes de nombres complexes conservant la norme du vecteur, qui doit rester à 1. Les portes quantiques non unitaires relient les qubits entre eux et en font évaluer les valeurs de manière conditionnelle. Sans ces différentes portes quantiques, on ne pourrait pas faire grand chose avec les qubits.
  • L’information dans les qubits qui est manipulée lors des calculs est “riche” avec une dimension de deux nombres réels, les angles theta et phi.
  • A part le cas où le qubit est initialisé à l’état |0> ou inversé à l’état |1>, celui-ci est en état de superposition. La représentation mathématique d’un qubit et son incarnation visuelle dans la sphère de Bloch montrent qu’un qubit peut-être dans une infinité d’états superposés différents.
  • Lorsque l’on lit la valeur du qubit, on retombe sur un 0 ou un 1 avec un retour probabiliste dépendant des paramètres du vecteur de l’état du qubit dans la sphère de Bloch.
  • Donc, nous avons un 0 en entrée, un 0 ou un 1 en sortie, et une infinité d’états entre les deux pendant les calculs !

Voici cela illustré dans le schéma ci-dessous :

Tout ça pour dire que la richesse mathématique du qubit intervient pendant les traitements et seulement pendant les traitements. Mais pas au départ ni à l’arrivée des traitements. La magie mathématique de l’ordinateur quantique est là !

Types de qubits

D’un point de vue physique, les qubits des calculateurs quantiques sont des dispositifs matériels qui intègrent des particules élémentaires qui ont deux états possibles que l’on peut initialiser, modifier avec des portes quantiques puis dont on peut évaluer l’état. Il s’agit parfois de particules élémentaires unitaires, comme avec les ions piégés ! Une seule à la fois ! Dans le cas des supraconducteurs ou de quantum dots à silicium, l’état quantique s’appuie sur un grande nombre de particules, ici, des électrons.

Voici les principaux types de qubits qui sont étudiés, expérimentés voir utilisés en production actuellement :

Les supraconducteurs : le qubit prend la forme de l’état d’un courant supraconducteur qui traverse une barrière très fine en s’appuyant sur l’effet Josephson. Il existe plusieurs types de qubits supraconducteurs, de flux, de phase et de charge. Ne rentrons pas dans les détails. Dans tous les cas, il s’agit de créer une superposition de deux états bien distincts d’un courant oscillant à haute fréquence et traversant la jonction Josephson dans une boucle supraconductrice. L’oscillation est rendue possible par le fait que la boucle intègre l’équivalent d’une inductance et d’une résistance. L’oscillation du courant est activée par l’application de micro-ondes de fréquences situées entre 5 et 10 GHz transmises par voie conductrice et physique. Ce ne sont pas des ondes émises par la voie “radio”. L’état du qubit est pour sa part mesuré avec un magnétomètre intégré dans le circuit.

C’est la technique la plus couramment employée aujourd’hui, notamment par IBM, Google et Intel pour des ordinateurs quantiques à circuits universels et avec les ordinateurs quantiques adiabatiques du Canadien D-Wave qui utilisent un autre arrangement de qubits de qualité moins bonne du côté du bruit et du taux d’erreurs. Elle est relativement facile à fabriquer car elle s’appuie sur les techniques de création de circuits CMOS même si certains des matériaux sont différents, comme le niobium qui est utilisé chez D-Wave. Cf Practical realization of Quantum Computation Supraconducting (36 slides). Notons, nous le verrons plus tard, que l’équipe de Daniel Estève au CEA de Saclay fait partie des précurseurs de la création de tels qubits supraconducteurs.

Les ions piégés : il s’agit d’ions d’atomes qui peuvent être du calcium et qui sont maintenus sous vide et suspendus par suspension électrostatique. Un pompage optique est réalisé pour leur initialisation. Un laser sert à la mesure et exploite le phénomène de fluorescence des ions excités par le laser. Le magnétisme est utilisé pour l’activation des portes quantiques. La startup IonQ issue de l’Université de Maryland planche là-dessus tout comme l’université d’Innsbruck en Autriche et sa spinoff AQT. Dans un tel système, plusieurs ions sont piégés de manière équidistante les uns des autres. Ils sont alignés en rang d’oignons. Il est difficile de faire “scaler” ce genre d’arrangement au-delà d’une centaine et quelques d’ions. Par contre, ils présentent l’avantage de bien pouvoir être intriqués les uns avec les autres, ce qui est moins le cas des qubits à base de supraconducteurs.

Les photons : leur état quantique est leur polarisation horizontale ou verticale. Cela fait partie du champ de l’optique linéaire. Il manipule des photons individuels. Les portes quantiques sont réalisées à l’aide de dispositifs optiques avec des filtres dichroïques ou polarisants. Il faut un grand nombre de lasers pour piloter l’ensemble. C’est donc pour l’instant assez embarrassant. L’avantage est que ces qubits fonctionnent à température ambiante. Mais ce type de qubit n’est pour l’instant utilisé qu’en laboratoire et à petite échelle et aucune société privée ne semble s’être embarquée dans ce chemin.

Le spin d’électron : haut ou bas, une sorte de sens de polarisation magnétique, que l’on retrouve dans les ordinateurs à base de quantum dots, notamment chez Intel  ou dans des prototypes de qubits réalisés au CEA à Grenoble. Ces qubits sont intégrés dans des circuits à base de semi-conducteurs CMOS. Ils bénéficient donc de la réutilisation de processus de fabrication de composants CMOS déjà bien maitrisés. Ces qubits sont cependant pour l’instant plutôt “bruyants”, même en les exploitant à des températures proches du zéro absolu. Le plan consiste à en aligner des batteries pour créer des qubits logiques, un concept que nous étudierons dans la prochaine partie.

Les centres NV pour Nitrogen Vacancy : ce sont des structures de diamant artificiel dans lesquelles un atome de carbone a été remplacé par un atome d’azote et à proximité duquel se situe une lacune d’atome de carbone. Les électrons des atomes de carbone adjacents à cette lacune occupent l’espace vide et avec un spin variable. Bref, l’état de la lacune est instable et quantique. Cet état est excité par laser et micro-ondes. La lecture de l’état du qubit est réalisée par une mesure de brillance de fluorescence. Seule la startup QDTI planche commercialement sur cette technique. Ce n’est visiblement pas encore au point.

Voici, ci-dessous, un schéma d’ensemble du mécanisme de contrôle de ces qubits vu dans Forefront engineering of nitrogen-vacancy centers in diamond for quantum technologies 2017 (235 pages).

Les fermions de Majorana : ce sont des anyons ou quasi-particules qui sont des états particuliers de nuages d’électrons organisées par paires. Pratiquement, ce sont des spins d’électrons aux deux bouts de fils supraconducteurs. On peut d’ailleurs considérer à ce titre là que c’est une variation des qubits supraconducteurs. De ce fait, ces ordinateurs quantiques doivent aussi être refroidis à une température voisine du zéro absolu, aux alentours de 10mK. C’est la voie choisie par Microsoft. Sachant que l’existence des fermions de Majorana est à peine prouvée. Et les quasi-particules sur lesquelles planchent les équipes de Microsoft ne seraient pas des fermions de Majorana. Il y a de quoi y perdre son latin ! Nous reviendrons dessus dans une prochaine partie de cette série.

Aucune de ces techniques n’est pour l’instant éprouvée à grande échelle. Elles ont toutes leurs avantages et inconvénients qui se situent dans plusieurs dimensions : la stabilité des qubits, la possibilité de les intriquer, le niveau d’erreurs dans les qubits et dans les portes qui les relient, la durée de cohérence des qubits, la température de fonctionnement, le niveau de miniaturisation, le processus de fabrication (qui est moins cher lorsque c’est du silicium en technologie CMOS). Nous verrons tout cela plus tard quand nous aborderons les différentes offres et projets d’ordinateurs quantiques.

En résumé, voici les principaux types de qubits et leurs caractéristiques clés avec la nature des qubits, la manière de stocker leur état, la nature des portes quantiques qui les modifient et de la mesure de l’état des qubits.

Leur niveau d’avancement est décrit dans cet excellent document de l’équivalent allemand de l’ANSSI, Entwicklungsstand Quantencomputer (état des lieux de l’informatique quantique, 231 pages, en anglais). Il évoque d’autres technologies que celles qui sont citées dans mon inventaire et que je ne cite pas car elles ont pour l’instant peu de chances d’aboutir et n’ont été adoptées par aucun industriel ou aucune startup.

______________________________

Voilà pour l’histoire des qubits. Mais celle-ci ne fait que commencer. Dans la partie suivante, nous allons aller plus loin et découvrir comment sont architecturés les ordinateurs quantiques dans la pratique avec leurs registres, leurs portes, leurs dispositifs de lecture, leurs systèmes de correction d’erreur ainsi que leur système de cryogénie pour ceux que cela concerne comme les ordinateurs à qubits supraconducteurs. Nous pourrons comprendre comment fonctionne un cycle complet de calcul dans un ordinateur quantique. La partie d’après donnera corps à l’ensemble avec un descriptif des principaux algorithmes connus pouvant être exploités sur les ordinateurs quantiques.

L’expert:

Olivier-Ezratty

Olivier Ezratty est consultant en nouvelles technologies et auteur d’Opinions Libres, un blog sur les médias numériques (TV numérique, cinéma numérique, photo numérique), et sur l’entrepreneuriat (innovation, marketing, politiques publiques…). Olivier est expert pour FrenchWeb.

Découvrez WE, le nouveau media d'intelligence économique consacré à l'innovation en europe. Retrouvez les informations de plus de 4500 startups et 600 fonds d'investissements Pour en savoir plus, cliquez ici

2 commentaires

  1. « L’état d’un qubit est représenté par un vecteur à deux dimensions, ce qui est une information bien plus riche qu’un 0 ou un 1 d’un simple bit ou même qu’une probabilité linéaire entre 0 et 1. Ce vecteur à deux dimensions comprend les deux composantes alpha et beta que nous venons de définir. »

    Le vecteur est plutôt me semble-t-il, en réalité, à trois dimensions puisque la sphère est une surface(2D) dans un ESPACE 3D. Le fait d’utiliser un coefficient complexe pour décrire une des composantes d’un vecteur à deux dimensions me semble juste une façon astucieuse de se projeter dans un espace de dimension supérieure, de la même façon que l’utilisation des nombres complexes(qui peuvent se paramétrer avec deux nombre réels) permet facilement de se projeter dans un espace de dimension 2 comparé à la simple utilisation des nombres réels(dimension 1).

    1. Une surface, c’est un objet à 2 dimensions; un de ses points peut se représenter par 2 réels, par exemple longitude et latitude, ici pour la surface sphérique. Un complexe, c’est bien 2 degrés de liberté, mais s’il est contraint à avoir 1 pour module, cela ne fait plus qu’un degré de liberté (par exemple son argument).

Bouton retour en haut de la page
Share This