Explicit Modelling of Spectral Bandshapes by a Mixed Quantum-Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distyrylbenzene.

The absorption and emission spectral shapes of a flexible organic probe, the distyrylbenzene (DSB) dye, are simulated accounting for the effect of different environments of increasing complexity, ranging from a homogeneous, low-molecular- weight solvent, to a long-chain alkane, and, eventually, a channel-forming organic matrix. Each embedding is treated explicitly, adopting a mixed quantum-classical approach, the Adiabatic Molecular Dynamics - generalized vertical Hessian (Ad-MD|gVH) model, which allows a direct simulation of the environment-induced constraining effects on the vibronic spectral shapes. In such a theoretical framework, the stiff modes of the dye are described at a quantum level within the harmonic approximation, including Duschinsky mixing effects, while flexible degrees of freedom of the solute (e. g. torsions) and those of the solvent are treated classically by means of molecular dynamics sampling. Such a setup is shown to reproduce the distinct effects exerted by the different environments in varied thermodynamic conditions. Besides allowing for a first-principles rationale on the supramolecular mechanism leading to the experimental spectral features, this result represents the first successful application of the Ad-MD|gVH method to complex embeddings and supports its potential application to other heterogeneous environments, such as for instance, pigment-protein complexes or organic dyes adsorbed into metal-organic frameworks.

Dynamic effects on the nonlinear optical properties of donor acceptor stenhouse adducts: insights from combined MD + QM simulations.

The second-order nonlinear optical (NLO) responses of a donor-acceptor stenhouse adduct (DASA) are investigated by using a computational approach combining molecular dynamics simulations and density functional theory (DFT) calculations. Specific force fields for the open and closed photoswitching forms are first parameterized and validated according to the Joyce protocol, in order to finely reproduce the geometrical features and potential energy surfaces of both isomers in chloroform solution. Then, DFT calculations are performed on structural snapshots extracted at regular time steps of the MD trajectories to address the influence of the thermalized conformational dynamics on the NLO responses related to hyper-Rayleigh scattering (HRS) experiments. We show that accounting for the structural dynamics largely enhances the HRS hyperpolarizability (ßHRS) compared to DFT calculations considering solely equilibrium geometries, and greatly improves the agreement with experimental measurements. Furthermore, we show that the NLO responses of the NLO-active open form are correlated with the bond order alternation along the triene bridge connecting the donor and acceptor moieties, which is rationalized using simple essential state models.

A Portrait of the Chromophore as a Young System-Quantum-Derived Force Field Unraveling Solvent Reorganization upon Optical Excitation of Cyclocurcumin Derivatives.

The study of fast non-equilibrium solvent relaxation in organic chromophores is still challenging for molecular modeling and simulation approaches, and is often overlooked, even in the case of non-adiabatic dynamics simulations. Yet, especially in the case of photoswitches, the interaction with the environment can strongly modulate the photophysical outcomes. To unravel such a delicate interplay, in the present contribution we resorted to a mixed quantum-classical approach, based on quantum mechanically derived force fields. The main task is to rationalize the solvent reorganization pathways in chromophores derived from cyclocurcumin, which are suitable for light-activated chemotherapy to destabilize cellular lipid membranes. The accurate and reliable decryption delivered by the quantum-derived force fields points to important differences in the solvent's reorganization, in terms of both structure and time scale evolution.

Solvent effects on catechol's binding affinity: investigating the role of the intra-molecular hydrogen bond through a multi-level computational approach.

