In this paper a novel approach to study the formation and relaxation of excited states in solution is presented within the integral equation formalism version of the polarizable continuum model. Such an approach uses the excited state relaxed density matrix to correct the time dependent density functional theory excitation energies and it introduces a state-specific solvent response, which can be further generalized within a time dependent formalism. This generalization is based on the use of a complex dielectric permittivity as a function of the frequency, ε̂(ω). The approach is here presented in its theoretical formulation and applied to the various steps involved in the formation and relaxation of electronic excited states in solvated molecules. In particular, vertical excitations (and emissions), as well as time dependent Stokes shift and complete relaxation from vertical excited states back to ground state, can be obtained as different applications of the same theory. Numerical results on two molecular systems are reported to better illustrate the features of the model.

We present a formal comparison between the two different approaches to the calculation of electronic excitation energies of molecules in solution within the continuum solvation model framework, taking also into account nonequilibrium effects. These two approaches, one based on the explicit evaluation of the excited state wave function of the solute and the other based on the linear response theory, are here proven to give formally different expressions for the excitation energies even when exact eigenstates are considered. Calculations performed for some illustrative examples show that this formal difference has sensible effects on absolute solvatochromic shifts (i.e., with respect to gas phase) while it has small effects on relative (i.e., nonpolar to polar solvent) solvatochromic shifts.

The study of molecular systems in a condensed phase with quantum-mechanical tools introduces problems not present in analogous studies on isolated molecules. Some of these problems are presented in the introductory part to justify a strategy to elaborate a theoretical model starting from an accurate continuum description of the solvent. This is here achieved through the polarizable continuum model (PCM), which also permits explicit consideration of solvent molecules. The decision of relying on the description of a large number of properties to refine the model is then justified. The consideration of complex molecular properties of a various nature has led to the elaboration of numerous additional features not present in the basic version of PCM, some among which are here briefly presented: local field and nonequilibrium effects, specific solute–solute and solute–solvent interactions, description of condensed phase with more complex structure. The paper ends with the presentation of some results, selected among a larger set of available examples, and limited to a few properties, to show the potentialities of the approach.

Understanding protein–inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein–surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein–surface interactions and discuss the successes and limitations of current approaches.