Noncovalent interactions in complex chemical systems are examined by considering model systems which capture the essential physics of the interactions and applying correlated electronic structure techniques to these systems. Noncovalent interactions are critical to understanding a host of energetic and structural properties in complex chemical systems, from base pair stacking in DNA to protein folding in organic solids. Complex chemical and biophysical systems, such as enzymes and proteins, are too large to be studied using computational techniques rigorous enough to capture the subtleties of noncovalent interactions. Thus, the larger chemical system must be truncated to a smaller model system to which rigorous methods can be applied in order to capture the essential physics of the interaction. Computational methodologies which can account for high levels of electron correlation, such as second-order perturbation theory and coupled-cluster theory, must be used. These computational techniques will be used to study several types (pi stacking, S/pi, and C-H/pi) of noncovalent interactions in two chemical contexts: biophysical systems and organic solids.

Disertacijoje nagrinėjamas galimas aplinkos poveikis organinių molekulių elektroninių būsenų savybėms. Tam tikslui yra naudojami kompiuterizuotieji kvantinės mechanikos metodai, kuriais remiantis nagrinėjamos įvairių molekulių savybės. Ištirtos 2-(N-metil-α-iminoethyl)-fenol ir N-triphenylmethylsalicylidene imine molekelulių, esančių poliniame tirpiklyje, struktūros pagrindinėje ir sužadintose elektroninėse būsenose. Pirmą kartą parodyta, kad, norint gauti teisingą kokybinį ir artimą kiekybiniam vidujmolekulinės protono pernašos potencinės energijos paviršių, būtina atsižvelgti į polinių tirpiklio molekulių kuriamą vandenilinių ryšių tinklą bei į nulinių svyravimų energijas. Pastarieji ir nulemia protono pernašos vyksmo kryptį bei efektyvumą. Parodyta, kad anilų klasės molekulių konformerų susiformavimas priklauso nuo tirpiklio poliškumo, o jų susidarymas savo ruožtu konkuruoja su klasterių iš tirpiklio molekulių susiformavimo galimybėmis. Pirmą kartą parodyta, kad dipolinio momento vertė bakteriorodopsine yra nulemta membranos paviršiuose esančių radikalų. Pademonstruota, kad stilbazolio molekulė deformuojasi ir sudaro naujus konformerus (pademonstruota dviejų formų atsiradimo galimybė) tik esant molekulėms tirpalo apsuptyje. Šis rezultatas paaiškino eksperimente stebimus skirtuminių spektrinių pokyčių evoliucijos prigimtį. Apskaičiuotos ir išanalizuotos karotinoidų - luteino, violaksantino ir zeaksantino molekulių - žemiausios sužadintos elektroninės būsenos. Parodyta... [toliau žr. visą tekstą] / To explore changes caused by the environment on the internal characteristics of an organic molecule is the objective of the thesis. For this purpose we investigate a variety of organic molecules. Using various methods of quantum mechanics calculations possible influence of a polar solvent on the ground and excited states of 2-(N-metil-α-iminoethyl)-fenol and N-triphenylmethylsalicylidene imine is considered. It is shown for the first time that in order to obtain the correct qualitative and quantitative interpretation of possible pathways of the intermolecular proton transfer the hydrogen network of the polar solvent molecules together with the zero point energy have to be taken into consideration. It is also shown that conformational variability of anil-type molecules in polar solvents is competing with clusters formation of solvent molecules. It is shown for the first time that the dipole moment of bacteriorhodopsin is mainly defined by cytoplasmic and extracellular coils on the surfaces of the membrane. It is also demonstrated that the stilbazole molecule experiences the deformation resulting in formation of new conformers (at least two forms are present) in the solvent surrounding. The experimental data of the transient spectroscopy were explained in the basis of these model calculations. The lowest excited states of carotinoids, such as lutein, zeaxanthin and violoxantin are calculated and analyzed. Sensitivity of the excited electronic state on the polar environment is... [to full text]

