Einzelmolekülspektroskopische und quantenchemische Untersuchungen zum elektronischen Energietransfer

In der vorliegenden Arbeit wurden Untersuchungen zum Mechanismus, der Dynamik und der Kontrolle des elektronischen Energietransfers in multichromophoren Modellsystemen durchgeführt. Als Untersuchungsmethoden wurden hauptsächlich die konfokale Einzelmolekülspektroskopie und die Quantenchemie eingesetzt. Der Aufbau des Einzelmolekülmikroskops wurde bezüglich der Anregungs- und Detektionskomponenten variiert, um die unterschiedlichen Experimente durchzuführen. Die quantenchemischen Rechnungen wurden auf Dichtefunktional- und Coupled-Cluster-Niveau durchgeführt. Die aus den Rechnungen erhaltenen zusätzlichen Informationen über experimentell zum Teil schwer zugängliche Eigenschaften der Farbstoffe unterstützten die Interpretation der experimentellen Befunde. rnIn früheren Untersuchungen der AG Basché wurden die Energietransfer-Raten von Donor-Akzeptor-Systemen gemessen, die erhebliche Abweichungen von nach der Förster-Theorie berechneten Raten zeigten. Daher war ein Ziel der vorliegenden Arbeit, diese Abweichungen zu erklären. Zu diesem Zweck wurde die Geometrie der Diaden experimentell untersucht, sowie die elektronische Kopplung zwischen den Chromophoren quantenchemisch berechnet. Die relative Orientierung der Chromophore in den Diaden wurde in einem Einzelmolekül-Experiment mit rotierender Anregungspolarisation abgefragt. Die erhaltenen Winkelverteilungen konnten schließlich eindeutig auf die Flexibilität der die Chromophore verbrückenden Oligophenyl-Einheiten zurückgeführt werden. Die Unterschiede der gemessenen Energietransfer-Raten zu den nach der Förster-Theorie ermittelten Werten konnten jedoch nicht über die molekulare Flexibilität der Systeme erklärt werden. Aufklärung über die Diskrepanzen zur Förster-Theorie ergaben die quantenchemischen Rechnungen. In Rahmen dieser Arbeit wurde zum ersten Mal die Coupled-Cluster-Theorie zur Berechnung der elektronischen Kopplung eingesetzt. Die Betrachtung der isolierten Chromophore reichte aber nicht aus, um die gemessenen Abweichungen von der Förster-Theorie zu erklären. Erst über die Berücksichtigung der molekularen Brücke konnten die gefunden Abweichungen erklärt werden. Die deutliche Verstärkung der elektronischen Kopplung ist auf die Polarisierbarkeit der Brücke zurückzuführen.rnNach diesen Betrachtungen stand die Kontrolle des Energietransfers im Fokus der weiteren Untersuchungen. In den durchgeführten Einzelmolekülexperimenten wurden die Chromophore der Donor-Akzeptor-Systeme selektiv mit zwei Laserpulsen unterschiedlicher Wellenlänge angeregt. Beim gleichzeitigen Anregen beider Chromophore wurde Singulett-Singulett-Annihilation (SSA) induziert, ein Energietransferprozess, bei dem die Anregungsenergie vom vorigen Akzeptor zum vorigen Donor übertragen wird. Da über SSA Fluoreszenzphotonen gelöscht wurden, konnte über den Abstand der Laserpulse die Fluoreszenzintensität des einzelnen Moleküls moduliert werden. Konzeptionell verwandte Einzelmolekülexperimente wurden an einem weiteren molekularen System durchgeführt, das aus einem Kern und einer Peripherie bestand. Fluoreszenzauszeiten des Gesamtsystems bei selektiver Anregung des Kerns wurden auf die Population eines Triplett-Zustandes zurückgeführt, der die Fluoreszenz der Peripherie löschte. rnAbschließend wurde der SSA-Prozess zwischen zwei gleichartigen Chromophoren untersucht. Es wurde eine Methode entwickelt, die es zum ersten Mal erlaubte, die SSA-Zeitkonstante individueller Moleküle zu bestimmen. Hierfür wurden die Daten der gemessenen Photonen-Koinzidenzhistogramme mittels eines im Rahmen dieser Arbeit hergeleiteten analytischen Zusammenhangs ausgewertet, der über Monte-Carlo-Simulationen bestätigt wurde.

