Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Dynamical simulation of inelastic quantum transport

A method for correlated quantum electron-ion dynamics is combined with a method for electronic open boundaries to simulate in real time the heating, and eventual equilibration at an elevated vibrational energy, of a quantum ion under current flow in an atomic wire, together with the response of the current to the ionic heating. The method can also be used to extract inelastic current voltage corrections under steady-state conditions. However, in its present form the open-boundary method contains an approximation that limits the resolution of current-voltage features. The results of the simulations are tested against analytical results from scattering theory. Directions for the improvement of the method are summarized at the end.

Inelastic quantum transport in nanostructures: The self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared to a steady-state method based on nonequilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits an enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations where differences begin to emerge.

Crystal structure and ab initio calculations of a cyano-carbamimidic acid ethyl ester

The structure of one tautomer (amine form) of cyano-carbamimidic acid ethyl ester or (amino-ethoxy-methylidene)aminoformonitrile (CAS: 13947-84-7) was determined by single crystal X-ray diffraction. Ab initio quantum chemical calculations at the B3LYP, MP2 and G3 levels were performed to investigate the stability and the formation of the different tautomers and conformers. The calculations indicate that the amine form is the more stable tautomer, showing a high degree of election conjugation. The most stable amine conformer located by the calculations corresponds to the crystallized structure. On the contrary, in the less stable imine form, the conjugation is separated by a N2-C2 single bond. (C) 2007 Elsevier B.V. All rights reserved.

Quantum chemistry studies on electronic properties of CN adsorbed on silver electrodes

The electronic properties of CN adsorbed on Ag electrodes at different potentials have been studied by using the method of self-consistent-charge discrete variational Xa (SCC-DV-Xa) cluster calculations. It is shown that the binding of NCAg is dominated by both electrostatic and polarization effects derived from the charge of CN, and that the transfer of s charge from CN to silver is the largest donation contribution. The electrode potential influences this charge transfer process and the interaction between CN adsorbate and silver electrode. The results of quantum chemistry calculations fit well with the experimental data of in situ spectroscopic studies on the CN/Ag electrode systems. © 1991.

Efficient simulations with electronic open boundaries

We present a reformulation of the hairy-probe method for introducing electronic open boundaries that is appropriate for steady-state calculations involving nonorthogonal atomic basis sets. As a check on the correctness of the method we investigate a perfect atomic wire of Cu atoms and a perfect nonorthogonal chain of H atoms. For both atom chains we find that the conductance has a value of exactly one quantum unit and that this is rather insensitive to the strength of coupling of the probes to the system, provided values of the coupling are of the same order as the mean interlevel spacing of the system without probes. For the Cu atom chain we find in addition that away from the regions with probes attached, the potential in the wire is uniform, while within them it follows a predicted exponential variation with position. We then apply the method to an initial investigation of the suitability of graphene as a contact material for molecular electronics. We perform calculations on a carbon nanoribbon to determine the correct coupling strength of the probes to the graphene and obtain a conductance of about two quantum units corresponding to two bands crossing the Fermi surface. We then compute the current through a benzene molecule attached to two graphene contacts and find only a very weak current because of the disruption of the π conjugation by the covalent bond between the benzene and the graphene. In all cases we find that very strong or weak probe couplings suppress the current.

Bismuth embrittlement of copper is an atomic size effect

Embrittlement by the segregation of impurity elements to grain boundaries is one of a small number of phenomena that can lead to metallurgical failure by fast fracture(1). Here we settle a question that has been debated for over a hundred years(2): how can minute traces of bismuth in copper cause this ductile metal to fail in a brittle manner? Three hypotheses for Bi embrittlement of Cu exist: two assign an electronic effect to either a strengthening(3) or weakening(4) of bonds, the third postulates a simple atomic size effect(5). Here we report first principles quantum mechanical calculations that allow us to reject the electronic hypotheses, while supporting a size effect. We show that upon segregation to the grain boundary, the large Bi atoms weaken the interatomic bonding by pushing apart the Cu atoms at the interface. The resolution of the mechanism underlying grain boundary weakening should be relevant for all cases of embrittlement by oversize impurities.

Tools for analysing configuration interaction wavefunctions

