Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Triplet excited electronic state switching induced by hydrogen bonding: A transient absorption spectroscopy and time-dependent DFT study

The solvent plays a decisive role in the photochemistry and photophysics of aromatic ketones. Xanthone (XT) is one such aromatic ketone and its triplet-triplet (T-T) absorption spectra show intriguing solvatochromic behavior. Also, the reactivity of XT towards H-atom abstraction shows an unprecedented decrease in protic solvents relative to aprotic solvents. Therefore, a comprehensive solvatochromic analysis of the triplet-triplet absorption spectra of XT was carried out in conjunction with time dependent density functional theory using the ad hoc explicit solvent model approach. A detailed solvatochromic analysis of the T-T absorption bands of XT suggests that the hydrogen bonding interactions are different in the corresponding triplet excited states. Furthermore, the contributions of non-specific and hydrogen bonding interactions towards differential solvation of the triplet states in protic solvents were found to be of equal magnitude. The frontier molecular orbital and electron density difference analysis of the T-1 and T-2 states of XT indicates that the charge redistribution in these states leads to intermolecular hydrogen bond strengthening and weakening, respectively, relative to the S-0 state. This is further supported by the vertical excitation energy calculations of the XT-methanol supra-molecular complex. The intermolecular hydrogen bonding potential energy curves obtained for this complex in the S-0, T-1, and T-2 states support the model. In summary, we propose that the different hydrogen bonding mechanisms exhibited by the two lowest triplet excited states of XT result in a decreasing role of the n pi* triplet state, and are thus responsible for its reduced reactivity towards H-atom abstraction in protic solvents. (C) 2016 AIP Publishing LLC.

Electronic spectroscopy of propofol, propofol homomers ad their hydratated clusters

538 p.

Vibrational predissociation spectroscopy of mass-selected ionic clusters

This thesis presents structural investigations of molecular ions and ionic clusters using vibrational predissociation spectroscopy. Experimentally, a pulsed beam of the mass-selected ion is crossed by a tunable infrared laser beam generated by a Nd:YAG pumped LiNbO_3 optical parametric oscillator. The resulting fragment ion is mass-analyzed and detected, with its intensity as a function of the laser wavelength being the "action" spectrum of the parent ion. In the case of SiH_7^+, we observed a vibrational band centered at 3866 cm^(-1) with clear P, Q, R branches, which is assigned as a perturbed H_2 stretch. The absence of a second H_2 band suggests that the ion forms a symmetric complex with a structure H_2•SiH_3^+•H_2 , in contrast to the species CH_7^+, which has the structure CH_5^+•H_2. The infrared spectra of NO_2^+(H_2O)_n clusters exhibit a marked change with cluster size, indicating that an intracluster reaction occurs with sufficient solvation. Specifically, in NO_2^+(H_2O)_n clusters where n≤3, H_2O binds to a nitronium ion core; but at n=4 the NO_2^+ reacts, transforming the cluster to a new structure of H_3O^+•(H_2O)_(n_2)•HNO_3. For protonated chlorine nitrate, we have observed two distinct isomers previously predicted by ab initio calculations: NO_2^+•(HOC1), the lowest energy isomer, and (ClO)(HO)NO^+, a covalently bonded isomer about 20 kcal/mol higher in energy. Both isomers decompose to NO_2^+ and HOCl upon photo-excitation. These results for HClONO_2^+ lend strong support to the involvement of an ionic mechanism in the reaction of ClONO_2 on polar stratospheric cloud surfaces, a critical step in the dramatic springtime depletion of ozone over Antarctica. Current research activities on halide-solvent clusters and metal-ligand complexes as well as technological improvements of the apparatus are also discussed.

Multiplexed DNA-mediated electrochemistry

The aromatic core of double helical DNA possesses the unique and remarkable ability to form a conduit for electrons to travel over exceptionally long molecular distances. This core of π-stacked nucleobases creates an efficient pathway for charge transfer to proceed that is exquisitely sensitive to even subtle perturbations. Ground state electrochemistry of DNA-modified electrodes has been one of the major techniques used both to investigate and to harness the property of DNA-mediated charge transfer. DNA-modified electrodes have been an essential tool for both gaining insights into the fundamental properties of DNA and, due to the exquisite specificity of DNA-mediated charge transfer for the integrity of the π-stack, for use in next generation diagnostic sensing. Here, multiplexed DNA-modified electrodes are used to (i) gain new insights on the electrochemical coupling of metalloproteins to the DNA π-stack with relevance to the fundaments of in vivo DNA-mediated charge transfer and (ii) enhance the overall sensitivity of DNA-mediated reduction for use in the detection of low abundance diagnostic targets.

