Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1 clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

Structural aspects of redox-mediated electron tunneling (Dedicated to the memory of professor Alexander M. Kuznetsov)

We investigated structural aspects of electron transfer (ET) in tunneling junctions (Au(1 1 1)vertical bar FcN vertical bar solution gap vertical bar Au STM tip) with four different redox-active N-thioalk(ano)ylferrocenes (FcN) embedded. The investigated molecules consist of a redox-active ferrocene (Fc) moiety connected via alkyl spacers with N = 4, 6, 8 and 11 carbon atoms to a thiol anchoring group. We found that for short FcNs (N = 4, 6,8) the redox-mediated ET response increases with the increase of the alkyl chain length, while no enhancement of the ET was observed for Fc1 1. The model of two-step ET with partial vibrational relaxation by Kuznetsov and Ulstrup was used to rationalize these results. The theoretical ET steps were assigned to two processes: (1) electron tunneling from the Fc group to the Au tip through the electrolyte layer and (2) electron transport from the Au(1 1 1) substrate to the Fc group through the organic adlayer. We argue that for the three short FcNs, the first process represents the rate-limiting step. The increase of the length of the alkyl chain leads to an approach of the Fc group to the STM tip, and consequently accelerates the first El' step. In case of the Fcl 1 junctions the rather high thickness of the organic layer leads to a decrease of the rate of the second ET step. In consequence, the contribution of the redox-mediated current enhancement to the total tunneling current appears to be insignificant. Our work demonstrates the importance of combined structural and transport approaches for the understanding of Er processes in electrochemical nanosystems. (C) 2010 Elsevier B.V. All rights reserved.

Femtosecond Rotational Raman Coherence Spectroscopy of Cyclohexane in a Pulsed Supersonic Jet

We combine the technique of femtosecond degenerate four-wave mixing (fs-DFWM) with a high repetition-rate pulsed supersonic jet source to obtain the rotational coherence spectrum (RCS) of cold cyclohexane (C(6)H(12)) with high signal/noise ratio. In the jet expansion, the near-parallel flow pattern combined with rapid translational cooling effectively eliminate dephasing collisions, giving near-constant RCS signal intensities over time delays up to 5 ns. The vibrational cooling in the jet eliminates the thermally populated vibrations that complicate the RCS coherences of cyclohexane at room temperature [Bragger, G.; et al. J. Phys. Chem. A 2011, 115, 9567]. The rotational cooling reduces the high-J rotational-state population, yielding the most accurate ground-state rotational constant to date, B(0) = 4305.859(9) MHz. Based on this B(0), a reanalysis of previous room-temperature gas-cell RCS measurements of cydohexane gives improved vibration rotation interaction constants for the v(32), v(6), v(16), and v(24) vibrational states. Combining the experimental B(0)(C(6)H(12)) with CCSD(T) calculations yields a very accurate semiexperimental equilibrium structure of the chair isomer of cyclohexane

Product Distributions Following Hyperthermal Collisions Between C2H4 And O(3P)

Transient Diode Laser Absorption Spectroscopy (TDLAS) was used to perform vibrational state population studies of the CO2 product from the hyperthermal reaction between C2H4 and O(3P) at room temperature using O3 as the O-atom precursor. Photodissociation of O3 using a frequency quadrupled Q-switch Nd:YAG laser pulse at 266 nm produced O(3P) atoms at high velocities which subsequently reacted with C2H4, producing several primary and secondary products including CO2. The CO2 product was detected using high-resolution TDLAS under five unique sets of reaction conditions. The vibrational distribution of the CO2 product did not follow a Boltzmann distribution at all five sets of conditions. The experiments showed a distribution in which there was a surprisingly high population in the (1000) (symmetric stretching) state compared with the other states probed, all of which contained bend excitation. In general, the CO2 population in the (1000) state was about 15-20% more populated than the Boltzmann distribution predicts. A possible explanation for this result may lie in the mechanism of CO2 evolution from the C2H4 + O(3P) reaction.

A Computational Analysis of the Collision-Induced Dissociation Mechanisms of the Dipeptide Glycine-Serine Under Mass Spectrometry

