Solvent Effects on the UV (n → π*) and NMR (13C and 17O) Spectra of Acetone in Aqueous Solution. An Integrated Car−Parrinello and DFT/PCM Approach

Accurate Thermochemical and Kinetic Parameters at Affordable Cost by Means of the Pisa Composite Scheme (PCS)

A new strategy for the computation at an affordable cost of geometrical structures, thermochemical parameters, and rate constants for medium-sized molecules in the gas phase is proposed. The most distinctive features of the new model are the systematic use of cc-pVnZ-F12 basis sets, the addition of MP2 core-valence correlation in geometry optimizations by a double-hybrid functional, the separate extrapolation of MP2 and post-MP2 contributions, and the inclusion of anharmonic contributions in zero-point energies and thermodynamic functions. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme outperforms the most well-known model chemistries without the need for any empirical parameter. Additional tests show that the computed zero-point energies and thermal contributions can be confidently used for obtaining accurate thermochemical and kinetic parameters. Since the whole computational workflow is translated in a black-box procedure, which can be followed with standard electronic structure codes, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.

Development, Validation, and Pilot Application of a Generalized Fluctuating Charge Model for Computational Spectroscopy in Solution

A general approach enforcing nonperiodic boundary conditions for the computation of spectroscopic properties in solution has been improved including an effective description of charge-transfer contributions and coordination number adjustment for explicit solvent molecules. Both contributions are obtained from a continuous description of intermolecular hydrogen bonds, which has been employed also for an effective clustering of molecular dynamics trajectories. Fine tuning of the model has been performed for several water clusters, and then its efficiency and reliability have been demonstrated by computing the absorption spectra of different creatinine tautomers in aqueous solution.

Machine Learning of Quasiparticle Energies in Molecules and Clusters

We present a Δ-machine learning approach for the prediction of GW quasiparticle energies (ΔMLQP) and photoelectron spectra of molecules and clusters, using orbital-sensitive representations (OSRs) based on molecular Cartesian coordinates in kernel ridge regression-based supervised learning. Coulomb matrix, bag-of-bond, and bond-angle-torsion representations are made orbital-sensitive by augmenting them with atom-centered orbital charges and Kohn-Sham orbital energies, both of which are readily available from baseline calculations at the level of density functional theory (DFT). We first illustrate the effects of different constructions of the OSRs on the prediction of frontier orbital energies of 22k molecules of the QM8 data set and show that it is possible to predict the full photoelectron spectrum of molecules within the data set using a single model with a mean absolute error below 0.1 eV. We further demonstrate that the OSR-based ΔMLQP captures the effects of intra- and intermolecular conformations in application to water monomers and dimers. Finally, we show that the approach can be embedded in multiscale simulation workflows, by studying the solvatochromic shifts of quasiparticle and electron-hole excitation energies of solvated acetone in a setup combining molecular dynamics, DFT, the GW approximation, and the Bethe-Salpeter equation. Our findings suggest that the ΔMLQP model allows us to predict quasiparticle energies and photoelectron spectra of molecules and clusters with GW accuracy at DFT cost.

Molecular and Supramolecular Interactions in Systems with Nitroxide-Based Radicals

Nitroxide-based radicals, having the advantage of firm chemical stability, are usable as probes in the detection of nanoscale details in the chemical environment of various multi-component systems, based on subtle variations in their electron paramagnetic resonance spectra. We propose a systematic walk through the vast area of problems and inquires that are implied by the rationalization of solvent effects on the spectral parameters, by first-principle methods of structural chemistry. Our approach consists of using state-of-the-art procedures, like Density Functional Theory (DFT), on properly designed systems, kept at the border of idealization and chemical realism. Thus, we investigate the case of real solvent molecules intervening in different configurations between two radical molecules, in comparison with radicals taken in vacuum or having the solvent that is treated by surrogate models, such as polarization continuum approximation. In this work, we selected the dichloromethane as solvent and the prototype radicals abbreviated TEMPO ((2,2,6,6-Tetramethylpiperidin-1-yl) oxyl). In another branch of the work, we check the interaction of radicals with large toroidal molecules, β-cyclodextrin, and cucurbit[6]uril, modeling the interaction energy profile at encapsulation. The drawn synoptic view offers valuable rationales for understanding spectroscopy and energetics of nitroxide radicals in various environments, which are specific to soft chemistry.

