Effect of the continuity of the pi conjugation on the conductance of ruthenium-octene-ruthenium molecular junctions

The conductance of a family of ruthenium-octene-ruthenium molecular junctions with different pi conjugation are investigated using a fully self-consistent ab initio approach which combines the nonequilibrium Green's function formalism with density functional theory. Our calculations demonstrate that the continuity of the pi conjugation in the contact region as well as along the molecular backbone affects the junction conductance significantly, showing the advantage of using the ruthenium-carbon double bond as the linkage of conjugated organic molecules.     

Implementation of molecular symmetry in valence bond calculation

A novel strategy for the construction of many-electron symmetry-adapted wave function is proposed for ab initio valence bond (VB) calculations and is implemented for valence bond self-consistent filed (VBSCF) and breathing orbital valence bond (BOVB) methods with various orbital optimization algorithms. Symmetry-adapted VB functions are constructed by the projection operator of symmetry group. The many-electron symmetry-adapted wave function is expressed in terms of symmetry-adapted VB functions, and thus the VB calculations can be performed with the molecular symmetry restriction. Test results show that molecular symmetry reduces the computational cost of both the iteration numbers and CPU time. Furthermore, excited states with specific symmetry can be conveniently obtained in VB calculations by using symmetry-adapted VB functions.     

Analysis on the contribution of molecular orbitals to the conductance of molecular electronic devices

We present a theoretical approach which allows one to extract the orbital contribution to the conductance of molecular electronic devices. This is achieved by calculating the scattering wave functions after the Hamiltonian matrix of the extended molecule is obtained from a self-consistent calculation that combines the nonequilibrium Green's function formalism with density functional theory employing a finite basis of local atomic orbitals. As an example, the contribution of molecular orbitals to the conductance of a model system consisting of a 4,4-bipyridine molecule connected to two semi-infinite gold monatomic chains is explored, illustrating the capability of our approach.     

Coexistence of 1,3-butadiene conformers in ionization energies and Dyson orbitals

The minimum-energy structures on the torsional potential-energy surface of 1,3-butadiene have been studied quantum mechanically using a range of models including ab initio Hartree-Fock and second-order M?ller-Plesset theories, outer valence Green's function, and density-functional theory with a hybrid functional and statistical average orbital potential model in order to understand the binding-energy (ionization energy) spectra and orbital cross sections observed by experiments. The unique full geometry optimization process locates the s-trans-1,3-butadiene as the global minimum structure and the s-gauche-1,3-butadiene as the local minimum structure. The latter possesses the dihedral angle of the central carbon bond of 32.81 degrees in agreement with the range of 30 degrees-41 degrees obtained by other theoretical models. Ionization energies in the outer valence space of the conformer pair have been obtained using Hartree-Fock, outer valence Green's function, and density-functional (statistical average orbital potentials) models, respectively. The Hartree-Fock results indicate that electron correlation (and orbital relaxation) effects become more significant towards the inner shell. The spectroscopic pole strengths calculated in the Green's function model are in the range of 0.85-0.91, suggesting that the independent particle picture is a good approximation in the present study. The binding energies from the density-functional (statisticaly averaged orbital potential) model are in good agreement with photoelectron spectroscopy, and the simulated Dyson orbitals in momentum space approximated by the density-functional orbitals using plane-wave impulse approximation agree well with those from experimental electron momentum spectroscopy. The coexistence of the conformer pair under the experimental conditions is supported by the approximated experimental binding-energy spectra due to the split conformer orbital energies, as well as the orbital momentum distributions of the mixed conformer pair observed in the orbital cross sections of electron momentum spectroscopy.     

Theory and range of modern semiempirical molecular orbital methods

Semiempirical molecular orbital methods have a long history. They serve to tackle large systems and complicated processes beyond the reach of ab initio or density functional methods. Although their setup is derived from Hartree–Fock theory, the design of approximate energy expressions and the empirical parameters are used to achieve higher accuracy than the underlying ab initio theory. In this way the effect of larger basis sets or correlation can be partially simulated. All widely used semiempirical methods establish their accuracy by error statistics for molecular properties with experimental and high-level ab initio or density functional theory calculations as a reference. Their computational efficiency makes them suitable for the study of biochemical systems and solid materials. The present review presents a variety of applications which demonstrate the need for and success of semiempirical methods.     

The ab initio model potential method and the optimized orbitals for the multiconfiguration self-consistent field-configuration interaction approach

A computational scheme is proposed for generating effective orbital sets for electron groups within the framework of the multiconfiguration self-consistent field-configuration interaction (MCSCF-CI) approach in the ab initio model potential method. © 1996 John Wiley & Sons, Inc.     