Collision-induced dissociation (CID) of peptides using tandem mass spectrometry (MS) has been used to determine the identity of peptides and other large biological molecules. Mass spectrometry (MS) is a useful tool for determining the identity of molecules based on their interaction with electromagnetic fields. If coupled with another method like infrared (IR) vibrational spectroscopy, MS can provide structural information, but in its own right, MS can only provide the mass-to-charge (m/z) ratio of the fragments produced, which may not be enough information to determine the mechanism of the collision-induced dissociation (CID) of the molecule. In this case, theoretical calculations provide a useful companion for MS data and yield clues about the energetics of the dissociation. In this study, negative ion electrospray tandem MS was used to study the CID of the deprotonated dipeptide glycine-serine (Gly-Ser). Though negative ion MS is not as popular a choice as positive ion MS, studies by Bowie et al. show that it yields unique clues about molecular structure which complement positive ion spectroscopy, such as characteristic fragmentations like the loss of formaldehyde from the serine residue.2 The increase in the collision energy in the mass spectrometer alters the flexibility of the dipeptide backbone, enabling isomerizations (reactions not resulting in a fragment loss) and dissociations to take place. The mechanism of the CID of Gly-Ser was studied using two computational methods, B3LYP/6-311+G* and M06-2X/6-311++G**. The main pathway for molecular dissociation was analyzed in 5 conformers in an attempt to verify the initial mechanism proposed by Dr. James Swan after examination of the MS data. The results suggest that the loss of formaldehyde from serine, which Bowie et al. indicates is a characteristic of the presence of serine in a protein residue, is an endothermic reaction that is made possible by the conversion of the translational energy of the ion into internal energy as the ion collides with the inert collision gas. It has also been determined that the M06-2X functional¿s improved description of medium and long-range correlation makes it more effective than the B3LYP functional at finding elusive transition states. M06-2X also more accurately predicts the energy of those transition states than does B3LYP. A second CID mechanism, which passes through intermediates with the same m/z ratio as the main pathway for molecular dissociation, but different structures, including a diketopiperazine intermediate, was also studied. This pathway for molecular dissociation was analyzed with 3 conformers and the M06-2X functional, due to its previously determined effectiveness. The results suggest that the latter pathway, which meets the same intermediate masses as the first mechanism, is lower in overall energy and therefore a more likely pathway of dissociation than the first mechanism.

Photodissociation dynamics of tert-butylnitrite following excitation to the S(1) and S(2) states. A study by velocity-map ion-imaging and 3D-REMPI spectroscopy

Excitation of tert-butylnitrite into the first and second UV absorption bands leads to efficient dissociation into the fragment radicals NO and tert-butoxy in their electronic ground states (2)Π and (2)E, respectively. Velocity distributions and angular anisotropies for the NO fragment in several hundred rotational and vibrational quantum states were obtained by velocity-map imaging and the recently developed 3D-REMPI method. Excitation into the well resolved vibronic progression bands (k = 0, 1, 2) of the NO stretch mode in the S(1) ← S(0) transition produces NO fragments mostly in the vibrational state with v = k, with smaller fractions in v = k - 1 and v = k - 2. It is concluded that dissociation occurs on the purely repulsive PES of S(1) without barrier. All velocity distributions from photolysis via the S(1)(nπ*) state are monomodal and show high negative anisotropy (β ≈ -1). The rotational distributions peak near j = 30.5 irrespective of the vibronic state S(1)(k) excited and the vibrational state v of the NO fragment. On average 46% of the excess energy is converted to kinetic energy, 23% and 31% remain as internal energy in the NO fragment and the t-BuO radical, respectively. Photolysis via excitation into the S(2) ← S(0) transition at 227 nm yields NO fragments with about equal populations in v = 0 and v = 1. The rotational distributions have a single maximum near j = 59.5. The velocity distributions are monomodal with positive anisotropy β ≈ 0.8. The average fractions of the excess energy distributed into translation, internal energy of NO, and internal energy of t-BuO are 39%, 23%, and 38%, respectively. In all cases ∼8500 cm(-1) of energy remain in the internal degrees of freedom of the t-BuO fragment. This is mostly assigned to rotational energy. An ab initio calculation of the dynamic reaction path shows that not only the NO fragment but also the t-BuO fragment gain large angular momentum during dissociation on the purely repulsive potential energy surface of S(2).

Probing chiral interfaces by infrared spectroscopic methods

Biological homochirality on earth and its tremendous consequences for pharmaceutical science and technology has led to an ever increasing interest in the selective production, the resolution and the detection of enantiomers of a chiral compound. Chiral surfaces and interfaces that can distinguish between enantiomers play a key role in this respect as enantioselective catalysts as well as for separation purposes. Despite the impressive progress in these areas in the last decade, molecular-level understanding of the interactions that are at the origin of enantiodiscrimination are lagging behind due to the lack of powerful experimental techniques to spot these interactions selectively with high sensitivity. In this article, techniques based on infrared spectroscopy are highlighted that are able to selectively target the chiral properties of interfaces. In particular, these methods are the combination of Attenuated Total Reflection InfraRed (ATR-IR) with Modulation Excitation Spectroscopy (MES) to probe enantiodiscriminating interactions at chiral solid-liquid interfaces and Vibrational Circular Dichroism (VCD), which is used to probe the structure of chirally-modified metal nanoparticles. The former technique aims at suppressing signals arising from non-selective interactions, which may completely hide the signals of interest due to enantiodiscriminating interactions. Recently, this method was successfully applied to investigate enantiodiscrimination at self-assembled monolayers of chiral thiols on gold surfaces. The nanometer size analogues of the latter--gold nanoparticles protected by a monolayer of a chiral thiol--are amenable to VCD spectroscopy. It is shown that this technique yields detailed structural information on the adsorption mode and the conformation of the adsorbed thiol. This may also turn out to be useful to clarify how chirality can be bestowed onto the metal core itself and the nature of the chirality of the latter, which is manifested in the metal-based circular dichroism activity of these nanoparticles.

Quantum transport in a single molecular junction

