Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Electron transfer rate modulation with mid-IR in butadiyne-bridged donor-bridge-acceptor compounds

Controlling electron transfer (ET) processes in donor-bridge-acceptor (DBA) compounds by mid-IR excitation can enhance our understanding of the ET dynamics and may find practical applications in molecular sensing and molecular-scale electronics. Alkyne moieties are attractive to serve as ET bridges, as they offer the possibility of fast ET and present convenient vibrational modes to perturb the ET dynamics. Yet, these bridges introduce complexity because of the strong torsion angle dependence of the ET rates and transition dipoles among electronic states and a shallow torsion barrier. In this study, we implemented ultrafast 3-pulse laser spectroscopy to investigate how the ET from the dimethyl aniline (D) electron donor to the N-isopropyl-1,8-napthalimide (NAP) electron acceptor can be altered by exciting the CC stretching mode (νCC) of the butadiyne bridge linking the donor and acceptor. The electron transfer was initiated by electronically exciting the acceptor moiety at 400 nm, followed by vibrational excitation of the alkyne, νCC, and detecting the changes in the absorption spectrum in the visible spectral region. The experiments were performed at different delay times t1 and t2, which are the delays between UV-mid-IR and mid-IR-Vis pulses, respectively. Two sets of torsion-angle conformers were identified, one featuring a very fast mean ET time of 0.63 ps (group A) and another featuring a slower mean ET time of 4.3 ps (group B), in the absence of the mid-IR excitation. TD-DFT calculations were performed to determine key torsion angle dependent molecular parameters, including the electronic and vibrational transition dipoles, transition frequencies, and electronic couplings. To describe the 3-pulse data, we developed a kinetic model that includes a locally excited, acceptor-based S2 state, a charge separated S1 state, and their vibrationally excited counterparts, with either excited νCC (denoted as S1Atr, S1Btr, S2Atr, and S2Btr, where tr stands for the excited triplet bond, νCC) or excited daughter modes of the νCC relaxation (S1Ah, S1Bh, S2Ah, and S2Bh, where h stands for vibrationally hot species). The kinetic model was solved analytically, and the species-associated spectra (SAS) were determined numerically using a matrix approach, treating first the experiments with longer t1 delays and then using the already determined SAS for modeling the experiments with shorter t1 delays. Strong vibronic coupling of νCC and of vibrationally hot states makes the analysis complicated. Nevertheless, the SAS were identified and the ET rates of the vibrationally excited species, S2Atr, S2Btr and S2Bh, were determined. The results show that the ET rate for the S2A species is ca. 1.2-fold slower when the νCC mode is excited. The ET rate for species S2B is slower by ca. 1.3-fold if the compound is vibrationally hot and is essentially unchanged when the νCC mode is excited. The SAS determined for the tr and h species resemble the SAS for their respective precursor species in the 2-pulse transient absorption experiments, which validates the procedure used and the results.

Electron-donating functional groups strengthen ligand-induced chiral imprinting on CsPbBr3 quantum dots

Chiral perovskite nanoparticles and films are promising for integration in emerging spintronic and optoelectronic technologies, yet few design rules exist to guide the development of chiral material properties. The chemical space of potential building blocks for these nanostructures is vast, and the mechanisms through which organic ligands can impart chirality to the inorganic perovskite lattice are not well understood. In this work, we investigate how the properties of chiral ammonium ligands, the most common organic ligand type used with perovskites, affect the circular dichroism of strongly quantum confined CsPbBr3 nanocrystals. We show that aromatic ammonium ligands with stronger electron-donating groups lead to higher-intensity circular dichroism associated with the lowest-energy excitonic transition of the perovskite nanocrystal. We argue that this behavior is best explained by a modulation of the exciton wavefunction overlap between the nanocrystal and the organic ligand, as the functional groups on the ligand can shift electron density toward the organic species-perovskite lattice interface to increase the imprinting.

Quantum Simulation of Polarized Light-Induced Electron Transfer with a Trapped-Ion Qutrit System

Electron transfer within and between molecules is crucial in chemistry, biochemistry, and energy science. This study describes a quantum simulation method that explores the influence of light polarization on electron transfer between two molecules. By implementing precise and coherent control among the quantum states of trapped atomic ions, we can induce quantum dynamics that mimic the electron-transfer dynamics in molecules. We use three-level systems (qutrits), rather than traditional two-level systems (qubits), to enhance the simulation efficiency and realize high-fidelity simulations of electron-transfer dynamics. We treat the quantum interference between the electron coupling pathways from a donor with two degenerate excited states to an acceptor and analyze the transfer efficiency. We also examine the potential error sources that enter the quantum simulations. The trapped-ion systems have favorable scalings with system size compared to those of classical computers, promising access to richer electron-transfer simulations.

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. We show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

Adenine Fine-Tunes DNA Photolyase's Repair Mechanism

The comparative study of DNA repair by mesophilic and extremophilic photolyases helps us understand the evolution of these enzymes and their role in preserving life on our changing planet. The mechanism of repair of cyclobutane pyrimidine dimer lesions in DNA by electron transfer from the flavin adenine dinucleotide cofactor is the subject of intense interest. The role of adenine in mediating this process remains unresolved. Using microsecond molecular dynamics simulations, we find that adenine mediates the electron transfer in both mesophile and extremophile DNA photolyases through a similar mechanism. In fact, in all photolyases studied, the molecular conformations with the largest electronic couplings between the enzyme cofactor and DNA show the presence of adenine in 10-20% of the strongest-coupling tunneling pathways between the atoms of the electron donor and acceptor. Our theoretical analysis finds that adenine serves the critical role of fine-tuning rather than maximizing the donor-acceptor coupling within the range appropriate for the repair function.

Interference between Molecular and Photon Field-Mediated Electron Transfer Coupling Pathways in Cavities