The subtle interplay between the inter-molecular interactions established by catechol with the surrounding solvent and the intra-molecular hydrogen bond (HB) characterizing its conformational dynamics is investigated through a multi-level computational approach. First, quantum mechanical (QM) calculations are employed to accurately characterize both large portions of the catechol's potential energy surface and the interaction energy with neighboring solvent molecules. The acquired information is thereafter exploited to develop a QM derived force-field (QMD-FF), in turn employed in molecular dynamics (MD) simulations based on classical mechanics. The reliability of the QMD-FF is further validated through a comparison with the outcomes of ab initio molecular dynamics, also purposely carried out in this work. In agreement with recent experimental findings, the MD results reveal remarkable differences in the conformational behavior of isolated and solvated catechol, as well as among the investigated solvents, namely water, acetonitrile or cyclohexane. The rather strong intramolecular HB, settled between the vicinal phenolic groups and maintained in the gas phase, loses stability when catechol is solvated in polar solvents, and is definitively lost in protic solvents such as water. In fact, the internal energy increase associated with the rotation of one hydroxyl group and the breaking of the internal HB is well compensated by the intermolecular HB network available when both phenolic hydrogens point toward the surrounding solvent. In such a case, catechol is stabilized in a chelating conformation, which in turn could be very effective in water removal and surface anchoring. Besides unraveling the role of the different contributors that govern catechol's conformational dynamics, the QMD-FF developed in this work could be in future employed to model larger catechol containing molecules, due to its accuracy to reliably model both internal flexibility and solvent effects, while exploiting MD computational benefits to include more complex players as for instance surfaces, ions or biomolecules.

Unraveling the contributions to the spectral shape of flexible dyes in solution: insights on the absorption spectrum of an oxyluciferin analogue.

We present a computational investigation of the absorption spectrum in water of 5,5-spirocyclopropyl-oxyluciferin (5,5-CprOxyLH), an analogue of the emitter compound responsible for the bioluminescence in fireflies. Several factors participate in determining the 5,5-CprOxyLH's spectral shape: (i) the contribution of the four close-energy excited states, which show significant non-adiabatic couplings, (ii) the flexible molecular structure and (iii) the specific interactions established with the surrounding environment, which strongly couple the protic solvent dynamics with the dye's spectral response. To tackle the challenge to capture and dissect the role of all these effects we preliminarily investigate the role of non-adiabatic couplings with quantum dynamics simulations and a linear vibronic coupling model in the gas phase. Then, we account for both the molecular flexibility and solvent interactions by resorting to a mixed quantum classical protocol, named Adiabatic Molecular Dynamics generalized Vertical Gradient (Ad-MD|gVG), which is built on a method recently proposed by some of us. It is rooted in the partition between stiff degrees of freedom of the dye, accounted for at the vibronic level within the harmonic approximation, and flexible degrees of freedom of the solute (and of the solvent), described classically through a sampling based on Molecular Dynamics (MD). Ad-MD|gVG avoids spurious effects arising in the excited state Hessians due to non-adiabatic couplings, and can therefore be applied to account for the contributions of the first four excited states to the 5,5-CprOxyLH absorption spectrum. The final simulated spectrum is in very good agreement with the experiment, especially when the MD is driven by a refined quantum-mechanically derived force-field. More importantly, the origin of each separate contribution to the spectral shape is appropriately accounted for, paving the way to future applications of the method to more complex systems or alternative spectroscopies, as emission or circular dichroism.

Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol.

We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.

Flexible Superlubricity Unveiled in Sidewinding Motion of Individual Polymeric Chains.

A combination of low temperature atomic force microcopy and molecular dynamic simulations is used to demonstrate that soft designer molecules realize a sidewinding motion when dragged over a gold surface. Exploiting their longitudinal flexibility, pyrenylene chains are indeed able to lower diffusion energy barriers via on-surface directional locking and molecular strain. The resulting ultralow friction reaches values on the order of tens of pN reported so far only for rigid chains sliding on an incommensurate surface. Therefore, we demonstrate how molecular flexibility can be harnessed to realize complex nanomotion while retaining a superlubric character. This is in contrast with the paradigm guiding the design of most superlubric nanocontacts (mismatched rigid contacting surfaces).

Complexes of Alkaline and Ammonium Cations with Dopamine and Eumelanin Precursors: Dissecting the Role of Noncovalent Cation-π and Cation-Lone Pair (σ-Type) Interactions.