First-principles calculation on the conductance of a single 1,4-diisocyanatobenzene molecule with single-walled carbon nanotubes as the electrodes

The conductance of a single 1,4-diisocyanatobenzene molecule sandwiched between two single-walled carbon nanotube (SWCNT) electrodes are studied using a fully self-consistent ab initio approach which combines nonequilibrium Green's function formalism with density functional theory calculations. Several metallic zigzag and armchair SWCNTs with different diameters are used as electrodes; dangling bonds at their open ends are terminated with hydrogen atoms. Within the energy range of a few eV of the Fermi energy, all the SWCNT electrodes couple strongly only with the frontier molecular orbitals that are related to nonlocal pi bonds. Although the chirality of SWCNT electrodes has significant influences on this coupling and thus the molecular conductance, the diameter of electrodes, the distance, and the torsion angle between electrodes have only minor influences on the conductance, showing the advantage of using SWCNTs as the electrodes for molecular electronic devices.     

Inter-ring rotational potential of biphenyl and its effect on poly(p-phenylene)

Detailed ab initio studies have been done on the inter-ring torsional states of the biphenyl molecule using self-consistent field molecular orbital method. The potential goes through a minimum at an angle of 38°. The height of the potential barrier for the coplanar state is 2.01 kcal/mol. When the phenyl rings are perpendicular to each other, this height increases to 2.37 kcal/mol. The role of correlation and polarization is found to be important. The shape of the potential suggests that polyparaphenylene may possibly exist as a super helix.     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions

The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.     

A self-consistent transport model for molecular conduction based on extended Hückel theory with full three-dimensional electrostatics

We present a transport model for molecular conduction involving an extended Hückel theoretical treatment of the molecular chemistry combined with a nonequilibrium Green's function treatment of quantum transport. The self-consistent potential is approximated by CNDO (complete neglect of differential overlap) method and the electrostatic effects of metallic leads (bias and image charges) are included through a three-dimensional finite element method. This allows us to capture spatial details of the electrostatic potential profile, including effects of charging, screening, and complicated electrode configurations employing only a single adjustable parameter to locate the Fermi energy. As this model is based on semiempirical methods it is computationally inexpensive and flexible compared to ab initio models, yet at the same time it is able to capture salient qualitative features as well as several relevant quantitative details of transport. We apply our model to investigate recent experimental data on alkane dithiol molecules obtained in a nanopore setup. We also present a comparison study of single molecule transistors and identify electronic properties that control their performance.     

On the construction of diabatic and adiabatic potential energy surfaces based on ab initio valence bond theory

A theoretical model is presented for deriving effective diabatic states based on ab initio valence bond self-consistent field (VBSCF) theory by reducing the multiconfigurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction. Note that the VDC-MOVB method has been described previously. Employing the symmetric S(N)2 reaction between NH(3) and CH(3)NH(3)(+) as a test system, we found that the results from ab initio VBSCF and from ab initio MOVB calculations using the same basis set are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy.     

Validation of intermolecular transfer integral and bandwidth calculations for organic molecular materials

We present an interpretation of the intermolecular transfer integral that is independent from the origin of the energy scale allowing convergence studies of this important parameter of organic molecular materials. We present extensive numerical studies by using an ethylene pi dimer to investigate the dependence of transfer integrals on the level of theory and intermolecular packing. Transfer integrals obtained from semiempirical calculations differ substantially from one another and from ab initio results. The ab initio results are consistent across all the levels used including Hartree-Fock, outer valence Green's function, and various forms of density functional theory (DFT). Validation of transfer integrals and bandwidths is performed by comparing the calculated values with the experimental values of tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ), bis[1,2,5]thiadiazolo-p-quinobis(1,3-dithiole), (BTQBT) K-TCNQ, and hexagonal graphite. DFT in one of its presently popular forms, such as Perdew-Wang functionals (PW91), in combination with sufficient basis sets provides reliable transfer integrals, and therefore can serve as a basis for energy band calculations for soft organic materials with van der Waals gaps.     

Ab initio studies of solvent effect on molecular conformation and vibrational spectra of diacetamide

The conformational equilibria and vibrational spectra of diacetamide have been investigated by ab initio molecular orbital studies using the basis sets 6-31g(d,p) and 6-31++g(d,p) at Hartree-Fock and MP2 levels. The vibrational spectra of diacetamide have been satisfactorily interpreted taking into consideration the agreement between the calculated harmonic vibrational frequencies, infrared and Raman band intensities and shifts in deuterated molecules with those observed. The solvent effects were investigated by the self-consistent reaction field (SCRF) theory. The effect of solvent on the conformational equilibria and vibrational spectra is discussed. The calculated changes in the geometry and vibrational spectra on going from the gas phase to the solvent medium are in accord with the increasing weight of the dipolar resonance structure of the amide group in more polar solvents.     