Cavity polaritonics creates novel opportunities to direct chemical reactions. Electron transfer (ET) reactions are among the simplest reactions, and they underpin energy conversion. New strategies to manipulate and direct electron flow at the nanoscale are of particular interest in biochemistry, energy science, bioinspired materials science, and chemistry. We show that optical cavities can modulate electron transfer pathway interferences and ET rates in donor-bridge-acceptor (DBA) systems. We derive the rate for DBA electron transfer when the molecules are coupled to cavity modes, emphasizing novel cavity-induced pathway interferences with the molecular electronic coupling pathways, as these interferences allow a new kind of ET rate tuning. The interference between the cavity-induced coupling pathways and the intrinsic molecular coupling pathway is dependent on the cavity properties. Thus, manipulating the interference between the cavity-induced DA coupling and the bridge-mediated coupling offers an approach to direct and manipulate charge flow at the nanoscale.

Regulating Singlet-Triplet Energy Gaps through Substituent-Driven Modulation of the Exchange and Coulomb Interactions

Control of the singlet-triplet energy gap (ΔE_ST) is central to realizing productive energy conversion reactions, photochemical reaction trajectories, and emergent applications that exploit molecular spin physics. Despite this, no systematic methods have been defined to tune ΔE_ST in simple molecular frameworks, let alone by an approach that also holds chromophore size and electronic structural parameters (such as the HOMO-LUMO gap) constant. Using a combination of molecular design, photophysical and potentiometric experiments, and quantum chemical analyses, we show that the degree of electron-electron repulsion in excited singlet and triplet states may be finely controlled through the substitution pattern of a simple porphyrin absorber, enabling regulation of relative electronically excited singlet and triplet state energies by the designed restriction of the electron-electron Coulomb (J) and exchange (K) interaction magnitudes. This approach modulates the ΔE_ST magnitude by controlling the densities of state in the occupied and virtual molecular orbital manifolds, natural transition orbital polarization, and the relative contributions of one electron transitions involving select natural transition orbital pairs. This road map, which regulates electron density overlaps in the occupied and virtual states that define the singlet and triplet wave functions of these chromophores, enables new approaches to preserve excitation energy despite intersystem crossing.

Mapping Simulated Two-Dimensional Spectra to Molecular Models Using Machine Learning

Two-dimensional (2D) spectroscopy encodes molecular properties and dynamics into expansive spectral data sets. Translating these data into meaningful chemical insights is challenging because of the many ways chemical properties can influence the spectra. To address the task of extracting chemical information from 2D spectroscopy, we study the capacity of simple feedforward neural networks (NNs) to map simulated 2D electronic spectra to underlying physical Hamiltonians. We examined hundreds of simulated 2D spectra corresponding to monomers and dimers with varied Franck-Condon active vibrations and monomer-monomer electronic couplings. We find the NNs are able to correctly characterize most Hamiltonian parameters in this study with an accuracy above 90%. Our results demonstrate that NNs can aid in interpreting 2D spectra, leading from spectroscopic features to underlying effective Hamiltonians.

Cofactor Dynamics Couples the Protein Surface to the Heme in Cytochrome c, Facilitating Electron Transfer

Electron transport through biomolecules and in biological transport networks is of great importance to bioenergetics and biocatalysis. More generally, it is of crucial importance to understand how the pathways that connect buried metallocofactors to other cofactors, and to protein surfaces, affect the biological chemistry of metalloproteins. In terms of electron transfer (ET), the strongest coupling pathways usually comprise covalent and hydrogen bonded networks, with a limited number of through-space contacts. Herein, we set out to determine the relative roles of hydrogen bonds involved in ET via an established heme-to-surface tunneling pathway in cytochrome (cyt) c (i.e., heme-W59-D60-E61-N62). A series of cyt c variants were produced where a ruthenium tris(diimine) photooxidant was placed at position 62 via covalent modification of the N62C residue. Surprisingly, variants where the H-bonding residues W59 and D60 were replaced (i.e., W59F and D60A) showed no change in ET rate from the ferrous heme to Ru(III). In contrast, changing the composition of an alternative tunneling pathway (i.e., heme-M64-N63-C62) with the M64L substitution shows a factor of 2 decrease in the rate of heme-to-Ru ET. This pathway involves a through-space tunneling step between the heme and M64 residue, and such steps are usually disfavored. To rationalize why the heme-M64-N63-C62 is preferred, molecular dynamics (MD) simulations and Pathways analysis were employed. These simulations show that the change in heme-Ru ET rates is attributed to different conformations with compressed donor-acceptor distances, by ∼2 Å in pathway distance, in the M64-containing protein as compared to the M64L protein. The change in distance is correlated with changes in the electronic coupling that are in accord with the experimentally observed heme-Ru ET rates. Remarkably, the M64L variation at the core of the protein translates to changes in cofactor dynamics at the protein surface. The surface changes identified by MD simulations include dynamic anion-π and dipole-dipole interactions. These interactions influence the strength of tunneling pathways and ET rates by facilitating decreases in through-space tunneling distances in key coupling pathways.

A Chirality-Based Quantum Leap

There is increasing interest in the study of chiral degrees of freedom occurring in matter and in electromagnetic fields. Opportunities in quantum sciences will likely exploit two main areas that are the focus of this Review: (1) recent observations of the chiral-induced spin selectivity (CISS) effect in chiral molecules and engineered nanomaterials and (2) rapidly evolving nanophotonic strategies designed to amplify chiral light-matter interactions. On the one hand, the CISS effect underpins the observation that charge transport through nanoscopic chiral structures favors a particular electronic spin orientation, resulting in large room-temperature spin polarizations. Observations of the CISS effect suggest opportunities for spin control and for the design and fabrication of room-temperature quantum devices from the bottom up, with atomic-scale precision and molecular modularity. On the other hand, chiral-optical effects that depend on both spin- and orbital-angular momentum of photons could offer key advantages in all-optical and quantum information technologies. In particular, amplification of these chiral light-matter interactions using rationally designed plasmonic and dielectric nanomaterials provide approaches to manipulate light intensity, polarization, and phase in confined nanoscale geometries. Any technology that relies on optimal charge transport, or optical control and readout, including quantum devices for logic, sensing, and storage, may benefit from chiral quantum properties. These properties can be theoretically and experimentally investigated from a quantum information perspective, which has not yet been fully developed. There are uncharted implications for the quantum sciences once chiral couplings can be engineered to control the storage, transduction, and manipulation of quantum information. This forward-looking Review provides a survey of the experimental and theoretical fundamentals of chiral-influenced quantum effects and presents a vision for their possible future roles in enabling room-temperature quantum technologies.