Cation-π interactions and their possible competition with other noncovalent interactions (NCI) might play a key role in both dopamine- and eumelanin-based bioinspired materials. In this contribution, to unravel the delicate interplay between cation-π interactions and other possible competing forces, the configurational space of noncovalent complexes formed by dopamine or eumelanin precursors (o-benzoquinone, DHI and a semiquinone dimer) and three different cations (Na+, K+, and NH4+) is sampled by means of accurate ab initio calculations. To this end, we resort to the mp2mod method, recently validated by us for benzene-, phenol-, and catechol-cation complexes, whose computational convenience allows for an extensive exploration of the cation-molecule interaction energy surface, by sampling a total of more than 104 arrangements. The mp2mod interaction energy landscapes reveal that, besides the expected cation-π driven arrangements, for all considered molecule-cation pairs the most stable complexes are found when the cation lies within the plane containing the six-membered ring, thus maximizing the σ-type interaction with the oxygen's lone pairs. Due to the loss of aromaticity, the σ-type/cation-π strength ratio is remarkably large in o-benzoquinone, where cation-π complexes seem unlikely to be formed. The above features are shared among all considered cations but are significantly larger when considering the smaller Na+. Besides delivering a deeper insight onto the NCI network established by the considered precursors in the presence of ions, the present results can serve as a reference database to validate or refine lower level methods, as, for instance, the force fields employed in classical simulations.

Sequential Bending and Twisting around C-C Single Bonds by Mechanical Lifting of a Pre-Adsorbed Polymer.

Bending and twisting around carbon-carbon single bonds are ubiquitous in natural and synthetic polymers. Force-induced changes were so far not measured at the single-monomer level, owing to limited ways to apply local forces. We quantified down to the submolecular level the mechanical response within individual poly-pyrenylene chains upon their detachment from a gold surface with an atomic force microscope at 5 K. Computer simulations based on a dedicated force field reproduce the experimental traces and reveal symmetry-broken bent and rotated conformations of the sliding physisorbed segment besides steric hindrance of the just lifted monomer. Our study also shows that the tip-molecule bond remains intact but remarkably soft and links force variations to complex but well-defined conformational changes.

The Role of Non-Covalent Interactions on Cluster Formation: Pentamer, Hexamers and Heptamer of Difluoromethane.

The role of non-covalent interactions (NCIs) has broadened with the inclusion of new types of interactions and a plethora of weak donor/acceptor partners. This work illustrates the potential of chirped-pulse Fourier transform microwave technique, which has revolutionized the field of rotational spectroscopy. In particular, it has been exploited to reveal the role of NCIs' in the molecular self-aggregation of difluoromethane where a pentamer, two hexamers and a heptamer were detected. The development of a new automated assignment program and a sophisticated computational screening protocol was essential for identifying the homoclusters in conditions of spectral congestion. The major role of dispersion forces leads to less directional interactions and more distorted structures than those found in polar clusters, although a detailed analysis demonstrates that the dominant interaction energy is the pairwise interaction. The tetramer cluster is identified as a structural unit in larger clusters, representing the maximum expression of bond between dimers.

The phenoxyl group-modulated interplay of cation-π and σ-type interactions in the alkali metal series.

