Theory of Tunneling Spectroscopy in a Mn12_{12} Single-Electron Transistor by Density-Functional Theory Methods

We consider tunneling transport through a Mn$_{12}$ molecular magnet using spin density functional theory. A tractable methodology for constructing many-body wavefunctions from Kohn-Sham orbitals allows for the determination of spin-dependent matrix elements for use in transport calculations. The tunneling conductance at finite bias is characterized by peaks representing transitions between spin multiplets, separated by an energy on the order of the magnetic anisotropy. The energy splitting of the spin multiplets and the spatial part of their many-body wave functions, describing the orbital degrees of freedom of the excess charge, strongly affect the electronic transport, and can lead to negative differential conductance.Comment: 4 pages, 3 figures, a revised version with minor change

Constraints on Fluctuations in Sparsely Characterized Biological Systems.

Biochemical processes are inherently stochastic, creating molecular fluctuations in otherwise identical cells. Such "noise" is widespread but has proven difficult to analyze because most systems are sparsely characterized at the single cell level and because nonlinear stochastic models are analytically intractable. Here, we exactly relate average abundances, lifetimes, step sizes, and covariances for any pair of components in complex stochastic reaction systems even when the dynamics of other components are left unspecified. Using basic mathematical inequalities, we then establish bounds for whole classes of systems. These bounds highlight fundamental trade-offs that show how efficient assembly processes must invariably exhibit large fluctuations in subunit levels and how eliminating fluctuations in one cellular component requires creating heterogeneity in another.The work was supported by grant 1137676 from the Division of Mathematical Sciences at the National Science Foundation, and grant GM081563 from the National Institutes of Health.This is the final version of the article. It first appeared from the American Physical Society via http://dx.doi.org/10.1103/PhysRevLett.116.05810

Kinetic Uncertainty Relations for the Control of Stochastic Reaction Networks.

Nonequilibrium stochastic reaction networks are commonly found in both biological and nonbiological systems, but have remained hard to analyze because small differences in rate functions or topology can change the dynamics drastically. Here, we conjecture exact quantitative inequalities that relate the extent of fluctuations in connected components, for various network topologies. Specifically, we find that regardless of how two components affect each other's production rates, it is impossible to suppress fluctuations below the uncontrolled equivalents for both components: one must increase its fluctuations for the other to be suppressed. For systems in which components control each other in ringlike structures, it appears that fluctuations can only be suppressed in one component if all other components instead increase fluctuations, compared to the case without control. Even the general N-component system-with arbitrary connections and parameters-must have at least one component with increased fluctuations to reduce fluctuations in others. In connected reaction networks it thus appears impossible to reduce the statistical uncertainty in all components, regardless of the control mechanisms or energy dissipation

The Smallest Molecular Switch

Ab-initio total energy calculations reveal benzene-dithiolate (BDT) molecules on a gold surface, contacted by a monoatomic gold STM tip to have two classes of low energy conformations with differing symmetries. Lateral motion of the tip or excitation of the molecule cause it to change from one conformation class to the other and to switch between a strongly and a weakly conducting state. Thus, surprisingly, despite their apparent simplicity these Au/BDT/Au nanowires are shown to be electrically bi-stable switches, the smallest two-terminal molecular switches to date. Experiments with a conventional or novel self-assembled STM are proposed to test these predictions.Comment: 8 pages, 3 figure

Inelastic Scattering in Metal-H2-Metal Junctions

We present first-principles calculations of the dI/dV characteristics of an H2 molecule sandwiched between Au and Pt electrodes in the presence of electron-phonon interactions. The conductance is found to decrease by a few percentage at threshold voltages corresponding to the excitation energy of longitudinal vibrations of the H2 molecule. In the case of Pt electrodes, the transverse vibrations can mediate transport through otherwise non-transmitting Pt $d$-channels leading to an increase in the differential conductance even though the hydrogen junction is characterized predominately by a single almost fully open transport channel. In the case of Au, the transverse modes do not affect the dI/dV because the Au d-states are too far below the Fermi level. A simple explanation of the first-principles results is given using scattering theory. Finally, we compare and discuss our results in relation to experimental data.Comment: Accepted in Phys. Rev.

Stochastic Simulations of the Repressilator Circuit

The genetic repressilator circuit consists of three transcription factors, or repressors, which negatively regulate each other in a cyclic manner. This circuit was synthetically constructed on plasmids in {\it Escherichia coli} and was found to exhibit oscillations in the concentrations of the three repressors. Since the repressors and their binding sites often appear in low copy numbers, the oscillations are noisy and irregular. Therefore, the repressilator circuit cannot be fully analyzed using deterministic methods such as rate-equations. Here we perform stochastic analysis of the repressilator circuit using the master equation and Monte Carlo simulations. It is found that fluctuations modify the range of conditions in which oscillations appear as well as their amplitude and period, compared to the deterministic equations. The deterministic and stochastic approaches coincide only in the limit in which all the relevant components, including free proteins, plasmids and bound proteins, appear in high copy numbers. We also find that subtle features such as cooperative binding and bound-repressor degradation strongly affect the existence and properties of the oscillations.Comment: Accepted to PR

Charging induced asymmetry in molecular conductors

We investigate the origin of asymmetry in various measured current-voltage (I-V) characteristics of molecules with no inherent spatial asymmetry, with particular focus on a recent break junction measurement. We argue that such asymmetry arises due to unequal coupling with the contacts and a consequent difference in charging effects, which can only be captured in a self-consistent model for molecular conduction. The direction of the asymmetry depends on the sign of the majority carriers in the molecule. For conduction through highest occupied molecular orbitals (i.e. HOMO or p-type conduction), the current is smaller for positive voltage on the stronger contact, while for conduction through lowest unoccupied molecular orbitals (i.e. LUMO or n-type conduction), the sense of the asymmetry is reversed. Within an extended Huckel description of the molecular chemistry and the contact microstructure (with two adjustable parameters, the position of the Fermi energy and the sulphur-gold bond length), an appropriate description of Poisson's equation, and a self-consistently coupled non-equilibrium Green's function (NEGF) description of transport, we achieve good agreement between theoretical and experimental I-V characteristics, both in shape as well as overall magnitude.Comment: length of the paper has been extended (4 pages to 6 pages), two new figures have been added (3 figures to 5 figures), has been accepted for PR

Creation and Reproduction of Model Cells with Semipermeable Membrane

A high activity of reactions can be confined in a model cell with a semipermeable membrane in the Schl\"ogl model. It is interpreted as a model of primitive metabolism in a cell. We study two generalized models to understand the creation of primitive cell systems conceptually from the view point of the nonlinear-nonequilibrium physics. In the first model, a single-cell system with a highly active state confined by a semipermeable membrane is spontaneously created from an inactive homogeneous state by a stochastic jump process. In the second model, many cell structures are reproduced from a single cell, and a multicellular system is created.Comment: 11 pages, 7 figure