The craze for faster and smaller electronic devices has never gone down and this has always kept researchers on their toes. Following Moore’s law, which states that the number of transistors in a single chip will double in every 18 months, today “30 million transistors can fit into the head of a 1.5 mm diameter pin”. But this miniaturization cannot continue indefinitely due to the ‘quantum leakage’ limit in the thickness of the insulating layer between the gate electrode and the current carrying channel. To bypass this limitation, scientists came up with the idea of using vastly available organic molecules as components in an electronic device. One of the primary challenges in this field was the ability to perform conductance measurements across single molecular junctions. Once that was achieved the focus shifted to a deeper understanding of the underlying physics behind the electron transport across these molecular scale devices. Our initial theoretical approach is based on the conventional Non-Equilibrium Green Function(NEGF) formulation, but the self-energy of the leads is modified to include a weighting factor that ensures negligible current in the absence of a molecular pathway as observed in a Mechanically Controlled Break Junction (MCBJ) experiment. The formulation is then made parameter free by a more careful estimation of the self-energy of the leads. The calculated conductance turns out to be atleast an order more than the experimental values which is probably due to a strong chemical bond at the metal-molecule junction unlike in the experiments. The focus is then shifted to a comparative study of charge transport in molecular wires of different lengths within the same formalism. The molecular wires, composed of a series of organic molecules, are sanwiched between two gold electrodes to make a two terminal device. The length of the wire is increased by sequentially increasing the number of molecules in the wire from 1 to 3. In the low bias regime all the molecular devices are found to exhibit Ohmic behavior. However, the magnitude of conductance decreases exponentially with increase in length of the wire. In the next study, the relative contribution of the ‘in-phase’ and the ‘out-of-phase’ components of the total electronic current under the influence of an external bias is estimated for the wires of three different lengths. In the low bias regime, the ‘out-of-phase’ contribution to the total current is minimal and the ‘in-phase’ elastic tunneling of the electrons is responsible for the net electronic current. This is true irrespective of the length of the molecular spacer. In this regime, the current-voltage characteristics follow Ohm’s law and the conductance of the wires is found to decrease exponentially with increase in length which is in agreement with experimental results. However, after a certain ‘off-set’ voltage, the current increases non-linearly with bias and the ‘out-of-phase’ tunneling of electrons reduces the net current substantially. Subsequently, the interaction of conduction electrons with the vibrational modes as a function of external bias in the three different oligomers is studied since they are one of the main sources of phase-breaking scattering. The number of vibrational modes that couple strongly with the frontier molecular orbitals are found to increase with length of the spacer and the external field. This is consistent with the existence of lowest ‘off-set’ voltage for the longest wire under study.

Fluidic oscillator design for water removal enhancement in a PEM fuel cell

This research initiative was triggered by the problems of water management of Polymer Electrolyte Membrane Fuel Cell (PEMFC). In low temperature fuel cells such as PEMFC, some of the water produced after the chemical reaction remains in its liquid state. Excess water produced by the fuel cell must be removed from the system to avoid flooding of the gas diffusion layers (GDL). The GDL is responsible for the transport of reactant gas to the active sites and remove the water produced from the sites. If the GDL is flooded, the supply gas will not be able to reach the reactive sites and the fuel cell fails. The choice of water removal method in this research is to exert a variable asymmetrical force on a liquid droplet. As the drop of liquid is subjected to an external vibrational force in the form of periodic wave, it will begin to oscillate. A fluidic oscillator is capable to produce a pulsating flow using simple balance of momentum fluxes between three impinging jets. By connecting the outputs of the oscillator to the gas channels of a fuel cell, a flow pulsation can be imposed on a water droplet formed within the gas channel during fuel cell operation. The lowest frequency produced by this design is approximately 202 Hz when a 20 inches feed-back port length was used and a supply pressure of 5 psig was introduced. This information was found by setting up a fluidic network with appropriate data acquisition. The components include a fluidic amplifier, valves and fittings, flow meters, a pressure gage, NI-DAQ system, Siglab®, Matlab software and four PCB microphones. The operating environment of the water droplet was reviewed, speed of the sound pressure which travels down the square channel was precisely estimated, and measurement devices were carefully selected. Applicable alternative measurement devices and its application to pressure wave measurement was considered. Methods for experimental setup and possible approaches were recommended, with some discussion of potential problems with implementation of this technique. Some computational fluid dynamic was also performed as an approach to oscillator design.

Looking at Chemistry in Hard to See Places: Understanding Energy Conversion at 1400° Celsius

Solid oxide fuel cells (SOFCs) are promising devices for stationary and portable power and heat generation, because they can use complex fuels such as hydro-carbons, CO, and alcohols. Extreme, non-equilibrium conditions and high tem-peratures (≥ 700 ˚C) required for SOFC operation hamper efforts to understand the mechanisms of component degradation in SOFCs. This talk focuses on new insights into SOFC chemistry and the conversion of carbon-containing fuels (both hydrocarbons and oxygenated) into electricity, carbon dioxide and water, gleaned from a combination of techniques including electrochemical impedance spectroscopy, voltammetry, and vibrational Raman scattering.