First, Methylene Blue (MB′) was covalently attached to DNA through a flexible C12 alkyl linker to yield a new redox reporter for DNA electrochemistry measurements with enhanced sensitivity. Tethered, intercalated MB′ was reduced through DNA-mediated charge transport. The redox signal intensity for MB′-dT-C12-DNA was found to be at least 3 fold larger than that of previously used Nile Blue (NB)-dT-DNA, which is coupled to the base stack via direct conjugation. The signal attenuation, due to an intervening mismatch, and therefore the degree of DNA-mediated reduction, does, however, depend on the DNA film morphology and the backfilling agent used to passivate the surface. These results highlight two possible mechanisms for the reduction of MB′ on the DNA-modified electrode that are distinguishable by their kinetics: reduction mediated by the DNA base pair stack and direct surface reduction of MB′ at the electrode. The extent of direct reduction at the surface can be minimized by overall DNA assembly conditions.

Next, a series of intercalation-based DNA-mediated electrochemical reporters were developed, using a flexible alkane linkage to validate and explore their DNA-mediated reduction. The general mechanism for the reduction of distally bound redox active species, covalently tethered to DNA through flexible alkyl linkages, was established to be an intraduplex DNA-mediated pathway. MB, NB, and anthraquinone were covalently tethered to DNA with three different covalent linkages. The extent of electronic coupling of the reporter was shown to correlate with the DNA binding affinity of the redox active species, supporting an intercalative mechanism. These electrochemical signals were shown to be exceptionally sensitive to a single intervening π-stack perturbation, an AC mismatch, in a densely packed DNA monolayer, which further supports that the reduction is DNA-mediated. Finally, this DNA-mediated reduction of MB occurs primarily via intra- rather than inter duplex intercalation, as probed through varying the proximity and integrity of the neighboring duplex DNA. Further gains to electrochemical sensitivity of our DNA-modified devices were then achieved through the application of electrocatalytic signal amplification using these solvent accessible intercalative reporters, MB-dT-C8, and hemoglobin as a novel electron sink. Electrocatalysis offers an excellent means of electrochemical signal amplification, yet in DNA based sensors, its application has been limited due to strict assembly conditions. We describe the use of hemoglobin as a robust and effective electron sink for electrocatalysis in DNA sensing on low density DNA films. Protein shielding of the heme redox center minimizes direct reduction at the electrode surface and permits assays on low density DNA films. Electrocatalysis of MB that is covalently tethered to the DNA by a flexible alkyl linkage allows for efficient interactions with both the base stack and hemoglobin. Consistent suppression of the redox signal upon incorporation of single CA mismatch in the DNA oligomer demonstrates that both the unamplified and the electrocatalytically amplified redox signals are generated through DNA-mediated charge transport. Electrocatalysis with hemoglobin is robust: it is stable to pH and temperature variations. The utility and applicability of electrocatalysis with hemoglobin is demonstrated through restriction enzyme detection, and an enhancement in sensitivity permits femtomole DNA sampling.

Finally, we expanded the application of our multiplexed DNA-modified electrodes to the electrochemical characterization of DNA-bound proteins containing [4Fe-4S] clusters. DNA-modified electrodes have become an essential tool for the characterization of the redox chemistry of DNA repair proteins that contain redox cofactors. Multiplexed analysis of EndonucleaseIII (EndoIII), a DNA repair protein containing a [4Fe-4S] cluster known to be accessible via DNA-mediated charge transport, elucidated subtle differences in the electrochemical behavior as a function of DNA morphology. DNA-bound EndoIII is seen to have two different electron transfer pathways for reduction, either through the DNA base stack or through direct surface reduction. Closely packed DNA films, where the protein has limited surface accessibility, produce electrochemical signals reflecting electron transfer that is DNA-mediated. The electrochemical comparison of EndoIII mutants, including a new family of mutations altering the electrostatics surrounding the [4Fe-4S] cluster, was able to be quantitatively performed. While little change in the midpoint potential was found for this family of mutants, significant variations in the efficiency of DNA-mediated electron transfer were apparent. Based on the stability of these proteins, examined by circular dichroism, we propose that the electron transfer pathway can be perturbed not only by the removal of aromatic residues, but also through changes in solvation near the cluster.

