On Ultrafast Time-Domain TeraHertz Spectroscopy in the Condensed Phase: Linear Spectroscopic Measurements of Hydrogen-Bond Dynamics of Astrochemical Ice Analogs and Nonlinear TeraHertz Kerr Effect Measurements of Vibrational Quantum Beats

Much of the chemistry that affects life on planet Earth occurs in the condensed phase. The TeraHertz (THz) or far-infrared (far-IR) region of the electromagnetic spectrum (from 0.1 THz to 10 THz, 3 cm^-1 to 300 cm^-1, or 3000 μm to 30 μm) has been shown to provide unique possibilities in the study of condensed-phase processes. The goal of this work is to expand the possibilities available in the THz region and undertake new investigations of fundamental interest to chemistry. Since we are fundamentally interested in condensed-phase processes, this thesis focuses on two areas where THz spectroscopy can provide new understanding: astrochemistry and solvation science. To advance these fields, we had to develop new instrumentation that would enable the experiments necessary to answer new questions in either astrochemistry or solvation science. We first developed a new experimental setup capable of studying astrochemical ice analogs in both the TeraHertz (THz), or far-Infrared (far-IR), region (0.3 - 7.5 THz; 10 - 250 cm^-1) and the mid-IR (400 - 4000 cm^-1). The importance of astrochemical ices lies in their key role in the formation of complex organic molecules, such as amino acids and sugars in space. Thus, the instruments are capable of performing variety of spectroscopic studies that can provide especially relevant laboratory data to support astronomical observations from telescopes such as the Herschel Space Telescope, the Stratospheric Observatory for Infrared Astronomy (SOFIA), and the Atacama Large Millimeter Array (ALMA). The experimental apparatus uses a THz time-domain spectrometer, with a 1750/875 nm plasma source and a GaP detector crystal, to cover the bandwidth mentioned above with ~10 GHz (~0.3 cm^-1) resolution.

Using the above instrumentation, experimental spectra of astrochemical ice analogs of water and carbon dioxide in pure, mixed, and layered ices were collected at different temperatures under high vacuum conditions with the goal of investigating the structure of the ice. We tentatively observe a new feature in both amorphous solid water and crystalline water at 33 cm^-1 (1 THz). In addition, our studies of mixed and layered ices show how it is possible to identify the location of carbon dioxide as it segregates within the ice by observing its effect on the THz spectrum of water ice. The THz spectra of mixed and layered ices are further analyzed by fitting their spectra features to those of pure amorphous solid water and crystalline water ice to quantify the effects of temperature changes on structure. From the results of this work, it appears that THz spectroscopy is potentially well suited to study thermal transformations within the ice.

To advance the study of liquids with THz spectroscopy, we developed a new ultrafast nonlinear THz spectroscopic technique: heterodyne-detected, ultrafast THz Kerr effect (TKE) spectroscopy. We implemented a heterodyne-detection scheme into a TKE spectrometer that uses a stilbazoiumbased THz emitter, 4-N,N-dimethylamino-4-N-methyl-stilbazolium 2,4,6-trimethylbenzenesulfonate (DSTMS), and high numerical aperture optics which generates THz electric field in excess of 300 kV/cm, in the sample. This allows us to report the first measurement of quantum beats at terahertz (THz) frequencies that result from vibrational coherences initiated by the nonlinear, dipolar interaction of a broadband, high-energy, (sub)picosecond THz pulse with the sample. Our instrument improves on both the frequency coverage, and sensitivity previously reported; it also ensures a backgroundless measurement of the THz Kerr effect in pure liquids. For liquid diiodomethane, we observe a quantum beat at 3.66 THz (122 cm^-1), in exact agreement with the fundamental transition frequency of the υ4 vibration of the molecule. This result provides new insight into dipolar vs. Raman selection rules at terahertz frequencies.

