Structural and Vibrational Properties of Carboxylates Intercalated into Layered Double Hydroxides: A Joint Computational and Experimental Study

Layered double hydroxides (LDHs) are fascinating clay-like materials that display versatile properties, making them an extremely fertile playground for diverse applications, ranging from bio-compatible materials to the pharmaceutical industry to catalysis and photocatalysis. When intercalating organic and bio-organic species between the inorganic layers, such materials are named hybrid LDHs. The structure-property relation in these systems is particularly relevant, since most of the properties of the materials may be fine-tuned if a comprehensive understanding of the microscopic structure in the interlamellar space is achieved, especially with respect to the reorganization under water uptake (swelling). In this work, we combined experiments and simulations to rationalize the behavior of LDHs intercalating three carboxylates, the general structure of which can be given as [Mg4Al2(OH)12]A2-·XH2O (with A2- = succinate, aspartate, or glutamate and X representing increasing water content). Following this strategy, we were able to provide an interpretation of the different shapes observed for the experimental water adsorption isotherms and for the evolution of the infrared carboxylate band of the anions. Apart from small differences, due to the different reorganization of the conformational space under confinement, the behavior of the two amino acids is very similar. However, such behavior is quite different in the case of succinate. We were able to describe the different response of the anions, which has a significant impact on the isotherm and on the size of the interlamellar region, in terms of a different interaction mechanism with the inorganic layer.

Hydration effects on the vibrational properties of carboxylates: From continuum models to QM/MM simulations

The presence of carboxyl groups in a molecule delivers an affinity to metal cations and a sensitivity to the chemical environment, especially for an environment that can give rise to intermolecular hydrogen bonds. Carboxylate groups can also induce intramolecular interactions, such as the formation of hydrogen bonds with donor groups, leading to an impact on the conformational space of biomolecules. In the latter case, the protonation state of the amino groups plays an important role. In order to provide an accurate description of the modifications induced in a carboxylated molecule by the formation of hydrogen bonds, one needs a compromise between a quantum chemical description of the system and the necessity to take into account explicit solvent molecules. In this work, we propose a bottom-up approach to study the conformational space and the carboxylate stretching band of (bio)organic anions. Starting from the anions in a continuum solvent, we then move to calculations using a microsolvation approach including one explicit water molecule per polar group, immersed in a continuum. Finally, we run QM/MM molecular dynamics simulations to analyze the solvation properties and to explore the anions conformational space. The results thus obtained are in good agreement with the description given by the microsolvation approach and they bring a more detailed description of the solvation shell and of the intermolecular hydrogen bonds.

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.

Optical properties of photodynamic therapy drugs in different environments: the paradigmatic case of temoporfin

Computational tools have been used to study the photophysical and photochemical features of photosensitizers in photodynamic therapy (PDT) - a minimally invasive, less aggressive alternative for cancer treatment. PDT is mainly based on the activation of molecular oxygen through the action of a photoexcited sensitizer (photosensitizer). Temoporfin, widely known as mTHPC, is a second-generation photosensitizer, which produces the cytotoxic singlet oxygen when irradiated with visible light and hence destroys tumor cells. However, the bioavailability of the mostly hydrophobic photosensitizer, and hence its incorporation into cells, is fundamental to achieve the desired effect on malignant tissues via PDT. In this study, we focus on the optical properties of the temoporfin chromophore in different environments -in vacuo, in solution, encapsulated in drug delivery agents, namely cyclodextrin, and interacting with a lipid bilayer.

Short- and Long-Range Solvation Effects on the Transient UV-Vis Absorption Spectra of a Ru(II)-Polypyridine Complex Disentangled by Nonequilibrium Molecular Dynamics

