Controlling Single Molecule Conductance by a Locally Induced Chemical Reaction on Individual Thiophene Units

Among the prerequisites for the progress of single‐molecule‐based electronic devices are a better understanding of the electronic properties at the individual molecular level and the development of methods to tune the charge transport through molecular junctions. Scanning tunneling microscopy (STM) is an ideal tool not only for the characterization, but also for the manipulation of single atoms and molecules on surfaces. The conductance through a single molecule can be measured by contacting the molecule with atomic precision and forming a molecular bridge between the metallic STM tip electrode and the metallic surface electrode. The parameters affecting the conductance are mainly related to their electronic structure and to the coupling to the metallic electrodes. Here, the experimental and theoretical analyses are focused on single tetracenothiophene molecules and demonstrate that an in situ‐induced direct desulfurization reaction of the thiophene moiety strongly improves the molecular anchoring by forming covalent bonds between molecular carbon and copper surface atoms. This bond formation leads to an increase of the conductance by about 50 % compared to the initial state.     

A self-consistent transport model for molecular conduction based on extended Hückel theory with full three-dimensional electrostatics

We present a transport model for molecular conduction involving an extended Hückel theoretical treatment of the molecular chemistry combined with a nonequilibrium Green's function treatment of quantum transport. The self-consistent potential is approximated by CNDO (complete neglect of differential overlap) method and the electrostatic effects of metallic leads (bias and image charges) are included through a three-dimensional finite element method. This allows us to capture spatial details of the electrostatic potential profile, including effects of charging, screening, and complicated electrode configurations employing only a single adjustable parameter to locate the Fermi energy. As this model is based on semiempirical methods it is computationally inexpensive and flexible compared to ab initio models, yet at the same time it is able to capture salient qualitative features as well as several relevant quantitative details of transport. We apply our model to investigate recent experimental data on alkane dithiol molecules obtained in a nanopore setup. We also present a comparison study of single molecule transistors and identify electronic properties that control their performance.     

Water at surfaces and interfaces: From molecules to ice and bulk liquid

The structure and growth of water films on surfaces is reviewed, starting from single molecules to two-dimensional wetting layers, and liquid interfaces. This progression follows the increase in temperature and vapor pressure from a few degrees Kelvin in ultra-high vacuum, where Scanning Tunneling and Atomic Force Microscopies (STM and AFM) provide crystallographic information at the molecular level, to ambient conditions where surface sensitive spectroscopic techniques provide electronic structure information. We show how single molecules bind to metal and non-metal surfaces, their diffusion and aggregation. We examine how water molecules can be manipulated by the STM tip via excitation of vibrational and electronic modes, which trigger molecular diffusion and dissociation. We review also the adsorption and structure of water on non-metal substrates including mica, alkali halides, and others under ambient humid conditions. We finally discuss recent progress in the exploration of the molecular level structure of solid-liquid interfaces, which impact our fundamental understanding of corrosion and electrochemical processes.     

Chirality- and diameter-dependent reactivity of NO2 on carbon nanotube walls

We report the density-functional calculations of NO2 adsorption on single-walled carbon nanotube walls. A single molecular adsorption was endothermic with an activation barrier, but a collective adsorption with several molecules became exothermic without an activation barrier. We find that NO2 adsorption is strongly electronic structure- and strain-dependent. The NO2 adsorption on metallic nanotubes was energetically more favorable than that on semiconducting nanotubes and furthermore the adsorption became less stable with increasing diameters of nanotubes. The adsorption barrier height shows similar dependence on the electronic structure and diameter to the adsorption energy. Our theoretical model can be a good guideline for the separation of nanotubes by electronic structures using various adsorbates.     

Using molecular similarity to construct accurate semiempirical electronic structure theories