Elektroninio sužadinimo evoliucija šviesai jautriose molekulėse yra reiškinys, kuriuo remiantis įmanoma nagrinėti daugelį natūralių ir dirbtinių procesų: augalų ir bakterijų fotosintezę, regos mechanizmą, optomechaninių bei optoelektroninių prietaisų (pavyzdžiui, organinių šviestukų) veikimą. Teoriškai šis reiškinys modeliuojamas sprendžiant laikinę Šriodingerio lygtį. Deja, toks sprendimas realiems, praktiškai panaudojamiems junginiams šiandien yra per sudėtingas uždavinys, todėl jį tenka keisti supaprastinant nagrinėjamų junginių modelius arba sprendimo metodiką. Šioje disertacijoje aprašomų tyrimų tikslas buvo elektroninės struktūros skaičiavimų metodais (t. y. sprendžiant paprastesnę nuostoviąją Šriodingerio lygtį) ištirti elektroninio sužadinimo sukeltus procesus fotoaktyviose molekulėse ir sudaryti sužadinimo relaksaciją apibūdinančius potencinės energijos paviršių modelius. Parodoma, jog ta pačia metodika atliekamų tyrimų rezultatai paaiškina įvairiuose junginiuose vykstančius reiškinius: bakteriorodopsino baltymo funkcinės grupės vykdomą protono pernašą poliniame tirpiklyje, indolo-benzoksazino junginio optomechaninį ciklą, našią fosforescenciją organiniame silicio polimere bei šviestukams naudojamo metaloorganinio komplekso su prijungtomis krūvininkų pernašos grupėmis ypatybes. / Evolution of the electronic excitation is a general process that can be used to explain many natural and artificial phenomena, such as photosynthesis in plants and bacteria, biological mechanism of vision, and operating principles of optomechanical and optoelectronic devices. This process is theoretically modeled by solving the time-dependent Schroedinger equation. However, such treatment is too computationally expensive to be used for practical molecular systems. Therefore, either models of the structure of the systems or the solving procedure itself must be simplified to get the desired results. The main goal of the research presented in this dissertation was to study processes caused by the electronic excitation in photoactive molecules using computational methods of electronic structure (i. e. solving the simpler time-independent Schroedinger equation) and to construct the potential energy surface models describing the energy relaxation in the investigated molecules. It is shown that the results of different investigations performed using the same procedure provide explanations of different phenomena in various compounds, such as: proton transfer in polar solvent, performed by a functional group of the bacteriorhodopsin protein; optomechanical cycle of the indolo-benzoxazine compound; efficient phosphorescence of the silicon-based organic polymer; and optical properties of organometallic emitter compound with additional charge-carrier groups.

Evolution of the electronic excitation is a general process that can be used to explain many natural and artificial phenomena, such as photosynthesis in plants and bacteria, biological mechanism of vision, and operating principles of optomechanical and optoelectronic devices. This process is theoretically modeled by solving the time-dependent Schroedinger equation. However, such treatment is too computationally expensive to be used for practical molecular systems. Therefore, either models of the structure of the systems or the solving procedure itself must be simplified to get the desired results. The main goal of the research presented in this dissertation was to study processes caused by the electronic excitation in photoactive molecules using computational methods of electronic structure (i. e. solving the simpler time-independent Schroedinger equation) and to construct the potential energy surface models describing the energy relaxation in the investigated molecules. It is shown that the results of different investigations performed using the same procedure provide explanations of different phenomena in various compounds, such as: proton transfer in polar solvent, performed by a functional group of the bacteriorhodopsin protein; optomechanical cycle of the indolo-benzoxazine compound; efficient phosphorescence of the silicon-based organic polymer; and optical properties of organometallic emitter compound with additional charge-carrier groups. / Elektroninio sužadinimo evoliucija šviesai jautriose molekulėse yra reiškinys, kuriuo remiantis įmanoma nagrinėti daugelį natūralių ir dirbtinių procesų: augalų ir bakterijų fotosintezę, regos mechanizmą, optomechaninių bei optoelektroninių prietaisų (pavyzdžiui, organinių šviestukų) veikimą. Teoriškai šis reiškinys modeliuojamas sprendžiant laikinę Šriodingerio lygtį. Deja, toks sprendimas realiems, praktiškai panaudojamiems junginiams šiandien yra per sudėtingas uždavinys, todėl jį tenka keisti supaprastinant nagrinėjamų junginių modelius arba sprendimo metodiką. Šioje disertacijoje aprašomų tyrimų tikslas buvo elektroninės struktūros skaičiavimų metodais (t. y. sprendžiant paprastesnę nuostoviąją Šriodingerio lygtį) ištirti elektroninio sužadinimo sukeltus procesus fotoaktyviose molekulėse ir sudaryti sužadinimo relaksaciją apibūdinančius potencinės energijos paviršių modelius. Parodoma, jog ta pačia metodika atliekamų tyrimų rezultatai paaiškina įvairiuose junginiuose vykstančius reiškinius: bakteriorodopsino baltymo funkcinės grupės vykdomą protono pernašą poliniame tirpiklyje, indolo-benzoksazino junginio optomechaninį ciklą, našią fosforescenciją organiniame silicio polimere bei šviestukams naudojamo metaloorganinio komplekso su prijungtomis krūvininkų pernašos grupėmis ypatybes.