Evidence of subtle effects in the dynamic reorganization of a protic solvent in its first- and farther-neighbor shells, in response to the sudden change in the solute's electronic distribution upon excitation, is unveiled by a multilevel computational approach. Through the combination of nonequilibrium molecular dynamics and quantum mechanical calculations, the experimental time evolution of the transient T₁ absorption spectra of a heteroleptic Ru(II)-polypyridine complex in ethanol or dimethyl sulfoxide solution is reproduced and rationalized in terms of both fast and slow solvent re-equilibration processes, which are found responsible for the red shift and broadening experimentally observed only in the protic medium. Solvent orientational correlation functions and a time-dependent analysis of the solvation structure confirm that the initial, fast observed red shift can be traced back to the destruction-formation of hydrogen bond networks in the first-neighbor shell, whereas the subsequent shift, evident in the [20-500] ps range and accompanied by a large broadening of the signal, is connected to a collective reorientation of the second and farther solvation shells, which significantly changes the electrostatic embedding felt by the excited solute.

Electronic Interactions in Iminophosphorane Superbase Complexes with Carbon Dioxide

Iminophosphoranes or phosphazenes are an important class of compounds with increasing use in synthetic organic chemistry as neutral organic superbases exhibiting low nucleophilicity. Their electronic structure and therefore their properties strongly depend on substitution, but there have been very few theoretical studies devoted to this topic, and more specifically to the formation of electron donor-acceptor complexes of iminophosphoranes with electrophiles. In this work, we have investigated the interaction with carbon dioxide at different ab initio levels. Carbon dioxide usually behaves as a Lewis acid and the reaction with iminiphosphoranes has been described as a nonconventional aza-Wittig process leading to isocyanates. The reaction can be conducted in supercritical CO₂ conditions (carbon dioxide acts as both solvent and reactant), which is a promising strategy in the context of green chemistry. Our calculations have been carried out at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level for model systems and at the M06-2X/6-611+G(d,p) level for a larger species used in experiments. The electronic interactions and the interaction energies are analyzed and discussed in detail using the natural bond orbital method. Proton affinities and gas-phase basicities are provided as well.

Modeling Solvation in Supercritical CO2

In recent decades, a microscopic understanding of solute-solvent intermolecular interactions has been key to advances in technologies based on supercritical carbon dioxide. In many cases, computational work has provided the impetus for new discoveries, shedding new light on important concepts such as the local structure around the solute in the supercritical medium, the influence of the peculiar properties of the latter on the molecular behavior of dissolved substances and, importantly, CO₂ -philicity. In this Review, the theoretical work that has been relevant to these developments is surveyed and, by presenting some crucial open questions, the possible routes to achieving further progress based on the interplay between theory and experiments is discussed.

Geometry-dependent distributed polarizability models for the water molecule

Geometry-dependent distributed polarizability models have been constructed by fits to ab initio calculations at the coupled cluster level of theory with up to noniterative triple excitations in an augmented triple-zeta quality basis set for the water molecule in the field of a point charge. The investigated models include (i) charge-flow polarizabilities between chemically bonded atoms, (ii) isotropic or anisotropic dipolar polarizabilities on oxygen atom or on all atoms, and (iii) combinations of models (i) and (ii). For each model, the polarizability parameters have been optimized to reproduce the induction energy of a water molecule polarized by a point charge successively occupying a grid of points surrounding the molecule. The quality of the models is ascertained by examining their ability to reproduce these induction energies as well as the molecular dipolar and quadrupolar polarizabilities. The geometry dependence of the distributed polarizability models has been explored by changing bond lengths and HOH angle to generate 125 molecular structures (reduced to 75 symmetry-unique ones). For each considered model, the distributed polarizability components have been fitted as a function of the geometry by a Taylor expansion in monomer coordinate displacements up to the sum of powers equal to 4.

Driving Forces Controlling Host-Guest Recognition in Supercritical Carbon Dioxide Solvent

