skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • The effect of solvent on determining highest occupied molecular orbital energies of semiconducting organic molecules: Insight from a combined computational approach
Title: The effect of solvent on determining highest occupied molecular orbital energies of semiconducting organic molecules: Insight from a combined computational approach
Abstract Although cyclic voltammetry (CV) measurements in solution have been widely used to determine the highest occupied molecular orbital energy (EHOMO) of semiconducting organic molecules, an understanding of the experimentally observed discrepancies due to the solvent used is lacking. To explain these differences, we investigate the solvent effects onEHOMOby combining density functional theory and molecular dynamics calculations for four donor molecules with a common backbone moiety. We compare the experimentalEHOMOvalues to the calculated values obtained from either implicit or first solvation shell theories. We find that the first solvation shell method can capture theEHOMOvariation arising from the functional groups in solution, unlike the implicit method. We further applied the first solvation shell method to other semiconducting organic molecules measured in solutions for different solvents. We find that theEHOMOobtained using an implicit method is insensitive to solvent choice. The first solvation shell, however, producesEHOMOvalues that are sensitive to solvent choices and agrees with published experimental results. The solvent sensitivity arises from a hierarchy of three effects: (1) the solute electronic state within a surrounding dielectric continuum, (2) ambient temperature or solvent atoms changing the solute geometry, and (3) electronic interactions between the solute and solvents. The implicit method, on the other hand, only captures the effect of a dielectric environment. Our findings suggest thatEHOMOobtained by CV measurements should account for the influence of solvent when the results are reported, interpreted, or compared to other molecules.  more » « less
Award ID(s):
1916612
PAR ID:
10389781
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Computational Chemistry
Volume:
44
Issue:
10
ISSN:
0192-8651
Page Range / eLocation ID:
p. 1064-1072
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; (, Journal of chemical theory and computation)
    null (Ed.)
    We develop a hybrid approach that combines the Monte Carlo (MC)method, a variational implicit-solvent model (VISM), and a binary level-set method forthe simulation of biomolecular binding in an aqueous solvent. The solvation free energy for the biomolecular complex is estimated by minimizing the VISM free-energy functional of all possible solute−solvent interfaces that are used as dielectric boundaries. This functional consists of the solute volumetric, solute−solvent interfacial, solute−solvent van der Waals interaction, and electrostatic free energy. A technique of shifting the dielectric boundary is used to accurately predict the electrostatic part of the solvation free energy.Minimizing such a functional in each MC move is made possible by our new and fast binary level-set method. This method is based on the approximation of surface area by the convolution of an indicator function with a compactly supported kernel and is implemented by simple flips of numerical grid cells locally around the solute−solvent interface. We apply our approach to the p53-MDM2 system for which the two molecules are approximated by rigid bodies. Our efficient approach captures some of the poses before the final bound state. All atom molecular dynamics simulations with most of such poses quickly reach the final bound state.Our work is a new step toward realistic simulations of biomolecular interactions. With further improvement of coarse graining and MC sampling, and combined with other models, our hybrid approach can be used to study the free-energy landscape and kinetic pathways of ligand binding to proteins. 
    more » « less
  2. ; (, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)
    The variational implicit-solvent model (VISM) is an efficient approach to biomolecular interactions, where electrostatic interactions are crucial. The total VISM free energy of a dielectric boundary (i.e. solute–solvent interface) consists of the interfacial energy, solute–solvent interaction energy and dielectric electrostatic energy. The last part is the maximum value of the classical and concave Poisson–Boltzmann (PB) energy functional of electrostatic potentials, with the maximizer being the equilibrium electrostatic potential governed by the PB equation. For the consistency of energy minimization and computational stability, here we propose alternatively to minimize the convex Legendre-transformed Poisson–Boltzmann (LTPB) electrostatic energy functional of all dielectric displacements constrained by Gauss’ Law in the solute region. Both integrable and discrete solute charge densities are treated, and the duality of the LTPB and PB functionals is established. A penalty method is designed for the constrained minimization of the LTPB functional. In application to biomolecular interactions, we minimize the total VISM free energy iteratively, while in each step of such iteration, minimize the LTPB energy. Convergence of such a min–min algorithm is shown. Our numerical results on the solvation of a single ion indicate that the LTPB performs better than the PB formulation, providing possibilities for efficient biomolecular simulations. 
    more » « less
  3. ; ; (, Journal of Computational Chemistry)
    Abstract Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self‐consistent continuum solvation (SCCS) and the soft‐sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation‐free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys.2019,150, 041710.) for the SCCS model, a single‐parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; (, Journal of molecular liquids)
    Non-toxic, chemically inert, organic polymers as polyethylene glycol (PEG) and polyoxymethylene (POM) have versatile applications in basic research, industry and pharmacy. In this work, we aim to characterize the hydration structure of PEG and POM oligomers by exploring how the solute disturbs the water structure compared to the bulk solvent and how the solute chain interacts with the solvent. We explore the effect of (i) the C–C–O (PEG) versus C-O (POM) constitution of the chain and (ii) chain length. To this end, MD simulations followed by clustering and topological analysis of the hydration network, as well as quantum mechanical calculations of atomic charges are used. We show that the hydration varies with chain conformation and length. The degree of folding of the chain impacts its degree of solvation, which is measurable by different parameters as for example the number of water molecules in the first solvation shell and the solvent accessible surface. Atomic charges calculated on the oligomers in gas phase are stable throughout conformation and chain length and seem not to determine solvation. Hydration however induces charge transfer from the solute molecule to the solvent, which depends on the degree of hydration. 
    more » « less
  5. ; ; ; ; ; ; (, The Journal of Chemical Physics)
    Computational simulation of biomolecules can provide important insights into protein design, protein-ligand binding interactions, and ab initio biomolecular folding, among other applications. Accurate treatment of the solvent environment is essential in such applications, but the use of explicit solvents can add considerable cost. Implicit treatment of solvent effects using a dielectric continuum model is an attractive alternative to explicit solvation since it is able to describe solvation effects without the inclusion of solvent degrees of freedom. Previously, we described the development and parameterization of implicit solvent models for small molecules. Here, we extend the parameterization of the generalized Kirkwood (GK) implicit solvent model for use with biomolecules described by the AMOEBA force field via the addition of corrections to the calculation of effective radii that account for interstitial spaces that arise within biomolecules. These include element-specific pairwise descreening scale factors, a short-range neck contribution to describe the solvent-excluded space between pairs of nearby atoms, and finally tanh-based rescaling of the overall descreening integral. We then apply the AMOEBA/GK implicit solvent to a set of ten proteins and achieve an average coordinate root mean square deviation for the experimental structures of 2.0 Å across 500 ns simulations. Overall, the continued development of implicit solvent models will help facilitate the simulation of biomolecules on mechanistically relevant timescales. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email