Este projeto tem dois objetivos. Primeiramente estudou-se a aplicabilidade dos splines tricúbicos para a construção de superfícies de energia potencial globais. Um dos obstáculos que este método tem de superar e a escolha de um sistema de coordenadas apropriado, que minimize a influência de pontos não físicos. Para isto, propôs-se o uso do sistema de coordenadas de Pekeris, nunca usado para este fim. Este procedimento foi realizado para três sistemas químicos bem descritos na literatura, [Cl,H2], [F,H,D] e [H,O,Cl], cujas superfícies de energia potencial e propriedades das reações foram usadas como referência. Com base nestes modelos, aplicamos o método proposto variando-se a quantidade e a disposição dos nós das interpolações, a fim de verificar sua influência na qualidade das superfícies interpoladas. Os resultados mostram que as superfícies construídas por este método reproduzem muito bem os cálculos de dinâmica química, tanto por métodos quânticos quanto por métodos clássicos. Para isto, os nós da interpolação devem cobrir as regiões mais importantes da superfície de energia potencial e os valores mais baixos das coordenadas de Pekeris devem ser priorizados. O segundo objetivo consiste na aplicação deste procedimento na construção da superfície de energia potencial [H,S,F]. Com esta superfície, diversas características deste sistema foram analisadas, tais como geometrias dos pontos estacionários, energias relativas e frequências vibracionais. Os valores obtidos estão de acordo com os dados descritos na literatura. A superfície construída também foi usada para a realização de cálculos de dinâmica para a reação F+HS → S+FH. Observamos a existência de dois tipos de mecanismos, um com a formação de um intermediário de longa duração e outro com a abstração direta do átomo de hidrogênio. / This project has two goals. First, we studied the applicability of the tricubic splines to construct global potential energy surfaces. One of the diculties this approach has to overcome is the choice of an appropriate coordinate system that minimises the in uence of non-physical points. For such, we proposed the use of the Pekeris coordinate system, never employed for this purpose. This procedure was carried out for three well described systems, [Cl,H2], [F,H,D] and [H,O,Cl], whose potential energy surfaces and reaction properties were taken as references. Based on these models, we applied the proposed method varying the amount and arrangement of the interpolation knots, to verify their influence on the quality of the interpolated surfaces. The results showed that surfaces constructed by this approach reproduce very well the chemical dynamics calculations, both for the quantum as well as for the classical methods, provided that the interpolation knots cover the most important regions of the potential energy surfaces, and the lower values of the Pekeris coordinates are prioritised. The second goal was the application of this procedure to the construction of the [H,S,F] potential energy surface. With this surface, several characteristics of this system were analysed, such as the geometry of the stationary points, relative energies and vibrational frequencies. The values obtained are in agreement with the data described in the literature. The constructed surface was also used for quantum dynamics calculations on the reaction F + HS → S + FH. We observed two kinds of mechanisms, one of them with the formation of a long-living intermediate and the other with the direct abstraction of the hydrogen atom.