Effects of self energy of the ions on the double layer structure and properties at the dielectric interface

Although numerous theoretical efforts have been put forth, a systematic, unified and predictive theoretical framework that is able to capture all the essential physics of the interfacial behaviors of ions, such as the Hofmeister series effect, Jones-Ray effect and the salt effect on the bubble coalescence remain an outstanding challenge. The most common approach to treating electrostatic interactions in the presence of salt ions is the Poisson-Boltzmann (PB) theory. However, there are many systems for which the PB theory fails to offer even a qualitative explanation of the behavior, especially for ions distributed in the vicinity of an interface with dielectric contrast between the two media (like the water-vapor/oil interface). A key factor missing in the PB theory is the self energy of the ion.

In this thesis, we develop a self-consistent theory that treats the electrostatic self energy (including both the short-range Born solvation energy and the long-range image charge interactions), the nonelectrostatic contribution of the self energy, the ion-ion correlation and the screening effect systematically in a single framework. By assuming a finite charge spread of the ion instead of using the point-charge model, the self energy obtained by our theory is free of the divergence problems and gives a continuous self energy across the interface. This continuous feature allows ions on the water side and the vapor/oil side of the interface to be treated in a unified framework. The theory involves a minimum set of parameters of the ion, such as the valency, radius, polarizability of the ions, and the dielectric constants of the medium, that are both intrinsic and readily available. The general theory is first applied to study the thermodynamic property of the bulk electrolyte solution, which shows good agreement with the experiment result for predicting the activity coefficient and osmotic coefficient.

Next, we address the effect of local Born solvation energy on the bulk thermodynamics and interfacial properties of electrolyte solution mixtures. We show that difference in the solvation energy between the cations and anions naturally gives rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The miscibility of the mixture can either increases or decreases depending on the competition between the solvation energy and translation entropy of the ions. The interfacial tension shows a non-monotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations, and decreases approximately as the square root of the salt concentration for dilute solutions, which is in agreement with the Jones-Ray effect observed in experiment.

Next, we investigate the image effects on the double layer structure and interfacial properties near a single charged plate. We show that the image charge repulsion creates a depletion boundary layer that cannot be captured by a regular perturbation approach. The correct weak-coupling theory must include the self-energy of the ion due to the image charge interaction. The image force qualitatively alters the double layer structure and properties, and gives rise to many non-PB effects, such as nonmonotonic dependence of the surface energy on concentration and charge inversion. The image charge effect is then studied for electrolyte solutions between two plates. For two neutral plates, we show that depletion of the salt ions by the image charge repulsion results in short-range attractive and long-range repulsive forces. If cations and anions are of different valency, the asymmetric depletion leads to the formation of an induced electrical double layer. For two charged plates, the competition between the surface charge and the image charge effect can give rise to like- charge attraction.

Then, we study the inhomogeneous screening effect near the dielectric interface due to the anisotropic and nonuniform ion distribution. We show that the double layer structure and interfacial properties is drastically affected by the inhomogeneous screening if the bulk Debye screening length is comparable or smaller than the Bjerrum length. The width of the depletion layer is characterized by the Bjerrum length, independent of the salt concentration. We predict that the negative adsorption of ions at the interface increases linearly with the salt concentration, which cannot be captured by either the bulk screening approximation or the WKB approximation. For asymmetric salt, the inhomogeneous screening enhances the charge separation in the induced double layer and significantly increases the value of the surface potential.

Finally, to account for the ion specificity, we study the self energy of a single ion across the dielectric interface. The ion is considered to be polarizable: its charge distribution can be self-adjusted to the local dielectric environment to minimize the self energy. Using intrinsic parameters of the ions, such as the valency, radius, and polarizability, we predict the specific ion effect on the interfacial affinity of halogen anions at the water/air interface, and the strong adsorption of hydrophobic ions at the water/oil interface, in agreement with experiments and atomistic simulations.

The theory developed in this work represents the most systematic theoretical technique for weak-coupling electrolytes. We expect the theory to be more useful for studying a wide range of structural and dynamic properties in physicochemical, colloidal, soft-matter and biophysical systems.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

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.