Dichtefunktionaltheorie für eine Flüssigkeit harter Scheiben auf Grafikkarten

Zweidimensionale Flüssigkeiten harter Scheiben sind in der Regel einfach zu simulieren, jedoch überraschend schwer theoretisch zu beschreiben. Trotz ihrer hohen Relevanz bleiben die meisten theoretischen Ansätze qualitativ. Hier wird eine Dichtefunktionaltheorie (DFT) vorgestellt, die erstmalig die Struktur solcher Flüssigkeiten bei hohen Dichten korrekt beschreibt und den Ansatz des Gefrierübergangs abbildet.rnEs wird gezeigt, dass der Ansatz der Fundamentalmaßtheorie zu einem solchen Funktional führt. Dabei werden sowohl Dichteverteilungen um ein Testteilchen als auch Zweiteilchen-Korrelationsfunktionen untersucht.rnGrafikkarten bieten sehr hohe Recheneffizienz und ihr Einsatz in der Wissenschaft nimmt stetig zu. In dieser Arbeit werden die Vor- und Nachteile der Grafikkarte für wissenschaftliche Berechnungen erörtert und es wird gezeigt, dass die Berechnung der DFT auf Grafikkarten effizient ausgeführt werden kann. Es wird ein Programm entwickelt, dass dies umsetzt. Dabei wird gezeigt, dass die Ergebnisse einfacher (bekannter) Funktionale mit denen von CPU-Berechnungen übereinstimmen, so dass durch die Nutzung der Grafikkarte keine systematischen Fehler zu erwarten sind.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Efficient and Accurate Characterization of the Bergman Cyclization for an Expanded Substructure of Esperamicin A1

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal·mol−1 (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.

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.

Conformationally Controlled Electron Delocalization in n-Type Rods: Synthesis, Structure, and Optical, Electrochemical, and Spectroelectrochemical Properties of Dicyanocyclophanes

A series of dicyanobiphenyl-cyclophanes 1-6 with various pi-backbone conformations and characteristic n-type semiconductor properties is presented. Their synthesis, optical, structural, electrochemical, spectroelectrochemical, and packing properties are investigated. The X-ray crystal structures of all n-type rods allow the systematic correlation of structural features with physical properties. In addition, the results are supported by quantum mechanical calculations based on density functional theory. A two-step reduction process is observed for all n-type rods, in which the first step is reversible. The potential gap between the reduction processes depends linearly on the cos(2) value of the torsion angle phi between the pi-systems. Similarly, optical absorption spectroscopy shows that the vertical excitation energy of the conjugation band correlates with the cos(2) value of the torsion angle phi. These correlations demonstrate that the fixed intramolecular torsion angle phi is the dominant factor determining the extent of electron delocalization in these model compounds, and that the angle phi measured in the solid-state structure is a good proxy for the molecular conformation in solution. Spectroelectrochemical investigations demonstrate that conformational rigidity is maintained even in the radical anion form. In particular, the absorption bands corresponding to the SOMO-LUMO+i transitions are shifted bathochromically, whereas the absorption bands corresponding to the HOMO-SOMO transition are shifted hypsochromically with increasing torsion angle phi.

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.

Effect of external electric field on hydrogen adsorption over activated carbon separated by dielectric materials

Energy crisis and worldwide environmental problem make hydrogen a prospective energy carrier. However, storage and transportation of hydrogen in large quantities at small volume is currently not practical. Lots of materials and devices have been developed for storage hydrogen, but to today none is able to meet the DOE targets. Activated carbon has been found to be a good hydrogen adsorbent due to its high surface area. However, the weak van der Waals force between hydrogen and the adsorbent has limited the adsorption capacity. Previous studies have found that enhanced adsorption can be obtained with applied electric field. Stronger interaction between the polarized hydrogen and the charged sorbents under high voltage is considered as the reason. This study was initiated to investigate if the adsorption can be further enhanced when the activated carbon particles are separated with a dielectric coating. Dielectric TiO2 nanoparticles were first utilized. Hydrogen adsorption measurements on the TiO2-coated carbon materials, with or without an external electric field, were made. The results showed that the adsorption capacity enhancement increased with the increasing amount of TiO2 nanoparticles with an applied electric field. Since the hydrogen adsorption capacity on TiO2 particles is very low and there is no hydrogen adsorption enhancement on TiO2 particles alone when electric field is applied, the effect of dielectric coating is demonstrated. Another set of experiments investigated the behavior of hydrogen adsorption over TiO2-coated activated carbon under various electric potentials. The results revealed that the hydrogen adsorption first increased and then decreased with the increase of electric field. The improved storage was due to a stronger interaction between charged carbon surface and polarized hydrogen molecule caused by field induced polarization of TiO2 coating. When the electric field was sufficient to cause considerable ionization of hydrogen, the decrease of hydrogen adsorption occurred. The current leak detected at 3000 V was a sign of ionization of hydrogen. Experiments were also carried out to examine the hydrogen adsorption performances over activated carbon separated by other dielectric materials, MgO, ZnO and BaTiO3, respectively. For the samples partitioned with MgO and ZnO, the measurements with and without an electric field indicated negligible differences. Electric field enhanced adsorption has been observed on the activated carbon separated with BaTiO3, a material with unusually high dielectric constant. Corresponding computational calculations using Density Functional Theory have been performed on hydrogen interaction with charged TiO2 molecule as well as TiO2 molecule, coronene and TiO2-doped coronene in the presence of an electric field. The simulated results were consistent with the observations from experiments, further confirming the proposed hypotheses.

Structure and properties of graphene nano disks (GND) with and without edge-dopants

Graphene is one of the most important materials. In this research, the structures and properties of graphene nano disks (GND) with a concentric shape were investigated by Density Functional Theory (DFT) calculations, in which the most effective DFT methods - B3lyp and Pw91pw91 were employed. It was found that there are two types of edges - Zigzag and Armchair in concentric graphene nano disks (GND). The bond length between armchair-edge carbons is much shorter than that between zigzag-edge carbons. For C24 GND that consists of 24 carbon atoms, only armchair edge with 12 atoms is formed. For a GND larger than the C24 GND, both armchair and zigzag edges co-exist. Furthermore, when the number of carbon atoms in armchair-edge are always 12, the number of zigzag-edge atoms increases with increasing the size of a GND. In addition, the stability of a GND is enhanced with increasing its size, because the ratio of edge-atoms to non-edge-atoms decreases. The size effect of a graphene nano disk on its HOMO-LUMO energy gap was evaluated. C6 and C24 GNDs possess HOMO-LUMO gaps of 1.7 and 2.1eV, respectively, indicating that they are semi-conductors. In contrast, C54 and C96 GNDs are organic metals, because their HOMO-LUMO gaps are as low as 0.3 eV. The effect of doping foreign atoms to the edges of GNDs on their structures, stabilities, and HOMO-LUMO energy gaps were also examined. When foreign atoms are attached to the edge of a GND, the original unsaturated carbon atoms become saturated. As a result, both of the C-C bonds lengths and the stability of a GND increase. Furthermore, the doping effect on the HOMO-LUMO energy gap is dependent on the type of doped atoms. The doping H, F, or OH into the edge of a GND increases its HOMO-LUMO energy gap. In contrast, a Li-doped GND has a lower HOMO-LUMO energy gap than that without doping. Therefore, Li-doping can increase the electrical conductance of a GND, whereas H, F, or OH-doping decreases its conductance.

Controlling electronic and magnetic properties of ultra narrow multilayered nanowires

Interest in the study of magnetic/non-magnetic multilayered structures took a giant leap since Grünberg and his group established that the interlayer exchange coupling (IEC) is a function of the non-magnetic spacer width. This interest was further fuelled by the discovery of the phenomenal Giant Magnetoresistance (GMR) effect. In fact, in 2007 Albert Fert and Peter Grünberg were awarded the Nobel Prize in Physics for their contribution to the discovery of GMR. GMR is the key property that is being used in the read-head of the present day computer hard drive as it requires a high sensitivity in the detection of magnetic field. The recent increase in demand for device miniaturization encouraged researchers to look for GMR in nanoscale multilayered structures. In this context, one dimensional(1-D) multilayerd nanowire structure has shown tremendous promise as a viable candidate for ultra sensitive read head sensors. In fact, the phenomenal giant magnetoresistance(GMR) effect, which is the novel feature of the currently used multilayered thin film, has already been observed in multilayered nanowire systems at ambient temperature. Geometrical confinement of the supper lattice along the 2-dimensions (2-D) to construct the 1-D multilayered nanowire prohibits the minimization of magnetic interaction- offering a rich variety of magnetic properties in nanowire that can be exploited for novel functionality. In addition, introduction of non-magnetic spacer between the magnetic layers presents additional advantage in controlling magnetic properties via tuning the interlayer magnetic interaction. Despite of a large volume of theoretical works devoted towards the understanding of GMR and IEC in super lattice structures, limited theoretical calculations are reported in 1-D multilayered systems. Thus to gauge their potential application in new generation magneto-electronic devices, in this thesis, I have discussed the usage of first principles density functional theory (DFT) in predicting the equilibrium structure, stability as well as electronic and magnetic properties of one dimensional multilayered nanowires. Particularly, I have focused on the electronic and magnetic properties of Fe/Pt multilayered nanowire structures and the role of non-magnetic Pt spacer in modulating the magnetic properties of the wire. It is found that the average magnetic moment per atom in the nanowire increases monotonically with an ~1/(N(Fe)) dependance, where N(Fe) is the number of iron layers in the nanowire. A simple model based upon the interfacial structure is given to explain the 1/(N(Fe)) trend in magnetic moment obtained from the first principle calculations. A new mechanism, based upon spin flip with in the layer and multistep electron transfer between the layers, is proposed to elucidate the enhancement of magnetic moment of Iron atom at the Platinum interface. The calculated IEC in the Fe/Pt multilayered nanowire is found to switch sign as the width of the non-magnetic spacer varies. The competition among short and long range direct exchange and the super exchange has been found to play a key role for the non-monotonous sign in IEC depending upon the width of the Platinum spacer layer. The calculated magnetoresistance from Julliere's model also exhibit similar switching behavior as that of IEC. The universality of the behavior of exchange coupling has also been looked into by introducing different non-magnetic spacers like Palladium, Copper, Silver, and Gold in between magnetic Iron layers. The nature of hybridization between Fe and other non-magnetic spacer is found to dictate the inter layer magnetic interaction. For example, in Fe/Pd nanowire the d-p hybridization in two spacer layer case favors anti-ferromagnetic (AFM) configuration over ferromagnetic (FM) configuration. However, the hybridization between half-filled Fe(d) and filled Cu(p) state in Fe/Cu nanowire favors FM coupling in the 2-spacer system.

First-principles studies of boron nanostructures

Boron is an 'electron deficient' element which has a rather fascinating chemical versatility. In the solid state, the elemental boron has neither a pure covalent nor a pure metallic character. As a result, its vast structural dimensionally and peculiar bonding features hold a unique place among other elements in the periodic table. In order to understand and properly describe these unusual bonding features, a detailed and systematic theoretical study is needed. In this work, I will show that some of the qualitative features of boron nanostructures, including clusters, sheets and nanotubes can easily be extracted from the results of first principles calculations based on density functional theory. Specifically, the size-dependent evolution of topological structures and bonding characteristics of boron clusters, Bn will be discussed. Based on the scenario observed in the boron clusters, the unique properties of boron sheets and boron nanotubes will be described. Moreover, the ballistic electron transport in single-walled carbon nanotubes will be considered. It is expected that the theoretical results obtained in the present thesis will initiate further studies on boron nanostructures, which will be helpful in understanding, designing and realizing boron-based nanoscale devices.

Development of a high performance parallel computing platform and its use in the study of nanostructures : clusters, sheets and tubes .

Small clusters of gallium oxide, technologically important high temperature ceramic, together with interaction of nucleic acid bases with graphene and small-diameter carbon nanotube are focus of first principles calculations in this work. A high performance parallel computing platform is also developed to perform these calculations at Michigan Tech. First principles calculations are based on density functional theory employing either local density or gradient-corrected approximation together with plane wave and gaussian basis sets. The bulk Ga2O3 is known to be a very good candidate for fabricating electronic devices that operate at high temperatures. To explore the properties of Ga2O3 at nonoscale, we have performed a systematic theoretical study on the small polyatomic gallium oxide clusters. The calculated results find that all lowest energy isomers of GamOn clusters are dominated by the Ga-O bonds over the metal-metal or the oxygen-oxygen bonds. Analysis of atomic charges suggest the clusters to be highly ionic similar to the case of bulk Ga2O3. In the study of sequential oxidation of these slusters starting from Ga2O, it is found that the most stable isomers display up to four different backbones of constituent atoms. Furthermore, the predicted configuration of the ground state of Ga2O is recently confirmed by the experimental result of Neumark's group. Guided by the results of calculations the study of gallium oxide clusters, performance related challenge of computational simulations, of producing high performance computers/platforms, has been addressed. Several engineering aspects were thoroughly studied during the design, development and implementation of the high performance parallel computing platform, rama, at Michigan Tech. In an attempt to stay true to the principles of Beowulf revolutioni, the rama cluster was extensively customized to make it easy to understand, and use - for administrators as well as end-users. Following the results of benchmark calculations and to keep up with the complexity of systems under study, rama has been expanded to a total of sixty four processors. Interest in the non-covalent intereaction of DNA with carbon nanotubes has steadily increased during past several years. This hybrid system, at the junction of the biological regime and the nanomaterials world, possesses features which make it very attractive for a wide range of applicatioins. Using the in-house computational power available, we have studied details of the interaction between nucleic acid bases with graphene sheet as well as high-curvature small-diameter carbon nanotube. The calculated trend in the binding energies strongly suggests that the polarizability of the base molecules determines the interaction strength of the nucleic acid bases with graphene. When comparing the results obtained here for physisorption on the small diameter nanotube considered with those from the study on graphene, it is observed that the interaction strength of nucleic acid bases is smaller for the tube. Thus, these results show that the effect of introducing curvature is to reduce the binding energy. The binding energies for the two extreme cases of negligible curvature (i.e. flat graphene sheet) and of very high curvature (i.e. small diameter nanotube) may be considered as upper and lower bounds. This finding represents an important step towards a better understanding of experimentally observed sequence-dependent interaction of DNA with Carbon nanotubes.

Single-Molecule Conductance of Functionalized Oligoynes: Length Dependence and Junction Evolution

We report a combined experimental and theoretical investigation of the length dependence and anchor group dependence of the electrical conductance of a series of oligoyne molecular wires in single-molecule junctions with gold contacts. Experimentally, we focus on the synthesis and properties of diaryloligoynes with n = 1, 2, and 4 triple bonds and the anchor dihydrobenzo[b]thiophene (BT). For comparison, we also explored the aurophilic anchor group cyano (CN), amino (NH2), thiol (SH), and 4-pyridyl (PY). Scanning tunneling microscopy break junction (STM-BJ) and mechanically controllable break junction (MCBJ) techniques are employed to investigate single-molecule conductance characteristics. The BT moiety is superior as compared to traditional anchoring groups investigated so far. BT-terminated oligoynes display a 100% probability of junction formation and possess conductance values which are the highest of the oligoynes studied and, moreover, are higher than other conjugated molecular wires of similar length. Density functional theory (DFT)-based calculations are reported for oligoynes with n = 1−4 triple bonds. Complete conductance traces and conductance distributions are computed for each family of molecules. The sliding of the anchor groups leads to oscillations in both the electrical conductance and the binding energies of the studied molecular wires. In agreement with experimental results, BT-terminated oligoynes are predicted to have a high electrical conductance. The experimental attenuation constants βH range between 1.7 nm−1 (CN) and 3.2 nm−1 (SH) and show the following trend: βH(CN) < βH(NH2) < βH(BT) < βH(PY) ≈ βH(SH). DFT-based calculations yield lower values, which range between 0.4 nm−1 (CN) and 2.2 nm−1 (PY).