Comparison of polarizable continuum model and quantum mechanics/molecular mechanics solute electronic polarization: Study of the optical and magnetic properties of diazines in water

Theoretical investigation of solvent and oxidation/deprotonation effects on the electronic structure of a mononuclear Ru-aqua-polypyridine complex in aqueous solution

Mononuclear polypyridine ruthenium (Ru) complexes can catalyze various reactions, including water splitting, and can also serve as photosensitizers in solar cells. Despite recent progress in their synthesis, accurately modeling their physicochemical properties, particularly in solution, remains challenging. Herein, we conduct a theoretical investigation of the structural and electronic properties of a mononuclear Ru-aqua polypyridine complex in aqueous solution, considering five of its possible oxidation/protonation states species: [RuII(H2O)(py)(bpy)2]2+, [RuII(OH)(py)(bpy)2]+, [RuIII(H2O)(py)(bpy)2]3+, [RuIII(OH)(py)(bpy)2]2+ and [RuIV(O)(py)(bpy)2]2+, where py = pyridine and bpy = 2,2'-bipyridine. At first, we investigate the impact of proton-coupled and non-coupled electron transfer reactions on the geometry and electronic structure of the complexes in vacuum and in solution, using an implicit solvent model. Then, using a sequential multiscale approach that combines quantum mechanics and molecular mechanics (S-QM/MM), we examine the explicit solvent effects on the electronic excitations of the complexes, and compare them with the experimental results. The complexes were synthesized, and their absorption spectra measured in aqueous solution. To accurately describe the QM interactions between the metal center and the aqueous ligand in the MM simulations, we developed new force field parameters for the Ru atom. We analyze the solvent structure around the complexes and account for its explicit influence on the polarization and electronic excitations of the complexes. Notably, accounting for the explicit solvent polarization effects of the first solvation shells is essential to correctly describe the energy of the electronic transitions, and the explicit treatment of the hydrogen bonds at the QM level in the excitation calculations improves the accuracy of the description of the metal-to-ligand charge-transfer bands. Transition density matrix analysis is used to characterize all electronic transitions in the visible and ultraviolet ranges according to their charge-transfer (CT) character. This study elucidates the electronic structure of those ruthenium polypyridyl complexes in aqueous solution and underscores the importance of precisely describing solvent effects, which can be achieved employing the S-QM/MM method.

The three pyridazines, three naphthyridines and two azoles: effect of the position of the second heteroatom on pK aH of their eight conjugate acids

In the present work, how the position of the second nitrogen in the conjugate acids of the three pyridazines viz., 1,2-pyridazine, 1,3-pyridazine (the pyrimidine) and 1,4-pyridazine (the pyrazine) and three naphthyridines viz., cinnoline, quinazoline and quinoxaline changes the pK aH systematically is taken up. They decrease nearly by a factor of half each time in the class of their own. In contrast there is an increase in the pK aH when we move from pyrazole to imidazole. The pK aH of pyrazole is less than imidazole by −4.45 units. Suitable explanations are given.

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the ¹H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (¹H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pK_as) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pK_a1 = 1.3/1.2, pK_a2 = 2.1/2.2 and pK_a5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pK_a3 = 3.7 and pK_a4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

Spectrophotometric studies of charge-transfer complexes formed with ions N,N'-alkyldiyl-bis(pyridinium) derivatives and iodide

The charge-transfer complexes (CTC) between electron donor iodide and pyridinium dimers and monomers as acceptors have been studied spectrophotometrically in acetonitrile, DMSO and water. The dimers were: N,N'-alkyldiyl-bis(4-cyanopyridinium) and N,N'-alkyldiyl-bis(2-bromopyridinium), and the monomers were: N-alkyl-4-cyanopyridinium and N-alkyl-2-bromopyridinium, bridged by n methylene spacers. The formation constant (K_CTC), molar absorptivity (ε_CTC), fluorescence-quenching constant (K_SV) were determined. The results indicate that the stoichiometry of the CTC for dimers and monomers is 1:1 (equimolar ratio), and its formation is favored for dimers more than for monomers, especially dimers with short methylene spacers forming a "sandwiched type-complex" in which the iodide is close to the two pyridinium rings. Solvents with low polarity favored the complex, which was destroyed by the presence of water. The experimental studies were complemented with theoretical studies with quantum mechanics (QM) calculations using density functional theory (DFT) and molecular mechanics (MM) simulations with Molecular Dynamics for identify the most stable conformers in acetonitrile solution. The electronic excitations were calculated with sequential QM/MM hybrid method, showing a charge-transfer wavelength in agreement with the UV-Vis absorption spectra obtained experimentally. It confirms the "sandwiched type-complex" conformers favoring the interaction of the iodide with one pyridinium rings and simultaneously forming one, or two, hydrogen bonds with the alkyl chain and additionally, for the case of dimers, with the other pyridinium ring.

Unveiling the relationship between structural and polarization effects on the first hyperpolarizability of a merocyanine dye

The nonlinear optical response, more specifically the Hyper-Rayleigh Scattering (HRS) response of the Brooker's merocyanine, has been calculated at the time-dependent density functional theory level and rationalized in terms of the structural changes and polarization effects induced by applied external electric fields. The structural change leads to large changes in the HRS response, while only slight variations were observed due to the polarization effects on the fixed quinoid form. Considering both structural and polarization contributions concurrently, the HRS response is dominated by cooperative behavior of those effects for weak and intermediate electric field strengths. At the same time, the competition between both effects was a crucial factor in the region of strong electric fields. The obtained results can lead to an easier understanding for upcoming studies considering more realistic models of solvents where it is not simple to disentangle these contributions.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.