The interaction potential energy surfaces (IPESs) of four alkaline metal cations (Na+, K+, Rb+ and Cs+) complexed with phenol and catechol were explored by accurate ab initio calculations to investigate the interplay of different noncovalent interactions and their behavior along the alkali metal series and upon -OH substitution. Selected one-dimensional interaction energy curves revealed two different minimum energy configurations for all phenol- and catechol-metal complexes, characterized either by cation-π or σ-type interactions. For each investigated complex several two-dimensional IPES maps were also computed, exploiting the computational advantages of the MP2mod approach. The size of the alkali cation was found to play a similar role in modulating both kinds of complexes, as the interaction strength always decreases along the metal series, from Na+ to Cs+. Conversely, the number of hydroxyl substituents markedly affected cation-π complexes vs. σ-type ones. As a most relevant finding, in catechol-metal complexes the strength of cation-π interactions is around half that of the σ-type ones. It is argued that the combined effect of cation dimensions and hydroxyl substitution in catechol-Na+ complexes makes σ-type configurations remarkably more stable and easily accessible than cation-π ones. Besides shedding new light on the origin of biological phenomena connected with underwater adhesion, the quantum mechanical interaction energy database provided herein may offer a useful reference for tuning accurate force fields, suitable for molecular dynamics simulations, where environmental effects might be also taken into account.

Benchmarking Cation-π Interactions: Assessment of Density Functional Theory and Möller-Plesset Second-Order Perturbation Theory Calculations with Optimized Basis Sets (mp2mod) for Complexes of Benzene, Phenol, and Catechol with Na+, K+, Rb+, and Cs.

The reliability of several density functional theory (DFT) functionals and of the Möller-Plesset second-order perturbation theory calculations with modified basis sets (mp2mod) approach in describing cation-π interactions is systematically investigated by benchmarking their performances with respect to high quality reference CCSD(T) calculations of the binding energies between alkaline cations of varying radius (Na+, K+, Rb+, and Cs+) and three aromatic species (benzene, phenol, and catechol). For this class of noncovalent interaction, mp2mod delivers, on average, results in very good agreement with the reference CCSD(T) data, yet at a very small computational cost, exploiting the reduced dimensions of the modified basis set. Conversely, the tested DFT functionals show a more erratic behavior, with different performances depending on both the investigated system and the combination of the employed functional and basis set. The mp2mod computational convenience is further exploited to extensively sample two-dimensional interaction energy surfaces of all investigated cation-π systems, which allow for a deeper insight on the effect of the increasing number of hydroxyl substituents, revealing the insurgence, upon substitution, of alternative minima, evident in particular for the smaller cations. The present results strongly support for further applications of the mp2mod method to study a larger variety of aromatic/metal cation species, relevant both in biological processes and in technological applications.

Automated parameterization of quantum-mechanically derived force-fields including explicit sigma holes: A pathway to energetic and structural features of halogen bonds in gas and condensed phase.

In classical molecular dynamics, general purpose atomistic force-fields (FFs) often deliver inaccurate results when dealing with halogen bonds (XBs), notwithstanding their crucial role in many fields of science, ranging from material design to drug development. Given the large dimensions of the systems of interest, it would be therefore desirable to increase the FF accuracy maintaining the simplicity of the standard Lennard-Jones (LJ) plus point charge description to avoid an excessive computational cost. A simple yet effective strategy consists in introducing a number of virtual sites able to mimic the so-called "explicit σ-hole." In this work, we present an automated FF parameterization strategy based on a global optimization of both LJ and charge parameters with respect to accurate quantum mechanical data, purposely computed for the system under investigation. As a test case, we report on two homologue series, characterized either by weak or strong XBs, namely, the di-halogenated methanes and the mono-, di-, and tri-substituted acetonitriles, taking into consideration Cl, Br, and I substituents. The resulting quantum mechanically derived FFs are validated for each compound in the gas and in the condensed phase by comparing them to general purpose and specific FFs without virtual sites and to highly accurate reference quantum mechanical data. The results strongly support the adoption of the specific FFs with virtual sites, which overcome the other investigated models in representing both gas phase energetics and the structural patterns of the liquid phase structure related to the presence of XBs.

Iron's Wake: The Performance of Quantum Mechanical-Derived Versus General-Purpose Force Fields Tested on a Luminescent Iron Complex.