Synthetic Control of Exciton Dynamics in Bioinspired Cofacial Porphyrin Dimers

Understanding how the complex interplay among excitonic interactions, vibronic couplings, and reorganization energy determines coherence-enabled transport mechanisms is a grand challenge with both foundational implications and potential payoffs for energy science. We use a combined experimental and theoretical approach to show how a modest change in structure may be used to modify the exciton delocalization, tune electronic and vibrational coherences, and alter the mechanism of exciton transfer in covalently linked cofacial Zn-porphyrin dimers (meso-beta linked ABm-β and meso-meso linked AAm-m). While both ABm-β and AAm-m feature zinc porphyrins linked by a 1,2-phenylene bridge, differences in the interporphyrin connectivity set the lateral shift between macrocycles, reducing electronic coupling in ABm-β and resulting in a localized exciton. Pump-probe experiments show that the exciton dynamics is faster by almost an order of magnitude in the strongly coupled AAm-m dimer, and two-dimensional electronic spectroscopy (2DES) identifies a vibronic coherence that is absent in ABm-β. Theoretical studies indicate how the interchromophore interactions in these structures, and their system-bath couplings, influence excitonic delocalization and vibronic coherence-enabled rapid exciton transport dynamics. Real-time path integral calculations reproduce the exciton transfer kinetics observed experimentally and find that the linking-modulated exciton delocalization strongly enhances the contribution of vibronic coherences to the exciton transfer mechanism, and that this coherence accelerates the exciton transfer dynamics. These benchmark molecular design, 2DES, and theoretical studies provide a foundation for directed explorations of nonclassical effects on exciton dynamics in multiporphyrin assemblies.

Discovery of the Xenon-Protein Interactome Using Large-Scale Measurements of Protein Folding and Stability

The intermolecular interactions of noble gases in biological systems are associated with numerous biochemical responses, including apoptosis, inflammation, anesthesia, analgesia, and neuroprotection. The molecular modes of action underlying these responses are largely unknown. This is in large part due to the limited experimental techniques to study protein-gas interactions. The few techniques that are amenable to such studies are relatively low-throughput and require large amounts of purified proteins. Thus, they do not enable the large-scale analyses that are useful for protein target discovery. Here, we report the application of stability of proteins from rates of oxidation (SPROX) and limited proteolysis (LiP) methodologies to detect protein-xenon interactions on the proteomic scale using protein folding stability measurements. Over 5000 methionine-containing peptides and over 5000 semi-tryptic peptides, mapping to ∼1500 and ∼950 proteins, respectively, in the yeast proteome, were assayed for Xe-interacting activity using the SPROX and LiP techniques. The SPROX and LiP analyses identified 31 and 60 Xe-interacting proteins, respectively, none of which were previously known to bind Xe. A bioinformatics analysis of the proteomic results revealed that these Xe-interacting proteins were enriched in those involved in ATP-driven processes. A fraction of the protein targets that were identified are tied to previously established modes of action related to xenon's anesthetic and organoprotective properties. These results enrich our knowledge and understanding of biologically relevant xenon interactions. The sample preparation protocols and analytical methodologies developed here for xenon are also generally applicable to the discovery of a wide range of other protein-gas interactions in complex biological mixtures, such as cell lysates.

Mechanism of Side Chain-Controlled Proton Conductivity in Bioinspired Peptidic Nanostructures

Bioinspired peptide assemblies are promising candidates for use as proton-conducting materials in electrochemical devices and other advanced technologies. Progress toward applications requires establishing foundational structure-function relationships for transport in these materials. This experimental-theoretical study sheds light on how the molecular structure and proton conduction are linked in three synthetic cyclic peptide nanotube assemblies that comprise the three canonical basic amino acids (lysine, arginine, and histidine). Experiments find an order of magnitude higher proton conductivity for lysine-containing peptide assemblies compared to histidine and arginine containing assemblies. The simulations indicate that, upon peptide assembly, the basic amino acid side chains are close enough to enable direct proton transfer. The proton transfer kinetics is determined in the simulations to be governed by the structure and flexibility of the side chains. Together, experiments and theory indicate that the proton mobility is the main determinant of proton conductivity, critical for the performance of peptide-based devices.

Excited-State Dynamics and Nonlinear Optical Properties of Hyperpolarizable Chromophores Based on Conjugated Bis(terpyridyl)Ru(II) and Palladium and Platinum Porphyrinic Components: Impact of Heavy Metals upon Supermolecular Electro-Optic Properties

