Computationally assisted vibrational spectroscopy of nucleic acid bases. 2. Thymine.

As in the case of cytosine [Phys. Chem. Chem. Phys. 2023, 25, 24121-24128], Raman and infrared (IR) spectra of aqueous thymine and its N-deuterated derivative, thymine-d2 have been computationally reproduced and interpreted with the use of the recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions [J. Phys. Chem. B, 2020, 124, 6664-6670]. A cluster model of a solute surrounded by 30 water molecules is shown to be sufficient to reproduce experimental vibrational frequencies and relative Raman intensities with the use of B3LYP-D3/def2-TZVP or B3LYP-D3/aug-cc-pVDZ simulations. Analogous PBE-D3 computations provided a less good, but still reasonably accurate, modeling of Raman spectra. It is shown that strong changes of frequencies and relative intensities of the Raman bands of thymine, caused by its hydration, can be interpreted mainly as a result of hydrogen bonding with 6 nearest water molecules. Non-negligible improvement of the quality of simulations for larger clusters comprising water molecules that do not have direct contacts with the solute, suggests that spectroscopic effects of hydration should be ascribed to the joined action of solute-solvent and solvent-solvent interactions. Nevertheless, the moderate number of water molecules required for successful simulations of the Raman spectra of aqueous thymine, suggests that the vibrational modes and derivatives of the polarizability of the solute are mainly locally influenced, while the effect of bulk water is rather modest.

Computational analysis of the vibrational spectra and structure of aqueous cytosine.

The recently developed efficient protocol for the explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is used to describe aqueous solutions of cytosine. The same cluster model of a solute surrounded by the first solvation shell of solvent molecules was shown to be sufficient to reproduce experimental vibrational frequencies and relative IR and Raman intensities. An equally good quality of Raman spectra was provided by B3LYP-D3/def2-TZVP and B3LYP-D3/aug-cc-pVDZ simulations. Computations using the PBE functional were sufficient for modeling of the IR spectra but failed in the simulations of Raman scattering. It is shown that strong changes of frequencies and relative intensities of Raman and IR bands of cytosine, caused by its hydration, cannot be completely assigned to the influence of hydrogen bonds (HBs) with 7 or 8 closest water molecules. They are rather ascribed to the combined effect of solute-solute and solute-solvent HBs with the participation of at least 30 water molecules separating cytosine from the bulk solvent. This suggests that the vibrational modes and derivatives of the polarizability and dipole moment of the solute are mainly locally influenced by its first hydration shell, while the influence of bulk water is rather modest.

Characterization of Conjugation Effects in the Series of Quinoxaline-2-ones by Means of Vibrational Raman Spectroscopy.

A broad series of quinoxalinone-based π-conjugated donor-acceptor fluoro- and NLO-phores is characterized by means of Raman spectroscopy and single-crystal X-ray analysis supported by quantum chemical computations. Intense Raman spectroscopic markers that allow the differentiation of even closely related structures are identified. The intensities of these bands are shown to be related to the conjugation of the different molecular moieties, and they can provide an estimation of its extent. The intensity redistribution between these markers serves as a source of auxiliary structural information capable of pointing to a distortion of the conjugation or to the influence of aggregation effects in the condensed state. A simple relation between the intensity of the marker and the position and oscillator strength of the lowest-energy electronic absorption band of quinoxalinones allows a linking of the Raman effect with the optical properties of these compounds, which can be used for the rational design of novel species with improved optical characteristics.

2,3-(Dibenzimidazol-2-yl)quinoxalines: Unexpected Dynamical Effect on Steady-State Electronic Absorption Spectra.

We report on the electronic absorption spectra, conformational behavior, and intra- and intermolecular hydrogen bonds of 2,3-(dibenzimidazol-2-yl)-quinoxaline (DBIQ). The experimentally found strong solvent dependence of the absorption spectra of DBIQ solutions cannot be assigned to electronic excitations of the equilibrium ground-state DBIQ structure. Extended consideration including the nonequilibrium structures within the framework of ab initio molecular dynamics (MD) revealed the importance of torsion molecular motions not covered by the static case. The strong impact of solute-solvent hydrogen bonding on stabilization of these nonequilibrium structures and on conformational composition of DBIQ was demonstrated. A presence of twisted nonplanar geometries along the whole MD trajectory was shown to drastically influence not only energies but also characters of electronic excitations, resulting in a change of local π-π* character in a solution of 1,2-dichloroethane to charge-transfer character in polar dimethylsulfoxide.

Benzimidazolylquinoxalines: novel fluorophores with tuneable sensitivity to solvent effects.