Recently synthetized iron complexes have achieved long-lived excited states and stabilities which are comparable, or even superior, to their ruthenium analogues, thus representing an eco-friendly and cheaper alternative to those materials based on rare metals. Most of computational tools which could help unravel the origin of this large efficiency rely on ab-initio methods which are not able, however, to capture the nanosecond time scale underlying these photophysical processes and the influence of their realistic environment. Therefore, it exists an urgent need of developing new low-cost, but still accurate enough, computational methodologies capable to deal with the steady-state and transient spectroscopy of transition metal complexes in solution. Following this idea, here we focus on the comparison between general-purpose transferable force-fields (FFs), directly available from existing databases, and specific quantum mechanical derived FFs (QMD-FFs), obtained in this work through the Joyce procedure. We have chosen a recently reported FeIII complex with nanosecond excited-state lifetime as a representative case. Our molecular dynamics (MD) simulations demonstrated that the QMD-FF nicely reproduces the structure and the dynamics of the complex and its chemical environment within the same precision as higher cost QM methods, whereas general-purpose FFs failed in this purpose. Although in this particular case the chemical environment plays a minor role on the photo physics of this system, these results highlight the potential of QMD-FFs to rationalize photophysical phenomena provided an accurate QM method to derive its parameters is chosen.

Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles.

The Spin Component Scaled (SCS) MP2 method using a reduced and optimized basis set (SCS-MP2mod) is employed to compute the interaction energies of nine homodimers, formed by aromatic heterocyclic molecules (pyrrole, furan, thiophene, oxazole, isoxazole, pyridine, pyridazine, pyrimidine, and pyrazine). The coefficients of the same-spin and opposite-spin correlation energies and the Gaussian type orbitals (GTO) polarization exponents of the 6-31G** basis set are simultaneously optimized in order to minimize the energy differences with respect to the coupled-cluster with single, double and perturbative triples excitations [CCSD(T)] reference interaction energies, extrapolated to a complete basis set. It is demonstrated that the optimization of the spin scale factors leads to a noticeable improvement of the accuracy with a root mean square deviation less than 0.1 kcal/mol and a largest unsigned deviation smaller than 0.25 kcal/mol. The pyrrole dimer provides an exception, with a slightly higher deviation from the reference data. Given the high benefit in terms of computational time with respect to the CCSD(T) technique and the small loss of accuracy, the SCS-MP2mod method appears to be particularly suitable for extensive sampling of intermolecular potential energy surfaces at a quantum mechanical level. Within this framework, a transferability test of the SCS-MP2mod parameters to a benchmark set of this class of molecules is very promising as the reference interaction energies of several heterocyclic aromatic heterodimers were reproduced with a standard deviation of 0.30 kcal/mol. The SCS-MP2mod remarkably outperforms the value of 1.95 kcal/mol obtained with standard MP2/6-31G**.

A General Treatment to Study Molecular Complexes Stabilized by Hydrogen-, Halogen-, and Carbon-Bond Networks: Experiment and Theory of (CH2 F2 )n ⋠⋠⋠(H2 O)m.

Rotational spectra of several difluoromethane-water adducts have been observed using two broadband chirped-pulse Fourier-transform microwave (CP-FTMW) spectrometers. The experimental structures of (CH2 F2 )⋅⋅⋅(H2 O)2 , (CH2 F2 )2 ⋅⋅⋅(H2 O), (CH2 F2 )⋅⋅⋅(H2 O)3 , and (CH2 F2 )2 ⋅⋅⋅(H2 O)2 were unambiguously identified with the aid of 18 isotopic substituted species. A subtle competition between hydrogen, halogen, and carbon bonds is observed and a detailed analysis was performed on the complex network of non-covalent interactions which stabilize each cluster. The study shows that the combination of stabilizing contact networks is able to reinforce the interaction strength through a cooperative effect, which can lead to large stable oligomers.

Toward a general mixed quantum/classical method for the calculation of the vibronic ECD of a flexible dye molecule with different stable conformers: Revisiting the case of 2,2,2-trifluoro-anthrylethanol.

