Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

Toward a numerically efficient description of bulk-solvated anionic states

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute-solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute-solvent and solvent-solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models.

Novel Approach for Predicting Vertical Electron Attachment Energies in Bulk-Solvated Molecules

When low-energy electrons interact with molecules, they can give rise to transient anion states commonly known as resonances. These states are formed through vertical electron attachment processes and have the potential to induce various forms of DNA lesions, including base damage, single- and double-strand breaks, cross-links, and clustered lesions that are challenging to repair. So far, most experimental and theoretical studies have investigated the formation of resonances of (bio)molecules in the gas phase or in microsolvated environments. Since cellular environments are mainly composed of water molecules, it is crucial to understand how bulk water affects the resonances of (bio)molecules. Given the existing gap in studies on resonances of bulk-solvated molecules, we propose a novel theoretical-computational approach to address this void. Our approach combines the multibasis-set (time-dependent-)density functional theory and self-consistent sequential quantum mechanics/molecular mechanics polarizable electrostatic embedding methods. We apply this combined methodology to predict the vertical electron attachment energies of 1-methyl-5-nitroimidazole (1M5NI), a well-known radiosensitizer model, in bulk water. In addition, we analyze the rapid mutual polarization between the resonances (both shape- and core-excited) of 1M5NI and the surrounding bulk water environment. For comparison, we also studied the isolated and microsolvated 1M5NI. Overall, while the polarization of the environment is clearly sensitive to the solute charge, causing a significant impact on the vertical electron affinity and consequently on the attachment electron energies, it does not have a significant impact on the excitation energies of the anion.

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 Importance of the Density Functional Theory Exchange-Correlation Hartree-Fock Term in Magnetic Resonance: Application to an Aqueous Environment

Within the framework of Density Functional Theory (DFT), the relevance of the term Hartree-Fock exchange (HFE) for a variety of molecular properties is a critical point. For this reason, we spend efforts to understand these relationships in nuclear magnetic resonance (NMR) parameters in a water solvent. This work takes advantage of the appropriate aug-cc-pVTZ-J basis set and the Minnesota family of DFT methods, which consider different portions of HFE contributions. With regard to solvent participation, the results are based on a sequential Monte Carlo/Quantum Mechanics procedure, which builds the structures of the liquid under realistic thermodynamic conditions. Compared to the accurate results of second-order polarization propagator approximation (SOPPA) and experimental data, all NMR parameters show a huge dependence on the size of the HFE contribution. For instance, the inclusion of this term in ¹J_OH and ²J_HH indirect spin-spin couplings does vary with 49.661 and 25.459 Hz, respectively. The M06-HF method accounts for 100% of HFE and better matches the σ_O and σ_H shielding constants. On the other hand, ¹J_OH and ²J_HH demand a medium contribution (54% of HFE), the best description being associated with the M06-2X method. Thus, the dependence varies regarding the phenomenology of the property in focus and the order for independent treatments. For elements that participate in hydrogen bonds simultaneously as donor and acceptor actors, the results indicate that explicit solvent molecules must be considered in the quantum mechanical calculations for better modeling of paramagnetic shielding constants.

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.

Breaking the Symmetry of Pyrimidine: Solvent Effects and Core-Excited State Dynamics

Symmetry and its breaking crucially define the chemical properties of molecules and their functionality. Resonant inelastic X-ray scattering is a local electronic structure probe reporting on molecular symmetry and its dynamical breaking within the femtosecond scattering duration. Here, we study pyrimidine, a system from the C_2v point group, in an aqueous solution environment, using scattering though its 2a₂ resonance. Despite the absence of clean parity selection rules for decay transitions from in-plane orbitals, scattering channels including decay from the 7b₂ and 11a₁ orbitals with nitrogen lone pair character are a direct probe for molecular symmetry. Computed spectra of explicitly solvated molecules sampled from a molecular dynamics simulation are combined with the results of a quantum dynamical description of the X-ray scattering process. We observe dominant signatures of core-excited Jahn-Teller induced symmetry breaking for resonant excitation. Solvent contributions are separable by shortening of the effective scattering duration through excitation energy detuning.

