Local Structure of Hydrogen-Bonded Liquids

Ordinary yet unique, water is the substance on which life is based. Water seems, at first sight, to be a very simple molecule, consisting of two hydrogen atoms attached to one oxygen. Its small size belies the complexity of its action and its numerous anomalies, central to a broad class of important phenomena, ranging from global current circulation, terrestrial water and CO2 cycles to corrosion and wetting. The explanation of this complex behavior comes from water's unique ability to form extensive three-dimensional networks of hydrogen-bonds, whose nature and structures, in spite of a great deal of efforts involving a plethora of experimental and theoretical techniques, still lacks a complete scientific understanding. This thesis is devoted to the study of the local structure of hydrogen-bonded liquids, with a particular emphasis on water, taking advantage of a combination of core-level spectroscopies and density functional theory spectra calculations. X-ray absorption, in particular, is found to be sensitive to the local hydrogen-bond environment, thus offering a very promising tool for spectroscopic identification of specific structural configurations in water, alcohols and aqueous solutions. More specifically, the characteristic spectroscopic signature of the broken hydrogen-bond at the hydrogen side is used to analyze the structure of bulk water, leading to the finding that most molecules are arranged in two hydrogen-bond configurations, in contrast to the picture provided by molecular dynamics simulations. At the liquid-vapor interface, an interplay of surface sensitive measurements and theoretical calculations enables us to distinguish a new interfacial species in equilibrium with the gas. In a similar approach the cluster form of the excess proton in highly concentrated acid solutions and the different coordination of methanol at the vacuum interface and in the bulk can also be clearly identified. Finally the ability of core-level spectroscopies, aided by sophisticated density functional theory calculations, to directly probe the valence electronic structure of a system is used to observe the nature of the interaction between water molecules and solvated ions in solution. Water around transition metal ions is found to interact with the solute via orbital mixing with the metal d-orbitals. The hydrogen-bond between water molecules is explained in terms of electrostatic interactions enhanced by charge rehybridization in which charge transfer between connecting molecules is shown to be fundamental.

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.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Theoretical insights on electron donor-acceptor interactions involving carbon dioxide

Electron donor-acceptor (EDA) interactions are widely involved in chemistry and their understanding is essential to design new technological applications in a variety of fields ranging from material sciences and chemical engineering to medicine. In this work, we study EDA complexes of carbon dioxide with ketones using several ab initio and Density Functional Theory methods. Energy contributions to the interaction energy have been analyzed in detail using both variational and perturbational treatments. Dispersion energy has been shown to play a key role in explaining the high stability of a non-conventional structure, which can roughly be described by a cooperative EDA interaction.

Polarized–unpolarized ground state of small polycyclic aromatic hydrocarbons

Do polyacenes, circumacenes, periacenes, nanographenes, and graphene nanoribbons show a spin polarized ground state? In this work, we present monodeterminantal (Hartree–Fock (HF) and density functional theory (DFT) types), and multideterminantal calculations (Møller–Plesset and Coupled Cluster), for several families of unsaturated organic molecules (n-Acenes, n-Periacenes and n-Circumacenes). All HF calculations and many DFT show a spin-polarized (antiferromagnetic) ground state, in agreement with previous calculations. Nevertheless, the multideterminantal calculations, carried out with perturbative and variational wavefunctions, show that the more stable state is obtained starting from the unpolarized HF wavefunction. The trend of the stabilization of wavefunctions (polarized or unpolarized) with respect to exchange and correlation potentials, and to the number of benzene rings, has been analyzed. A study of the spin (〈Ŝ2〉) and the spin density on the carbon atoms has also been carried out.

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. All rights reserved.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

Density-functional description of superconducting and magnetic proximity effects across a tunneling barrier

A density-functional formalism for superconductivity and magnetism is presented. The resulting relations unify previously derived Kohn-Sham equations for superconductors and for noncollinear magnetism. The formalism, which discriminates Cooper-pair singlets from triplets, is applied to two quantum liquids coupled by tunneling through a barrier. An exact expression is derived, relating the eigenstates and eigenvalues of the Kohn-Sham equations, unperturbed by tunneling, on one side of the barrier to the proximity-induced ordering potential on the other.

A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

Comparative study of the adsorption and dissociation of vinylacetic acid and acrylic acid on silicon (001)

In this work, we employ the state of the art pseudopotential method, within a generalized gradient approximation to the density functional theory, to investigate the adsorption process of acrylic acid (AAc) and vinylacetic acid (VAA) on the silicon surface. Our total energy calculations support the proposed experimental process, as it indicates that the chemisorption of the molecule is as follows: The gas phase VAA (AAc) adsorbs molecularly to the electrophilic surface Si atom and then dissociates into H(2)C = CH - COO and H, bonded to the electrophilic and nucleophilic surface silicon dimer atoms, respectively. The activation energy for both processes correspond to thermal activations that are smaller than the usual growth temperature. In addition, the electronic structure, calculated vibrational modes, and theoretical scanning tunneling microscopy images are discussed, with a view to contribute to further experimental investigations.

First-principles study of the adsorption of atomic and molecular hydrogen on BC(2)N nanotubes

The adsorption of atomic and molecular hydrogen on armchair and zigzag boron carbonitride nanotubes is investigated within the ab initio density functional theory. The adsorption of atomic H on the BC(2)N nanotubes presents properties which are promising for nanoelectronic applications. Depending on the adsorption site for the H, the Fermi energy moves toward the bottom of the conduction band or toward the top of the valence band, leading the system to exhibit donor or acceptor characteristics, respectively. The H(2) molecules are physisorbed on the BC(2)N surface for both chiralities. The binding energies for the H(2) molecules are slightly dependent on the adsorption site, and they are near to the range to work as a hydrogen storage medium.

Control of the persistent currents in two interacting quantum rings through the Coulomb interaction and interring tunneling

The persistent current in two vertically coupled quantum rings containing few electrons is studied. We find that the Coulomb interaction between the rings in the absence of tunneling affects the persistent current in each ring and the ground-state configurations. Quantum tunneling between the rings alters significantly the ground state and the persistent current in the system.

Ground-state functional for quantum spin chains with bond defects

We have developed a nonlocal functional of the exchange interaction for the ground-state energy of quantum spin chains described by the Heisenberg Hamiltonian. An alternating chain is used to obtain the correlation energy and a local unit-cell approximation is defined in the context of the density-functional theory. The agreement with our exact numerical data, for small chains, is significantly better than a previous formulation, even for chains with several ferromagnetic or antiferromagnetic bond defects. The results can be particularly relevant in the study of finite spin-1/2 Heisenberg chains, with exchange couplings changing, magnitude, or even sign, from bond-to-bond.