Molecular structure does not easily identify the intricate non-covalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect non-covalent interactions in real space, based on the electron density and its derivatives. Our approach reveals underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of non-bonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.

Non-covalent interactions hold the key to understanding many chemical, biological, and technological problems. Describing these non-covalent interactions accurately, including their positions in real space, constitutes a first step in the process of decoupling the complex balance of forces that define non-covalent interactions. Because of the size of macromolecules, the most common approach has been to assign van der Waals interactions (vdW), steric clashes (SC), and hydrogen bonds (HBs) based on pairwise distances between atoms according to their van der Waals radii. We recently developed an alternative perspective, derived from the electronic density: the Non-Covalent Interactions (NCI) index [J. Am. Chem. Soc. 2010, 132, 6498]. This index has the dual advantages of being generally transferable to diverse chemical applications and being very fast to compute, since it can be calculated from promolecular densities. Thus, NCI analysis is applicable to large systems, including proteins and DNA, where analysis of non-covalent interactions is of great potential value. Here, we describe the NCI computational algorithms and their implementation for the analysis and visualization of weak interactions, using both selfconsistent fully quantum-mechanical, as well as promolecular, densities. A wide range of options for tuning the range of interactions to be plotted is also presented. To demonstrate the capabilities of our approach, several examples are given from organic, inorganic, solid state, and macromolecular chemistry, including cases where NCI analysis gives insight into unconventional chemical bonding. The NCI code and its manual are available for download at

The effects of structural fluctuations on charge transfer in double-stranded DNA and peptide nucleic acid (PNA) are investigated. A palindromic sequence with two guanine bases that play the roles of hole donor and acceptor, separated by a bridge of two adenine bases, was analyzed using combined molecular dynamics (MD) and quantum-chemical methods. Surprisingly, electronic structure calculations on individual MD snapshots show significant frontier orbital electronic population on the bridge in approximately 10% of the structures. Electron-density delocalization to the bridge is found to be gated by fluctuations of the covalent conjugated bond structure of the aromatic rings of the nucleic bases. It is concluded, therefore, that both thermal hopping and superexchange should contribute significantly to charge transfer even in short DNA/PNA fragments. PNA is found to be more flexible than DNA, and this flexibility is predicted to produce larger rates of charge transfer.

ConspectusElectron transfer (ET) reactions provide a nexus among chemistry, biochemistry, and physics. These reactions underpin the "power plants" and "power grids" of bioenergetics, and they challenge us to understand how evolution manipulates structure to control ET kinetics. Ball-and-stick models for the machinery of electron transfer, however, fail to capture the rich electronic and nuclear dynamics of ET molecules: these static representations disguise, for example, the range of thermally accessible molecular conformations. The influence of structural fluctuations on electron-transfer kinetics is amplified by the exponential decay of electron tunneling probabilities with distance, as well as the delicate interference among coupling pathways. Fluctuations in the surrounding medium can also switch transport between coherent and incoherent ET mechanisms-and may gate ET so that its kinetics is limited by conformational interconversion times, rather than by the intrinsic ET time scale. Moreover, preparation of a charge-polarized donor state, or of a donor state with linear or angular momentum, can have profound dynamical and kinetic consequences. In this Account, we establish a vocabulary to describe how the conformational ensemble and the prepared donor state influence ET kinetics in macromolecules. This framework is helping to unravel the richness of functional biological ET pathways, which have evolved in within fluctuating macromolecular structures.The conceptual framework for describing nonadiabatic ET seems disarmingly simple: compute the ensemble averaged (mean-squared) donor-acceptor (DA) tunneling interaction, , and the Franck-Condon weighted density of states, ρ FC , to describe the rate: (2π/ħ) ρ FC . Modern descriptions of the thermally averaged electronic coupling and of the Franck-Condon factor establish a useful predictive framework in biology, chemistry, and nanoscience. Describing the influence of geometric and energetic fluctuations on ET allows us to address a rich array of mechanistic and kinetic puzzles. How strongly is a protein's fold imprinted on the ET kinetics, and might thermal fluctuations "wash out" signatures of structure? What is the influence of thermal fluctuations on ET kinetics beyond averaging of the tunneling barrier structure? Do electronic coupling mechanisms change as donor and acceptor reposition in a protein, and what are the consequences for the ET kinetics? Do fluctuations access minority species that dominate tunneling? Can energy exchanges between the electron and bridge vibrations generate vibronic signatures that label some of the D-to-A pathway traversed by the electron, thus eliminating unmarked pathways that would otherwise contribute to the DA coupling (as in other "which way" or double-slit experiments)? Might medium fluctuations drive tunneling-hopping mechanistic transitions? How does the donor-state preparation -in particular its polarization toward the acceptor and its momentum characteristics (which may introduce complex rather than pur...