The formation of supramolecular host-guest complexes is a very useful and widely employed tool in chemistry. However, supramolecular chemistry in non-conventional solvents such as supercritical carbon dioxide (scCO2 ), one of the most promising sustainable solvents, is still in its infancy. In this work, we explored a successful route to the development of green processes in supercritical CO2 by combining a theoretical approach with experiments. We were able to synthesize and characterize an inclusion complex between a polar aromatic molecule (benzoic acid) and peracetylated-β-cyclodextrin, which is soluble in the supercritical medium. This finding opens the way to wide, environmental friendly, applications of scCO2 in many areas of chemistry, including supramolecular synthesis, reactivity and catalysis, micro and nano-particle formation, molecular recognition, as well as enhanced extraction processes with increased selectivity.

Water interactions with hydrophobic groups: assessment and recalibration of semiempirical molecular orbital methods

In this work, we present a study of the ability of different semiempirical methods to describe intermolecular interactions in water solution. In particular, we focus on methods based on the Neglect of Diatomic Differential Overlap approximation. Significant improvements of these methods have been reported in the literature in the past years regarding the description of non-covalent interactions. In particular, a broad range of methodologies has been developed to deal with the properties of hydrogen-bonded systems, with varying degrees of success. In contrast, the interactions between water and a molecule containing hydrophobic groups have been little analyzed. Indeed, by considering the potential energy surfaces obtained using different semiempirical Hamiltonians for the intermolecular interactions of model systems, we found that none of the available methods provides an entirely satisfactory description of both hydrophobic and hydrophilic interactions in water. In addition, a vibrational analysis carried out in a model system for these interactions, a methane clathrate cluster, showed that some recent methods cannot be used to carry out studies of vibrational properties. Following a procedure established in our group [M. I. Bernal-Uruchurtu, M. T. C. Martins-Costa, C. Millot, and M. F. Ruiz-López, J. Comput. Chem. 21, 572 (2000); W. Harb, M. I. Bernal-Uruchurtu, and M. F. Ruiz-López, Theor. Chem. Acc. 112, 204 (2004)], we developed new parameters for the core-core interaction terms based on fitting potential energy curves obtained at the MP2 level for our model system. We investigated the transferability of the new parameters to describe a system, having both hydrophilic and hydrophobic groups, interacting with water. We found that only by introducing two different sets of parameters for hydrophilic and hydrophobic hydrogen atom types we are able to match the features of the ab initio calculated properties. Once this assumption is made, a good agreement with the MP2 reference is achieved. The results reported in this work provide therefore a direction for future developments of semiempirical approaches that are still required to investigate chemical processes in biomolecules and in large disordered systems.

Flexible AlN flags for efficient wind energy harvesting at ultralow cut-in wind speed

Flexible AlN for better harvesting of wind power at low wind speed.

Hydration Effect on Amide I Infrared Bands in Water: An Interpretation Based on an Interaction Energy Decomposition Scheme

The sensitivity of some infrared bands to the local environment can be exploited to shed light on the structure and the dynamics of biological systems. In particular, the amide I band, which is specifically related to vibrations within the peptide bonds, can give information on the ternary structure of proteins, and can be used as a probe of energy transfer. In this work, we propose a model to quantitatively interpret the frequency shift on the amide I band of a model peptide induced by the formation of hydrogen bonds in the first solvation shell. This method allows us to analyze to what extent the electrostatic interaction, electronic polarization and charge transfer affect the position of the amide I band. The impact of the anharmoniticy of the pontential energy surface on the hydration induced shift is elucidated as well.

Cavity closure dynamics of peracetylated β-cyclodextrins in supercritical carbon dioxide

Structural properties of peracetylated β-cyclodextrin in supercritical carbon dioxide were investigated by means of molecular dynamics simulations. The study indicated a strong reduction of the cavity accessibility to guest molecules, compared to native β-cyclodextrin in water. Indeed, the cavity is self-closed during the largest part of the simulation, which agrees well with suggestions made on the basis on high-pressure NMR experiments. Self-closure happens because one glucose unit undergoes a main conformational change (from chair to skew) that brings one of the acetyl groups in the wide rim of the cyclodextrin to the cavity interior. This arrangement turns out to be quite favorable, persisting for several nanoseconds. In addition to the wide rim self-closure, a narrow rim self-closure may also occur, though it is less likely and exhibits short duration (<1 ns). Therefore, the number of solvent molecules reaching the cavity interior is much smaller than that found in the case of native β-cyclodextrin in water after correction to account for different molar densities. These findings support the weak tendency of the macromolecule to form host-guest complexes in this nonconventional medium, as reported by some experiments. Finally, Lewis acid/base interactions between the acetyl carbonyl groups and the solvent CO(2) molecules were analyzed through ab initio calculations that revealed the existence of a quite favorable four-member ring structure not yet reported. The ensemble of these results can contribute to establish general thermodynamic principles controlling the formation of inclusion complexes in supercritical CO(2), where the hydrophilicity/hydrophobicity balance is not applicable.