We extend a recently proposed mixed quantum/classical method for computing the vibronic electronic circular dichroism (ECD) spectrum of molecules with different conformers, to cases where more than one hindered rotation is present. The method generalizes the standard procedure, based on the simple Boltzmann average of the vibronic spectra of the stable conformers, and includes the contribution of structures that sample all the accessible conformational space. It is applied to the simulation of the ECD spectrum of (S)-2,2,2-trifluoroanthrylethanol, a molecule with easily interconvertible conformers, whose spectrum exhibits a pattern of alternating positive and negative vibronic peaks. Results are in very good agreement with experiment and show that spectra averaged over all the sampled conformational space can deviate significantly from the simple average of the contributions of the stable conformers. The present mixed quantum/classical method is able to capture the effect of the nonlinear dependence of the rotatory strength on the molecular structure and of the anharmonic couplings among the modes responsible for molecular flexibility. Despite its computational cost, the procedure is still affordable and promises to be useful in all cases where the ECD shape arises from a subtle balance between vibronic effects and conformational variety.

Quantitative prediction and interpretation of spin energy gaps in polyradicals: the virtual magnetic balance.

Open-shell organic molecules possessing more than two unpaired electrons and sufficient stability even at room temperature are very unusual, but few were recently synthesized that promise a number of fascinating applications. Unfortunately, reliable structural information is not available and only lower limits can be estimated for energy splittings between the different spin states. On these grounds, we introduce here an effective 'virtual magnetic balance', a robust and user-friendly tool purposely tailored for polyradicals and devised to be used in parallel with experimental studies. The main objective of this tool is to provide reliable structures and quantitative splittings of spin states of large, complex molecules. We achieved this objective with reasonable computation times and in a theoretical framework that allows disentanglement of different stereo-electronic effects contributing to the overall experimental result. A recently synthesized tetraradical with remarkable chemical stability was used as a case study.

Magnetic gaps in organic tri-radicals: From a simple model to accurate estimates.

The calculation of the energy gap between the magnetic states of organic poly-radicals still represents a challenging playground for quantum chemistry, and high-level techniques are required to obtain accurate estimates. On these grounds, the aim of the present study is twofold. From the one side, it shows that, thanks to recent algorithmic and technical improvements, we are able to compute reliable quantum mechanical results for the systems of current fundamental and technological interest. From the other side, proper parameterization of a simple Hubbard Hamiltonian allows for a sound rationalization of magnetic gaps in terms of basic physical effects, unraveling the role played by electron delocalization, Coulomb repulsion, and effective exchange in tuning the magnetic character of the ground state. As case studies, we have chosen three prototypical organic tri-radicals, namely, 1,3,5-trimethylenebenzene, 1,3,5-tridehydrobenzene, and 1,2,3-tridehydrobenzene, which differ either for geometric or electronic structure. After discussing the differences among the three species and their consequences on the magnetic properties in terms of the simple model mentioned above, accurate and reliable values for the energy gap between the lowest quartet and doublet states are computed by means of the so-called difference dedicated configuration interaction (DDCI) technique, and the final results are discussed and compared to both available experimental and computational estimates.

Perturbative multireference configuration interaction (CI-MRPT2) calculations in a focused dynamical approach: a computational study of solvatochromism in pyrimidine.

We have investigated solvatochromic effects over a solvent series of increasing polarity on the prototype molecule pyrimidine as a solute species. The line shape profiles, obtained by a time-dependent approach based on quantum mechanical calculations performed over frames sampled from classical molecular dynamics trajectories, were directly compared to the available experimental bands. The multireference configuration interaction second-order perturbation (CI-MRPT2) calculations are in quantitative agreement with the experiment. The results also confirm how nonprotic solvents can be confidently modeled by continuous solvation models as the polarizable continuum model, whereas protic solvents, as water, require the inclusion of explicit solvent molecules to account for the effects of hydrogen bonds.