Ab initio electronic structure methods give accurate results for small systems, but do not scale well to large systems. Chemical insight tells us that molecular functional groups will behave approximately the same way in all molecules, large or small. This molecular similarity is exploited in semiempirical methods, which couple simple electronic structure theories with parameters for the transferable characteristics of functional groups. We propose that high-level calculations on small molecules provide a rich source of parametrization data. In principle, we can select a functional group, generate a large amount of ab initio data on the group in various small-molecule environments, and "mine" this data to build a sophisticated model for the group's behavior in large environments. This work details such a model for electron correlation: a semiempirical, subsystem-based correlation functional that predicts a subsystem's two-electron density matrix as a functional of its one-electron density matrix. This model is demonstrated on two small systems: chains of linear, minimal-basis (H-H)(5), treated as a sum of four overlapping (H-H)(2) subsystems; and the aldehyde group of a set of HOC-R molecules. The results provide an initial demonstration of the feasibility of the approach.     

Structure and dynamics of the solvation of bovine pancreatic trypsin inhibitor in explicit water: a comparative study of the effects of solvent and protein polarizability

To isolate the effects of the inclusion of polarizability in the force field model on the structure and dynamics of the solvating water in differing electrostatic environments of proteins, we present the results of molecular dynamics simulations of the bovine pancreatic trypsin inhibitor (BPTI) in water with force fields that explicitly include polarization for both the protein and the water. We use three model potentials for water and two model potentials for the protein. Two of the water models and one of the protein models are polarizable. A total of six systems were simulated representing all combinations of these polarizable and nonpolarizable protein and water force fields. We find that all six systems behave in a similar manner in regions of the protein that are weakly electrostatic (either hydrophobic or weakly hydrophilic). However, in the vicinity of regions of the protein with relatively strong electrostatic fields (near positively or negatively charged residues), we observe that the water structure and dynamics are dependent on both the model of the protein and the model of the water. We find that a large part of the dynamical dependence can be described by small changes in the local environments of each region that limit the local density of non-hydrogen-bonded waters, precisely the water molecules that facilitate the dynamical relaxation of the water-water hydrogen bonds. We introduce a simple method for rescaling for this effect. When this is done, we are able to effectively isolate the influence of polarizability on the dynamics. We find that the solvating water's relaxation is most affected when both the protein and the water models are polarizable. However, when only one model (or neither) is polarizable, the relaxation is similar regardless of the models used.     

Dynamics on the electronically excited state surface of the bioluminescent firefly luciferase-oxyluciferin system

Dynamics of the firefly luciferase-oxyluciferin complex in its electronic ground and excited states are studied using various theoretical approaches. By mimicking the physiological conditions with realistic models of the chromophore oxyluciferin, the enzyme luciferase, and solvating water molecules and by performing real time simulations with a molecular dynamics technique on the model surfaces, we reveal that the local chromophore-surrounding interaction patterns differ rather severely in the two states. Because of the presence of protein, the solvation dynamics of water around the chromophore is also peculiar and shows widely different time scales on the two terminal oxygen atoms. In addition, simulations of the emission with the quantum-mechanics/molecular-mechanics approach show a close relationship between the emission color variation and the environmental dynamics, mostly through electrostatic effects from the chromophore-surrounding interaction. We also discuss the importance of considering the time scales of the luminescence and the dynamics of the interaction.     

Structure and bonding of ethoxy species adsorbed on transition metal surfaces

The interaction of the ethoxy radical with Cu(111), Ag(111), Pd(111) and Au(111) has been studied using a periodic density functional approach. The most stable adsorption site is the fcc with adsorption energies in the 1.1–2.2 eV interval. All analyses consistently indicate that ethoxy becomes negatively charged, that the presence of ethoxy slightly perturbs the electronic structure of the metallic surface, that the interaction is essentially of electrostatic character and not directional predicting a rather mobile species. The calculated adsorption energies are found to correlate almost linearly with the total net charge on the ethoxy moiety thus confirming that the electrostatic interactions dominate the bonding between this organic species and the underlying metallic surfaces.     

Adsorption of water on sodium chloride surfaces: electrostatics – guided ab initio studies