Pathways for H2O bend vibrational relaxation in liquid water

The mechanism of the H2O bend vibrational relaxation in liquid water has been examined via classical MD simulations and an analysis of work and power contributions. The relaxation is found to be dominated by energy flow to the hindered rotation of the bend excited water molecule. This energy transfer, representing approximately 2/3 of the transferred energy, is due to a 2:1 Fermi resonance for the centrifugal coupling between the water bend and rotation. The remaining energy flow (approximately 1/3) from the excited water bend is dominated by transfer to the excited water molecule's first four water neighbors, i.e., the first hydration shell, and is itself dominated by energy flow to the two water molecules hydrogen (H)-bonded to the hydrogens of the central H2O. The energy flow from the produced rotationally excited central molecule is less local in character, with approximately half of its rotational kinetic energy being transferred to water molecules outside of the first hydration shell, whereas the remaining half is preferentially transferred to the two first hydration shell water molecules donating H-bonds to the central water oxygen. The overall energy flow is well described by an approximate kinetic scheme.

Electron Transfer in a Radical Ion Pair: Quantum Calculations of the Solvent Reorganization Energy

Results are presented for an investigation of intermolecular electron transfer (ET) in solution by means of quantum calculations. The two molecules that are involved in the ET reaction form a solvent-separated radical ion pair. The solvent plays an important role in the ET between the two molecules. In particular, it can give rise to specific solute-solvent interactions with the solutes. An example of specific interactions is the formation of a hydrogen bond between a protic solvent and one of the molecules involved in the ET. We address the study of this system by means of quantum calculations on the solutes immersed in a continuum solvent. However, when the solvent can give rise to hydrogen bond formation with the negatively charged ion after ET, we explicitly consider solvent molecules in the solute cavity, determining the hydrogen bond energetic contribution to the overall interaction energy. Solute-solvent pair distribution functions, showing the different arrangement of solvent molecules before and after ET in the first solvation shell, are reported. We provide results of the solvent reorganization energy from quantum calculations for both the two isolated fragments and the ion pair in solution. Results are in agreement with available experimental data.

Solvation Dynamics by Computer Simulation: Coumarin C153 in 1,4-Dioxane

Computer simulation results are presented for an atomistic pair potential model of 1,4-dioxane which takes into account molecular flexibility. The model has been conceived to be applied to the study of solvatochromism and solvation dynamics in the presence of the polar probe coumarin C153. Computer simulations on the pure liquid have produced thermodynamical, structural, and dynamical data in good agreement with available experimental measures. This constitutes a valuable test of the 1,4-dioxane all-atom model employed. The study of solute-solvent interactions for C153 in 1,4-dioxane has been motivated by the aim of casting light, through simulations, on the interesting experimental findings according to which such a solvent behaves as a "polar" solvent with respect to dynamic solvation properties. Molecular dynamics is particularly suitable to model the process and provides an interpretation of the so-called "dioxane anomaly". An investigation of the structure of the solvation shell and of the dynamics of solvation is presented and discussed. In particular, the satisfactory accordance between simulated and experimental solvation response implies that the simulations give a reliable description of both solute and solvent at a molecular level and reinforces the idea that the explicit inclusion of discrete solvent molecules is needed for a realistic treatment of solvation phenomena in which the local structure of the liquid plays a key role.