Development of a molecular mechanics (MM2) force field for α-chlorosilanes

Molecular mechanics (MM2) parameters for silanes which have a Si-C-Cl fragment have been developed based on available experimental data and ab initio molecular orbital (MO) calculations. Molecular properties, mainly rotational barriers and geometries, of α-chlorosilanes have been studied using our new MM2 parameter set. Changes in the Si-C bond lengths and several bond angles of α-chlorosilanes due to the additional attachment of polar atom(s) have been investigated utilizing ab initio calculations. An electronegativity correction to both bond lengths and angles helps MM2 to reproduce results from ab initio calculations. The new force field has been applied to the conformational analysis of l-(chloromethyl)-1,2-dimethylsilacyclopentane, a model used in our studies of rearrangements of α-halosilanes.     

An empirical potential for the O-H·O Hydrogen bond: Part I. Comparison with ab initio molecular orbital results

An empirical potential energy function has been devised for the O-H·O hydrogen bond, for use with the MMI force field. The energy of the hydrogen bond is described as the sum of van der Waals, electrostatic and Morse components. The function has been used to calculate the potential energy hypersurface of the water dimer, and the results are compared with published ab initio molecular orbital studies. Satisfactory agreement is obtained except for orientations involving very short H·H contacts. The geometry and hydrogen bond energy of the equilibrium linear form of (H₂O)₂ are calculated to be r(O·O) = 2.84 Å, θ = 36°, ΔE = ?5.35 kcal mol^?1, which are close to the values obtained by experiment, and from molecular orbital calculations. The relative importance of the electrostatic component of the empirical hydrogen bond energy is consistent with molecular orbital energy decomposition studies. The empirical function has also been used to calculate the energy of the water trimer in orientations which serve as models for the crystallographic bifurcated hydrogen bond. The results indicate that, in these orientations, the trimer is typically 0–3 kcal mol^?1 more stable than the dimer, a result which is consistent with ab initio calculations.     

Functionality in single-molecule devices: model calculations and applications of the inelastic electron tunneling signal in molecular junctions

We analyze how functionality could be obtained within single-molecule devices by using a combination of non-equilibrium Green's functions and ab initio calculations to study the inelastic transport properties of single-molecule junctions. First, we apply a full non-equilibrium Green's function technique to a model system with electron-vibration coupling. We show that the features in the inelastic electron tunneling spectra (IETS) of the molecular junctions are virtually independent of the nature of the molecule-lead contacts. Since the contacts are not easily reproducible from one device to another, this is a very useful property. The IETS signal is much more robust versus modifications at the contacts and hence can be used to build functional nanodevices. Second, we consider a realistic model of a organic conjugated molecule. We use ab initio calculations to study how the vibronic properties of the molecule can be controlled by an external electric field which acts as a gate voltage. The control, through the gate voltage, of the vibron frequencies and (more importantly) of the electron-vibron coupling enables the construction of functionality: nonlinear amplification and/or switching is obtained from the IETS signal within a single-molecule device.     

Determination of molecular vibrational state energies using the ab initio semiclassical initial value representation: application to formaldehyde

We have demonstrated the use of ab initio molecular dynamics (AIMD) trajectories to compute the vibrational energy levels of molecular systems in the context of the semiclassical initial value representation (SC-IVR). A relatively low level of electronic structure theory (HF/3-21G) was used in this proof-of-principle study. Formaldehyde was used as a test case for the determination of accurate excited vibrational states. The AIMD-SC-IVR vibrational energies have been compared to those from curvilinear and rectilinear vibrational self-consistent field/vibrational configuration interaction with perturbation selected interactions-second-order perturbation theory (VSCF/VCIPSI-PT2) and correlation-corrected vibrational self-consistent field (cc-VSCF) methods. The survival amplitudes were obtained from selecting different reference wavefunctions using only a single set of molecular dynamics trajectories. We conclude that our approach is a further step in making the SC-IVR method a practical tool for first-principles quantum dynamics simulations.     

STOP: A slater-type orbital package for molecular electronic structure determination

A general ab initio package using Slater-type atomic orbitals is presented. This package, called STOP, uses the one-center two-range expansion method to evaluate the multicenter electronic integrals. Thoroughly optimized numerical techniques, in particular, convergence accelerators and suitable Gauss quadratures, are used in the algorithms which provide accurate numerical values for all these integrals. STOP thus provides wavefunctions for general molecular structures at the self-consistent field level for the first time over a Slater-type orbital basis. © 1996 John Wiley & Sons, Inc.