The configuration interaction (CI) approach to quantum chemical calculations is a well-established means of calculating accurately the solution to the Schrodinger equation for many-electron systems. It represents the many-body electron wavefunction as a sum of spin-projected Slater determinants of orthogonal one-body spin-orbitals. The CI wavefunction becomes the exact solution of the Schrodinger equation as the length of the expansion becomes infinite, however, it is a difficult quantity to visualise and analyse for many-electron problems. We describe a method for efficiently calculating the spin-averaged one- and two-body reduced density matrices rho(psi)((r) over bar; (r) over bar' ) and Gamma(psi)((r) over bar (1), (r) over bar (2); (r) over bar'(1), (r) over bar'(2)) of an arbitrary CI wavefunction Psi. These low-dimensional functions are helpful tools for analysing many-body wavefunctions; we illustrate this for the case of the electron-electron cusp. From rho and Gamma one can calculate the matrix elements of any one- or two-body spin-free operator (O) over cap. For example, if (O) over cap is an applied electric field, this field can be included into the CI Hamiltonian and polarisation or gating effects may be studied for finite electron systems. (C) 2003 Elsevier B.V. All rights reserved.

On the origin of the forward peak and backward oscillations in the F + H2(v = 0) to HF(v' = 2) + H reaction

We present a semiclassical complex angular momentum (CAM) analysis of the forward scattering peak which occurs at a translational collision energy around 32 meV in the quantum mechanical calculations for the F + H2(v = 0, j = 0) ? HF(v' = 2, j' = 0) + H reaction on the Stark–Werner potential energy surface. The semiclassical CAM theory is modified to cover the forward and backward scattering angles. The peak is shown to result from constructive/destructive interference of the two Regge states associated with two resonances, one in the transition state region and the other in the exit channel van der Waals well. In addition, we demonstrate that the oscillations in the energy dependence of the backward differential cross section are caused by the interference between the direct backward scattering and the decay of the two resonance complexes returning to the backward direction after one full rotation.

Metastable dissociation of anions formed by electron attachment

A comprehensive analysis of metastable dissociation of 2, 4-dinitrotoluene (DNT) parent anions formed by attachment of electrons of controlled energy is presented. We characterize the energy dependence and kinetic energy release of the reaction which competes with autodetachment. A surprising finding is a highly exothermic metastable reaction triggered by the attachment of thermal electrons which we relate to the well-known electrostatic ignition hazards of DNT and other explosives. Quantum chemical calculations are performed for dinitrobenzene in order to elucidate the process of NO abstraction.

Mutual Solubilities of Water and Hydrophobic Ionic Liquids

The ionic nature of ionic liquids (ILs) results in a unique combination of intrinsic properties that produces increasing interest in the research of these fluids as environmentally friendly "neoteric" solvents. One of the main research fields is their exploitation as solvents for liquid-liquid extractions, but although ILs cannot vaporize leading to air pollution, they present non-negligible miscibility with water that may be the cause of some environmental aquatic risks. It is thus important to know the mutual solubilities between ILs and water before their industrial applications. In this work, the mutual solubilities of hydrophobic yet hygroscopic imidazolium-, pyridinium-, pyrrolidinium-, and piperidinium-based ILs in combination with the anions bis(trifluoromethylsulfonyl)imide, hexafluorophosphate, and tricyanomethane with water were measured between 288.15 and 318.15 K. The effect of the ILs structural combinations, as well as the influence of several factors, namely cation side alkyl chain length, the number of cation substitutions, the cation family, and the anion identity in these mutual solubilities are analyzed and discussed. The hydrophobicity of the anions increases in the order [C(CN)3] <[PF6] <[Tf2N] while the hydrophobicity of the cations increases from [Cnmim] <[Cnmpy] [Cnmpyr] <[Cnmpip] and with the alkyl chain length increase. From experimental measurements of the temperature dependence of ionic liquid solubilities in water, the thermodynamic molar functions of solution, such as Gibbs energy, enthalpy, and entropy at infinite dilution were determined, showing that the solubility of these ILs in water is entropically driven and that the anion solvation at the IL-rich phase controls their solubilities in water. The COSMO-RS, a predictive method based on unimolecular quantum chemistry calculations, was also evaluated for the description of the water-IL binary systems studied, where it showed to be capable of providing an acceptable qualitative agreement with the experimental data.

Mutual Solubilities of Water and the [Cnmim][Tf2N] Hydrophobic Ionic Liquids

Ionic liquids (ILs) have recently garnered increased attention because of their potential environmental benefits as "green" replacements over conventional volatile organic solvents. While ILs cannot significantly volatilize and contribute to air pollution, even the most hydrophobic ones present some miscibility with water posing environmental risks to the aquatic ecosystems. Thus, the knowledge of ILs toxicity and their water solubility must be assessed before an accurate judgment of their environmental benefits and prior to their industrial applications. In this work, the mutual solubilities for [C2-C8mim][Tf2N] (n-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) and water between 288.15 and 318.15 K at atmospheric pressure were measured. Although these are among the most hydrophobic ionic liquids known, the solubility of water in these compounds is surprisingly large, ranging from 0.17 to 0.36 in mole fraction, while the solubility of these ILs in water is much lower ranging from 3.2 × 10-5 to 1.1 × 10-3 in mole fraction, in the temperature and pressure conditions studied. From the experimental data, the molar thermodynamic functions of solution and solvation such as Gibbs energy, enthalpy, and entropy at infinite dilution were estimated, showing that the solubility of these ILs in water is entropically driven. The predictive capability of COSMO-RS, a model based on unimolecular quantum chemistry calculations, was evaluated for the description of the binary systems investigated providing an acceptable agreement between the model predictions and the experimental data both with the temperature dependence and with the ILs structural variations.

Nitrile-Functionalized Pyridinium, Pyrrolidinium, and Piperidinium Ionic Liquids

Two series of 1-alkylpyridinium and N-alkyl-N-methylpiperidinium ionic liquids fiinctionalized with a nitrile group at the end of the alkyl chain have been synthesized. Structural modifications include a change of the alkyl spacer length between the nitrile group and the heterocycle of the cationic core, as well as adding methyl or ethyl substituents on different positions of the pyridinium ring. The anions are the bromide and the bis(trifluoromethylsulfonyl)imide ion. All the bis(trifluoromethylsulfonyl)imide salts as well as the bromide salts with a long alkyl spacer were obtained as viscous liquids at room temperature, but some turned out to be supercooled liquids. In addition, pyrrolidinium and piperidinium ionic liquids with two nitrile functions attached to the heterocyclic core have been prepared. The crystal structures of seven pyridinium bis(trifluoromethylsulfonyl)imide salts are reported. Quantum chemical calculations have been performed on model cations and ion pairs with the bis(trifluoromethylsulfonyl)imide anion. A continuum model has been used to take solvation effects into account. These calculations show that the natural partial charge on the nitrogen atom of the nitrile group becomes more negative when the length of the alkyl spacer between the nitrile functional group and the heterocyclic core of the cation is increased. Methyl or methoxy substituents on the pyridinium ring slightly increase the negative charge on the nitrile nitrogen atom due to their electron-donating abilities. The position of the substituent (ortho, meta, or para) has only a very minor effect on the charge of the nitrogen atom. The N-15 NMR spectra of the bis(trifluoromethylsulfonyl)imide ionic liquids were recorded with the nitrogen-15 nucleus at its natural abundance. The chemical shift of the N-15 nucleus of the nitrile nitrogen atom could be correlated with the calculated negative partial charge on the nitrogen atom.

Dispersion interactions in room-temperature ionic liquids: Results from a non-empirical density functional

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Roman-Perez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by approximate to 7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Popolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields. (C) 2011 American Institute of Physics. [doi:10.1063/1.3652897]

Stellar Surface Magneto-Convection as a Source of Astrophysical Noise. I. Multi-component Parameterisation of Absorption Line Profiles

We outline our techniques to characterise photospheric granulation as an astrophysical noise source. A four component parameterisation of granulation is developed that can be used to reconstruct stellar line asymmetries and radial velocity shifts due to photospheric convective motions. The four components are made up of absorption line profiles calculated for granules, magnetic intergranular lanes, non-magnetic intergranular lanes, and magnetic bright points at disc centre. These components are constructed by averaging Fe I $6302 \mathrm{\AA}$ magnetically sensitive absorption line profiles output from detailed radiative transport calculations of the solar photosphere. Each of the four categories adopted are based on magnetic field and continuum intensity limits determined from examining three-dimensional magnetohydrodynamic simulations with an average magnetic flux of $200 \mathrm{G}$. Using these four component line profiles we accurately reconstruct granulation profiles, produced from modelling 12 x 12 Mm$^2$ areas on the solar surface, to within $\sim \pm$ 20 cm s$^{-1}$ on a $\sim$ 100 m s$^{-1}$ granulation signal. We have also successfully reconstructed granulation profiles from a $50 \mathrm{G}$ simulation using the parameterised line profiles from the $200 \mathrm{G}$ average magnetic field simulation. This test demonstrates applicability of the characterisation to a range of magnetic stellar activity levels.