Quantum Mechanics Studies of Fuel Cell Catalysts and Proton Conducting Ceramics with Validation by Experiment

We carried out quantum mechanics (QM) studies aimed at improving the performance of hydrogen fuel cells. This led to predictions of improved materials, some of which were subsequently validated with experiments by our collaborators.

In part I, the challenge was to find a replacement for the Pt cathode that would lead to improved performance for the Oxygen Reduction Reaction (ORR) while remaining stable under operational conditions and decreasing cost. Our design strategy was to find an alloy with composition Pt3M that would lead to surface segregation such that the top layer would be pure Pt, with the second and subsequent layers richer in M. Under operating conditions we expect the surface to have significant O and/or OH chemisorbed on the surface, and hence we searched for M that would remain segregated under these conditions. Using QM we examined surface segregation for 28 Pt₃M alloys, where M is a transition metal. We found that only Pt₃Os and Pt₃Ir showed significant surface segregation when O and OH are chemisorbed on the catalyst surfaces. This result indicates that Pt₃Os and Pt3Ir favor formation of a Pt-skin surface layer structure that would resist the acidic electrolyte corrosion during fuel cell operation environments. We chose to focus on Os because the phase diagram for Pt-Ir indicated that Pt-Ir could not form a homogeneous alloy at lower temperature. To determine the performance for ORR, we used QM to examine all intermediates, reaction pathways, and reaction barriers involved in the processes for which protons from the anode reactions react with O₂ to form H₂O. These QM calculations used our Poisson-Boltzmann implicit solvation model include the effects of the solvent (water with dielectric constant 78 with pH 7 at 298K). We found that the rate determination step (RDS) was the O_ad hydration reaction (O_ad + H₂O_ad -> OH_ad + OH_ad) in both cases, but that the barrier for pure Pt of 0.50 eV is reduced to 0.48 eV for Pt₃Os, which at 80 degrees C would increase the rate by 218%. We collaborated with the Pu-Wei Wu’s group to carry out experiments, where we found that the dealloying process-treated Pt2Os catalyst showed two-fold higher activity at 25 degrees C than pure Pt and that the alloy had 272% improved stability, validating our theoretical predictions.

We also carried out similar QM studies followed by experimental validation for the Os/Pt core-shell catalyst fabricated by the underpotential deposition (UPD) method. The QM results indicated that the RDS for ORR is a compromise between the OOH formation step (0.37 eV for Pt, 0.23 eV for Pt_2ML/Os core-shell) and H₂O formation steps (0.32 eV for Pt, 0.22 eV for Pt_2ML/Os core-shell). We found that Pt_2ML/Os has the highest activity (compared to pure Pt and to the Pt₃Os alloy) because the 0.37 eV barrier decreases to 0.23 eV. To understand what aspects of the core shell structure lead to this improved performance, we considered the effect on ORR of compressing the alloy slab to the dimensions of pure Pt. However this had little effect, with the same RDS barrier 0.37 eV. This shows that the ligand effect (the electronic structure modification resulting from the Os substrate) plays a more important role than the strain effect, and is responsible for the improved activity of the core- shell catalyst. Experimental materials characterization proves the core-shell feature of our catalyst. The electrochemical experiment for Pt_2ML/Os/C showed 3.5 to 5 times better ORR activity at 0.9V (vs. NHE) in 0.1M HClO₄ solution at 25 degrees C as compared to those of commercially available Pt/C. The excellent correlation between experimental half potential and the OH binding energies and RDS barriers validate the feasibility of predicting catalyst activity using QM calculation and a simple Langmuir–Hinshelwood model.

In part II, we used QM calculations to study methane stream reforming on a Ni-alloy catalyst surfaces for solid oxide fuel cell (SOFC) application. SOFC has wide fuel adaptability but the coking and sulfur poisoning will reduce its stability. Experimental results suggested that the Ni4Fe alloy improves both its activity and stability compared to pure Ni. To understand the atomistic origin of this, we carried out QM calculations on surface segregation and found that the most stable configuration for Ni₄Fe has a Fe atom distribution of (0%, 50%, 25%, 25%, 0%) starting at the bottom layer. We calculated that the binding of C atoms on the Ni4Fe surface is 142.9 Kcal/mol, which is about 10 Kcal/mol weaker compared to the pure Ni surface. This weaker C binding energy is expected to make coke formation less favorable, explaining why Ni₄Fe has better coking resistance. This result confirms the experimental observation. The reaction energy barriers for CHx decomposition and C binding on various alloy surface, Ni₄X (X=Fe, Co, Mn, and Mo), showed Ni₄Fe, Ni₄Co, and Fe₄Mn all have better coking resistance than pure Ni, but that only Ni₄Fe and Fe₄Mn have (slightly) improved activity compared to pure Ni.

