Theoretical modeling of spectroscopic properties of molecules in solution: toward an effective dynamical discrete/continuum approach

An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning

Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.

Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

Direct or Indirect ESPT Mechanism in CFP psamFP488? A Theoretical-Computational Investigation

Fluorescent Proteins are widely studied for their multiple applications in technological and biotechnological fields. Despite this, they continue to represent a challenge in terms of a complete understanding of all the non-equilibrium photo-induced processes that rule their properties. In this context, a theoretical-computational approach can support experimental results in unveiling and understanding the processes taking place after electronic excitation. A non-standard cyan fluorescent protein, psamFP488, is characterized by an absorption maximum that is blue-shifted in comparison to other cyan fluorescent proteins. This protein is characterized by an extended Stokes shift and an ultrafast (170 fs) excited state proton transfer. In this work, a theoretical-computational study, including excited state ab initio dynamics, is performed to help understanding the reaction mechanism and propose new hypotheses on the role of the residues surrounding the chromophore. Our results suggest that the proton transfer could be indirect toward the acceptor (Glu167) and involves other residues surrounding the chromophore, despite the ultrafast kinetics.

Direct observation of the solvent organization and nuclear vibrations of [Ru(dcbpy)2(NCS)2]4-, [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], via ab initio molecular dynamics

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy)₂(NCS)₂]^4- [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C₂ symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies.

Energy, Structures, and Response Properties with a Fully Coupled QM/AMOEBA/ddCOSMO Implementation

We present the implementation of a fully coupled polarizable QM/MM/continuum model based on the AMOEBA polarizable force field and the domain decomposition implementation of the conductor-like screening model. Energies, response properties, and analytical gradients with respect to both QM and MM nuclear positions are available, and a generic, atomistic cavity can be employed. The model is linear scaling in memory requirements and computational cost with respect to the number of classical atoms and is therefore suited to model large, complex systems. Using three variants of the green-fluorescent protein, we investigate the overall computational cost of such calculations and the effect of the continuum model on the convergence of the computed properties with respect to the size of the embedding. We also demonstrate the fundamental role of polarization effects by comparing polarizable and nonpolarizable embeddings to fully QM ones.

Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy

The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.

Modeling Excited-State Proton Transfer to Solvent: A Dynamics Study of a Super Photoacid with a Hybrid Implicit/Explicit Solvent Model

The rapid growth of time-resolved spectroscopies and the theoretical advances in ab initio molecular dynamics (AIMD) pave the way to look at the real-time molecular motion following the electronic excitation. Here, we exploited the capabilities of AIMD combined with a hybrid implicit/explicit model of solvation to investigate the ultrafast excited-state proton transfer (ESPT) reaction of a super photoacid, known as QCy9, in water solution. QCy9 transfers a proton to a water solvent molecule within 100 fs upon the electronic excitation in aqueous solution, and it is the strongest photoacid reported in the literature so far. Because of the ultrafast kinetics, it has been experimentally hypothesized that the ESPT escapes the solvent dynamics control (Huppert et al., J. Photochem. Photobiol. A2014,277, 90). The sampling of the solvent configuration space on the ground electronic state is the first key step toward the simulation of the ESPT event. Therefore, several configurations in the Franck–Condon region, describing an average solvation, were chosen as starting points for the excited-state dynamics. In all cases, the excited-state evolution spontaneously leads to the proton transfer event, whose rate is strongly dependent on the hydrogen bond network around the proton acceptor solvent molecule. Our study revealed that the explicit representation at least of three solvation shells around the proton acceptor molecule is necessary to stabilize the excess proton. Furthermore, the analysis of the solvent molecule motions in proximity of the reaction site suggested that even in the case of the strongest photoacid, the ESPT is actually assisted by the solvation dynamics of the first and second solvation shells of the water accepting molecule.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Computational Spectroscopy in Solution by Integration of Variational and Perturbative Approaches on Top of Clusterized Molecular Dynamics

Multiscale QM/MM approaches have become the most suitable and effective methods for the investigation of spectroscopic properties of medium- or large-size chromophores in condensed phases. On these grounds, we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available tools. In the present paper, we report the latest developments of such an approach with specific reference to a general workplan starting with the addition of acetonitrile to the panel of solvents already available in the General Liquid Optimized Boundary (GLOB) model enforcing nonperiodic boundary conditions (NPBC). Next, the solvatochromic shifts induced by acetonitrile on both rigid (uracil and thymine) and flexible (thyrosine) chromophores have been studied introducing in our software a number of new features ranging from rigid-geometry NPBC molecular dynamics based on the quaternion formalism to a full integration of variational (ONIOM) and perturbative (perturbed matrix method (PMM)) approaches for describing different solute–solvent topologies and local fluctuations, respectively. Finally, thymine and uracil have been studied also in methanol to point out the generality of the computational strategy. While further developments are surely needed, the strengths of our integrated approach even in its present version are demonstrated by the accuracy of the results obtained by an unsupervised approach and coupled to a computational cost strongly reduced with respect to that of conventional QM/MM models without any appreciable accuracy deterioration.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

The Virtual Multifrequency Spectrometer: a new paradigm for spectroscopy

On going developments of hardware and software are changing computational spectroscopy from a strongly specialized research area to a general tool in the inventory of most researchers. Increased interactions between experimentally-oriented users and theoretically-oriented developers of new methods and models would result in more robust, flexible and reliable tools and studies for the systems of increasing complexity, which are of current scientific and technological interest. This is the philosophy behind this review, which presents the development of a so-called virtual multi-frequency spectrometer (VMS) including state-of-the-art approaches in a user-friendly frame. The current status of the VMS tool will be illustrated by a number of case studies with special reference to infrared and UV-vis regions of the electro-magnetic spectrum including also chiral spectroscopies. Only the basic theoretical background will be provided avoiding explicit equations as far as possible, and pointing out the most recent advancements beyond the standard rigid-rotor harmonic-oscillator model coupled to vertical electronic excitation energies.

Development and Validation of the B3LYP/N07D Computational Model for Structural Parameter and Magnetic Tensors of Large Free Radicals

Extensive calculations on a large set of free radicals containing atoms of the second and third row show that the B3LYP/N07D computational model provides remarkably accurate structural parameters and magnetic tensors at reasonable computational costs. The key of this success is the optimization of core-valence s functions for hyperfine coupling constants, while retaining (and even improving) the good performances of the parent 6-31+G(d,p) basis set for valence properties through reoptimization of polarization and diffuse p functions.

Combining the Fluctuating Charge Method, Non-Periodic Boundary Conditions and Meta-Dynamics: Aqua Ions as case studies

We present the current status of development of our code for performing Molecular Dynamics (MD) simulations exploiting a polarizable force field based on the Fluctuating Charge (FQ) method and non-Periodic Boundary Conditions (NPBC). Continuing on the path set in a previous work, we increased the capabilities of the code by implementing a number of new features, including: a non-iterative algorithm for rigid trigonal molecule simulations; two additional temperature coupling schemes; a meta-dynamics based approach for effective free energy evaluations. Although these are well known algorithms, each present in one or more widely used MD packages, they have now been tested, for the first time, in the context of the FQ model coupled with NPBC. As case studies, we considered three aqueous ions of increasing charge, namely Na⁺, Ca²⁺ and La³⁺, at infinite dilution. In particular, by exploiting a computational approach recently proposed by our group and based on the metadynamics technique, we focused on the important role played by solvent polarization on ionic hydration structures, also investigating the free energy landscapes of ion coordination and the water exchange rates. Such an approach, previously tested with standard non-polarizable models, was applied here to evaluate the effects of explicit polarization on water exchange barriers between different solvent coordination structures. Moreover, we have analyzed and discussed in some detail non-linear electrostatic effects arising from solvent polarization while going from a mono- to a di- and trivalent ion.

Vibrational analysis of x-ray absorption fine structure thermal factors by ab initio molecular dynamics: the Zn(II) ion in aqueous solution as a case study

In this work, we consider a new combination of vibrational analysis and normal-like mode decomposition of Debye-Waller factors of solvated ions entirely based on molecular dynamics data. Such a novel time-dependent analysis procedure provides a direct link between x-ray absorption fine structure parameters and normal mode contributions for an ion-solvent system. The potentialities of such a methodology rely on two fundamental aspects which distinguish it from already available tools. First, a general vibrational analysis that does not require any Gaussian or harmonic model for describing atomic fluctuations in liquids. Second, a very accurate sampling of the short range motions around the structural probe via the recently developed atom centered density matrix propagation/general liquid optimized boundary method. This novel molecular dynamics methodology is based on an integrated ab initio/classical potential using localized basis functions and nonperiodic boundary conditions. As a case study we have chosen the Zn(II) ion in aqueous solution. The consistency of our results and the observed good agreement with experiments show how the key support to advanced structural techniques from molecular dynamics can be further expanded and investigated.

Predicting the ¹H and ¹³C NMR spectra of paramagnetic Ru(III) complexes by DFT

Nuclear shieldings, including the Fermi-contact and pseudocontact terms, have been calculated with density functional theory (DFT) (nonrelativistic and relativistic) methods in several Ru(III) complexes, thereby predicting (1)H and (13)C paramagnetic shifts. A fair agreement with experimental values is observed. Structural, magnetic and dynamic parameters have also been input to the Solomon-Bloembergen equation in order to predict signal lineshapes. It is shown that DFT-predicted paramagnetic shifts can greatly aid in obtaining and understanding NMR spectra of paramagnetic Ru(III) complexes.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

Fully Integrated Approach to Compute Vibrationally Resolved Optical Spectra: From Small Molecules to Macrosystems

A general and effective time-independent approach to compute vibrationally resolved electronic spectra from first principles has been integrated into the Gaussian computational chemistry package. This computational tool offers a simple and easy-to-use way to compute theoretical spectra starting from geometry optimization and frequency calculations for each electronic state. It is shown that in such a way it is straightforward to combine calculation of Franck-Condon integrals with any electronic computational model. The given examples illustrate the calculation of absorption and emission spectra, all in the UV-vis region, of various systems from small molecules to large ones, in gas as well as in condensed phases. The computational models applied range from fully quantum mechanical descriptions to discrete/continuum quantum mechanical/molecular mechanical/polarizable continuum models.

Model systems for poly(acrylic acid) main-chain radicals based on the Kemp's triacid framework

Kemp's triacid (KTA) and cyclohexane tricarboxylic acid (CTA) are small-molecule model systems for acrylic acid polymers, having the same functionalities and stereoregularities as isotactic poly(methacrylic acid) (PMAA) and poly(acrylic acid) (PAA), respectively. As part of an ongoing investigation of radicals produced by photolysis of acrylic polymers, the photochemistry and free radicals from the model systems have been studied using time-resolved EPR spectroscopy as a function of temperature and pH. Radicals are created by direct photolysis of the acids at 248 nm or by sensitized photo-oxidation using quinone triplet states at 308 nm. The two methods of radical production lead to different chemically induced electron spin polarization (CIDEP) patterns in the ensuing radicals, which are simulated and discussed. Well-resolved spectra are obtained at all temperatures for the model system radicals, which are determined to be in the slow motion condition. DFT calculations of the model system radicals are presented and discussed in support of the experimental data.

Quantum mechanical computations and spectroscopy: from small rigid molecules in the gas phase to large flexible molecules in solution

Interpretation of structural properties and dynamic behavior of molecules in solution is of fundamental importance to understand their stability, chemical reactivity, and catalytic action. While information can be gained, in principle, by a variety of spectroscopic techniques, the interpretation of the rich indirect information that can be inferred from the analysis of experimental spectra is seldom straightforward because of the subtle interplay of several different effects, whose specific role is not easy to separate and evaluate. In such a complex scenario, theoretical studies can be very helpful at two different levels: (i) supporting and complementing experimental results to determine the structure of the target molecule starting from its spectral properties; (ii) dissecting and evaluating the role of different effects in determining the observed spectroscopic properties. This is the reason why computational spectroscopy is rapidly evolving from a highly specialized research field into a versatile and widespread tool for the assignment of experimental spectra and their interpretation in terms of chemical physical effects. In such a situation, it becomes important that both computationally and experimentally oriented chemists are aware that new methodological advances and integrated computational strategies are available, providing reliable estimates of fundamental spectral parameters not only for relatively small molecules in the gas phase but also for large and flexible molecules in condensed phases. In this Account, we review the most significant methodological contributions from our research group in this field, and by exploiting some recent results of their application to the computation of IR, UV-vis, NMR, and EPR spectral parameters, we discuss the microscopic mechanisms underlying solvent and vibrational effects on the spectral parameters. After reporting some recent achievements for the study of excited states by first principle quantum mechanical approaches, we focus on the treatment of environmental effects by means of mixed discrete-continuum solvent models and on effective methods for computing vibronic contributions to the spectra. We then discuss some new developments, mainly based on time-dependent approaches, allowing us to go beyond the determination of spectroscopic parameters toward the simulation of line widths and shapes. Although further developments are surely needed to improve the accuracy and effectiveness of several items in the proposed approach, we try to show that the first important steps toward a direct comparison between the results obtained in vitro and those obtained in silico have been made, making easier fruitful crossovers among experiments, computations and theoretical models, which would be decisive for a deeper understanding of the spectral behavior associated with complex systems and processes.

Unraveling Solvent-Driven Equilibria between α- and 310-Helices through an Integrated Spin Labeling and Computational Approach

In this work we present an effective and flexible computational approach, which is the result of an ongoing development in our groups, allowing the complete a priori simulation of the ESR spectra of complex systems in solution. The usefulness and reliability of the method are demonstrated on the very demanding playground represented by the tuning of the equilibrium between 3(10)- and alpha-helices of polypeptides by different solvents. The starting point is the good agreement between computed and X-ray diffraction structures for the 3(10)-helix adopted by the double spin-labelled heptapeptide Fmoc-(Aib-Aib-TOAC)2-Aib-OMe. Next, density functional computations, including dispersion interactions and bulk solvent effects, suggest another energy minimum corresponding to an alpha-helix in polar solvents, which, eventually, becomes the most stable structure. Computation of magnetic and diffusion tensors provides the basic ingredients for the building of complete spectra by methods rooted in the Stochastic Liouville Equation (SLE). The remarkable agreement between computed and experimental spectra at different temperatures allowed us to identify helical structures in the various solvents. The generality of the computational strategy and its implementation in effective and user-friendly computer codes pave the route toward systematic applications in the field of biomolecules and other complex systems.