Solvent Dependence of n−π* Absorption in Acetone

Method to Compute the Solute-Solvent Dispersion Contribution to the Electronic Excitation Energy in Solution

A method formulated within the polarizable continuum model of the solvent and a quantum Monte Carlo treatment of the electronic states of the solute molecule is presented for the calculation of the solute–solvent dispersion contribution to the electronic excitation energy in solution. Variational quantum Monte Carlo is exploited to measure the fluctuations of the electronic electric field of the solute molecule to compute the London’s dispersion forces with the solvent. The method previously applied to the ground state of the solute is here extended to excited states. To perform the Casimir–Polder integration, we introduce a positive parameter Ω whose value is properly chosen for this purpose. We derive a general expression that for Ω = 0 reduces to that already proposed for the ground state. For an excited state, Ω must be less than the first transition electronic energy of the solvent molecule but greater than the transition energy from the ground to excited electronic state of the solute molecule. Benchmark calculations were performed on the n → π* transition for formaldehyde, acrolein, and acetone in six solvents, including water, ethanol, cyclohexane, chloroform, carbon tetrachloride, and toluene, and the π → π* transition of acrolein in cyclohexane. Solvents are characterized by their ionization potential and the refractive index at frequency Ω. In all cases, we found that the dispersion solute–solvent interaction stabilizes the excited state of the solutes leading to red (negative) solvatochromic shifts.

Multiphoton Phosphorescence of Simple Ketones by Visible-light Excitation and Its Consideration for Active Sensing in Space

Acetone and butanone were seen to emit blue light around 450 nm when excited in the green by a high intensity pulsed laser. The pathway of this anti-Stokes emission is believed to be multiphoton absorption followed by phosphorescence, with emission being observed in the samples at cryogenic temperatures below their melting point and not seen from either ketone in their cold liquid state. Given the widespread nature of these simple ketones in off-world bodies and their potential importance as an organic resource for Space Resource Utilization, signals which enable the identification and tracing of these materials are of use in applications from remote sensing and mapping to monitoring during extraction processes. While the excitation process has a low efficiency, the ability to use visible light for sensing of these targets has advantages over UV sources, such as the wider availability of high-powered lasers which could be utilized.

Evaluation of an Efficient 3D-RISM-SCF Implementation as a Tool for Computational Spectroscopy in Solution

The 3D-RISM-SCF solvent-model implementation of Gusarov et al. [ J. Phys. Chem. A 2006, 110, 6083-6090] in the Amsterdam density functional program has been improved and extended. In particular, an accurate yet efficient representation of the solute electrostatic potential is provided. The Coulomb-potential fitting of many DFT codes can be used advantageously in this context. The extra effort compared to a point-charge representation is small for a given SCF cycle and compensated by faster SCF convergence. This allows applications to large solutes, as demonstrated by evaluation of the solvatochromism of Reichardt's dye. In general, TDDFT applications to excitation energies in solution stand out and are highlighted. Applications to the ¹⁷O NMR chemical shifts of N-methylformamide in different solvents also demonstrate the distinct advantages of 3D-RISM over continuum solvents. Limitations are observed in this case for water solvent, where the solvent shielding is overestimated. This shortcoming applies also to the ¹⁷O gas-to-liquid shift of water, where we used localized molecular orbital analyses for a deeper understanding. For such cases of extremely strong solute-solvent interactions, couplings between solute and solvent orbitals induced by the magnetic perturbation are relevant. These clearly require a quantum-mechanical treatment of the most closely bound solvent molecules. Except for such extreme cases, 3D-RISM-SCF is very well suited to treat solvent effects on NMR parameters. More serious limitations pertain to the treatment of vibrational spectra, where the absence of the coupling between solute and solvent vibrational modes limits the accuracy of applications of 3D-RISM-SCF. The reported extended, efficient, and numerically accurate 3D-RISM-SCF implementation should provide a useful tool to study chemical and spectroscopic properties of molecules of appreciable size in a realistic solvent environment.

Solvent effects for vertical absorption and emission processes in solution using a self-consistent state specific method based on constrained equilibrium thermodynamics

A self-consistent state specific (SS) method in the framework of TDDFT is presented to account for solvent effects on absorption and emission processes for molecules in solution. In these processes, the initial state is an equilibrium state, while the polarization of the solvent is in nonequilibrium with the electron density of the solute in the final state. Nonequilibrium solvation free energy is calculated based on a novel nonequilibrium solvation model with constrained equilibrium manipulation. The bulk solvent effects are considered using the polarizable continuum method (PCM), where the solvent-solute interaction is described with a reaction field. Molecular orbitals and orbital energies in the presence of the reaction field corresponding to the excited state are employed and the response of the solvent is not included in the TDDFT calculations. A self-consistent procedure is designed to obtain the excited state reaction field. The equations based on this new nonequilibrium solvation model in the framework of the self-consistent SS-PCM/TDDFT method for calculation of vertical absorption and emission energies are presented and implemented in the Q-Chem package. Vertical absorption and emission energies for several small molecules in solution using the newly developed code are calculated and compared with available experimental data and the results of other theoretical studies. Solvent shifts of absorption and emission energies are reasonably reproduced with this approach. The new model is a promising approach to study nonequilibrium absorption and emission processes in solution.