In part III, we used QM to examine the proton transport in doped perovskite-ceramics. Here we used a 2x2x2 supercell of perovskite with composition Ba₈X₇M₁(OH)₁O₂₃ where X=Ce or Zr and M=Y, Gd, or Dy. Thus in each case a 4⁺ X is replace by a 3⁺ M plus a proton on one O. Here we predicted the barriers for proton diffusion allowing both includes intra-octahedron and inter-octahedra proton transfer. Without any restriction, we only observed the inter-octahedra proton transfer with similar energy barrier as previous computational work but 0.2 eV higher than experimental result for Y doped zirconate. For one restriction in our calculations is that the O_donor-O_acceptor atoms were kept at fixed distances, we found that the barrier difference between cerates/zirconates with various dopants are only 0.02~0.03 eV. To fully address performance one would need to examine proton transfer at grain boundaries, which will require larger scale ReaxFF reactive dynamics for systems with millions of atoms. The QM calculations used here will be used to train the ReaxFF force field.

Investigations of DNA-mediated protein oxidation

DNA charge transport (CT) involves the efficient transfer of electrons or electron holes through the DNA π-stack over long molecular distances of at least 100 base-pairs. Despite this shallow distance dependence, DNA CT is sensitive to mismatches or lesions that disrupt π-stacking and is critically dependent on proper electronic coupling of the donor and acceptor moieties into the base stack. Favorable DNA CT is very rapid, occurring on the picosecond timescale. Because of this speed, electron holes equilibrate along the DNA π-stack, forming a characteristic pattern of DNA damage at low oxidation potential guanine multiplets. Furthermore, DNA CT may be used in a biological context. DNA processing enzymes with 4Fe4S clusters can perform DNA-mediated electron transfer (ET) self-exchange reactions with other 4Fe4S cluster proteins, even if the proteins are quite dissimilar, as long as the DNA-bound [4Fe4S]^3+/2+ redox potentials are conserved. This mechanism would allow low copy number DNA repair proteins to find their lesions efficiently within the cell. DNA CT may also be used biologically for the long-range, selective activation of redox-active transcription factors. Within this work, we pursue other proteins that may utilize DNA CT within the cell and further elucidate aspects of the DNA-mediated ET self-exchange reaction of 4Fe4S cluster proteins.

Dps proteins, bacterial mini-ferritins that protect DNA from oxidative stress, are implicated in the survival and virulence of pathogenic bacteria. One aspect of their protection involves ferroxidase activity, whereby ferrous iron is bound and oxidized selectively by hydrogen peroxide, thereby preventing formation of damaging hydroxyl radicals via Fenton chemistry. Understanding the specific mechanism by which Dps proteins protect the bacterial genome could inform the development of new antibiotics. We investigate whether DNA-binding E. coli Dps can utilize DNA CT to protect the genome from a distance. An intercalating ruthenium photooxidant was employed to generate oxidative DNA damage via the flash-quench technique, which localizes to a low potential guanine triplet. We find that Dps loaded with ferrous iron, in contrast to Apo-Dps and ferric iron-loaded Dps which lack available reducing equivalents, significantly attenuates the yield of oxidative DNA damage at the guanine triplet. These data demonstrate that ferrous iron-loaded Dps is selectively oxidized to fill guanine radical holes, thereby restoring the integrity of the DNA. Luminescence studies indicate no direct interaction between the ruthenium photooxidant and Dps, supporting the DNA-mediated oxidation of ferrous iron-loaded Dps. Thus DNA CT may be a mechanism by which Dps efficiently protects the genome of pathogenic bacteria from a distance.