L’étude des mécanismes de formation et de destruction des anions moléculaires est devenu un champ d’intérêt prononcé après la détection récente de six anions moléculaires (C4H-, C6H-, C8H-, CN-, C3N-, C5N-) dans le milieu interstellaire.Dans les environnements interstellaires où la densité d’électrons est relativement importante, le canal principal de formation de ces anions devrait être l’attachement électronique radiatif. Mais il manque aujourd’hui des données expérimentales et théoriques permettant d’évaluer cette hypothèse. D’autre part, le photodétachement est la principale cause de destruction de ces anions dans les nuages diffus et les régions de photodissociation. Une approche basée sur un développement monocentrique est appliquée à l’étude de ces deux processus opposés que sont le photodétachement et l’attachement électronique radiatif. Les résultats obtenus avec la présente méthode sont comparés à des données expérimentales et théoriques précédemment rapportées et montrent un bon accord. Cette méthode est ensuite utilisée pour déterminer les constantes de vitesse nécessaires pour confirmer si ces mécanismes sont cruciaux pour la chimie d’anions interstellaires. En plus des constantes de vitesse de formation et de destruction des anions, les constantes de vitesse d’excitation collisionnelle sont nécessaires pour modéliser les abondances observées des anions. Nous avons choisi de porter notre effort sur le calcul des constantes de vitesse de transition entre états rotationels de la molécule C3N- dans son état vibrationnel fondamental lors des collisions avec H2 et He en utilisant de nouvelles surfaces d’énergie potentielles. / The mechanisms of formation and destruction of molecular anions have become a field of special interest after the recent detection of six molecular anions (C4H-,C6H-, C8H-, CN-, C3N-, C5N-) in the interstellar medium. The main channel of formation of these anions is expected to be radiative electron attachment in environments where the density of electron is relatively important. There is however at themoment a lack of experimental and theoretical data allowing to assess this hypothesis. Photodetachment, on the other hand, is the main source of destruction of the anions in diffuse clouds and photodissociation regions. A single center expansion approach is applied to the study of both processes: photodetachment and radiative electron attachment. The results obtained with the present method are compared to previously reported experimental and theoretical data and show a good agreement.This method is then employed to determine the rate constants which are needed to confirm whether or not these mechanisms are crucial for the chemistry of the interstellar anions. Along with the formation and destruction rates, rotational excitation rate coefficients are needed to accurately model the observed anions abundances.We focus on the calculation of state-to-state rotational transitions rate coefficients of the C3N- molecule in its ground vibrational state in collisions with H2 and He using new potential energy surfaces

Disertacijoje nagrinėjamas galimas aplinkos poveikis organinių molekulių elektroninių būsenų savybėms. Tam tikslui yra naudojami kompiuterizuotieji kvantinės mechanikos metodai, kuriais remiantis nagrinėjamos įvairių molekulių savybės. Ištirtos 2-(N-metil-α-iminoethyl)-fenol ir N-triphenylmethylsalicylidene imine molekelulių, esančių poliniame tirpiklyje, struktūros pagrindinėje ir sužadintose elektroninėse būsenose. Pirmą kartą parodyta, kad, norint gauti teisingą kokybinį ir artimą kiekybiniam vidujmolekulinės protono pernašos potencinės energijos paviršių, būtina atsižvelgti į polinių tirpiklio molekulių kuriamą vandenilinių ryšių tinklą bei į nulinių svyravimų energijas. Pastarieji ir nulemia protono pernašos vyksmo kryptį bei efektyvumą. Parodyta, kad anilų klasės molekulių konformerų susiformavimas priklauso nuo tirpiklio poliškumo, o jų susidarymas savo ruožtu konkuruoja su klasterių iš tirpiklio molekulių susiformavimo galimybėmis. Pirmą kartą parodyta, kad dipolinio momento vertė bakteriorodopsine yra nulemta membranos paviršiuose esančių radikalų. Pademonstruota, kad stilbazolio molekulė deformuojasi ir sudaro naujus konformerus (pademonstruota dviejų formų atsiradimo galimybė) tik esant molekulėms tirpalo apsuptyje. Šis rezultatas paaiškino eksperimente stebimus skirtuminių spektrinių pokyčių evoliucijos prigimtį. Apskaičiuotos ir išanalizuotos karotinoidų - luteino, violaksantino ir zeaksantino molekulių - žemiausios sužadintos elektroninės būsenos. Parodyta... [toliau žr. visą tekstą] / To explore changes caused by the environment on the internal characteristics of an organic molecule is the objective of the thesis. For this purpose we investigate a variety of organic molecules. Using various methods of quantum mechanics calculations possible influence of a polar solvent on the ground and excited states of 2-(N-metil-α-iminoethyl)-fenol and N-triphenylmethylsalicylidene imine is considered. It is shown for the first time that in order to obtain the correct qualitative and quantitative interpretation of possible pathways of the intermolecular proton transfer the hydrogen network of the polar solvent molecules together with the zero point energy have to be taken into consideration. It is also shown that conformational variability of anil-type molecules in polar solvents is competing with clusters formation of solvent molecules. It is shown for the first time that the dipole moment of bacteriorhodopsin is mainly defined by cytoplasmic and extracellular coils on the surfaces of the membrane. It is also demonstrated that the stilbazole molecule experiences the deformation resulting in formation of new conformers (at least two forms are present) in the solvent surrounding. The experimental data of the transient spectroscopy were explained in the basis of these model calculations. The lowest excited states of carotinoids, such as lutein, zeaxanthin and violoxantin are calculated and analyzed. Sensitivity of the excited electronic state on the polar environment is... [to full text]