Accurate Prediction of Absorption Spectral Shifts of Proteorhodopsin Using a Fragment-Based Quantum Mechanical Method

Many experiments have been carried out to display different colors of Proteorhodopsin (PR) and its mutants, but the mechanism of color tuning of PR was not fully elucidated. In this study, we applied the Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps (EE-GMFCC) method to the prediction of excitation energies of PRs. Excitation energies of 10 variants of Blue Proteorhodopsin (BPR-PR105Q) in residue 105GLN were calculated with the EE-GMFCC method at the TD-B3LYP/6-31G* level. The calculated results show good correlation with the experimental values of absorption wavelengths, although the experimental wavelength range among these systems is less than 50 nm. The ensemble-averaged electric fields along the polyene chain of retinal correlated well with EE-GMFCC calculated excitation energies for these 10 PRs, suggesting that electrostatic interactions from nearby residues are responsible for the color tuning. We also utilized the GMFCC method to decompose the excitation energy contribution per residue surrounding the chromophore. Our results show that residues ASP97 and ASP227 have the largest contribution to the absorption spectral shift of PR among the nearby residues of retinal. This work demonstrates that the EE-GMFCC method can be applied to accurately predict the absorption spectral shifts for biomacromolecules.

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.

First principles theoretical spectroscopy of methylene blue: Between limitations of time-dependent density functional theory approximations and its realistic description in the solvent

Methylene blue [3,7-Bis(di-methylamino) phenothiazin-5-ium chloride] is a phenothiazine dye with applications as a sensitizer for photodynamic therapy, photoantimicrobials, and dye-sensitized solar cells. Time-dependent density functional theory (TDDFT), based on (semi)local and global hybrid exchange-correlation functionals, fails to correctly describe its spectral features due to known limitations for describing optical excitations of π-conjugated systems. Here, we use TDDFT with a non-empirical optimally tuned range-separated hybrid functional to explore the optical excitations of gas phase and solvated methylene blue. We compute solvated configurations using molecular dynamics and an iterative procedure to account for explicit solute polarization. We rationalize and validate that by extrapolating the optimized range separation parameter to an infinite amount of solvating molecules, the optical gap of methylene blue is well described. Moreover, this method allows us to resolve contributions from solvent-solute intermolecular interactions and dielectric screening. We validate our results by comparing them to first-principles calculations based on the GW+Bethe-Salpeter equation approach and experiment. Vibronic calculations using TDDFT and the generating function method account for the spectra's subbands and bring the computed transition energies to within 0.15 eV of the experimental data. This methodology is expected to perform equivalently well for describing solvated spectra of π-conjugated systems.

A new interpretation of the absorption and the dual fluorescence of Prodan in solution

Remarkable interest is associated with the interpretation of the Prodan fluorescent spectrum. A sequential hybrid Quantum Mechanics/Molecular Mechanics method was used to establish that the fluorescent emission occurs from two different excited states, resulting in a broad asymmetric emission spectrum. The absorption spectra in several solvents were measured and calculated using different theoretical models presenting excellent agreement. All theoretical models [semiempirical, time dependent density functional theory and and second-order multiconfigurational perturbation theory] agree that the first observed band at the absorption spectrum in solution is composed of three electronic excitations very close in energy. Then, the electronic excitation around 340 nm-360 nm may populate the first three excited states (π-π^*L_b, n-π^*, and π-π^*L_a). The ground state S₀ and the first three excited states were analyzed using multi-configurational calculations. The corresponding equilibrium geometries are all planar in vacuum. Considering the solvent effects in the electronic structure of the solute and in the solvent relaxation around the solute, it was identified that these three excited states can change the relative order depending on the solvent polarity, and following the minimum path energy, internal conversions may occur. A consistent explanation of the experimental data is obtained with the conclusive interpretation that the two bands observed in the fluorescent spectrum of Prodan, in several solvents, are due to the emission from two independent states. Our results indicate that these are the n-π^* S₂ state with a small dipole moment at a lower emission energy and the π-π^*L_b S₁ state with large dipole moment at a higher emission energy.