Solvent effects on excitation energies obtained using the state-specific TD-DFT method with a polarizable continuum model based on constrained equilibrium thermodynamics

Nonequilibrium solvation effects need to be treated properly in the study of electronic absorption processes of solutes since solvent polarization is not in equilibrium with the excited-state charge density of the solute. In this work, we developed a state specific (SS) method based on the novel nonequilibrium solvation model with constrained equilibrium manipulation to account for solvation effects in electronic absorption processes. Time-dependent density functional theory (TD-DFT) is adopted to calculate electronic excitation energies and a polarizable continuum model is employed in the treatment of bulk solvent effects on both the ground and excited electronic states. The equations based on this novel nonequilibrium solvation model in the framework of TDDFT to calculate vertical excitation energy are presented and implemented in the Q-Chem package. The implementation is validated by comparing reorganization energies for charge transfer excitations between two atoms obtained from Q-Chem and those obtained using a two-sphere model. Solvent effects on electronic transitions of coumarin 153 (C153), acetone, pyridine, (2E)-3-(3,4-dimethoxyphenyl)-1-(2-hydroxyphenyl)prop-2-en-1-one (DMHP), and uracil in different solvents are investigated using the newly developed code. Our results show that the obtained vertical excitation energies as well as spectral shifts generally agree better with the available experimental values than those obtained using the traditional nonequlibrium solvation model. This new model is thus appropriate to study nonequilibrium excitation processes in solution.

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.

Pick and Choose the Spectroscopic Method to Calibrate the Local Electric Field inside Proteins

Electrostatic interactions in proteins play a crucial role in determining the structure-function relation in biomolecules. In recent years, fluorescent probes have been extensively employed to interrogate the polarity in biological cavities through dielectric constants or semiempirical polarity scales. A choice of multiple spectroscopic methods, not limited by fluorophores, along with a molecular level description of electrostatics involving solute-solvent interactions, would allow more flexibility to pick and choose the experimental technique to determine the local electrostatics within protein interiors. In this work we report that ultraviolet/visible-absorption, infrared-absorption, or (13)C NMR can be used to calibrate the local electric field in both hydrogen bonded and non-hydrogen bonded protein environments. The local electric field at the binding site of a serum protein has been determined using the absorption wavelength as well as the carbonyl stretching frequency of its natural steroid substrate, testosterone. Excellent agreement is observed in the results obtained from two independent spectroscopic techniques.

Electrostatic Interactions Are Key to C═O n-π* Shifts: An Experimental Proof

Carbonyl n-π* transitions are known to undergo blue shift in polar and hydrogen-bonding solvents. Using semiempirical expressions, previous studies hypothesized several factors like change in dipole moment and hydrogen-bond strength upon excitation to cause the blue shift. Theoretically, ground-state electrostatics has been predicted to be the key to the observed shifts, however, an experimental proof has been lacking. Our experimental results demonstrate a consistent linear correlation between IR (ground-state phenomenon) and n-π* frequency shifts (involves both ground and excited electronic-states) of carbonyls in hydrogen-bonded and non-hydrogen-bonded environments. The carbonyl hydrogen-bonding status is experimentally verified from deviation in n-π*/fluorescence correlation. The IR/n-π* correlation validates the key role of electrostatic stabilization of the ground state toward n-π* shifts and demonstrates the electrostatic nature of carbonyl hydrogen bonds. n-π* shifts show linear sensitivity to calculated electrostatic fields on carbonyls. Our results portray the potential for n-π* absorption to estimate local polarity in biomolecules and to probe chemical reactions involving carbonyl activation/stabilization.

Frozen-density embedding theory with average solvent charge densities from explicit atomistic simulations

Besides molecular electron densities obtained within the Born–Oppenheimer approximation (ρ_B(r)) to represent the environment, the ensemble averaged density (〈ρ_B〉(r)) is also admissible in frozen-density embedding theory (FDET) [Wesolowski, Phys. Rev. A, 2008, 77, 11444].

Solvent Induced Transformations of n-π* Absorption in Formaldehyde, Acetaldehyde, and Acetone