Decomposition of molecular properties

We review recent work on property decomposition techniques using quantum chemical methods and discuss some topical applications in terms of quantum mechanics-molecular mechanics calculations and the constructing of properties of large molecules and clusters. Starting out from the so-called LoProp decomposition scheme [Gagliardi et al., J. Chem. Phys., 2004, 121, 4994] for extracting atomic and inter-atomic contributions to molecular properties we show how this method can be generalized to localized frequency-dependent polarizabilities, to localized hyperpolarizabilities and to localized dispersion coefficients. Some applications of the generalized decomposition technique are reviewed - calculations of frequency-dependent polarizabilities, Rayleigh scattering of large clusters, and calculations of hyperpolarizabilities of proteins.

Stability of melamine-exfoliated graphene in aqueous media: quantum-mechanical insights at the nanoscale

In recent experiments, melamine (1,3,5-triazine-2,4,6-triamine) has been proposed as an effective exfoliating agent to obtain high quality graphene from graphite. After washing out the melamine in excess, small amounts (ppm) are still needed to stabilize the dispersion of graphene flakes in aqueous media. To understand the origin of this behaviour, we investigated the melamine-graphene-water system and the fundamental interactions that determine its structure and energetics. To disentangle the subtle interplay of hydrogen-bonding and dispersive forces we used state-of-the-art ab initio calculations based on density functional theory. First, we focused on the case of water molecules interacting with melamine-graphene assemblies at different melamine coverages. We found that water-melamine interactions provide the driving force for washing off the melamine from graphene. Then, we addressed the interaction of single and double layers of water molecules with the graphene surface in the presence of an adsorbed melamine molecule. We found that this melamine acts as a non-covalent anchor for keeping a number of water molecules conveniently close to the graphene surface, thus helping its stabilization in aqueous media. Our analysis helps understanding how competing weak forces can lead to a stable graphene water suspension thanks to small amounts of adsorbed melamine. From our results, we derive simple indications on how the water-graphene interfacial properties can be tuned via non-covalent adsorption of small functional molecules with H-bond donor/acceptor groups. These new hints can be helpful to prepare stable graphene dispersions in water and so to unlock graphene potential in aqueous environments.

Tuning of dye optical properties by environmental effects: a QM/MM and experimental study

The present work is aimed at a deeper investigation of two recently synthesized heteroaromatic fluorophores by means of a computational multilayer approach, integrating quantum mechanics (QM) and molecular mechanics (MM).

The present work is aimed at a deeper investigation of two recently synthesized heteroaromatic fluorophores by means of a computational multilayer approach, integrating quantum mechanics (QM) and molecular mechanics (MM). In particular, dispersion of the title dyes in a polymer matrix is studied in connection with potential applications as photoactive species in luminescent solar concentrators (LSCs). Molecular dynamics simulations, based on accurate QM-derived force fields, reveal increased stiffness of these organic dyes when going from CHCl₃ solution to the polymer matrix. QM/MM computations of UV spectra for snapshots extracted from MD simulations show that this different flexibility permits explaining the different spectral shapes obtained experimentally for the two different environments. Moreover, the general spectroscopic trends are reproduced well by static computations employing a polarizable continuum description of environmental effects.

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.

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.

Essential dynamics for the study of microstructures in liquids