Averaged Condensed Phase Model for Simulating Molecules in Complex Environments

The need for configurational sampling dramatically increases the cost of combined quantum mechanics/molecular mechanics (QM/MM) simulations of chemical processes in solution. We developed an averaged condensed phase environment (ACPE) model that constructs an effective polarizable environment directly from explicitly sampled molecular dynamics configurations via the K-means++ algorithm and a mathematically rigorous translation of the molecular mechanics parameters. The model captures detailed heterogeneous features in the environment that may be difficult to describe using a conventional polarizable continuum model. Instead of performing repeated QM/MM calculations for each new configuration of the environment, the ACPE approach allows one to perform a single QM calculation on an averaged configuration. Here, we demonstrate the model by computing electronic excitation energies for several small molecules in solution. The ACPE model predicts the excitation energies in excellent agreement with conventional configurational averaging yet with orders of magnitude of reduction in the computational cost.

The Effect of Molecular Conformation on the Accuracy of Theoretical (1)H and (13)C Chemical Shifts Calculated by Ab Initio Methods for Metabolic Mixture Analysis

NMR spectroscopy is a powerful method for analyzing metabolic mixtures. The information obtained from an NMR spectrum is in the form of physical parameters, such as chemical shifts, and construction of databases for many metabolites will be useful for data interpretation. To increase the accuracy of theoretical chemical shifts for development of a database for a variety of metabolites, the effects of sets of conformations (structural ensembles) and the levels of theory on computations of theoretical chemical shifts were systematically investigated for a set of 29 small molecules in the present study. For each of the 29 compounds, 101 structures were generated by classical molecular dynamics at 298.15 K, and then theoretical chemical shifts for 164 (1)H and 123 (13)C atoms were calculated by ab initio quantum chemical methods. Six levels of theory were used by pairing Hartree-Fock, B3LYP (density functional theory), or second order Møller-Plesset perturbation with 6-31G or aug-cc-pVDZ basis set. The six average fluctuations in the (1)H chemical shift were ±0.63, ± 0.59, ± 0.70, ± 0.62, ± 0.75, and ±0.66 ppm for the structural ensembles, and the six average errors were ±0.34, ± 0.27, ± 0.32, ± 0.25, ± 0.32, and ±0.25 ppm. The results showed that chemical shift fluctuations with changes in the conformation because of molecular motion were larger than the differences between computed and experimental chemical shifts for all six levels of theory. In conclusion, selection of an appropriate structural ensemble should be performed before theoretical chemical shift calculations for development of an accurate database for a variety of metabolites.

Isotropic magnetic shielding constants of retinal derivatives in aprotic and protic solvents

We investigate the nuclear isotropic shielding constants σ((13)C) and σ((17)O) of isomers of retinoic acid and retinal in gas-phase and in chloroform, acetonitrile, methanol, and water solutions via Monte Carlo simulation and quantum mechanics calculations using the GIAO-B3LYP∕6-311++G(2d,2p) approach. Electronic solute polarization effects due to protic and aprotic solvents are included iteratively and play an important role in the quantitative determination of oxygen shielding constants. Our MP2∕6-31G+(d) results show substantial increases of the dipole moment of both retinal derivatives in solution as compared with the gas-phase results (between 22% and 26% in chloroform and between 55% and 99% in water). For the oxygen atoms the influence of the solute polarization is mild for σ((17)O) of hydroxyl group, even in protic solvents, but it is particularly important for σ((17)O) of carbonyl group. For the latter, there is a sizable increase in the magnitude with increasing solvent polarity. For the carbon atoms, the solvent effects on the σ((13)C) values are in general small, being more appreciable in carbon atoms of the polyene chain than in the carbon atoms of the β-ionone ring and methyl groups. The results also show that isomeric changes on the backbones of the polyene chains have marked influence on the (13)C chemical shifts of carbon atoms near to the structural distortions, in good agreement with the experimental results measured in solution.