First-principles study of rectification in bis-2-(5-ethynylthienyl)ethyne molecular junctions

Using density functional theory (DFT) combined with the first-principles nonequilibrium Green's function (NEGF), we investigated the electron-transport properties and rectifying behaviors of several molecular junctions based on the bis-2-(5-ethynylthienyl)ethyne (BETE) molecule. To examine the roles of different rectification factors, asymmetric electrode-molecule contacts and donor-acceptor substituent groups were introduced into the BETE-based molecular junction. The asymmetric current-voltage characteristics were obtained for the molecular junctions containing asymmetric contacts and donor-acceptor groups. In our models, the computed rectification ratios show that the mode of electrode-molecule contacts plays a crucial role in rectification and that the rectifying effect is not enhanced significantly by introducing the additional donor-acceptor components for the molecular rectifier with asymmetric electrode-molecule contacts. The current-voltage characteristics and rectifying behaviors are discussed in terms of transmission spectra, molecular projected self-consistent Hamiltonian (MPSH) states, and energy levels of MPSH states.     

Correlation effects in molecular conductors

The source-sink potential (SSP) model introduced previously [F. Goyer, M. Ernzerhof, and M. Zhuang, J. Chem. Phys. 126, 144104 (2007)] enables one to eliminate the semi-infinite contacts in molecular electronic devices (MEDs) in favor of complex potentials. SSP has originally been derived for independent electrons and extended to interacting two-electron systems subsequently [A. Goker, F. Goyer, and M. Ernzerhof, J. Chem. Phys. 129, 194901 (2008)]. Here we generalize SSP to N-electron systems and consider the impact of electron correlation on the transmission probability. In our correlated method for molecular conductors, the molecular part of the Hu?ckel Hamiltonian of the original SSP is replaced by the Hubbard Hamiltonian. For the contacts, however, the single-electron picture is retained and they are assumed to be spin polarized. Using our method, we study electron transmission in molecular wires, cross-conjugated chains, as well as aromatic systems. We find that, for realistic values of the electron-electron repulsion parameter, correlation effects modify the transmission probability quantitatively, the qualitative features remain. However, we find subtle new effects in correlated MEDs, such as Coulomb drag, that are absent in uncorrelated systems.     

Density functional theory of complex transition densities

We present an extension of Hohenberg-Kohn-Sham density functional theory to the domain of complex local potentials and complex electron densities. The approach is applicable to resonance (Siegert) [Phys. Rev. 56, 750 (1939)] states and other scattering and transport problems that can be described by a normalized state of a Hamiltonian containing a complex local potential. Such Hamiltonians are non-Hermitian and their eigenvalues are in general complex, the imaginary part being inversely proportional to the lifetime of the system. The one-to-one correspondence between complex local potentials nu and complex electron densities rho is established provided that the complex variables are sufficiently close to real local potentials and densities of nondegenerate ground states. We show that the exchange-correlation functionals, contributing to the complex energy, are determined through analytic continuation of their ground-state-theory counterparts. This implies that the exchange-correlation effects on the lifetime of a resonance are, under appropriate conditions, already determined by the functionals of the ground-state theory.     

A simple model of molecular electronic devices and its analytical solution

I present an analytically solvable model for molecular electronic devices (MEDs). The model contains all the essential components: semi-infinite contacts, molecule-contact interface, and of course, the molecule. I obtain explicit expression for the reflection coefficient r(E), as well as for the transmission probability T(E) [T(E)=1-mid R:r(E)mid R:(2)]. r(E) exhibits a surprisingly simple structure if studied as a function of a complex energy variable E. In this case, r(E) can be expressed in terms of a finite number of eigenvalues of a non-Hermitian Hamiltonian. This Hamiltonian also yields the molecular part of the MED wave function. Considering various MEDs, it is illustrated that the theory presented allows for a transparent interpretation of molecular conductance in terms of discrete eigenstates.     

A model electronic Hamiltonian to study low-lying electronic states of [Fe(bpy)3]2+ in aqueous solution

A simple model electronic Hamiltonian to describe the potential energy surfaces of several low-lying d-d states of the [Fe(bpy)(3)](2+) complex is developed for use in molecular dynamics (MD) simulation studies. On the basis of a method proposed previously for first-row transition metal ions in aqueous solution, the model Hamiltonian is constructed using density functional theory calculations for the lowest singlet and quintet states. MD simulations are then carried out for the two spin states in aqueous solution in order to examine the performance of the model Hamiltonian. The simulation results indicate that the present model electronic Hamiltonian reasonably describes the potential energy surfaces of the two spin states of the aqueous [Fe(bpy)(3)](2+) system, while retaining sufficient simplicity for application in simulation studies on excited state dynamics.     

PyCDFT: A Python package for constrained density functional theory

We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097 .     

Excited states of OsO4: A comprehensive time‐dependent relativistic density functional theory study

A large number of scalar as well as spinor excited states of OsO₄, in the experimentally accessible energy range of 3–11 eV, have been captured by time‐dependent relativistic density functional linear response theory based on an exact two‐component Hamiltonian resulting from the symmetrized elimination of the small component. The results are grossly in good agreement with those by the singles and doubles coupled‐cluster linear response theory in conjunction with relativistic effective core potentials. The simulated‐excitation spectrum is also in line with the available experiment. Furthermore, combined with detailed analysis of the excited states, the nature of the observed optical transitions is clearly elucidated. It is found that a few scalar states of ³T₁ and ³T₂ symmetries are split significantly by the spin‐orbit coupling. The possible source for the substantial spin‐orbit splittings of ligand molecular orbitals is carefully examined, leading to a new interpretation on the primary valence photoelectron ionization spectrum of OsO₄. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010