In molecular electronics, the linker group, which attaches the functional molecular core to the electrode, plays a crucial role in determining the overall conductivity of the molecular junction. While much focus has been placed on optimizing molecular core conductivity, there have been relatively few attempts at designing optimal linker groups to metallic or semiconducting electrodes. The vast majority of molecular electronic studies use thiol linker groups; work probing alternative amine linker systems has only recently been explored. Here, we probe single-molecule conductances in phenylene-ethynylene molecules terminated with thiol and carbodithioate linkers, experimentally using STM break-junction methods and theoretically using a nonequilibrium Green's function approach. Experimental studies demonstrate that the carbodithioate linker augments electronic coupling to the metal electrode and lowers the effective barrier for charge transport relative to the conventional thiol linker, thus enhancing the conductance of the linker-phenylene-ethynylene-linker unit; these data underscore that phenylene-ethynylene-based structures are more highly conductive than originally appreciated in molecular electronics applications. The theoretical analysis shows that the nature of sulfur hybridization in these species is responsible for the order-of-magnitude increased conductance in carbodithioate-terminated systems relative to identical conjugated structures that feature classic thiol linkers, independent of the mechanism of charge transport. Interestingly, in these systems, the tunneling current is not dominated by the frontier molecular orbitals. While barriers >k(B)T to produce the low beta values seen in our experiments. Taken together, these experimental and theoretical studies indicate a promising role for carbodithioate-based connectivity in molecular-scale electronics applications involving metallic and semiconducting electrodes.

Theoretical studies of charge transport in deoxyribonucleic acid (DNA) and peptide nucleic acid (PNA) indicate that structure and dynamics modulate the charge transfer rates, and that different members of a structural ensemble support different charge transport mechanisms. Here, we review the influences of nucleobase geometry, electronic structure, solvent environment, and thermal conformational fluctuations on the charge transfer mechanism. We describe an emerging framework for understanding the diversity of charge transport mechanisms seen in nucleic acids.

Self-assembled monolayers of single-stranded (ss) peptide nucleic acids (PNAs) containing seven nucleotides (TTTXTTT), a C-terminus cysteine, and an N-terminus ferrocene redox group were formed on gold electrodes. The PNA monomer group (X) was selected to be either cytosine (C), thymine (T), adenine (A), guanine (G), or a methyl group (Bk). The charge transfer rate through the oligonucleotides was found to correlate with the oxidation potential of X. Kinetic measurements and computational studies of the ss-PNA fragments show that a nucleobase mediated charge transport mechanism in the deep tunneling superexchange regime can explain the observed dependence of the kinetics of charge transfer on the PNA sequence. Theoretical analysis suggests that the charge transport is dominantly hole-mediated and takes place through the filled bridge orbitals. The strongest contribution to conductance comes from the highest filled orbitals (HOMO, HOMO-1, and HOMO-2) of individual bases, with a rapid drop off in contributions from lower lying filled orbitals. Our studies further suggest that the linear correlation observed between the experimental charge transfer rates and the oxidation potential of base X arises from weak average interbase couplings and similar stacking geometries for the four TTTXTTT systems.

Nonlinear optical chromophores containing (porphyrinato)Zn(II), proquinoid, and (terpyridyl)metal(II) building blocks were optimized in a library containing approximately 10(6) structures using the linear combination of atomic potentials (LCAP) methodology. We report here the library design and molecular property optimizations. Two basic structural types of large beta(0) chromophores were examined: linear and T-shaped motifs. These T-shaped geometries suggest a promising NLO chromophoric architecture for experimental investigation and further support the value of performing LCAP searches in large chemical spaces.