To conclude we discuss future directions for the nonlinear THz spectroscopy in the Blake lab. We report the first results from an experiment using a plasma-based THz source for nonlinear spectroscopy that has the potential to enable nonlinear THz spectra with a sub-100 fs temporal resolution, and how the optics involved in the plasma mechanism can enable THz pulse shaping. Finally, we discuss how a single-shot THz detection scheme could improve the acquisition of THz data and how such a scheme could be implemented in the Blake lab. The instruments developed herein will hopefully remain a part of the groups core competencies and serve as building blocks for the next generation of THz instrumentation that pushes the frontiers of both chemistry and the scientific enterprise as a whole.

On the solvation of L-aspartic acid

We use molecular statics and dynamics to study the stability of L-aspartic acid both in vacuo and solvated by polar and non-polar molecules using density functional theory in the generalized gradient approximation. We find that structures stable in vacuo are unstable in aqueous solution and vice versa. From our simulations we are able to come to some conclusions about the mechanism of stabilisation of zwitterions by polar protic solvents, water and methanol.

Neutron diffraction, NMR and molecular dynamics study of glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate

beta-D-glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate in a 6 : 1 molar ratio (ionic liquid : glucose) has been studied by neutron scattering, NMR and molecular dynamics simulations. Good agreement was found between simulated neutron scattering profiles generated for isotopically substituted liquid systems and those experimentally determined as well as between simulated and experimental diffusion coefficients obtained by Pulsed Field Gradient NMR spectroscopy. The overriding glucose-ionic liquid interactions in the liquid are hydrogen-bonding between acetate oxygens and sugar hydroxyl groups. The ionic liquid cation was found to play only a minor role in the solvation of the sugar and does not participate in hydrogen-bonding with the sugar to any significant degree. NOESY experiments lend further evidence that there is no direct interaction between sugar hydroxyl groups and acidic hydrogens on the ionic liquid cation.

A Neutron Diffraction and Molecular Dynamics Investigation of Acetate-based Ionic Liquids as Solvents for Glucose

The liquid state structure of the ionic liquid, 1-ethyl-3-methylimidazolium acetate, and the solute/solvent structure of glucose dissolved in the ionic liquid at a 1: 6 molar ratio have been investigated at 323 K by molecular dynamics simulations and neutron diffraction experiments using H/D isotopically substituted materials. Interactions between hydrogen-bond donating cation sites and polar, directional hydrogen-bond accepting acetate anions are examined. Ion-ion radial distribution functions for the neat ionic liquid, calculated from both MD and derived from the empirical potential structure refinement model to the experimental data, show the alternating shell-structure of anions around the cation, as anticipated. Spatial probability distributions reveal the main anion-to-cation features as in-plane interactions of anions with imidazolium ring hydrogens and cation-cation planar stacking. Interestingly, the presence of the polarised hydrogen-bond acceptor anion leads to increased anion-anion tail-tail structuring within each anion shell, indicating the onset of hydrophobic regions within the anion regions of the liquid.

Elementary steps in aqueous proton transfer reactions : a first principles molecular dynamics study