Further work focused on spectroscopic characterization of the DNA-mediated oxidation of ferrous iron-loaded Dps. X-band EPR was used to monitor the oxidation of DNA-bound Dps after DNA photooxidation via the flash-quench technique. Upon irradiation with poly(dGdC)₂, a signal arises with g = 4.3, consistent with the formation of mononuclear high-spin Fe(III) sites of low symmetry, the expected oxidation product of Dps with one iron bound at each ferroxidase site. When poly(dGdC)₂ is substituted with poly(dAdT)₂, the yield of Dps oxidation is decreased significantly, indicating that guanine radicals facilitate Dps oxidation. The more favorable oxidation of Dps by guanine radicals supports the feasibility of a long-distance protection mechanism via DNA CT where Dps is oxidized to fill guanine radical holes in the bacterial genome produced by reactive oxygen species.

We have also explored possible electron transfer intermediates in the DNA-mediated oxidation of ferrous iron-loaded Dps. Dps proteins contain a conserved tryptophan residue in close proximity to the ferroxidase site (W52 in E. coli Dps). In comparison to WT Dps, in EPR studies of the oxidation of ferrous iron-loaded Dps following DNA photooxidation, W52Y and W52A mutants were deficient in forming the characteristic EPR signal at g = 4.3, with a larger deficiency for W52A compared to W52Y. In addition to EPR, we also probed the role of W52 Dps in cells using a hydrogen peroxide survival assay. Bacteria containing W52Y Dps survived the hydrogen peroxide challenge more similarly to those containing WT Dps, whereas cells with W52A Dps died off as quickly as cells without Dps. Overall, these results suggest the possibility of W52 as a CT hopping intermediate.

DNA-modified electrodes have become an essential tool for the study of the redox chemistry of DNA processing enzymes with 4Fe4S clusters. In many cases, it is necessary to investigate different complex samples and substrates in parallel in order to elucidate this chemistry. Therefore, we optimized and characterized a multiplexed electrochemical platform with the 4Fe4S cluster base excision repair glycosylase Endonuclease III (EndoIII). Closely packed DNA films, where the protein has limited surface accessibility, produce EndoIII electrochemical signals sensitive to an intervening mismatch, indicating a DNA-mediated process. Multiplexed analysis allowed more robust characterization of the CT-deficient Y82A EndoIII mutant, as well as comparison of a new family of mutations altering the electrostatics surrounding the 4Fe4S cluster in an effort to shift the reduction potential of the cluster. While little change in the DNA-bound midpoint potential was found for this family of mutants, likely indicating the dominant effect of DNA-binding on establishing the protein redox potential, significant variations in the efficiency of DNA-mediated electron transfer were apparent. On the basis of the stability of these proteins, examined by circular dichroism, we proposed that the electron transfer pathway in EndoIII can be perturbed not only by the removal of aromatic residues but also through changes in solvation near the cluster.

While the 4Fe4S cluster of EndoIII is relatively insensitive to oxidation and reduction in solution, we have found that upon DNA binding, the reduction potential of the [4Fe4S]^3+/2+ couple shifts negatively by approximately 200 mV, bringing this couple into a physiologically relevant range. Demonstrated using electrochemistry experiments in the presence and absence of DNA, these studies do not provide direct molecular evidence for the species being observed. Sulfur K-edge X-ray absorbance spectroscopy (XAS) can be used to probe directly the covalency of iron-sulfur clusters, which is correlated to their reduction potential. We have shown that the Fe-S covalency of the 4Fe4S cluster of EndoIII increases upon DNA binding, stabilizing the oxidized [4Fe4S]³⁺ cluster, consistent with a negative shift in reduction potential. The 7% increase in Fe-S covalency corresponds to an approximately 150 mV shift, remarkably similar to DNA electrochemistry results. Therefore we have obtained direct molecular evidence for the shift in 4Fe4S reduction potential of EndoIII upon DNA binding, supporting the feasibility of our model whereby these proteins can utilize DNA CT to cooperate in order to efficiently find DNA lesions inside cells.

In conclusion, in this work we have explored the biological applications of DNA CT. We discovered that the DNA-binding bacterial ferritin Dps can protect the bacterial genome from a distance via DNA CT, perhaps contributing to pathogen survival and virulence. Furthermore, we optimized a multiplexed electrochemical platform for the study of the redox chemistry of DNA-bound 4Fe4S cluster proteins. Finally, we have used sulfur K-edge XAS to obtain direct molecular evidence for the negative shift in 4Fe4S cluster reduction potential of EndoIII upon DNA binding. These studies contribute to the understanding of DNA-mediated protein oxidation within cells.