Water adsorption is studied on medium-sized clusters of sodium chloride representing (100) and (110) surfaces at the ab initio level. Topographical features of molecular electrostatic potential (MESP) have been employed for predicting the potent sites for binding of one to four water molecules on these surfaces. Such guess geometries are initially optimized using an electrostatics-based model, electrostatic potential for intermolecular complexation (EPIC) and further at the Hartree–Fock and B3LYP/6-31G(d, p) levels. The corresponding interaction energies are examined for assessing co-operative binding effects. The geometry and interaction energy of four water molecules adsorbed on NaCl(100) clearly brings out the co-operative binding among the water molecules. Further, water binding to (110) surface is stronger than that with (100) surface. This is also in confirmation with the electrostatic properties of (110) surface. Many-body decomposition analysis brings out the stronger interaction between NaCl clusters with water molecules vis-a-vis water–water interaction.     

Water-mediated electron transfer between protein redox centers

Recent experimental and theoretical investigations show that water molecules between or near redox partners can significantly affect their electron-transfer (ET) properties. Here we study the effects of intervening water molecules on the electron self-exchange reaction of azurin (Az), by performing a conformational sampling on the water medium and by using a newly developed ab initio method to calculate transfer integrals between molecular redox sites. We show that the insertion of water molecules at the interface between the copper active sites of Az dimers slightly increases the overall ET rate, while some favorable water conformations can considerably enhance the ET kinetics. These features are traced back to the interplay of two competing factors: the electrostatic interaction between the water and protein subsystems (mainly opposing the ET process for the water arrangements drawn from MD simulations) and the effectiveness of water in mediating ET coupling pathways. Such an interplay provides a physical basis for the found absence of correlation between the electronic couplings derived through ab initio electronic structure calculations and the related quantities obtained through the Empirical Pathways (EP) method. In fact, the latter does not account for electrostatic effects on the transfer integrals. Thus, we conclude that the water-mediated electron tunneling is not controlled by the geometry of a single physical pathway. We discuss the results in terms of the interplay between different ET pathways controlled by the conformational changes of one of the water molecules via its electrostatic influence. Finally, we examine the dynamical effects of the interfacial water and check the validity of the Condon approximation.     

Electrostatic field-adapted molecular fractionation with conjugated caps for energy calculations of charged biomolecules

An electrostatic field-adapted molecular fractionation with conjugated caps (EFA-MFCC) approach is implemented for treating macromolecules with several charge centers. The molecular fragmentation is performed in an "electrostatic field," which is described by putting point charges on charge centers, directly affecting the Hamiltonians of both fragments and conjugated caps. So the present method does not need truncation during the calculation of electrostatic interactions. Our test calculations on a series of charged model systems and biological macromolecules using the HF and B3LYP methods have demonstrated that this approach is capable of describing the electronic structure with accuracy comparable to other fragment-based methods. The EFA-MFCC approach is an alternative way for predicting the total energies of charged macromolecules with acyclic, loop, and intersectional loop structures and interaction energies between two molecules.     

Switching in molecular transport junctions: polarization response

We discuss several proposed explanations for the switching and negative differential resistance (NDR) behavior seen in some molecular junctions. Several theoretical models are discussed, and we present results of electronic structure calculations on a series of substituted oligo(phenylene ethynylene) molecules. It is shown that a previously proposed polaron model is successful in predicting NDR behavior, and the model is elaborated with image charge effects and parameters from electronic structure calculations. This model now incorporates substituent effects and includes the effects of conformational change, charging, and image charge stabilization.     

Continuum treatment of electronic polarization effect

A continuum treatment of electronic polarization has been explored for in molecular mechanics simulations in implicit solvents. The dielectric constant for molecule interior is the only parameter in the continuum polarizable model. A value of 4 is found to yield optimal agreement with high-level ab initio quantum mechanical calculations for the tested molecular systems. Interestingly, its performance is not sensitive to the definition of molecular volume, in which the continuum electronic polarization is defined. In this model, quantum mechanical electrostatic field in different dielectric environments from vacuum, low-dielectric organic solvent, and water can be used simultaneously in atomic charge fitting to achieve consistent treatment of electrostatic interactions. The tests show that a single set of atomic charges can be used consistently in different dielectric environments and different molecular conformations, and the atomic charges transfer well from training monomers to tested dimers. The preliminary study gives us the hope of developing a continuum polarizable force field for more consistent simulations of proteins and nucleic acids in implicit solvents.     