La nature des acides dans un environnement aqueux est primordiale dans de nombreux aspects de la chimie et de la biologie. La caractéristique principale d'un acide est sa capacité à transférer un proton vers une molécule d'eau ou vers n'importe quelle base, mais ce procédé n'est pas aussi simple qu'il y paraît. Il peut au contraire être extrêmement complexe et dépendre de manière cruciale de la solvatation des différents intermédiaires de réaction impliqués. Cette thèse décrit les études computationnelles basées sur des simulations de dynamique moléculaire ab initio qui ont pour but d'obtenir une description à l'échelle moléculaire des divers procédés de transferts de proton entre acide et bases dans un milieu aqueux. Pour cela, nous avons étudié une serie de système, dont l'acide hydrofluorique aqueux, l'acide trifluoroacétique aqueux, et un système modèle constitué d'un phénol et d'une entité carboxylate reliés entre eux par une molécule d'eau en solution aqueuse. Deux états intermédiaires ont été identifiés pour le transfert d'un proton depuis un acide. Ces intermédiaires apparaissent stabilisés par un motif local de solvatation via des ponts H. Leurs signatures spectroscopiques ont été caractérisées au moyen de la spectroscopie infrarouge, en utilisant le formalisme de la dynamique moléculaire ab initio, qui inclut l'effet quantique nucléaire de manière explicite. Cette étude a aussi identifié trois chemins de réaction élémentaire, qui sont responsable pour le transfert d'un proton d'un acide à une base, ainsi que leurs échelles de temps caractéristiques. Les conclusions tirées de ces études sont discutées dans les détails, au niveau moléculaire, avec une emphase sur les comparaisons entre les résultats théoriques et les mesures expérimentales obtenues dans a littérature ou via des collaborateurs.

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Function of the hydration layer around an antifreeze protein revealed by atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations are used to investigate the mechanism by which the antifreeze protein from the spruce budworm, Choristoneura fumiferana, binds to ice. Comparison of structural and dynamic properties of the water around the three faces of the triangular prism-shaped protein in aqueous solution reveals that at low temperature the water structure is ordered and the dynamics slowed down around the ice-binding face of the protein, with a disordering effect observed around the other two faces. These results suggest a dual role for the solvation water around the protein. The preconfigured solvation shell around the ice-binding face is involved in the initial recognition and binding of the antifreeze protein to ice by lowering the barrier for binding and consolidation of the protein:ice interaction surface. Thus, the antifreeze protein can bind to the molecularly rough ice surface by becoming actively involved in the formation of its own binding site. Also, the disruption of water structure around the rest of the protein helps prevent the adsorbed protein becoming covered by further ice growth.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.

Picosecond Dynamics of Proton Transfer of a 7-Hydroxyflavylium Salt in Aqueous-Organic Solvent Mixtures

The intermediacy of the geminate base proton pair (A*center dot center dot center dot H(+)) in excited-state proton-transfer (ESPT) reactions (two-step mechanism) has been investigated employing the synthetic flavylium salt 7-hydroxy-4-methyl-flavylium chloride (HMF). In aqueous solution, the ESPT mechanism involves solely the excited acid AH* and base A* forms of HMF as indicated by the fluorescence spectra and double-exponential fluorescence decays (two species, two decay times). However, upon addition of either 1,4-dioxane or 1,2-propylene glycol, the decays become triple-exponential with a term consistent with the presence of the geminate base proton pair A*center dot center dot center dot H(+). The geminate pair becomes detectable because of the increase in the recombination rate constant, k(rec), of (A*center dot center dot center dot H(+)) with increasing the mole fraction of added organic cosolvent. Because the two-step ESPT mechanism splits the intrinsic prototropic reaction rates (deprotonation of AH(+)*, k(d), and recombination, k(rec) of A*center dot center dot center dot H(+)) from the diffusion controlled rates (dissociation, k(diss) and formation, k(diff)[H(+)], of A*center dot center dot center dot H+), the experimental detection of the geminate pair provides a wealth of information on the proton-transfer reaction (k(d) and k(rec)) as well as on proton diffusion/migration (k(diss) and k(diff)).

Correlation between the mobility of spin-labeled peptide chains and resin solvation: An approach to optimize the synthesis of aggregating sequences

Resin solvation properties affect the efficiency of the coupling reactions in solid-phase peptide synthesis. Here we report a novel approach to evaluate resin solvation properties, making use of spin label electron paramagnetic resonance (EPR) spectroscopy. The aggregating VVLGAAIV and ING sequences were assembled in benzhydrylamine-resin with different amino group contents (up to 2.6 mmol/g) to examine the extent of chain association within the beads. These model peptidyl-resins were first labeled at their N-terminus with the amino acid spin label 2,2,6,6-tetramethylpiperidine-N-oxyl-4-amino-4-carboxylic acid (Toac). Their solvation properties in different solvents were estimated, either by bead swelling measurement or by assessing the dynamics of their polymeric matrixes through the analysis of Toac EPR spectra, and were correlated with the yield of the acylation reaction. In most cases the coupling rate was found to depend on bead swelling. Comparatively, the EPR approach was more effective. Line shape analysis allowed the detection of more than one peptide chain population, which influenced the reaction. The results demonstrated the unique potential of EPR spectroscopy not only for improving the yield of peptide synthesis, even in challenging conditions, but also for other relevant polymer-supported methodologies in chemistry and biology.

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Chiroptical properties of bioactive molecules: sensitivity to conformation and solvation