Absorption spectra of formaldehyde (FA), acetaldehyde (AA), and acetone are compared in the vapor phase, nonpolar, and polar solutions at 295 K. The vibronic n-π* transition of carbonyl chromophore is mainly composed of the overtones of >C═O stretching vibration. A new phenomenon is observed in liquid solutions, consisting of a relative increase of Franck-Condon factors for the second and third harmonics in FA, and the second to fourth replica in AA, with respect to the gas phase. In AA and acetone with poorly resolved vibronic structure, the redistribution of intensities produces a false "solvent shift" of the band maximum between the vapor and nonpolar liquid phase by -250 ± 50 cm(-1). Modification in vibronic coupling can also explain unusual narrowing of the band contour in the solution, reported earlier for acetone (Renge , I. J. Phys. Chem. A 2009, 113, 10678). No detectable shift occurs as a function of solvent polarizability (refractive index function (n(2) - 1)/(n(2) + 2)) in n-alkanes for FA, AA, and acetone, as well as for cyclopentanone and camphor. Incidentally, the bathochromic dispersive shift is almost exactly compensated by a hypsochromic induction shift. The latter is due to the diminishing dipole moment in the excited state of the carbonyl chromophore. Differences in polarizability α and dipole moments μ were estimated for FA (Δα = 0.33 ± 0.1 Å(3)), AA (Δμ = -1.05 ± 0.2 D, Δα = 0.5 ± 0.2 Å(3)), and acetone (Δμ = -1.3 ± 0.2 D, Δα = 0.65 ± 0.2 Å(3)). The increase of α by ∼10% upon excitation is plausible for a weak n-π* transition. By contrast, near doubling of α in the upper state has been reported recently for several ketones, with Δα reaching 10 Å(3) (Catalán, J.; Catalán, J. P. Phys. Chem. Chem. Phys. 2011, 13, 4072). Empirical partitioning of solvent shifts into repulsive-dispersive, induction, dipole-dipole, and hydrogen bonding contributions was proposed to serve as a benchmark in computer chemical calculations.

Spectral shift mechanisms of chlorophylls in liquids and proteins

Origins of non-excitonic spectral shifts of chlorophylls that can reach -1,000 cm(-1) in pigment-protein complexes are actively debated in literature. We investigate possible shift mechanisms, basing on absorption and fluorescence measurements in large number of liquids. Transition wavelength in solvent-free state was estimated (±2 nm) for chlorophyll a (Chl a, 647 nm), Chl b (624 nm), bacteriochlorophyll a (BChl a, 752 nm), and pheophytines. The dispersive-repulsive shift is a predominating mechanism. It depends on polarizability difference between the ground and the excited state Δα and the Lorenz-Lorentz function of refractive index of solvent (n). The approximate (± 2Å(3)) increase of polarizability Δα is close to 15Å(3) for S(1) bands of Chl a, BChl a, and BPheo a, slightly larger for Chl b (18Å(3)), and less for Pheo a (11Å(3)). The effect of solvent polarity, expressed in terms of static dielectric permittivity (ε) is relatively minor, but characteristic for different pigments and transitions. Remarkably, maximum influence of ε on S(1) band of BChl a is less (-20 ± 10 cm(-1)) than that for Chl a (-50 ± 10 cm(-1)), and not correlated with dipole moment changes on excitation Δμ (∼2D and 0.1 ± 0.1D, respectively). Hydrogen bonding in protic solvents produces red shifts in Chl a (-60 cm(-1)) and BChl a (-100 cm(-1)), but not in Chl b. Second axial ligand of BChl a has no influence on the S(1) band, whereas the S(2) transition suffers a -400 to -600 cm(-1) down shift. Aromatic character of solvent is responsible for a ∼-100 cm(-1) red shift of both Q transitions in BChl a. The S(1) bands in chlorophylls are relatively insensitive with respect to dielectric properties and specific solvation. Therefore, nontrivial mechanisms, yielding large site-energy shifts are expected in photosynthetic chlorophyll-proteins.

Solvent effects on the steady state photophysics of estrone and 17β-estradiol

Absorption and emission yields for estrone and 17β-estradiol were measured in a variety of room temperature solvents. Molar extinction coefficients were found to not vary as a function of solvent, while fluorescence yields were found to be significantly affected by the polarity and hydrogen-bond accepting ability of the solvent, with the yield for 17β-estradiol being highest in nonpolar, hydrogen-bond donating solvents, and lowest in the nonpolar, hydrogen-bond accepting solvent ethyl acetate. Estrone's emission yield was found to be a factor of ten smaller than 17β-estradiol's. Strong solvent and excitation wavelength dependences were found for the relative amounts of emission between estrone's two emission bands, with increased relative emission occurring in nonpolar aprotic solvents, and under higher excitation energies. These results are interpreted with the aid of vertical excitation energies from time-dependent density functional calculations using both explicit and implicit solvation models.

Solvation Effects on Electronic Transitions: Exploring the Performance of Advanced Solvent Potentials in Polarizable Embedding Calculations

The polarizable embedding (PE) approach, which combines quantum mechanics (QM) and molecular mechanics (MM), is applied to predict solvatochromic effects on excitation energies of several representative molecules in aqueous, methanol, acetonitrile, and carbon tetrachloride solutions. Good agreement with experimental results for excitation energies and for solvatochromic shifts is demonstrated on the basis of either density functional theory or coupled cluster methods. Solvent-dependent trends are fully reproduced in this diverse set of solvents. Furthermore, it is shown that the inclusion of higher order multipole moments and anisotropic polarizabilities in the electrostatic embedding potentials leads to a faster convergence with respect to a full QM treatment (within about 0.1 eV of estimated full QM treatments). It is thereby illustrated that the use of advanced solvent potentials can provide higher accuracy compared to various simpler approaches for the prediction of solvent shifts and do so in a computationally competitive manner.