Triisopropylsilylethynyl‐Functionalized Graphene‐Like Fragment Semiconductors: Synthesis,Crystal Packing,and Density Functional Theory Calculations

Tri‐isopropylsilylethynyl (TIPS)‐functionalized polycyclic aromatic hydrocarbon (PAH) molecules incorporate structural components of graphene nanoribbons and represent a family of model molecules that form organic crystal semiconductors for electronic devices. Here, we report a series of TIPS‐functionalized PAHs and discuss their electronic properties and crystal packing features. We observe that these soluble compounds easily form one‐dimensional (1 D) packing arrangements and allow a direct evolution of the π stacking by varying the geometric shape. We find that the aspect ratio between length and width plays an important role on crystal packing. Our result indicates that when the PAH molecules have zigzag edges, these can provide enough volume for the molecules to rotate and reorient, alleviating the unfavorable electrostatic interactions found in perfectly cofacial π–π stacking. Density functional theory calculations were carried out to provide insights into how the molecular geometric shape influences the electronic structure and transport properties. The calculations indicate that, among the compounds studied here, “brick‐layer” stacks provide the highest hole mobility.     

Electronic g-tensors of solvated molecules using the polarizable continuum model

We present the implementation of density functional response theory combined with the polarizable continuum model (PCM), enabling first principles calculations of molecular g-tensors of solvated molecules. The calculated g-tensor shifts are compared with experimental g-tensor shifts obtained from electron paramagnetic resonance spectra for a few solvated species. The results indicate qualitative agreement between the calculations and the experimental data for aprotic solvents, whereas PCM fails to reproduce the electronic g-tensor behavior for protic solvents. This failure of PCM for protic solvents can be resolved by including into the model those solvent molecules which are involved in hydrogen bonding with the solute. The results for the protic solvents show that the explicit inclusion of the solvent molecules of the first solvation sphere is not sufficient in order to reproduce the behavior of the electronic g-tensor in protic solvents, and that better agreement with experimental data can be obtained by including the long-range electrostatic effects accounted for by the PCM approach on top of the explicit hydrogen-bonded complexes.     

Modeling intermolecular interactions of physisorbed organic molecules using pair potential calculations

The understanding and control of epitaxial growth of organic thin films is of crucial importance in order to optimize the performance of future electronic devices. In particular, the start of the submonolayer growth plays an important role since it often determines the structure of the first layer and subsequently of the entire molecular film. We have investigated the structure formation of 3,4,9,10-perylene-tetracarboxylic dianhydride and copper-phthalocyanine molecules on Au(111) using pair-potential calculations based on van der Waals and electrostatic intermolecular interactions. The results are compared with the fundamental lateral structures known from experiment and an excellent agreement was found for these weakly interacting systems. Furthermore, the calculations are even suitable for chemisorptive adsorption as demonstrated for copper-phthalocyanine/Cu(111), if the influence of charge transfer between substrate and molecules is known and the corresponding charge redistribution in the molecules can be estimated. The calculations are of general applicability for molecular adsorbate systems which are dominated by electrostatic and van der Waals interaction.     

Partial atomic multipoles for internally consistent microelectrostatic calculations

An extension of the extant microelectrostatic methodologies, based on the concept of distributed generalized polarizability matrix derived from the Coupled Perturbed Hartree–Fock (CPHF) equations, is proposed for self‐consistent calculation of charge carrier and charge‐transfer (CT) state electrostatic energies in molecular solids, including the doped, defected and disordered ones. The CPHF equations are solved only once and the generalized molecular polarizability they yield enables low cost iterations that mutually adjust the molecular electronic distributions and the local electric field in which the molecules are immersed. The approach offers a precise picture of molecular charge densities, accounting for atomic partial multipoles up to order 2, which allows one to reproduce the recently reported large charge‐quadrupole contributions to CT state energies in low‐symmetry local environments. It is particularly well suited for repetitive calculations for large clusters (up to 300,000 atoms), and may potentially be useful for describing electrostatic solvent effects. © 2017 Wiley Periodicals, Inc.