We report on the photophysical properties, conjugation, conformational behavior, intra- and intermolecular hydrogen bonds (HBs) of a series of novel fluorophores, consisting of 3-arylquinoxaline and benzimidazole moieties linked by a single CC bond. Computations employing density functional theory (DFT) reveal that conjugation between these moieties stabilizes syn-conformers with two HB centers located on the same side of the molecule. Anti-conformers form stronger intermolecular HBs with DMSO and DMF than syn-conformers, and this influences the energy gap between syn- and anti-forms, especially upon excitation of the molecules to the S1 state. Substituents introduced in various positions of the molecules modify their conformational behavior, and mutual disposition of excited singlet states relative to the ground states. Various substitution patterns produce very different effects on relative quantum yield of luminescence: from a moderate increase in polar DMSO and DMF relative to 1,2-dichloroethane solutions to complete quenching of emission which is observable in polar media. The observed behavior is understood with the aid of DFT and time-dependent DFT calculations. The tuneability of the spectroscopic range of the luminescence and especially of its sensitivity to environmental effects allows rational design of the novel fluorophores of this family for various applications.

Conjugation in and optical properties of 1-R-1,2-diphospholes and 1-R-phospholes.

The strength of conjugation between the diene moieties of 1-R-1,2-diphospholes and 1-R-phospholes and exocyclic phenyl groups of these P-heteroles has been quantitatively characterized by the use of Raman activities of the bands of the phenyl substituents. It is shown that conjugation in both types of phospholes is very similar to the conjugation of phenyl groups with the diene system of cyclopentadiene. Introduction of substituents (-OMe, -C(═O)H, -NO2, -NMe2, and -CH═CH2) in the para-position of the phenyl groups of 1-R-1,2-diphospholes extends π-delocalization of exocyclic groups into the electronic system of the 1,2-diphosphole ring, producing bathochromic shifts of the absorption bands up to 63 nm. In contrast, hypsochromic shifts up to 40 nm can be achieved by introduction of SnMe3 or SiMe3 groups at the phosphorus(III) atom of the 1,2-diphosphole and concomitant increase of aromaticity of the P-heterole. Conjugation shifts the "centre of gravity" of the whole electronic absorption spectrum, whereas positions of separate absorption bands are not simply dependent on conjugation lengths.

Application of time-dependent density functional theory and optical spectroscopy toward the rational design of novel 3,4,5-triaryl-1-R-1,2-diphospholes.

Twenty 3,4,5-triaryl-1-R-1,2-diphospholes were studied within the framework of density functional theory (DFT) and experimentally by UV/vis spectroscopy to check their suitability for opto-electronic applications. Time-dependent DFT (TD-DFT) calculations employing the PBE0 hybrid density functional combined with moderately sized def-TZVP basis set were shown to excellently reproduce the experimental absorption spectra of various 1,2-diphospholes. Frontier molecular orbital analysis reveals that HOMO and LUMO are mainly localized on the diphosphole ring and, to some extent, on the aryl moieties. The HOMO-LUMO energy gap can be easily tuned by variation of substituents introduced in para-positions of the aryl moieties and, to a lesser extent, by modification of the R group at phosphorus atom. As a result, both position and intensity of the absorption bands with highest wavelength are strongly influenced by the above structural changes. The UV-spectra simulations reveal that the introduction of donor groups like para-OMe, para-NMe2, and para-N(H)Ph, which are in n-π conjugation with the aryl moieties, should result in absorption of visible light by the corresponding 1,2-diphospholes, thus making them promising candidates for new functional materials.

Is there a simple way to reliable simulations of infrared spectra of organic compounds?

To assess the ability of the quantum-chemical computations to reproduce the experimental relative intensities in the infrared (IR) spectra of both the gas- and condensed-phase systems, the hybrid DFT functional B3LYP has been applied to simulation of IR spectra for species containing from three to twelve first- or second-row atoms, both in the gas phase and in CCl4 solutions. The results demonstrate that B3LYP, combined with the highly compact double-ζ basis set 6-31+G* and "scaled quantum mechanics" techniques, offers excellent quantitative performance in the calculations of relative IR intensities and frequencies (ν ≤ 2200 cm(-1)) for the bands of vibrations of medium-size isolated molecules, whereas it produces unsatisfactory results for the solutions of the same species. Neither larger basis sets nor implicit treatment of the media effects improve the agreement of the simulated spectra with the condensed-phase experiment. At the same time, some preliminary results suggest that explicit modeling of media effects could offer better quality of the IR spectral simulations for the condensed-phase systems.