Ce protocole est la porte d’entrée d’un débutant dans le traitement, l’ajustement et l’interprétation des spectres d’absorption transitoire. L’objectif de ce protocole est la préparation d’ensembles de données et l’ajustement en utilisant à la fois la cinétique à longueur d’onde unique et l’analyse globale de la durée de vie. Les défis associés aux données d’absorption transitoire et à leur ajustement sont discutés.
La spectroscopie d’absorption transitoire (TA) est une puissante méthode spectroscopique résolue en temps utilisée pour suivre l’évolution des processus à l’état excité à travers les changements dans le spectre d’absorption du système. Les premières mises en œuvre de l’AT étaient confinées à des laboratoires spécialisés, mais l’évolution des systèmes commerciaux clés en main a rendu la technique de plus en plus accessible aux groupes de recherche du monde entier. Les systèmes modernes d’AT sont capables de produire de grands ensembles de données avec une résolution énergétique et temporelle élevée et riches en informations photophysiques. Cependant, le traitement, l’ajustement et l’interprétation des spectres TA peuvent être difficiles en raison du grand nombre de caractéristiques à l’état excité et d’artefacts instrumentaux. De nombreux facteurs doivent être soigneusement pris en compte lors de la collecte, du traitement et de l’ajustement des données d’AT afin de réduire l’incertitude quant au modèle ou à l’ensemble de paramètres d’ajustement qui décrit le mieux les données. L’objectif de la préparation et de l’ajustement des données est de réduire le plus grand nombre possible de ces facteurs étrangers tout en préservant les données pour l’analyse. Dans cette méthode, les débutants reçoivent un protocole de traitement et de préparation des données TA ainsi qu’une brève introduction aux procédures et modèles d’ajustement sélectionnés, en particulier l’ajustement à longueur d’onde unique et l’analyse globale de la durée de vie. Des commentaires sur un certain nombre de défis courants en matière de préparation des données et sur les méthodes pour les résoudre sont fournis, suivis d’une discussion sur les défis et les limites de ces méthodes d’ajustement simples.
La spectroscopie d’absorption transitoire (TA) est une technique spectroscopique résolue en temps qui surveille l’évolution des espèces photoexcitées par des changements temporels de leur spectre d’absorption après une excitation avec une impulsion lumineuse. Comme l’AT est une technique absorbante, les signaux spectroscopiques provenant d’états qui subissent à la fois des transitions radiatives (c’est-à-dire des états qui émettent généralement un photon) et des transitions non radiatives (états qui sont généralement non fluorescents et subissent une conversion interne, un croisement intersystème ou participent à des photoréactions) peuvent être identifiés et leur évolution suivie 1,2. Selon les spécificités de la source d’excitation et de la méthode de détection, l’AT permet d’accéder à la cinétique de la femtoseconde à la microseconde au-delà et de l’UV à l’IR lointain, ce qui en fait un outil spectroscopique polyvalent. La commercialisation des spectromètres TA a considérablement progressé au cours des dernières décennies, ce qui a permis à un plus grand nombre de laboratoires et d’installations d’avoir accès à cette technique puissante2.
Les systèmes d’AT modernes sont capables de produire de grands ensembles de données avec une résolution énergétique et temporelle élevée. Les jeux de données prennent généralement la forme d’une matrice 2D de valeurs de différence de transmittance ou d’absorbance en fonction de la longueur d’onde et du retard temporel par rapport à l’impulsion d’excitation. Cet ensemble de données peut être considéré comme une carte thermique bidimensionnelle ou une carte topographique tridimensionnelle. L’interprétation de ces données est devenue plus complexe car les chercheurs s’efforcent d’inclure l’ensemble des données lors de la génération d’ajustements qui décrivent le mieux leur système d’intérêt3.
Bien que l’AT puisse couvrir une large gamme de longueurs d’onde et d’échelles de temps, ce protocole se concentre sur l’une de ses formes les plus largement accessibles4 : la spectroscopie à large bande dans la région UV-visible pilotée par un laser pulsé femtoseconde. Un schéma 5,6 d’un tel instrument est fourni à la figure 1. L’expérience commence par prendre une impulsion du laser et la diviser en deux copies. Une copie de l’impulsion, appelée « pompe », est utilisée pour exciter l’échantillon. Un dispositif tel qu’un amplificateur paramétrique optique (OPA) est généralement utilisé pour transformer l’impulsion de la pompe à la longueur d’onde d’excitation souhaitée 5,7. La deuxième copie de l’impulsion, appelée « sonde », entre dans une étape de retard mécanique, qui peut faire varier le délai entre les impulsions de la pompe et de la sonde en faisant varier la distance parcourue par l’impulsion. L’impulsion de la sonde à longueur d’onde unique est ensuite transformée en un continuum de lumière blanche à l’aide d’un cristal de saphir ou de fluorure de calcium (CaF2)8. L’impulsion de lumière blanche traverse l’échantillon et son spectre est mesuré à l’aide d’un détecteur à large bande tel qu’une caméra CCD (dispositif à couplage de charge). En mesurant les changements dans le spectre de l’impulsion de lumière blanche avec et sans la pompe, les changements dans le spectre d’absorption de l’échantillon induits par la pompe, ΔA(T), peuvent être mesurés. Les lecteurs intéressés sont invités à consulter cette revue utile9 pour plus d’informations sur le processus de détection.
Dans toutes les formes de spectroscopie TA, les spectres ΔA(t) sont calculés en prenant la différence entre l’absorption à l’état fondamental, unesonde, et celle de l’état excité, unepompe + sonde, à un délai donné, t, entre les deux impulsions 2,5,9,10.
(1)
Notez que lasonde A est équivalente au spectre d’absorption à l’état d’équilibre de l’échantillon et est indépendante du temps ; la résolution temporelle de l’expérience provient du délai entre la pompe et la sonde capturé dans Apompe+sonde(t). Une simulation de ces données est présentée à la figure 2A.
Contrairement aux spectres d’absorption à l’état stationnaire, les spectres TA peuvent avoir des caractéristiques positives et négatives en raison de la différence prise dans l’équation 1. Les caractéristiques positives sont le résultat de nouvelles espèces absorbantes créées par l’impulsion de la pompe et peuvent représenter des états chromophores excités, des états triplets, des réarrangements géométriques, des effets de solvatation ou des photoproduits à état excité3. Des lignes directrices générales pour identifier ces caractéristiques et les attribuer aux espèces chimiques seront présentées dans la discussion. Des caractéristiques négatives peuvent provenir de l’eau de Javel à l’état fondamental (GSB) ou de l’émission stimulée (SE) (Figure 2B). La GSB est due à la perte de population à l’état fondamental suite à l’absorption de l’impulsion de la pompe. Les molécules promues à l’état excité n’absorbent plus dans la même région que leur état fondamental ; par conséquent, moins d’impulsion de la sonde est absorbée et la différence dans l’équation 1 peut être négative dans cette région. Le GSB se caractérise par la même forme spectrale que celle de l’absorption à l’état fondamental mais avec un signe opposé. Les signaux SE résultent de l’émission d’une espèce à l’état excité stimulée par l’impulsion de la sonde3. L’émission de ces espèces entraîne une plus grande quantité de lumière atteignant le détecteur, ce qui équivaut à une absorption moindre à ces longueurs d’onde. Le signal SE aura une forme spectrale similaire à celle du spectre d’émission spontanée de l’espèce, mais avec un signe négatif et une pondération fréquentielle différente10.
En plus des informations sur les espèces à l’état excité, les spectres TA peuvent contenir un certain nombre d’artefacts et de caractéristiques étrangères qui peuvent déformer la dynamique sous-jacente et obscurcir l’attribution des bandes d’absorption11. Un traitement inapproprié de ces artefacts dans la préparation et l’analyse des données peut conduire à l’application de modèles photophysiques inappropriés aux données et, par conséquent, à des conclusions trompeuses11. Par conséquent, la première partie de ce protocole se concentre sur la façon de traiter correctement les ensembles de données d’assistance technique après leur collecte. L’objectif de cette section est de fournir aux chercheurs novices en TA un ensemble de lignes directrices qui les aideront à développer une intuition et une appréciation de la préparation et du traitement rigoureux de leurs données.
Une fois qu’un ensemble de données est traité, une pléthore d’outils et de modèles sont disponibles pour ajuster et interpréter les spectres avec différents niveaux de complexité et de rigueur10. L’objectif de la deuxième section de ce protocole est de préparer le lecteur à appliquer l’ajustement d’une seule longueur d’onde et l’analyse globale aux données et de fournir des conseils sur le moment où ces modèles sont appropriés pour décrire leurs données. Des logiciels commerciaux sont maintenant facilement disponibles pour la préparation et le traitement des données TA, tels que Surface Xplorer12,13 à partir de systèmes ultrarapides (téléchargement et utilisation gratuits, voir Tableau des matériaux). D’autres alternatives gratuites ont été publiées par des chercheurs universitaires, comme Glotaran14. Glotaran est un logiciel libre développé pour l’analyse globale et cible de données de spectroscopie et de microscopie résolues en temps. Il sert d’interface utilisateur graphique (GUI) pour le package R TIMP14. De plus, les utilisateurs peuvent utiliser des langages de programmation comme Python pour écrire leurs propres codes qui effectuent l’analyse. Chacun de ces logiciels et solutions de programmation a des caractéristiques positives qui en font des contributions importantes. Pour les besoins de cette étude, nous ne pouvons présenter qu’un seul logiciel pour le volet visuel de cette activité. Une discussion approfondie de chaque logiciel d’ajustement dépasse le cadre de cet article.
Cet article fournit une procédure étape par étape pour (1) traiter les données TA, (2) ajuster les données TA à l’aide de la cinétique à longueur d’onde unique et de l’analyse globale, et (3) extraire les données et les ajuster à d’autres modèles. Un ensemble de données représentatives sur l’assistance technique est inclus pour que le lecteur puisse l’utiliser en pratique (fichier supplémentaire 1 et fichier supplémentaire 2). Les données sont une mesure d’un échantillon de 165 μM de 1,4-bis(5-phényloxazole-2-yl)benzène (POPOP) dans de l’éthanol excité à 330 nm et prélevé sur une plage de −5 ps à 5,5 ns. De plus, un échantillon « blanc » contenant uniquement de l’éthanol et aucun échantillon n’a été prélevé dans les mêmes conditions expérimentales sur une plage de −5 ps à 5 ps, qui est utilisée pour préparer les données pour l’ajustement (étape 1). Les spectres ont été collectés à l’aide d’un spectromètre d’absorption transitoire ultrarapide. L’échantillon a été contenu dans une cuvette de 2 mm de longueur et soumis à une agitation constante. La procédure de traitement et d’ajustement décrite est basée sur le logiciel Surface Xplorer qui ajuste les données au format *.ufs, et qui sera ci-après appelé le « programme d’ajustement ». Des programmes permettant de convertir des ensembles de données dans d’autres formats en fichiers *.ufs sont disponibles15. Bien que les détails de ce protocole soient spécifiques à Surface Xplorer, les étapes qui suivent sont généralisables à tout logiciel, commercial ou maison. En outre, les résultats du traitement des données peuvent être extraits et ajustés à l’aide de ces autres logiciels. Un dossier d’information à l’appui (dossier supplémentaire 3) fournit des conseils supplémentaires sur l’ajustement.
Considérations générales pour la préparation des données
L’ajustement des données d’AT peut sembler à première vue relativement simple, et on pourrait s’attendre à ce qu’une « réponse » correcte claire en résulte pour un ensemble de données donné. Cependant, comme le souligne le protocole, de nombreux facteurs dans l’acquisition, la préparation et l’analyse des données doivent être soigneusement pris en compte et peuvent entraîner une incertitude quant au modèle ou à l’ensemble de paramètres d’ajustement qui décrit le mieux les données. L’objectif de la préparation et de l’ajustement des données est de réduire autant que possible ces facteurs superflus, tout en préservant les données pour l’analyse. La tâche à accomplir peut sembler intimidante pour un débutant, car il y a beaucoup à considérer. Pour développer l’intuition sur le processus d’ajustement, le débutant est encouragé à essayer de préparer les mêmes données plusieurs fois à partir de zéro de manière légèrement différente pour vérifier l’impact des étapes de préparation des données sur le meilleur ajustement. De plus, deux chercheurs différents peuvent préparer et ajuster les mêmes données et comparer les résultats. Ce processus peut prendre du temps les premières fois, mais cela permettra au débutant de développer une intuition sur la façon de préparer systématiquement les données pour les échantillons futurs. Comme toute compétence, cette préparation et cet ajustement des données prendront du temps à développer, et le débutant est encouragé à être patient et discipliné lors de l’expérimentation et de l’apprentissage du processus. L’ensemble de données utilisé dans cette étude est fourni pour donner au débutant la possibilité de s’adapter directement au tutoriel et de comparer directement les résultats avec ceux produits dans le tutoriel.
Les données peuvent contenir des caractéristiques de fond qui sont présentes à tout moment (figure supplémentaire 2 et figure supplémentaire 3) telles que la diffusion du faisceau de la pompe et l’émission spontanée de l’échantillon. Ces caractéristiques indésirables doivent être supprimées afin d’isoler le signal d’absorption transitoire de l’espèce d’intérêt11. La suppression de ces caractéristiques se fait en choisissant, en faisant la moyenne et en supprimant la contribution d’un certain nombre de spectres de différence de temps négative. Lors de la sélection des spectres de fond, il est important de s’assurer qu’aucune caractéristique qui pourrait faire partie du processus d’intérêt n’est incluse pour la suppression. Des caractéristiques de fond provenant du solvant, telles que l’absorption des impuretés ou du solvant lui-même, peuvent également être observées dans les données d’AT. Lorsque le solvant produit un signal, un ensemble de données « vierges » contenant uniquement le solvant analysé dans les mêmes conditions expérimentales que l’échantillon devra être soustrait de l’ensemble de données de l’échantillon. Les détails de cette procédure sont inclus dans le dossier supplémentaire 3.
La correction du chirp est un autre facteur à considérer attentivement. Le bip se produit lorsque l’impulsion de la sonde se déplace vers l’échantillon et s’élargit en raison d’imperfections dans les rétroviseurs de direction ou en passant à travers des optiques dispersives telles que des lentilles ou des filtres. Le résultat final est que les photons de plus faible énergie dans l’impulsion de la sonde (c’est-à-dire le côté rouge du spectre de la sonde) arrivent à l’échantillon avant les photons de plus haute énergie (c’est-à-dire le côté bleu du spectre de la sonde). Il en résulte que le « temps zéro » des spectres TA est étalé sur plusieurs femtosecondes ou picosecondes18, ce qui se manifeste par une courbe distincte dans l’ensemble de données brutes commençant dans les longueurs d’onde bleues puis s’aplatissant à mesure qu’il s’approche du rouge (Figure supplémentaire 7). Le chirp est plus perceptible sur des échelles de temps plus courtes telles que celles auxquelles accède l’AT ultrarapide. Ce temps zéro dépendant de la longueur d’onde peut être corrigé comme décrit dans le protocole, mais l’application de ce processus peut être délicate et subjective. Le fait de disposer d’un échantillon « vierge » ou d’une mesure de la réponse Kerr du solvant peut minimiser la nature subjective des points de prélèvement manuel pour la correction du chirp nécessaire pour générer l’ajustement polynomial utilisé pour ajuster et corriger le chirp. Le but de la correction de chirp est de supprimer la « courbe » distincte du temps zéro. Il peut falloir plusieurs tentatives d’ajustement du chirp pour obtenir les meilleures données corrigées par chirp. Les données peuvent être ajustées plusieurs fois avec différentes corrections de chirp appliquées afin de comprendre l’impact de la correction de chirp sur les valeurs des courtes durées de vie des TA.
Artefacts qui apparaissent au « temps zéro »
Plusieurs artefacts peuvent être observés près du « temps zéro » dans les données TA, notamment la diffusion de Rayleigh, la diffusion Raman stimulée et la modulation de phase croisée. La diffusion de Rayleigh du faisceau de la pompe est une diffusion élastique qui se produit sans changement d’énergie. Cette fonction apparaîtra à la même longueur d’onde que l’impulsion de la pompe. Une diffusion Raman stimulée peut accompagner le signal de diffusion de la pompe19. La diffusion Raman, qui résulte de la diffusion inélastique d’un photon de pompe, produit des pics à la fois à une énergie plus élevée (anti-Stokes) et plus faible (Stokes) que l’énergie incidente de la pompe. Dans les données TA, une diffusion Raman stimulée est observée en raison de l’irradiation simultanée de l’échantillon avec les faisceaux de la pompe et de la sonde. Lorsque le faisceau de la sonde interagit avec l’échantillon en même temps que le faisceau de la pompe, il stimule le processus Raman. Par conséquent, la diffusion Raman stimulée se produit autour du temps zéro et entraîne des pics supplémentaires dans les spectres dans les premières centaines de femtosecondes (Figure 6, observée dans le spectre bleu plus foncé dans la région mise en évidence et Figure supplémentaire 17). La modulation de phase croisée provient de la modulation de l’indice de réfraction du solvant à partir de l’interaction avec le champ électrique intense d’une impulsion.
La diffusion Raman stimulée peut être distinguée de la modulation de phase croisée car les pics Raman apparaissent à des fréquences spécifiques qui correspondent aux modes vibratoires du solvant. Comme il s’agit d’un processus Raman, on peut observer des lignes de Stokes et d’anti-Stokes de part et d’autre de l’excitation. Les solvants chlorés comme le chlorure de méthylène présentent des bandes Raman très importantes en raison de la grande polarisabilité du chlore. Les signatures spectrales de la modulation de phase croisée sont uniques à un solvant mais ne sont pas aussi faciles à prédire que les caractéristiques de diffusion Raman.
Selon la cinétique de l’échantillon mesuré, la diffusion de Rayleigh, la diffusion Raman et la modulation de phase croisée peuvent chevaucher les premières caractéristiques des données TA et peuvent être difficiles à supprimer des données. En principe, ces caractéristiques peuvent être observées dans une mesure nette du solvant et, soustraites des données, les programmes d’analyse de données peuvent avoir des fonctions d’ajustement pour tenir compte de ces caractéristiques, mais en pratique, cela peut être difficile. Lorsqu’il est trop difficile de soustraire ces artefacts sans compromettre les données de l’échantillon, il peut être préférable de recadrer les spectres compromis vers le temps zéro pour éliminer les artefacts. Cela aura pour effet secondaire malheureux de supprimer les 300 premiers fs de données, mais rendra l’ajustement plus fiable par la suite. Au cours de l’analyse de plusieurs ensembles de données du même échantillon et d’échantillons différents, le débutant acquerra l’intuition pour atteindre cet équilibre entre la soustraction de la surface d’arrière-plan et le recadrage des données initiales de 100 à 200 fs.
Un recadrage général peut être nécessaire pour les parties des spectres qui contiennent un faible rapport signal/bruit. L’instabilité du faisceau de la sonde à certaines régions, la faible intensité de la lumière de la sonde, les concentrations d’échantillon trop élevées (bloquant ainsi une grande partie de la sonde incidente), la faible intensité de la pompe et la section efficace d’absorption de l’échantillon sont des responsables typiques d’un faible rapport signal/bruit qui peut rendre l’ajustement des données difficile. Dans ces cas, le recadrage de l’ensemble de données de chaque côté de la fenêtre optique afin d’atteindre le niveau souhaité de rapport signal/bruit peut faciliter le processus d’ajustement.
Un ensemble de données est prêt à être analysé une fois qu’il a été suffisamment recadré pour éliminer les sections médiocres de l’ensemble de données, que le chirp a été corrigé et que les spectres de fond ont été moyennés et soustraits. Cette procédure doit aboutir à des données qui ne contiennent que les parties les plus pertinentes pour la photophysique et la photochimie d’intérêt. En effet, il est clair qu’il y a un certain degré de subjectivité dans ce processus. L’objectif de la préparation des données est de trouver un équilibre entre la suppression des artefacts afin qu’ils ne perturbent pas l’ajustement, mais pas au point de compromettre l’intégrité de l’ensemble de données, entravant ainsi son interprétation. Trouver cet équilibre prend du temps et de l’expérience pour construire l’intuition de ce qui est un artefact et de ce qui est une donnée. Ajuster (et réajuster) le même ensemble de données à plusieurs jours différents, ou demander à deux chercheurs d’ajuster les mêmes données, peut être un moyen de minimiser l’erreur humaine et la subjectivité de la préparation et de l’analyse des données.
Considérations générales pour l’ajustement et l’interprétation
Une fois les spectres TA bruts traités, ils doivent être interprétés et modélisés pour extraire des informations sur les espèces et la dynamique présente dans le système d’intérêt. Ce processus peut être décrit comme une procédure en trois étapes qui comprend l’interprétation spectrale initiale, la modélisation quantitative et l’attribution de l’interprétation spectrale au modèle ou à l’ajustement.
Interprétation spectrale initiale : Dans l’étape d’interprétation spectrale, l’objectif est d’attribuer des caractéristiques présentes dans les spectres TA à des états électroniques accessibles dans l’évolution photophysique ou photochimique du système. Pour commencer, il convient d’identifier différents États. Dans ce travail, les états font référence à des états électroniques uniques qui font partie de l’évolution photophysique ou photochimique du système. Un état, représenté par exemple par une courbe d’énergie potentielle spécifique (PEC), possède un ensemble de pics caractéristiques représentant son spectre d’absorption. Un changement qui se produit dans un seul état est appelé un processus. Un processus photophysique peut apparaître dans les spectres TA sous la forme d’un décalage de pic ou d’une modification de la largeur du spectre. L’aspect clé d’un processus est que la population de l’État reste la même (c’est-à-dire que le processus se déroule dans un PEC donné) ; C’est la distribution de l’énergie au sein de l’État qui change. Un changement dans la population d’un État sera appelé transition. Au cours d’une transition, le système évolue vers un autre PEC (c’est-à-dire un état électronique). Les transitions peuvent inclure la conversion interne (IC), le croisement intersystème (ISC), le transfert de charge, le transfert d’énergie, la formation de nouveaux produits ou le retour à l’état fondamental. Les directives pour l’affectation d’états, de processus et de transitions sont abordées dans les paragraphes suivants.
Affectation d’états
La première étape de ce processus consiste à attribuer des caractéristiques spectrales à des espèces ou des états chimiques spécifiques. L’état S1 en TA doit montrer une durée de vie qui correspond à la durée de vie de la fluorescence prise à l’aide de la spectroscopie d’émission résolue en temps. Un état triplet peut être vérifié si sa durée de vie est éteinte par l’oxygène. Si un anion ou un cation radical est suspecté dans l’évolution photophysique, une spectroélectrochimie ou une oxydation/réduction chimique peut être effectuée pour générer l’espèce radicalaire, et un spectre d’absorption de cette espèce peut être obtenu et comparé à la forme de bande TA. La spectroscopie par résonance de spin électronique (ESR) peut être réalisée pour vérifier la présence de radicaux libres. Un excellent tutoriel organisé par la division de chimie inorganique de l’ACS donne un aperçu de l’AT et de ces considérations dans l’attribution des caractéristiques20. Une fois que les bandes ont été attribuées aux espèces, l’étape suivante de l’interprétation des spectres TA consiste à décrire qualitativement les processus dynamiques qui se produisent dans le système. Cette étape est essentielle car elle donne au chercheur une idée des modèles appropriés pour décrire son système et lui donnera une base de référence à laquelle comparer les paramètres d’ajustement.
Changements au sein d’un état
Le refroidissement vibratoire, le réarrangement géométrique ou la solvatation sont des processus extrêmement rapides (inférieurs à 10 ps) qui peuvent être observés avec l’AT. Le refroidissement vibratoire est observé comme un décalage rapide vers le bleu du spectre TA sur une échelle de temps de plusieurs picosecondes 21,22,23. Un réarrangement géométrique peut se produire sur l’échelle de temps de 10 ps. La dynamique de solvatation est observée sous la forme d’un décalage vers le rouge et d’un rétrécissement du spectre sur plusieurs picosecondes dans les liquides dipolaires conventionnels, mais les solvants à haute viscosité tels que le glycérol, le polyéthylène glycol (PEG), les liquides ioniques et les solvants eutectiques profonds peuvent présenter une dynamique de solvatation se produisant sur plusieurs nanosecondes 24,25,26.
Changements dans la population d’un État
Les réactions sont caractérisées par un changement d’intensité d’une bande, où une diminution d’intensité est associée à une diminution de la concentration de ses espèces chimiques et inversement pour une augmentation. Dans certains cas, les espèces de réactif et de produit sont visibles dans les spectres, tandis que dans d’autres, les états du produit sont trop courts ou trop décalés vers le rouge pour être observés. Souvent, les transitions d’état à état peuvent être observées par la présence d’un point isobestique dans les spectres.
Modélisation/ajustement quantitatif : Un modèle doit ensuite être ajusté aux données afin d’extraire des informations quantitatives sur la dynamique du système. Comme décrit précédemment dans l’introduction, il existe une vaste gamme de modèles à utiliser. Ce protocole se concentre sur deux des méthodes les plus courantes : l’ajustement d’une seule longueur d’onde et l’analyse globale. La méthode à longueur d’onde unique consiste à ajuster les traces de longueur d’onde individuelles des spectres à une forme fonctionnelle, généralement une somme d’exponentielles :
(2)
où ΔA(t) est le signal TA à une longueur d’onde choisie, n est le nombre de composantes exponentielles, et aiest l’amplitude de la composante exponentielle, i, avec une constante de temps τi. Plusieurs composants peuvent être ajoutés jusqu’à ce que l’ajustement reproduise les données expérimentales. L’objectif de tout processus d’ajustement est de modéliser les données en utilisant suffisamment de durées de vie pour bien reproduire les données, mais pas de surajuster les données en incluant trop de composants. Par conséquent, des paramètres pondérés de qualité de l’ajustement tels que , sont utilisés pour aider à déterminer quand les données sont ajustées dans les incertitudes expérimentales5.
Une fois la désintégration ajustée de manière satisfaisante, les paramètres du modèle peuvent être utilisés pour caractériser la dynamique du système. Les constantes de temps résultantes peuvent ensuite être extraites et interprétées. Malheureusement, le grand nombre de caractéristiques qui se chevauchent dans les spectres TA signifie qu’une seule longueur d’onde dans le spectre peut contenir des dynamiques correspondant à différentes espèces dont les signatures spectrales se chevauchent, ce qui signifie que les constantes de temps extraites d’un ajustement de longueur d’onde unique peuvent représenter un composite de plusieurs processus coïncidants. De plus, tout changement de forme et de position de la bande influencera également les amplitudes et les constantes de temps extraites de l’ajustement d’une seule longueur d’onde. Ces problèmes peuvent être contournés dans certains cas par une méthode d’ajustement appelée « analyse de la forme de bande », où l’on détermine ou suppose une forme fonctionnelle pour les bandes TA de chaque espèce absorbante dans le système. Ces formes sont ensuite pondérées par des amplitudes dépendantes du temps et additionnées afin de reproduire le spectre observé. Cette procédure est couramment utilisée dans l’analyse des spectres de fluorescence résolus en temps, mais les formes plus compliquées et les composantes qui se chevauchent des bandes TA rendent cette méthode défendable dans quelques cas simples seulement, comme détaillé ailleurs10.
Un autre inconvénient de l’ajustement à longueur d’onde unique est qu’il ne tire pas intrinsèquement parti de la large gamme spectrale offerte par les expériences TA modernes. On pourrait, en principe, ajuster méthodiquement chaque longueur d’onde individuelle des spectres, mais une telle analyse est lourde, longue et coûteuse en calcul. Pour relever ce défi, une méthode appelée « analyse globale » peut être utilisée pour ajuster simultanément un ensemble complet de spectres TA à un ensemble de paramètres dynamiques partagés4. L’analyse globale, et une méthode étroitement liée appelée analyse cible, sont des méthodes efficaces et largement utilisées, mais elles présentent également leur propre ensemble d’inconvénients et de limites. Comme pour tout modèle, il est impératif de comprendre les hypothèses qui sont utilisées pour le créer ainsi que les limites qu’elles présentent.
Dans l’analyse globale, les spectres TA sont représentés par une matrice m par n, où m représente le nombre de longueurs d’onde mesurées dans chaque spectre et n représente le nombre de points temporels collectés. Cette matrice est alors supposée décomposable en produit de deux autres matrices :
(3)
où C(t) est une matrice n par k et S(λ) est une matrice m par k. La valeur k représente le nombre de composantes spectrales distinctes utilisées pour reproduire les spectres. Chacun de ces composants représente une espèce absorbante avec une signature spectrale et une dynamique uniques. La matrice S(λ) représente les spectres TA des k composantes et C(t) leurs concentrations en fonction du temps. Dans l’implémentation la plus simple et la plus courante de l’analyse globale, chaque composante est supposée avoir une cinétique exponentielle unique (i = 1 dans l’équation 2, chaque composante se voyant attribuer sa propre constante de temps). En résumé, le spectre TA complet peut être représenté par la somme des composantes de k spectres, chacune ayant son propre spectre d’absorption caractéristique et sa décroissance exponentielle unique.
Lorsque les spectres TA sont ajustés, l’utilisateur devine combien de composantes (c’est-à-dire une valeur pour k) sont nécessaires et devine la constante de temps associée à une décroissance exponentielle unique de ces espèces. L’ajusteur génère ensuite Cguess(t) et résout l’équation 3 pour Sfit(t). Ensuite, Sfit(λ) et Cguess(t) sont multipliés comme dans l’équation 3 pour créer les spectres ajustés, ΔA(λ,t)fit. Enfin, les résidus, ΔA(λ,t)exp − A(λ,t)fit, sont minimisés et les constantes optimalesd’ajustement S (λ) et de temps sont renvoyées. La relative simplicité de l’analyse globale, représentant un ensemble complet de spectres à l’aide d’une poignée de constantes de temps et de composantes spectrales fixes, en fait une méthode attrayante (et efficace) pour démêler les formes de bande et la dynamique complexes rencontrées en spectroscopie TA. Cependant, il faut veiller à ce que l’analyse globale soit un modèle approprié pour le système en question.
Une hypothèse clé de l’analyse globale, illustrée dans l’équation 3, est la séparabilité complète des parties longueur d’onde et temps de la dynamique, une propriété appelée « bilinéarité ». Cette hypothèse exige que les formes de bande des composantes soient indépendantes du temps (c’est-à-dire qu’elles ont une forme spectrale fixe qui ne varie pas ou ne se déplace pas avec le temps). La seule chose qui change au cours de l’expérience est la population relative de chaque composant, représentée par C(t). Sur de longues échelles de temps, ~1 ns environ, cette hypothèse est généralement valable et l’analyse globale peut être utilisée sans trop de soucis. D’autre part, les processus à l’état excité tels que le refroidissement vibrationnel et la dynamique de solvatation, importants sur les échelles de temps ultrarapides accessibles par l’AT femtoseconde, entraînent des changements dans la signature spectrale d’une espèce en fonction du temps et une rupture de la bilinéarité. Cela ne signifie pas que l’analyse globale ne peut pas reproduire un ensemble de données, en fait, elle peut toujours produire un ajustement satisfaisant à condition qu’un nombre suffisant de composants soient utilisés. Le problème réside alors dans l’interprétation des spectres des composantes et l’attribution des constantes de temps à des processus particuliers à l’état excité, car les composantes peuvent ne plus correspondre à des espèces absorbantes distinctes. Par conséquent, il faut toujours faire preuve de prudence lorsqu’on applique une analyse globale à des situations où la bilinéarité ne peut être supposée.
Attribution de l’interprétation spectrale au modèle/ajustement : Une fois qu’un ajustement est obtenu, l’interprétation spectrale doit être mise en correspondance avec les durées de vie obtenues dans l’ajustement. Les durées de vie à partir de l’ajustement sont attribuées aux processus et aux réactions qui ont été identifiés lors de l’interprétation initiale des spectres. Cependant, l’évaluation initiale à partir des spectres et le nombre de durées de vie ajustées obtenues par le modèle peuvent ne pas correspondre immédiatement les uns aux autres. Dans cette situation (courante !), l’ajusteur doit revenir en arrière et évaluer l’interprétation initiale. Peut-être y avait-il un refroidissement vibratoire ou un autre processus qui a été omis dans l’évaluation initiale, mais qui a été identifié dans le processus de modélisation et d’ajustement. Ou peut-être que deux ensembles différents de paramètres d’ajustement pourraient bien reproduire les données et que l’interprétation initiale pourrait guider l’ensemble de paramètres d’ajustement choisi. Dans cette dernière étape, l’ajusteur doit faire des allers-retours entre l’interprétation et l’ajustement pour trouver une description qui mène à une attribution photophysique plausible de l’espèce et de la dynamique du système. D’autres programmes d’ajustement qui incluent des modèles d’ajustement séquentiel, tels que l’analyse des cibles, peuvent également être explorés pour compléter les ajustements produits par l’analyse globale et le logiciel d’ajustement présenté dans cet article4.
En résumé, ce protocole traite de la préparation et de l’ajustement des données d’absorption transitoire. Son objectif est de mettre en évidence les défis associés au processus et de commenter les moyens d’éviter ou d’atténuer ces défis de manière pratique. L’ajustement des données d’AT, comme celui de la plupart des données rencontrées dans les domaines techniques, peut être délicat et, parfois, subjectif. Par conséquent, il est essentiel d’être conscient du processus et des limites des données, de la préparation des données et des outils mathématiques utilisés pour modéliser et attribuer un sens aux données. Les scientifiques doivent aborder les données et la modélisation avec un œil critique.
On peut tenter d’atténuer la subjectivité de leurs crises. Par exemple, les données peuvent être préparées et ajustées à partir de différents points de départ et à différents jours pour s’assurer que le même ajustement est produit. Les données prises à des jours différents avec différentes préparations d’échantillons peuvent être comparées. Plusieurs chercheurs peuvent ajuster les mêmes données et comparer leurs résultats. Au fil du temps, les chercheurs peuvent se forger une intuition sur les données qu’ils obtiennent (en fonction des spécificités de leur dispositif expérimental et de leurs paramètres expérimentaux) qui leur permettra d’être plus confiants dans leurs ajustements.
Il y a beaucoup à apprendre sur l’ajustement des données TA et les détails des modèles abordés dans cet article. Plusieurs excellents articles de synthèse sont recommandés avec enthousiasme qui approfondissent ce sujet 4,10,27. Ce protocole est destiné à être la porte d’entrée d’un débutant dans le processus d’analyse et d’ajustement qui suscite l’intérêt pour une compréhension plus approfondie du processus.
The authors have nothing to disclose.
Ce travail a été rendu possible grâce au programme d’instrumentation de recherche majeure de la NSF qui a établi l’installation laser multi-utilisateurs pour l’absorption transitoire (CHE-1428633). Ce matériel est basé sur des travaux soutenus par la National Science Foundation dans le cadre de la subvention No. CHE-2313290.
EtOH 200% Proof | Decon Laboratories Inc | CAS 64-17-5 | Solvent used to prepare Sample |
Helios transient absorption spectrometer | Ultrafast systems | https://ultrafast.systems/products/spectrometers-accessories/helios/ | Transient absorption spectrometer |
POPOP 1,4-Bis[2-(5-phenyloxazolyl)]benzene | Tokyo Chemical Industry | CAS 1806-34-4 | Sample used for Examples |
Surface Xplorer | Ultrafast systems | https://ultrafast.systems/products/spectrometers-accessories/surface-xplorer/ | Fitting program |