Chiroptical spectroscopies play a fundamental role in pharmaceutical analysis for the stereochemical characterisation of bioactive molecules, due to the close relationship between chirality and optical activity and the increasing evidence of stereoselectivity in the pharmacological and toxicological profiles of chiral drugs. The correlation between chiroptical properties and absolute stereochemistry, however, requires the development of accurate and reliable theoretical models. The present thesis will report the application of theoretical chiroptical spectroscopies in the field of drug analysis, with particular emphasis on the huge influence of conformational flexibility and solvation on chiroptical properties and on the main computational strategies available to describe their effects by means of electronic circular dichroism (ECD) spectroscopy and time-dependent density functional theory (TD-DFT) calculations. The combination of experimental chiroptical spectroscopies with state-of-the-art computational methods proved to be very efficient at predicting the absolute configuration of a wide range of bioactive molecules (fluorinated 2-arylpropionic acids, β-lactam derivatives, difenoconazole, fenoterol, mycoleptones, austdiol). The results obtained for the investigated systems showed that great care must be taken in describing the molecular system in the most accurate fashion, since chiroptical properties are very sensitive to small electronic and conformational perturbations. In the future, the improvement of theoretical models and methods, such as ab initio molecular dynamics, will benefit pharmaceutical analysis in the investigation of non-trivial effects on the chiroptical properties of solvated systems and in the characterisation of the stereochemistry of complex chiral drugs.

Solvation-Driven Charge Transfer and Localization in Metal Complexes

In any physicochemical process in liquids, the dynamical response of the solvent to the solutes out of equilibrium plays a crucial role in the rates and products: the solvent molecules react to the changes in volume and electron density of the solutes to minimize the free energy of the solution, thus modulating the activation barriers and stabilizing (or destabilizing) intermediate states. In charge transfer (CT) processes in polar solvents, the response of the solvent always assists the formation of charge separation states by stabilizing the energy of the localized charges. A deep understanding of the solvation mechanisms and time scales is therefore essential for a correct description of any photochemical process in dense phase and for designing molecular devices based on photosensitizers with CT excited states. In the last two decades, with the advent of ultrafast time-resolved spectroscopies, microscopic models describing the relevant case of polar solvation (where both the solvent and the solute molecules have a permanent electric dipole and the mutual interaction is mainly dipole−dipole) have dramatically progressed. Regardless of the details of each model, they all assume that the effect of the electrostatic fields of the solvent molecules on the internal electronic dynamics of the solute are perturbative and that the solvent−solute coupling is mainly an electrostatic interaction between the constant permanent dipoles of the solute and the solvent molecules. This well-established picture has proven to quantitatively rationalize spectroscopic effects of environmental and electric dynamics (time-resolved Stokes shifts, inhomogeneous broadening, etc.). However, recent computational and experimental studies, including ours, have shown that further improvement is required. Indeed, in the last years we investigated several molecular complexes exhibiting photoexcited CT states, and we found that the current description of the formation and stabilization of CT states in an important group of molecules such as transition metal complexes is inaccurate. In particular, we proved that the solvent molecules are not just spectators of intramolecular electron density redistribution but significantly modulate it. Our results solicit further development of quantum mechanics computational methods to treat the solute and (at least) the closest solvent molecules including the nonperturbative treatment of the effects of local electrostatics and direct solvent−solute interactions to describe the dynamical changes of the solute excited states during the solvent response.