A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Moller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr_2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

/pdf/advances-in-molecular-quantum-chemistry-contained-in-the-q-3j7dwbtpjv.pdf

Advances in molecular quantum chemistry contained in the Q-Chem 4 program package

The Atomic Simulation Environment (ASE) is a software package written in the Python programming language with the aim of setting up, steering, and analyzing atomistic simula- tions. In ASE, tasks are fully scripted in Python. The powerful syntax of Python combined with the NumPy array library make it possible to perform very complex simulation tasks. For example, a sequence of calculations may be performed with the use of a simple "for-loop" construction. Calculations of energy, forces, stresses and other quantities are performed through interfaces to many external electronic structure codes or force fields using a uniform interface. On top of this calculator interface, ASE provides modules for performing many standard simulation tasks such as structure optimization, molecular dynamics, handling of constraints and performing nudged elastic band calculations.

/pdf/the-atomic-simulation-environment-a-python-library-for-48jb1pkf0e.pdf

The Atomic Simulation Environment - A Python library for working with atoms

The electronic ground state of a periodic system is usually described in terms of extended Bloch orbitals, but an alternative representation in terms of localized "Wannier functions" was introduced by Gregory Wannier in 1937. The connection between the Bloch and Wannier representations is realized by families of transformations in a continuous space of unitary matrices, carrying a large degree of arbitrariness. Since 1997, methods have been developed that allow one to iteratively transform the extended Bloch orbitals of a first-principles calculation into a unique set of maximally localized Wannier functions, accomplishing the solid-state equivalent of constructing localized molecular orbitals, or "Boys orbitals" as previously known from the chemistry literature. These developments are reviewed here, and a survey of the applications of these methods is presented. This latter includes a description of their use in analyzing the nature of chemical bonding, or as a local probe of phenomena related to electric polarization and orbital magnetization. Wannier interpolation schemes are also reviewed, by which quantities computed on a coarse reciprocal-space mesh can be used to interpolate onto much finer meshes at low cost, and applications in which Wannier functions are used as efficient basis functions are discussed. Finally the construction and use of Wannier functions outside the context of electronic-structure theory is presented, for cases that include phonon excitations, photonic crystals, and cold-atom optical lattices.

Maximally-localized Wannier Functions: Theory and Applications

Molecular conductance junctions are structures in which single molecules or small groups of molecules conduct electrical current between two electrodes. In such junctions, the connection between the molecule and the electrodes greatly affects the current-voltage characteristics. Despite several experimental and theoretical advances, including the understanding of simple systems, there is still limited correspondence between experimental and theoretical studies of these systems.

/pdf/electron-transport-in-molecular-wire-junctions-2m2v92jao5.pdf

Electron transport in molecular wire junctions.

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

RESEARCH ARTICLE Advances in molecular quantum chemistry contained in the Q-Chem 4 program package

We describe an ab initio method for calculating the electronic structure, electronic transport, and forces acting on the atoms, for atomic scale systems connected to semi-infinite electrodes and with an applied voltage bias. Our method is based on the density-functional theory (DFT) as implemented in the well tested SIESTA approach (which uses nonlocal norm-conserving pseudopotentials to describe the effect of the core electrons, and linear combination of finite-range numerical atomic orbitals to describe the valence states). We fully deal with the atomistic structure of the whole system, treating both the contact and the electrodes on the same footing. The effect of the finite bias (including self-consistency and the solution of the electrostatic problem) is taken into account using nonequilibrium Green's functions. We relate the nonequilibrium Green's function expressions to the more transparent scheme involving the scattering states. As an illustration, the method is applied to three systems where we are able to compare our results to earlier ab initio DFT calculations or experiments, and we point out differences between this method and existing schemes. The systems considered are: (i) single atom carbon wires connected to aluminum electrodes with extended or finite cross section, (ii) single atom gold wires, and finally (iii) large carbon nanotube systems with point defects.

/pdf/density-functional-method-for-nonequilibrium-electron-2m3kqra8o3.pdf

Density-functional method for nonequilibrium electron transport

Theoretical study of the nonlinear conductance of di-thiol benzene coupled to Au(111) surfaces via thiol and thiolate bonds

We report first-principles studies of electronic transport and rectification in molecular wires attached to gold electrodes. Our ab initio calculation gives an accurate description of the voltage drop as well as the broadening and alignment of the molecular levels in the metal-molecule-metal complex. We find that the operation range and rectification in such strongly chemisorbed molecules is limited by the width of the transmission resonances and their proximity to the Fermi level.

/pdf/theory-of-rectification-in-tour-wires-the-role-of-electrode-dq3f6e0vcx.pdf

Theory of rectification in tour wires: the role of electrode coupling.

We present state-of-the-art first principles calculations for the IV characteristics of a donor-insulator-acceptor (DsigmaA) type molecular diode anchored with thiolate bonds to two gold electrodes. We find very poor diode characteristics of the device, and the origin of this is analyzed in terms of the bias-dependent electronic structure. At zero bias, the highest occupied molecular orbital (HOMO) is confined to the D part, and the lowest unoccupied molecular orbital (LUMO) is confined to the A part, while at 3.8 V the two states align, and this gives rise to an increasing current. The latter is a potential mechanism for rectification and may in some cases lead to favorable diode characteristics. We identify the origin of the vanishing rectification for the investigated molecule, and on the basis of this we suggest parameters which are important for successful chemical engineering of DsigmaA rectifiers.

Do Aviram-Ratner diodes rectify?

Our recently developed method, TranSIESTA, enables modelling of molecular electronic devices under operation conditions. The method is based on density functional theory, and calculates the self-consistent electronic structure of a nanostructure coupled to three-dimensional electrodes with different electrochemical potentials. It uses a full atomistic ab initio description of both the electrodes and the nanoscale device. The calculations reveal information about the scattering states, transmission coefficients, electron current, and non-equilibrium forces in the systems. In this paper we use the method to investigate the electrical properties of three ring phenyl-ethynylene oligomers (OPE). We present results for the electrical effect of side groups and molecular conformations of the molecules. The calculations indicate that molecular switching and negative differential conductance (NDC) are related to rotations of the middle phenyl ring.

Jeremy Taylor

Papers

Density-functional method for nonequilibrium electron transport

Theoretical study of the nonlinear conductance of di-thiol benzene coupled to Au(111) surfaces via thiol and thiolate bonds

Theory of rectification in tour wires: the role of electrode coupling.

Do Aviram-Ratner diodes rectify?

TranSIESTA: a spice for molecular electronics