Disertacijoje nagrinėjamas galimas aplinkos poveikis organinių molekulių elektroninių būsenų savybėms. Tam tikslui yra naudojami kompiuterizuotieji kvantinės mechanikos metodai, kuriais remiantis nagrinėjamos įvairių molekulių savybės. Ištirtos 2-(N-metil-α-iminoethyl)-fenol ir N-triphenylmethylsalicylidene imine molekelulių, esančių poliniame tirpiklyje, struktūros pagrindinėje ir sužadintose elektroninėse būsenose. Pirmą kartą parodyta, kad, norint gauti teisingą kokybinį ir artimą kiekybiniam vidujmolekulinės protono pernašos potencinės energijos paviršių, būtina atsižvelgti į polinių tirpiklio molekulių kuriamą vandenilinių ryšių tinklą bei į nulinių svyravimų energijas. Pastarieji ir nulemia protono pernašos vyksmo kryptį bei efektyvumą. Parodyta, kad anilų klasės molekulių konformerų susiformavimas priklauso nuo tirpiklio poliškumo, o jų susidarymas savo ruožtu konkuruoja su klasterių iš tirpiklio molekulių susiformavimo galimybėmis. Pirmą kartą parodyta, kad dipolinio momento vertė bakteriorodopsine yra nulemta membranos paviršiuose esančių radikalų. Pademonstruota, kad stilbazolio molekulė deformuojasi ir sudaro naujus konformerus (pademonstruota dviejų formų atsiradimo galimybė) tik esant molekulėms tirpalo apsuptyje. Šis rezultatas paaiškino eksperimente stebimus skirtuminių spektrinių pokyčių evoliucijos prigimtį. Apskaičiuotos ir išanalizuotos karotinoidų - luteino, violaksantino ir zeaksantino molekulių - žemiausios sužadintos elektroninės būsenos. Parodyta... [toliau žr. visą tekstą] / To explore changes caused by the environment on the internal characteristics of an organic molecule is the objective of the thesis. For this purpose we investigate a variety of organic molecules. Using various methods of quantum mechanics calculations possible influence of a polar solvent on the ground and excited states of 2-(N-metil-α-iminoethyl)-fenol and N-triphenylmethylsalicylidene imine is considered. It is shown for the first time that in order to obtain the correct qualitative and quantitative interpretation of possible pathways of the intermolecular proton transfer the hydrogen network of the polar solvent molecules together with the zero point energy have to be taken into consideration. It is also shown that conformational variability of anil-type molecules in polar solvents is competing with clusters formation of solvent molecules. It is shown for the first time that the dipole moment of bacteriorhodopsin is mainly defined by cytoplasmic and extracellular coils on the surfaces of the membrane. It is also demonstrated that the stilbazole molecule experiences the deformation resulting in formation of new conformers (at least two forms are present) in the solvent surrounding. The experimental data of the transient spectroscopy were explained in the basis of these model calculations. The lowest excited states of carotinoids, such as lutein, zeaxanthin and violoxantin are calculated and analyzed. Sensitivity of the excited electronic state on the polar environment is... [to full text]