Effect of the interplay between protein and surface on the properties of adsorbed protein layers.

Although protein adsorption to surface is a common phenomenon, investigation of the process is challenging due to the complexity of the interplay between external factors, protein and surface properties. Therefore experimental approaches have to measure the properties of adsorbed protein layers with high accuracy in order to achieve a comprehensive description of the process. To this end, we used a combination of two biosensing techniques, dual polarization interferometry and quartz crystal microbalance with dissipation. From this, we are able to extract surface coverage values, layer structural parameters, water content and viscoelastic properties to examine the properties of protein layers formed at the liquid/solid interface. Layer parameters were examined upon adsorption of proteins of varying size and structural properties, on surfaces with opposite polarity. We show that "soft" proteins such as unfolded α-synuclein and high molecular weight albumin are highly influenced by the surface polarity, as they form a highly diffuse and hydrated layer on the hydrophilic silica surface as opposed to the denser, less hydrated layer formed on a hydrophobic methylated surface. These layer properties are a result of different orientations and packing of the proteins. By contrast, lysozyme is barely influenced by the surface polarity due to its intrinsic structural stability. Interestingly, we show that for a similar molecular weight, the unfolded α-synuclein forms a layer with the highest percentage of solvation not related to surface coverage but resulting from the highest water content trapped within the protein. Together, these data reveal a trend in layer properties highlighting the importance of the interplay between protein and surface for the design of biomaterials.

Effect of the interplay between protein and surface on the properties of adsorbed protein layers

Although protein adsorption to surface is a common phenomenon, investigation of the process is challenging due to the complexity of the interplay between external factors, protein and surface properties. Therefore experimental approaches have to measure the properties of adsorbed protein layers with high accuracy in order to achieve a comprehensive description of the process. To this end, we used a combination of two biosensing techniques, dual polarization interferometry and quartz crystal microbalance with dissipation. From this, we are able to extract surface coverage values, layer structural parameters, water content and viscoelastic properties to examine the properties of protein layers formed at the liquid/solid interface. Layer parameters were examined upon adsorption of proteins of varying size and structural properties, on surfaces with opposite polarity. We show that "soft" proteins such as unfolded α-synuclein and high molecular weight albumin are highly influenced by the surface polarity, as they form a highly diffuse and hydrated layer on the hydrophilic silica surface as opposed to the denser, less hydrated layer formed on a hydrophobic methylated surface. These layer properties are a result of different orientations and packing of the proteins. By contrast, lysozyme is barely influenced by the surface polarity due to its intrinsic structural stability. Interestingly, we show that for a similar molecular weight, the unfolded α-synuclein forms a layer with the highest percentage of solvation not related to surface coverage but resulting from the highest water content trapped within the protein. Together, these data reveal a trend in layer properties highlighting the importance of the interplay between protein and surface for the design of biomaterials. © 2014 The Authors.

Polymerization of ionic liquid-based microemulsions: A versatile method for the synthesis of polymer electrolytes

电子邮箱fyan@suda.edu.cn

Effects of phase-breaking on long-range charge transfer in DNA: Partially-coherent-tunneling model study

The mechanism of hole charge transfer in DNA of various lengths and sequences is investigated based on a partially coherent tunneling theory (Zhang et al., J Chem Phys 117:4578, 2002), where the effects of phase-breaking in adenine-thymine and guanine-cytosine base pairs are treated on equal foot. This work aims at providing a self-consistent microscopic interpretation for rate experiments on various DNA systems. We will also clarify the condition under which the simple superexchange-mediated-hopping picture is valid, and make some comments on the further development of present theory.

Scattering matrix approach to electronic dephasing in long-range electron transfer

Based on the Buttiker dephasing model, we propose an analytical scattering matrix approach to the long-range electron transfer phenomena. The present efficient scheme smoothly interpolates between the superexchange and the sequential hopping mechanisms. Various properties such as the drastic dephasing-assisted enhancement and turnover behaviors are demonstrated in good agreement with those obtained via the dynamical reduced density-matrix methods. These properties are further elucidated as results of the interplay among the dephasing strength, the tunneling parameter, and the bridge length of the electron transfer system. (C) 2001 American Institute of Physics.