A new series of strongly coupled oscillators based upon (porphinato)Pd, (porphinato)Pt, and bis(terpyridyl)ruthenium(II) building blocks is described. These RuPPd, RuPPt, RuPPdRu, and RuPPtRu chromophores feature bis(terpyridyl)Ru(II) moieties connected to the (porphinato)metal unit via an ethyne linker that bridges the 4'-terpyridyl and porphyrin macrocycle meso-carbon positions. Pump-probe transient optical data demonstrate sub-picosecond excited singlet-to-triplet-state relaxation. The relaxed lowest-energy triplet (T₁) excited states of these chromophores feature absorption manifolds that span the 800-1200 nm spectral region, microsecond triplet-state lifetimes, and large absorptive extinction coefficients [ε(T₁ → T_n) > 4 × 10⁴ M^-1 cm^-1]. Dynamic hyperpolarizability (β_λ) values were determined from hyper-Rayleigh light scattering (HRS) measurements carried out at several incident irradiation wavelengths over the 800-1500 nm spectral region. Relative to benchmark RuPZn and RuPZnRu chromophores which showed large β_HRS values over the 1200-1600 nm range, RuPPd, RuPPt, RuPPdRu, and RuPPtRu displayed large β_HRS values over the 850-1200 nm region. Generalized Thomas-Kuhn sum (TKS) rules and experimental hyperpolarizability values were utilized to determine excited state-to-excited state transition dipole terms from experimental electronic absorption data and thus assessed frequency-dependent β_λ values, including two- and three-level contributions for both β_zzz and β_xzx tensor components to the RuPPd, RuPPt, RuPPdRu, and RuPPtRu hyperpolarizability spectra. These analyses qualitatively rationalize how the β_zzz and β_xzx tensor elements influence the observed irradiation wavelength-dependent hyperpolarizability magnitudes. The TKS analysis suggests that supermolecules related to RuPPd, RuPPt, RuPPdRu, and RuPPtRu will likely feature intricate dependences of experimentally determined β_HRS values as a function of irradiation wavelength that derive from substantial singlet-triplet mixing, and complex interactions among multiple different β tensor components that modulate the long wavelength regime of the nonlinear optical response.

Electron ratcheting in self-assembled soft matter

Ratcheted multi-step hopping electron transfer systems can plausibly produce directional charge transport over very large distances without requiring a source-drain voltage bias. We examine molecular strategies to realize ratcheted charge transport based on multi-step charge hopping, and we illustrate two ratcheting mechanisms with examples based on DNA structures. The charge transport times and currents that may be generated in these assemblies are also estimated using kinetic simulations. The first ratcheting mechanism described for nanoscale systems requires local electric fields on the 10⁹ V/m scale to realize nearly 100% population transport. The second ratcheting mechanism for even larger systems, based on electrochemical gating, is estimated to generate currents as large as 0.1 pA for DNA structures that are a few μm in length with a gate voltage of about 5 V, a magnitude comparable to currents measured in DNA wires at the nanoscale when a source-drain voltage bias of similar magnitude is applied, suggesting an approach to considerably extend the distance range over which DNA charge transport devices may operate.

Why Do Most Aromatics Fail to Support Hole Hopping in the Cytochrome c Peroxidase-Cytochrome c Complex?

Electron transport through aromatic species (especially tryptophan and tyrosine) plays a central role in water splitting, redox signaling, oxidative damage protection, and bioenergetics. The cytochrome c peroxidase (CcP)-cytochrome c (Cc) complex (CcP:Cc) is used widely to study interprotein electron transfer (ET) mechanisms. Tryptophan 191 (Trp191) of CcP supports hole hopping charge recombination in the CcP:Cc complex. Experimental studies find that when Trp191 is substituted by tyrosine, phenylalanine, or redox-active aniline derivatives bound in the W191G cavity, enzymatic activity and charge recombination rates both decrease. Theoretical analysis of these CcP:Cc complexes finds that the ET kinetics depend strongly on the chemistry of the modified Trp site. The computed electronic couplings in the W191F and W191G species are orders of magnitude smaller than in the native protein, due largely to the absence of a hopping intermediate and the large tunneling distance. Small molecules bound in the W191G cavity are weakly coupled electronically to the Cc heme, and the structural disorder of the guest molecule in the binding pocket may contribute further to the lack of enzymatic activity. The couplings in W191Y are not substantially weakened compared to the native species, but the redox potential difference for tyrosine vs tryptophan oxidation accounts for the slower rate in the Tyr mutant. Thus, theoretical analysis explains why only the native Trp supports rapid hole hopping in the CcP:Cc complex. Favorable free energies and electronic couplings are essential for establishing an efficient hole hopping relay in this protein-protein complex.

Oxalate decarboxylase uses electron hole hopping for catalysis

The hexameric low-pH stress response enzyme oxalate decarboxylase catalyzes the decarboxylation of the oxalate mono-anion in the soil bacterium Bacillus subtilis. A single protein subunit contains two Mn-binding cupin domains, and catalysis depends on Mn(III) at the N-terminal site. The present study suggests a mechanistic function for the C-terminal Mn as an electron hole donor for the N-terminal Mn. The resulting spatial separation of the radical intermediates directs the chemistry toward decarboxylation of the substrate. A π-stacked tryptophan pair (W96/W274) links two neighboring protein subunits together, thus reducing the Mn-to-Mn distance from 25.9 Å (intrasubunit) to 21.5 Å (intersubunit). Here, we used theoretical analysis of electron hole-hopping paths through redox-active sites in the enzyme combined with site-directed mutagenesis and X-ray crystallography to demonstrate that this tryptophan pair supports effective electron hole hopping between the C-terminal Mn of one subunit and the N-terminal Mn of the other subunit through two short hops of ∼8.5 Å. Replacement of W96, W274, or both with phenylalanine led to a large reduction in catalytic efficiency, whereas replacement with tyrosine led to recovery of most of this activity. W96F and W96Y mutants share the wildtype tertiary structure. Two additional hole-hopping networks were identified leading from the Mn ions to the protein surface, potentially protecting the enzyme from high Mn oxidation states during turnover. Our findings strongly suggest that multistep hole-hopping transport between the two Mn ions is required for enzymatic function, adding to the growing examples of proteins that employ aromatic residues as hopping stations.

Temperature Dependence of Charge and Spin Transfer in Azurin

The steady-state charge and spin transfer yields were measured for three different Ru-modified azurin derivatives in protein films on silver electrodes. While the charge-transfer yields exhibit weak temperature dependences, consistent with operation of a near activation-less mechanism, the spin selectivity of the electron transfer improves as temperature increases. This enhancement of spin selectivity with temperature is explained by a vibrationally induced spin exchange interaction between the Cu(II) and its chiral ligands. These results indicate that distinct mechanisms control charge and spin transfer within proteins. As with electron charge transfer, proteins deliver polarized electron spins with a yield that depends on the protein's structure. This finding suggests a new role for protein structure in biochemical redox processes.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G₀, where G₀ is the quantum of conductance, for 5 nm duplexes.

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

Quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyl)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI₂ (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI₂^•- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI₂^•-. Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the ^3*ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

Correlation between Charge Transport and Base Excision Repair in the MutY-DNA Glycosylase

Experimental evidence suggests that DNA-mediated redox signaling between high-potential [Fe₄S₄] proteins is relevant to DNA replication and repair processes, and protein-mediated charge transfer (CT) between [Fe₄S₄] clusters and nucleic acids is a fundamental process of the signaling and repair mechanisms. We analyzed the dominant CT pathways in the base excision repair glycosylase MutY using molecular dynamics simulations and hole hopping pathway analysis. We find that the adenine nucleobase of the mismatched A·oxoG DNA base pair facilitates [Fe₄S₄]-DNA CT prior to adenine excision by MutY. We also find that the R153L mutation in MutY (linked to colorectal adenomatous polyposis) influences the dominant [Fe4S4]-DNA CT pathways and appreciably decreases their effective CT rates.

Electrostatic Field-Induced Oscillator Strength Focusing in Molecules

The development of light-harvesting devices based on molecular materials depends critically on the ability to focus the electronic oscillator strength of molecules into the UV-vis spectral window. Typical molecular chromophores have only about 1% of their total electronic oscillator strength in this spectral region and thus perform at a small fraction of their possible effectiveness. This theoretical study finds that the electronic oscillator strength of polyenes in the UV-vis region may be enhanced by 1 order of magnitude using electrostatic fields, motivating specific experimental studies of oscillator strength focusing. We find scaling relationships between the polyene length, the intensity of the applied field, and the field-induced increase in oscillator strength that are useful for the implementation of light-harvesting strategies based on polyenes.

Conductance and configuration of molecular gold-water-gold junctions under electric fields

Water molecules can mediate charge transfer in biological and chemical reactions by forming electronic coupling pathways. Understanding the mechanism requires a molecular-level electrical characterization of water. Here, we describe the measurement of single water molecular conductance at room temperature, characterize the structure of water molecules using infrared spectroscopy, and perform theoretical studies to assist in the interpretation of the experimental data. The study reveals two distinct states of water, corresponding to a parallel and perpendicular orientation of the molecules. Water molecules switch from parallel to perpendicular orientations on applying an electric field, producing switching from high to low conductance states, thus enabling the determination of single water molecular dipole moments. The work further shows that water-water interactions affect the atomic scale configuration and conductance of water molecules. These findings demonstrate the importance of the discrete nature of water molecules in electron transfer and set limits on water-mediated electron transfer rates.

A Robust Bioderived Wavelength-Specific Photosensor Based on BLUF Proteins

Photosensitive proteins are naturally evolved photosensors that often respond to light signals of specific wavelengths. However, their poor stability under ambient conditions hinders their applications in non-biological settings. In this proof-of-principle study, we grafted the blue light using flavin (BLUF) protein reconstructed with flavin adenine dinucleotide (FAD) or roseoflavin (RoF) onto pristine graphene, and achieved selective sensitivity at 450 nm or 500 nm, respectively. We improved the thermal and operational stability substantially via structure-guided cross-linking, achieving 6-month stability under ambient condition and normal operation at temperatures up to 200 °C. Furthermore, the device exhibited rare negative photoconductivity behavior. The origins of this negative photoconductivity behavior were elucidated via a combination of experimental and theoretical analysis. In the photoelectric conversion studies, holes from photoexcited flavin migrated to graphene and recombined with electrons. The device allows facile modulation and detection of charge transfer, and provides a versatile platform for future studies of photoinduced charge transfer in biosensors as well as the development of stable wavelength-selective biophotosensors.

Revisiting the Hole Size in Double Helical DNA with Localized Orbital Scaling Corrections

The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier's wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine-cytosine (GC) base pairs and among 3 adenine-thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π-π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.

Mutation effects on charge transport through the p58c iron-sulfur protein

Growing experimental evidence indicates that iron–sulfur proteins play key roles in DNA repair and replication. In particular, charge transport between [Fe₄S₄] clusters, mediated by proteins and DNA, may convey signals to coordinate enzyme action. Human primase is a well studied [Fe₄S₄] protein, and its p58c domain (which contains an [Fe₄S₄] cluster) plays a role in the initiation of DNA replication. The Y345C mutation in p58c is linked to gastric tumors and may influence the protein-mediated charge transport. The complexity of protein–DNA systems, and the intricate electronic structure of [Fe₄S₄] clusters, have impeded progress into understanding functional charge transport in these systems. In this study, we built force fields to describe the high potential [Fe₄S₄] cluster in both oxidation states. The parameterization is compatible with AMBER force fields and enabled well-balanced molecular dynamics simulations of the p58c–RNA/DNA complex relevant to the initiation of DNA replication. Using the molecular mechanics Poisson–Boltzmann and surface area solvation method on the molecular dynamics trajectories, we find that the p58c mutation induces a modest change in the p58c–duplex binding free energy in agreement with recent experiments. Through kinetic modeling and analysis, we identify key features of the main charge transport pathways in p58c. In particular, we find that the Y345C mutation partially changes the composition and frequency of the most efficient (and potentially relevant to the biological function) charge transport pathways between the [Fe₄S₄] cluster and the duplex. Moreover, our approach sets the stage for a deeper understanding of functional charge transfer in [Fe₄S₄] protein–DNA complexes.

Functional electron transfer between the [Fe₄S₄] cluster and the nucleic acid is impacted by a Y345C mutation in the p58c subunit of human primase.

Symmetry controlled photo-selection and charge separation in butadiyne-bridged donor–bridge–acceptor compounds

Electron transfer (ET) in donor–bridge–acceptor (DBA) compounds featuring alkyne bridges depends strongly on the torsion angle between the donor and acceptor.

Orientational Dependence of Cofacial Porphyrin-Quinone Electronic Interactions within the Strong Coupling Regime

We examine the relative magnitudes of electronic coupling H_DA in two face-to-face rigid and diastereomeric (porphinato)zinc(II)-quinone (PZn-Q) assemblies, 1β-ZnA and 1β-ZnB, in which the six quinonyl carbon atoms lie in virtually identical arrangements relative to the PZn plane at sub-van der Waals donor-acceptor (D-A) interplanar separations. Steady-state and time-resolved transient optical data and computational studies show that minor differences in relative D-A cofacial orientation give rise to disparate H_DA magnitudes for both photoinduced charge separation (CS) and thermal charge recombination (CR). Time-dependent density functional theory (TDDFT) computations illuminate the nature of direct charge transfer states and the electronic structural factors that give rise to these differential H_DAs. These data show more extensive mixing of locally excited (LE) and CS states in 1β-ZnA relative to 1β-ZnB and that these H_DA differences track the magnitudes of electronic coupling matrix elements determined from steady-state electronic spectral data and thermal CR rate constants measured via pump-probe spectroscopy. Collectively, this work shows that electron transfer dynamics may be manipulated in cofacial D-A systems, even at sub-van der Waals contact, provided that conformational rigidity precludes structural fluctuations that modulate D-A interactions on the charge transfer time scale.

Assessing Possible Mechanisms of Micrometer-Scale Electron Transfer in Heme-Free Geobacter sulfurreducens Pili

The electrically conductive pili of Geobacter sulfurreducens are of both fundamental and practical interest. They facilitate extracellular and interspecies electron transfer (ET) and also provide an electrical interface between living and nonliving systems. We examine the possible mechanisms of G. sulfurreducens electron transfer in regimes ranging from incoherent to coherent transport. For plausible ET parameters, electron transfer in G. sulfurreducens bacterial nanowires mediated only by the protein is predicted to be dominated by incoherent hopping between phenylalanine (Phe) and tyrosine (Tyr) residues that are 3 to 4 Å apart, where Phe residues in the hopping pathways may create delocalized "islands." This mechanism could be accessible in the presence of strong oxidants that are capable of oxidizing Phe and Tyr residues. We also examine the physical requirements needed to sustain biological respiration via nanowires. We find that the hopping regimes with ET rates on the order of 10⁸ s^-1 between Phe islands and Tyr residues, and conductivities on the order of mS/cm, can support ET fluxes that are compatible with cellular respiration rates, although sustaining this delocalization in the heterogeneous protein environment may be challenging. Computed values of fully coherent electron fluxes through the pili are orders of magnitude too low to support microbial respiration. We suggest experimental probes of the transport mechanism based on mutant studies to examine the roles of aromatic amino acids and yet to be identified redox cofactors.

Why Are DNA and Protein Electron Transfer So Different?

The corpus of electron transfer (ET) theory provides considerable power to describe the kinetics and dynamics of electron flow at the nanoscale. How is it, then, that nucleic acid (NA) ET continues to surprise, while protein-mediated ET is relatively free of mechanistic bombshells? I suggest that this difference originates in the distinct electronic energy landscapes for the two classes of reactions. In proteins, the donor/acceptor-to-bridge energy gap is typically several-fold larger than in NAs. NA ET can access tunneling, hopping, and resonant transport among the bases, and fluctuations can enable switching among mechanisms; protein ET is restricted to tunneling among redox active cofactors and, under strongly oxidizing conditions, a few privileged amino acid side chains. This review aims to provide conceptual unity to DNA and protein ET reaction mechanisms. The establishment of a unified mechanistic framework enabled the successful design of NA experiments that switch electronic coherence effects on and off for ET processes on a length scale of multiple nanometers and promises to provide inroads to directing and detecting charge flow in soft-wet matter.

A single AT-GC exchange can modulate charge transfer-induced p53-DNA dissociation

Using molecular dynamics simulations and electronic structure theory, we shed light on the charge dynamics that causes the differential interaction of tumor suppressor protein p53 with the p21 and Gadd45 genes in response to oxidative stress. We show that the sequence dependence of this selectivity results from competing charge transfer to the protein and through the DNA, with implications on the use of genome editing tools to influence the p53 regulatory function.

Electron bifurcation: progress and grand challenges

Electron bifurcation moves electrons from a two-electron donor to reduce two spatially separated one-electron acceptors.

Electron transfer characteristics of 2'-deoxy-2'-fluoro-arabinonucleic acid, a nucleic acid with enhanced chemical stability

The non-biological nucleic acid 2'-deoxy-2'-fluoro-arabinonucleic acid (2'F-ANA) may be of use because of its higher chemical stability than DNA in terms of resistance to hydrolysis and nuclease degradation. In order to investigate the charge transfer characteristics of 2'F-ANA, of relevance to applications in nucleic acid-based biosensors and chip technologies, we compare the electronic couplings for hole transfer between stacked nucleobase pairs in DNA and 2'F-ANA by carrying out density functional theory (DFT) calculations on geometries taken from molecular dynamics simulations. We find similar averages and distribution widths of the base-pair couplings in the two systems. On the basis of this result, 2'F-ANA is expected to have charge transfer properties similar to those of DNA, while offering the advantage of enhanced chemical stability. As such, 2'F-ANA may serve as a possible alternative to DNA for use in a broad range of nanobiotechnological applications. Furthermore, we show that the (experimentally observed) enhanced chemical stability resulting from the backbone modifications does not cause reduced fluctuations of the base-pair electronic couplings around the values found for "ideal" B-DNA (with standard step parameter values). Our study also supports the use of a DFT implementation, with the M11 functional, of the wave function overlap method to compute effective electronic couplings in nucleic acid systems.

Directing Charge Transfer in Quantum Dot Assemblies

The optical and electronic properties of semiconductor quantum dots (QDs) make them attractive candidates for applications in photovoltaics, spintronics, photocatalysis, and optoelectronics. Understanding how to control the flow of charge in QD assemblies is essential for realizing novel applications. This Account explores some unique characteristics of charge transport in QD dyads, triads, and their assemblies. The emerging features of these assemblies that provide new opportunities to manipulate charge flow at the nanoscale are (1) cascading energy landscapes and band offsets to inhibit charge recombination, (2) electrostatic fields that direct charge flow through QD-QD and QD-conjugated polymer junctions, and (3) QD chirality and chiral imprinting that promotes vectorial electron and spin selective transport. Charge flow kinetics is determined by a combination of familiar electron transfer parameters (reaction free energy, reorganization energy, and electronic coupling), donor and acceptor electronic densities of states, and internal electric fields. Electron transfer and electronic structure theory, combined with kinetic modeling, place the measured kinetics of QD electron transfer donor-acceptor assemblies into a unified conceptual context. The experimental transfer rates measured in these systems depend upon structure and the internal electric fields that are present in the assemblies. A negatively charged donor and positively charged acceptor, for example, facilitates (inhibits) electron (hole) transfer, while an electric field of opposite orientation (reversal of charges) inhibits (promotes) electron (hole) transfer. These and other emerging rules that govern charge flow in NP assemblies provide a strategy to design the directionality and yield of interfacial charge transport. Chirality at the nanoscale can induce spin selective charge transport, providing new ways to direct charge (and spin) flow in QD assemblies. Magnetoresistance and magnetic conductive probe atomic force microscopy experiments show spin selective electron transport for chirally imprinted QD assemblies. Photoinduced electron transfer from achiral donor-QDs to chiral acceptor-QDs depends on the electron spin and chiroptical properties of the acceptor-QDs. These assemblies show transport characteristics that correlate with features of the QDs' circular dichroism spectra, presenting intriguing challenges to theory, and indicating that spectroscopic signatures may assist in the design and diagnosis of functional molecular assemblies. Theoretical and experimental studies of charge transport in well-defined QD assemblies are establishing design principles for vectorial charge transport and are also refining questions surrounding the mechanism and control of these processes. These intensified efforts are forging links between fundamental discoveries regarding mechanism and practical applications for these novel assembled nanostructures.

Energy Transduction in Nitrogenase

Nitrogenase is a complicated two-component enzyme system that uses ATP binding and hydrolysis energy to achieve one of the most difficult chemical reactions in nature, the reduction of N₂ to NH₃. One component of the Mo-based nitrogenase system, Fe protein, delivers electrons one at a time to the second component, the catalytic MoFe protein. This process occurs through a series of synchronized events collectively called the "Fe protein cycle". Elucidating details of the events associated with this cycle has constituted an important challenge in understanding the nitrogenase mechanism. Electron delivery is a multistep process involving three metal clusters with intra- and interprotein events. It is proposed that the first electron transfer event is a gated intraprotein transfer of one electron from the MoFe protein P-cluster to the FeMo cofactor. Measurement of the effect of osmotic pressure on the rate of this electron transfer process revealed that it is gated by protein conformational changes. This first electron transfer is activated by binding of the Fe protein containing two bound ATP molecules. The mechanism of how this protein-protein association triggers electron transfer remains unknown. The second electron transfer event is proposed to be a rapid interprotein "backfill" with transfer of one electron from the reduced Fe protein 4Fe-4S cluster to the oxidized P-cluster. In this way, electron delivery can be viewed as a case of "deficit spending". Such a deficit-spending electron transfer process can be envisioned as a way to achieve one-direction electron flow, limiting the potential for back electron flow. Hydrolysis of two ATP molecules associated with the Fe protein occurs after the electron transfer and therefore is not used to directly drive the electron transfer. Rather, ATP hydrolysis is proposed to contribute to relaxation of the "activated" conformational state associated with the ATP form of the complex, with the free energy from ATP hydrolysis being used to pay back energy associated with component protein association and electron transfer. Release of inorganic phosphate (Pi) and protein-protein dissociation follow electron transfer and ATP hydrolysis. The rate-limiting step for the Fe protein cycle is not dissociation of the two proteins, as previously believed, but rather is release of Pi after ATP hydrolysis, which is then followed by rapid protein-protein complex dissociation. Nitrogenase is composed of two catalytic halves that do not function independently but rather exhibit anticooperative nuclear motion in which electron transfer in one-half of the complex partially inhibits electron transfer and ATP hydrolysis in the other half. Calculations indicated the existence of anticooperative interactions across the entire nitrogenase complex, suggesting a mechanism for the control of events on opposite ends of this large complex. The mechanistic necessity for this anticooperative process remains unknown. This Account presents a working model for how all of these processes work together in the nitrogenase "machine" to transduce the energy from ATP binding and hydrolysis to drive N₂ reduction.

Control of electron transfer in nitrogenase

The bacterial enzyme nitrogenase achieves the reduction of dinitrogen (N₂) to ammonia (NH₃) utilizing electrons, protons, and energy from the hydrolysis of ATP. Building on earlier foundational knowledge, recent studies provide molecular-level details on how the energy of ATP hydrolysis is utilized, the sequencing of multiple electron transfer events, and the nature of energy transduction across this large protein complex. Here, we review the state of knowledge about energy transduction in nitrogenase.

A new era for electron bifurcation

Electron bifurcation, or the coupling of exergonic and endergonic oxidation-reduction reactions, was discovered by Peter Mitchell and provides an elegant mechanism to rationalize and understand the logic that underpins the Q cycle of the respiratory chain. Thought to be a unique reaction of respiratory complex III for nearly 40 years, about a decade ago Wolfgang Buckel and Rudolf Thauer discovered that flavin-based electron bifurcation is also an important component of anaerobic microbial metabolism. Their discovery spawned a surge of research activity, providing a basis to understand flavin-based bifurcation, forging fundamental parallels with Mitchell's Q cycle and leading to the proposal of metal-based bifurcating enzymes. New insights into the mechanism of electron bifurcation provide a foundation to establish the unifying principles and essential elements of this fascinating biochemical phenomenon.

How can infra-red excitation both accelerate and slow charge transfer in the same molecule?

A UV-IR-Vis 3-pulse study of infra-red induced changes to electron transfer (ET) rates in a donor-bridge-acceptor species finds that charge-separation rates are slowed, while charge-recombination rates are accelerated as a result of IR excitation during the reaction. We explore the underpinning mechanisms for this behavior, studying IR-induced changes to the donor-acceptor coupling, to the validity of the Condon approximation, and to the reaction coordinate distribution. We find that the dominant IR-induced rate effects in the species studied arise from changes to the density of states in the Marcus curve crossing region. That is, IR perturbation changes the probability of accessing the activated complex for the ET reactions. IR excitation diminishes the population of the activated complex for forward (activationless) ET, thus slowing the rate. However, IR excitation increases the population of the activated complex for (highly activated) charge recombination ET, thus accelerating the charge recombination rate.

Charge splitters and charge transport junctions based on guanine quadruplexes

Self-assembling circuit elements, such as current splitters or combiners at the molecular scale, require the design of building blocks with three or more terminals. A promising material for such building blocks is DNA, wherein multiple strands can self-assemble into multi-ended junctions, and nucleobase stacks can transport charge over long distances. However, nucleobase stacking is often disrupted at junction points, hindering electric charge transport between the two terminals of the junction. Here, we show that a guanine-quadruplex (G4) motif can be used as a connector element for a multi-ended DNA junction. By attaching specific terminal groups to the motif, we demonstrate that charges can enter the structure from one terminal at one end of a three-way G4 motif, and can exit from one of two terminals at the other end with minimal carrier transport attenuation. Moreover, we study four-way G4 junction structures by performing theoretical calculations to assist in the design and optimization of these connectors.

Determinants of Photolyase's DNA Repair Mechanism in Mesophiles and Extremophiles

Light-driven DNA repair by extremophilic photolyases is of tremendous importance for understanding the early development of life on Earth. The mechanism for flavin adenine dinucleotide repair of DNA lesions is the subject of debate and has been studied mainly in mesophilic species. In particular, the role of adenine in the repair process is poorly understood. Using molecular docking, molecular dynamics simulations, electronic structure calculations, and electron tunneling pathways analysis, we examined adenine's role in DNA repair in four photolyases that thrive at different temperatures. Our results indicate that the contribution of adenine to the electronic coupling between the flavin and the cyclobutane pyrimidine dimer lesion to be repaired is significant in three (one mesophilic and two extremophilic) of the four enzymes studied. Our analysis suggests that thermophilic and hyperthermophilic photolyases have evolved structurally to preserve the functional position (and thus the catalytic function) of adenine at their high temperatures of operation. Water molecules can compete with adenine in establishing the strongest coupling pathway for the electron transfer repair process, but the adenine contribution remains substantial. The present study also reconciles prior seemingly contradictory conclusions on the role of adenine in mesophile electron transfer repair reactions, showing how adenine-mediated superexchange is conformationally gated.

On the nature of organic and inorganic centers that bifurcate electrons, coupling exergonic and endergonic oxidation–reduction reactions

Bifurcating electrons to couple endergonic and exergonic electron-transfer reactions has been shown to have a key role in energy conserving redox enzymes.

Electron Bifurcation: Thermodynamics and Kinetics of Two-Electron Brokering in Biological Redox Chemistry

How can proteins drive two electrons from a redox active donor onto two acceptors at very different potentials and distances? And how can this transaction be conducted without dissipating very much energy or violating the laws of thermodynamics? Nature appears to have addressed these challenges by coupling thermodynamically uphill and downhill electron transfer reactions, using two-electron donor cofactors that have very different potentials for the removal of the first and second electron. Although electron bifurcation is carried out with near perfection from the standpoint of energy conservation and electron delivery yields, it is a biological energy transduction paradigm that has only come into focus recently. This Account provides an exegesis of the biophysical principles that underpin electron bifurcation. Remarkably, bifurcating electron transfer (ET) proteins typically send one electron uphill and one electron downhill by similar energies, such that the overall reaction is spontaneous, but not profligate. Electron bifurcation in the NADH-dependent reduced ferredoxin: NADP⁺ oxidoreductase I (Nfn) is explored in detail here. Recent experimental progress in understanding the structure and function of Nfn allows us to dissect its workings in the framework of modern ET theory. The first electron that leaves the two-electron donor flavin (L-FAD) executes a positive free energy "uphill" reaction, and the departure of this electron switches on a second thermodynamically spontaneous ET reaction from the flavin along a second pathway that moves electrons in the opposite direction and at a very different potential. The singly reduced ET products formed from the bifurcating flavin are more than two nanometers distant from each other. In Nfn, the second electron to leave the flavin is much more reducing than the first: the potentials are said to be "crossed." The eventually reduced cofactors, NADH and ferredoxin in the case of Nfn, perform crucial downstream redox processes of their own. We dissect the thermodynamics and kinetics of electron bifurcation in Nfn and find that the key features of electron bifurcation are (1) spatially separated transfer pathways that diverge from a two-electron donor, (2) one thermodynamically uphill and one downhill redox pathway, with a large negative shift in the donor's reduction potential after departure of the first electron, and (3) electron tunneling and activation factors that enable bifurcation, producing a 1:1 partitioning of electrons onto the two pathways. Electron bifurcation is found in the CO₂ reducing pathways of methanogenic archaea, in the hydrogen pathways of hydrogenases, in the nitrogen fixing pathway of Fix, and in the mitochondrial charge transfer chain of complex III, cytochrome bc₁. While crossed potentials may offer the biological advantage of producing tightly regulated high energy reactive species, neither kinetic nor thermodynamic considerations mandate crossed potentials to generate successful electron bifurcation. Taken together, the theoretical framework established here, focusing on the underpinning electron tunneling barriers and activation free energies, explains the logic of electron bifurcation that enables energy conversion and conservation in Nfn, points toward bioinspired schemes to execute multielectron redox chemistry, and establishes a roadmap for examining novel electron bifurcation networks in nature.