Solvation Dynamics of C153 in Supercritical Fluoroform: A Simulation Study Based on Two-Site and Five-Site Models of the Solvent

Molecular dynamics (MD) simulations of a probe solute (coumarin C153) in supercritical fluoroform are used to study time-dependent solute-solvent interactions. We study the dynamics of solvent reorganization in response to electronic excitation of C153 at a temperature of 1.03 T(c) (the critical temperature) and a series of densities above and below the critical density. Simulations of a two-site and five-site models of fluoroform are presented and compared. The time-dependent solvation response after solute electronic excitation is studied in the two cases, and the five-site results present an earlier onset of exponential decay that is closer to what is expected to be the experimental response. This is confirmed by comparison to experiment. In addition to obtaining the solvation response from nonequilibrium MD trajectories, approximate solvation responses were obtained from equilibrium time correlations of the fluctuations in the solvation energy change in the presence of ground- and excited-state solutes. For the five-site model, the equilibrium excited-state response shows stronger density dependence than the ground-state one. The nonequilibrium response appears to have an intermediate decay rate between the two equilibrium functions. The solute-partial-charge-solvent-induced-dipole interaction was also taken into account by means of a perturbative approach, which improved the agreement with experimental measurements available at densities corresponding to 1.4-1.6 rho(c) (where rho(c) the critical density). From the comparison between the two models, it is possible to conclude that an atomistic description is necessary for correctly representing the portion of solvation dynamics that is related to reorientation. This consideration is supported by providing results for orientational time correlation functions and by comparing the correlation times with the experimental ones.

Formation and relaxation of excited states in solution: A new time dependent polarizable continuum model based on time dependent density functional theory

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, epsilonomega. 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.

Solvation of Coumarin 153 in Supercritical Fluoroform

We present a study of local density augmentation around an attractive solute (i.e., giving rise to more attractive interaction with the solvent than solvent-solvent interactions) in supercritical fluoroform. This work is based on molecular dynamics simulations of coumarin 153 in supercritical fluoroform at densities both above and below the critical density, ranging from dilute gas-like to liquid-like, at a reduced temperature (T/T(c)) of 1.03. We focused on studying the structure of the solvation shell and the variation of the solute electronic absorption and emission shifts with density. Quantum calculations at the density functional theory (DFT) level were run on the solute in the ground state, and time-dependent DFT calculations were performed in the solute excited state in order to determine the solute-solvent potential parameters. The results obtained for the Stokes shift are in agreement with the experimental measurements. To evaluate local density augmentation from simulations, we used two different definitions, one based on the solvation number and the other derived from solvatochromic shifts. In the former case, the agreement with experimental results is good, while, in the latter case, better agreement is achieved by perturbatively including the induced-dipole contribution to the solvation energy.

A time-dependent polarizable continuum model: Theory and application

This work presents an extention of the polarizable continuum model to explicitly describe the time-dependent response of the solvent to a change in the solute charge distribution. Starting from an initial situation in which solute and solvent are in equilibrium, we are interested in modeling the time-dependent evolution of the solvent response, and consequently of the solute-solvent interaction, after a perturbation in this equilibrium situation has been switched on. The model introduces an explicit time-dependent treatment of the polarization by means of the linear-response theory. Two strategies are tested to account for this time dependence: the first one employs the Debye model for the dielectric relaxation, which assumes an exponential decay of the solvent polarization; the second one is based on a fitting of the experimental data of the solvent complex dielectric permittivity. The first approach is simpler and possibly less accurate but allows one to write an analytic expression of the equations. By contrast, the second approach is closer to the experimental evidence but it is limited to the availability of experimental data. The model is applied to the ionization process of N,N-dimethyl-aniline in both acetonitrile and water. The nonequilibrium free-energy profile is studied both as a function of the solvent relaxation coordinate and as a function of time. The solvent reorganization energy is evaluated as well.

Solvation Dynamics in Acetonitrile: A Study Incorporating Solute Electronic Response and Nuclear Relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.