Ce travail de thèse a pour objectif d'étudier les propriétés des molécules de métaux alcalins M2+ immergées dans des agrégats de néon, par le biais de simulations numériques. Nous développons une approche globale dans laquelle la détermination de la structure électronique se réduit à un problème à un électron. Ce dernier évolue dans un potentiel modélisé par des pseudo-potentiels semi-locaux à cœur polarisable. Nous les avons paramétrés après avoir calculé les courbes de potentiel des dimères MNe et M+Ne de manière ab initio. Nous effectuons une dynamique moléculaire classique, en y incorporant un traitement des couplages non-adiabatiques grâce à un algorithme de saut de surface. Nous avons trouvé les géométries d'équilibre des systèmes M2+Nen jusqu'à la première couche de solvatation de la molécule. Nous en avons déduit les propriétés statiques de ces systèmes en examinant les énergies de liaison, les distances d'équilibre, et les spectres optiques d'absorption. Nous avons ensuite étudié la dynamique de ces systèmes placés initialement sur un état excité. Nous avons établi que le taux de photodissociation dépend fortement du nombre d'atomes de néon et de l'ordre des transitions électroniques. Nous avons observé un effet de cage pour les systèmes Li2+Nen à partir de 18 atomes de néon. Nous avons également effectué des analyses sur la distribution des fragments produits, sur les états moléculaires stabilisés, et sur la localisation de la charge dans les systèmes asymétriques.