Essential Dynamics (ED) is a powerful tool for analyzing molecular dynamics (MD) simulations and it is widely adopted for conformational analysis of large molecular systems such as, for example, proteins and nucleic acids. In this study, we extend the use of ED to the study of clusters of arbitrary size constituted by weakly interacting particles, for example, atomic clusters and supramolecular systems. The key feature of the method we present is the identification of the relevant atomic-molecular clusters to be analyzed by ED for extracting the information of interest. The application of this computational approach allows a straightforward and unbiased conformational study of the local microstructures in liquids, as emerged from semiclassical MD simulations. The good performance of the method is demonstrated by calculating typical observables of liquid water, that is, NMR, NEXAFS O1s, and IR spectra, known to be rather sensitive both to the presence and to the conformational features of hydrogen-bonded clusters.

π-Stacking effects on the EPR parameters of a prototypical DNA spin label

The character and value of spin labels for probing environments like double-stranded DNA depend on the degree of change in the spin Hamiltonian parameters of the spin label induced by the environment. Herein we provide a systematic theoretical investigation of this issue, based on a density functional theory method applied to a spin labeled DNA model system, focusing on the dependence of the EPR properties of the spin label on the π stacking and hydrogen bonding that occur upon incorporating the spin label into the selected base pair inside DNA. It is found that the EPR spin Hamiltonian parameters of the spin label are only negligibly affected by its incorporation into DNA, when compared to its free form. This result gives a theoretical ground for the common empirical assumption regarding the behaviour of spin Hamiltonian parameters made in EPR based measurements of the distance between spin labels incorporated into DNA.

SINGLE-PARTICLE AND PAIR DYNAMICAL PROPERTIES OF ACETONE–METHANOL MIXTURES CONTAINING CHARGED AND NEUTRAL SOLUTES: A MOLECULAR DYNAMICS STUDY

The dynamical properties of acetone–methanol mixtures containing either an ionic or a neutral hydrophobic solute are investigated by means of a series of molecular dynamics simulations. The primary goal has been to study how the solute and solvent dynamical properties change with variation of composition of the mixture ranging from pure acetone to pure methanol. The variations of structure and energetics of the mixture with composition are also calculated. The diffusion coefficients of both ionic and neutral solutes are found to show nonlinear variation with composition of the mixture, although the extent of nonlinearity in the diffusion of the neutral solute is much weaker. Calculations of appropriate solute-solvent distribution functions reveal the extent and nature of selective solvation of these solute species which play a role in determining the nonideal dynamical characteristics of these solutes. The free energies of solvation of the ionic solutes are also calculated and the results are discussed in the context of their dynamical behavior. The hydrogen bond statistics and dynamics of these mixtures are also calculated over their entire composition range. The energies and lifetimes of hydrogen bonds between an acetone and a methanol molecule or between two methanol molecules are found to increase with increase of acetone mole fraction of the mixture. Residence times of methanol molecules in solvation shells of acetone and methanol are also found to follow the same trend as relaxation times. However, these pair dynamical properties show essentially linear dependence on composition, thus behave almost ideally with respect to changes in composition of the mixture.

UV-vis spectra of the anticancer camptothecin family drugs in aqueous solution: specific spectroscopic signatures unraveled by a combined computational and experimental study

The ultraviolet-visible absorption spectrum of camptothecin (CPT) has been been recorded in aqueous solution at pH 5.3, where the equilibrium among the different CPT forms is shifted toward the lactonic one. Time-dependent density functional theory (TD-DFT) computations lead to a remarkable reproduction of the experimental spectrum only upon addition of explicit water molecules in interaction with specific moieties of the camptothecin molecule. Molecular dynamics (MD) simulations enforcing boundary periodic conditions for CPT embedded with 865 water molecules, with a force field derived from DFT computations, show that the experimental spectrum is due to the contributions of CPT molecules with different solvation patterns. A similar solvent effect is observed for several CPT derivatives, including the clinically relevant SN-38 and topotecan drugs. The quantitative agreement between TD-DFT/MD computations and experimental data allow us to identify specific spectroscopic signatures diagnostic of the drug environment and to develop procedures that can be used to monitor the drug-DNA/protein interaction.