Thermodynamic stabilities of eighteen cyclic C5H6 isomers were explored computationally both on singlet and triplet state potential energy surfaces (PES). All isomers have singlet ground states except for bicyclo[2.1.0]pent-5-ylidene (B5) having no stable geometry on the singlet C5H6 PES. Cyclopenta-1,3-diene (M1) is the most stable cyclic C5H6 isomer while cyclopent-1,4-diylidene is the least stable one among all. Cyclopenta-1,2-diene (M2) and cyclopentyne (M3) have biradical characters of 46.9 and 21.5%, respectively. Seven unimolecular isomerization reactions occurring among several of these molecules were investigated by DFT and ab initio methods. The conversion of bicyclo[2.1.0]pent-2-ene (B1) and tricyclo[2.1.0.02,5]-pentane (T1) into 1,3-cyclopentadiene (M1) are shown to be concerted processes whose reaction paths pass through TSs with a high degree of biradical character. The reaction enthalpies (DH0) are predicted to be -47.7 kcal/mol for B1 and -63.8 kcal/mol for T1 at UB3LYP/6-31G(d) level. The activation enthalpy (DH0&sup1 / ) for the ring opening of B1 was calculated by the CR-CCSD(T) method to be 25.2 kcal/mol, in good agreement with experiment. Furthermore, the DH0&sup1 / for the ring opening of T1 was obtained by the CR-CCSD(T) method to be 48.2 kcal/mol. The self-conversion of M1 via 1,5-hydrogen shift is a facile and concerted reaction with aromatic TS. The DH0&sup1 / estimations of B3LYP and CC methods are 25.24 and 28.78 kcal/mol, respectively. For 1,2-hydrogen shift reactions of cyclopent-3-enylidene (M4) and cyclopenten-2-ylidene (M5), the single point CC calculations predicted the DH0&sup1 / values of 3.13 and 10.12 kcal/mol, as well as, the DH0 values of -71.28 and -64.05 kcal/mol, respectively. The reason of M5 being more stable than M4 is due to the conjugation of the carbene carbon and the double bond in M5. The reaction path of cyclobutylidene methylene to cyclopentyne rearrangement is found to be rather shallow. The DH0&sup1 / and DH0 values predicted by the RCCSD(T) method to be 3.65 and -5.72 kcal/mol, respectively. Finally, triplet state isomerization of bicyclo[2.1.0]pent-5-ylidene to cyclopenta-1,2-diene, as well as, its parent reaction, cyclopropylidene to 1,2-propadiene were investigated at several levels of theory including DFT, CASSCF and CC methods. The UCCSD(T) method estimated a moderate barrier whose value is 8.12 kcal/mol for the isomerization of 3B5 with the reaction enthalpy of -44.63 kcal/mol.

Este projeto tem dois objetivos. Primeiramente estudou-se a aplicabilidade dos splines tricúbicos para a construção de superfícies de energia potencial globais. Um dos obstáculos que este método tem de superar e a escolha de um sistema de coordenadas apropriado, que minimize a influência de pontos não físicos. Para isto, propôs-se o uso do sistema de coordenadas de Pekeris, nunca usado para este fim. Este procedimento foi realizado para três sistemas químicos bem descritos na literatura, [Cl,H2], [F,H,D] e [H,O,Cl], cujas superfícies de energia potencial e propriedades das reações foram usadas como referência. Com base nestes modelos, aplicamos o método proposto variando-se a quantidade e a disposição dos nós das interpolações, a fim de verificar sua influência na qualidade das superfícies interpoladas. Os resultados mostram que as superfícies construídas por este método reproduzem muito bem os cálculos de dinâmica química, tanto por métodos quânticos quanto por métodos clássicos. Para isto, os nós da interpolação devem cobrir as regiões mais importantes da superfície de energia potencial e os valores mais baixos das coordenadas de Pekeris devem ser priorizados. O segundo objetivo consiste na aplicação deste procedimento na construção da superfície de energia potencial [H,S,F]. Com esta superfície, diversas características deste sistema foram analisadas, tais como geometrias dos pontos estacionários, energias relativas e frequências vibracionais. Os valores obtidos estão de acordo com os dados descritos na literatura. A superfície construída também foi usada para a realização de cálculos de dinâmica para a reação F+HS → S+FH. Observamos a existência de dois tipos de mecanismos, um com a formação de um intermediário de longa duração e outro com a abstração direta do átomo de hidrogênio. / This project has two goals. First, we studied the applicability of the tricubic splines to construct global potential energy surfaces. One of the diculties this approach has to overcome is the choice of an appropriate coordinate system that minimises the in uence of non-physical points. For such, we proposed the use of the Pekeris coordinate system, never employed for this purpose. This procedure was carried out for three well described systems, [Cl,H2], [F,H,D] and [H,O,Cl], whose potential energy surfaces and reaction properties were taken as references. Based on these models, we applied the proposed method varying the amount and arrangement of the interpolation knots, to verify their influence on the quality of the interpolated surfaces. The results showed that surfaces constructed by this approach reproduce very well the chemical dynamics calculations, both for the quantum as well as for the classical methods, provided that the interpolation knots cover the most important regions of the potential energy surfaces, and the lower values of the Pekeris coordinates are prioritised. The second goal was the application of this procedure to the construction of the [H,S,F] potential energy surface. With this surface, several characteristics of this system were analysed, such as the geometry of the stationary points, relative energies and vibrational frequencies. The values obtained are in agreement with the data described in the literature. The constructed surface was also used for quantum dynamics calculations on the reaction F + HS → S + FH. We observed two kinds of mechanisms, one of them with the formation of a long-living intermediate and the other with the direct abstraction of the hydrogen atom.