Switching mechanism of photochromic diarylethene derivatives molecular junctions

The electronic transport properties and switching mechanism of single photochromic diarylethene derivatives sandwiched between two gold surfaces with closed and open configurations are investigated by a fully self-consistent nonequilibrium Green's function method combined with density functional theory. The calculated transmission spectra of two configurations are strikingly distinctive. The open form lacks any significant transmission peak within a wide energy window, while the closed structure has two significant transmission peaks on both sides of the Fermi level. The electronic transport properties of the molecular junction with closed structure under a small bias voltage are mainly determined by the tail of the transmission peak contributed unusually by the perturbed lowest perturbed unoccupied molecular orbital. The calculated on-off ratio of currents between the closed and open configurations is about two orders of magnitude, which reproduces the essential features of the experimental measured results. Moreover, we find that the switching behavior within a wide bias voltage window is extremely robust to both substituting F or S for H or O and varying end anchoring atoms from S to Se and Te.     

Nonequilibrium electronic transport of 4,4'-bipyridine molecular junction

The electronic transport properties of a 4,4'-bipyridine molecule sandwiched between two Au(111) surfaces are studied with a fully self-consistent nonequilibrium Green's-function method combined with the density-functional theory. The 4,4'-bipyridine molecule prefers to adsorb near the hollow site of the Au(111) surface and distorts slightly. The modifications on the electronic structure of the molecule due to the presence of the electrodes are described by the renormalized molecular orbitals, which correspond well to the calculated transmission peaks. The average Fermi level lies close to the lowest unoccupied renormalized molecular orbital, which determines the electronic transport property of the molecular junction under a small bias voltage. The total transmission is contributed by a single channel. The transmission peaks shift with the applied bias voltage, and this behavior depends on the spatial distribution of the renormalized molecular orbitals and the voltage drop along the molecular junction. The shape of the calculated conductance curve of the equilibrium geometric configuration reproduces the main feature of the experimental results, but the value is larger than the measured data by about 6 times. Good agreement with the experimental measurements can be obtained by elongating the molecular junction. The electronic transport behaviors depend strongly on the interface configuration.     

The nature of singlet excitons in oligoacene molecular crystals

A theory for polarized absorption in crystalline oligoacenes is presented, which includes Frenkel exciton coupling, the coupling between Frenkel and charge-transfer (CT) excitons, and the coupling of all neutral and ionic excited states to the dominant ring-breathing vibrational mode. For tetracene, spectra calculated using all Frenkel couplings among the five lowest energy molecular singlet states predict a Davydov splitting (DS) of the lowest energy (0-0) vibronic band of only -32 cm(-1), far smaller than the measured value of 631 cm(-1) and of the wrong sign-a negative sign indicating that the polarizations of the lower and upper Davydov components are reversed from experiment. Inclusion of Frenkel-CT coupling dramatically improves the agreement with experiment, yielding a 0-0 DS of 601 cm(-1) and a nearly quantitative reproduction of the relative spectral intensities of the 0-n vibronic components. Our analysis also shows that CT mixing increases with the size of the oligoacenes. We discuss the implications of these results on exciton dissociation and transport.     

Electronic Transport Properties of Spin-Crossover Magnet Fe(Ⅱ)-N4S2 Complexes

Spin-crossover (SCO) magnets can act as one of the most possible building blocks in molecular spintronics due to their magnetic bistability between the high-spin (HS) and low-spin (LS) states. Here, the electronic structures and transport properties through SCO magnet Fe(Ⅱ)-N₄S₂ complexes sandwiched between gold electrodes are explored by performing extensive density functional theory calculations combined with non-equilibrium Green's function formalism. The optimized Fe-N and Fe-S distances and predicted magnetic moment of the SCO magnet Fe(Ⅱ)-N₄S₂ complexes agree well with the experimental results. The reversed spin transition between the HS and LS states can be realized by visible light irradiation according to the estimated SCO energy barriers. Based on the obtained transport results, we observe nearly perfect spin-filtering effect in this SCO magnet Fe(Ⅱ)-N₄S₂ junction with the HS state, and the corresponding current under small bias voltage is mainly contributed by the spin-down electrons, which is obviously larger than that of the LS case. Clearly, these theoretical findings suggest that SCO magnet Fe(Ⅱ)-N₄S₂ complexes hold potential applications in molecular spintronics.     

Huge Tunneling Magnetoresistance in Magnetic Tunnel Junction with Heusler Alloy Co2MnSi Electrodes

Magnetic tunnel junction with a large tunneling magnetoresistance has attracted great attention due to its importance in the spintronics applications. By performing extensive density functional theory calculations combined with the nonequilibrium Green's function method, we explore the spin-dependent transport properties of a magnetic tunnel junction, in which a non-polar SrTiO$_{3}$ barrier layer is sandwiched between two Heusler alloy Co$_{2}$MnSi electrodes. Theoretical results clearly reveal that the near perfect spin-filtering effect appears in the parallel magnetization configuration. The transmission coefficient in the parallel magnetization configuration at the Fermi level is several orders of magnitude larger than that in the antiparallel magnetization configuration, resulting in a huge tunneling magnetoresistance (i.e. $>10^6$), which originates from the coherent spin-polarized tunneling, due to the half-metallic nature of Co$_{2}$MnSi electrodes and the significant spin-polarization of the interfacial Ti 3d orbital.     

Single quintuple bond [PhCrCrPh] molecule as a possible molecular switch

The electronic transport properties of a single quintuple bond [PhCrCrPh] molecule sandwiched between two Au(111) surfaces with the trans-bent and linear configurations are studied by a fully self-consistent nonequilibrium Green's function method combined with density functional theory. The calculated transmission spectra of two chemical isomers are remarkably distinctive. Theoretical results suggest that the current through the trans-bent configuration is significantly larger than the corresponding linear one. The predicted on-off ratio of currents ranging from around 50 to 200 in the applied bias window [-1.5 V, 1.5 V] suggests that multiple bond compounds have attractive potential in molecular switch technology.     

Molecular tunnel junctions based on π-conjugated oligoacene thiols and dithiols between Ag, Au, and Pt contacts: effect of surface linking group and metal work function

The tunneling resistance and electronic structure of metal-molecule-metal junctions based on oligoacene (benzene, naphthalene, anthracene, and tetracene) thiol and dithiol molecules were measured and correlated using conducting probe atomic force microscopy (CP-AFM) in conjunction with ultraviolet photoelectron spectroscopy (UPS). Nanoscopic tunnel junctions (~10 nm(2)) were formed by contacting oligoacene self-assembled monolayers (SAMs) on flat Ag, Au, or Pt substrates with metalized AFM tips (Ag, Au, or Pt). The low bias (<0.2 V) junction resistance (R) increased exponentially with molecular length (s), i.e., R = R(0) exp(βs), where R(0) is the contact resistance and β is the tunneling attenuation factor. The R(0) values for oligoacene dithiols were 2 orders of magnitude less than those of oligoacene thiols. Likewise, the β value was 0.5 per ring (0.2 ?(-1)) for the dithiol series and 1.0 per ring (0.5 ?(-1)) for the monothiol series, demonstrating that β is not simply a characteristic of the molecular backbone but is strongly affected by the number of chemical (metal-S) contacts. R(0) decreased strongly as the contact work function (Φ) increased for both monothiol and dithiol junctions, whereas β was independent of Φ within error. This divergent behavior was explained in terms of the metal-S bond dipoles and the electronic structure of the junction; namely, β is independent of contact type because of weak Fermi level pinning (UPS revealed E(F) - E(HOMO) varied only weakly with Φ), but R(0) varies strongly with contact type because of the strong metal-S bond dipoles that are responsible for the Fermi level pinning. A previously published triple barrier model for molecular junctions was invoked to rationalize these results in which R(0) is determined by the contact barriers, which are proportional to the size of the interfacial bond dipoles, and β is determined by the bridge barrier, E(F) - E(HOMO). Current-voltage (I-V) characteristics obtained over a larger voltage range 0-1 V revealed a characteristic transition voltage V(trans) at which the current increased more sharply with voltage. V(trans) values were generally >0.5 V and were well correlated with the bridge barrier E(F) - E(HOMO). Overall, the combination of electronic structure determination by UPS with length- and work function-dependent transport measurements provides a remarkably comprehensive picture of tunneling transport in molecular junctions based on oligoacenes.     

Spin-dependent transport characteristics of Fe met-cars

Using non-equilibrium Green’s function and first-principles calculations we study structural, electronic, and transport properties of Fe₈C₁₂ met-car cluster sandwiched between two Au (1 0 0) electrodes. Several orientations were considered for the cluster attached to the gold surface and full structural optimization has been performed for the whole two-probe system. It was found a large current value for the present device and the molecular orientation plays an important role in the conducting behavior of the system. In energetically favorable case the I–V characteristic remains almost linear at low bias voltage (up to 1.5 V). This finding can be attributed to this fact that the transmission coefficient is almost flat around the gold Fermi level since the transmission is dominated by several broad molecular orbitals. We show that the electronic transmission is significantly spin-polarized while its size is large for the C atoms linkage. We also observe and discuss the NDR behavior of this novel molecular device in the range of 1.0–1.5 V for the energetically favorable configuration. The results are rationalized by analyzing the device transmission coefficient and density of states spectra.     

Au-MgB2-Au纳米结点的负微分电阻效应

运用密度泛函理论结合非平衡格林函数的方法对MgB2直线原子链与两半无限Au（100）电极构成纳米结点的电子输运特性进行了第一性原理计算．在模拟Au-MgB2-Au纳米结点的拉伸过程中,计算了结点在不同距离下的结合能与电导．结果发现结点的Au-B键长为1．90A,B-Mg键长为2．22A时,结合能最大,结构最稳定,此时结点平衡电导为0．51G0（G0=2e^0／h）．通过计算投影态密度发现电子通过结点时主要是通过B、Mg原子的px和py电子轨道形成的π键进行传输的．在-1．5～1．5V的电压范围内,结点的电流-电压近似为线性关系,表现出类似金属的导电性质,而当正负电压高于15V时,电流对称性逐渐减小,即存在负微分电阻效应,从不同电压下透射谱的变化对负微分电阻现象进行了分析与讨论．     

噻吩低聚物分子器件电输运特性的理论研究

Based on the first-principles computational method and the elastic scattering Green's func-tion theory, we have investigated the electronic transport properties of different oligothio-phene molecular junctions theoretically. The numerical results show that the difference of geometric symmetries of the oligothiophene molecules leads to the difference of the contact configurations between the molecule and the electrodes, which results in the difference of the coupling parameters between the molecules and electrodes as well as the delocalization properties of the molecular orbitals. Hence, the series of oligothiophene molecular junctions display unusual conductive properties on the length dependence.     

Theoretical investigation on electron transport through an organic molecule: effect of the contact structure

Knowing how the contact geometry influences the conductance of a molecular wire junction requires both a precise determination of the molecule/metallic-electrode interface structure and an evaluation of the conductance for different contact geometries with a fair accuracy. With a greatly improved method to solve the Lippmann-Schwinger equation, we are able to include at least one atomic layer of each electrode into the extended molecule. The artificial effect of the jellium model used for the electrodes is therefore significantly reduced. Our first-principles calculations on the transport properties of a single benzene dithiolate molecule sandwiched between Au(111) surfaces show that the transmission of the bridge site contact, which is the most stable adsorption configuration in equilibrium, displays different features from those of other configurations, and that the inclusion of the surface layers of Au electrodes into the extended molecule shifts and broadens the transmission peaks due to a stronger and more realistic S-Au bonding. We discuss the geometry dependence of the transport properties by analyzing the density of states of the molecular orbitals.     

Accurate calculation of transport properties for organic molecular semiconductors with spin-component scaled MP2 and modern density functional theory methods

At ambient temperatures, intermolecular hopping of charge carriers dominates the field effect mobility and thus the performance of organic molecular semiconductors for organic-based electronic devices. We have used a wide variety of modern and accurate computational methods to calculate the main parameters associated with charge transport, taking oligoacenes, and its derivatives as the exemplary organic materials. We tackle the problem from a combined inter- and intramolecular approach, in which the parameters are calculated for an isolated single molecule concomitantly with the stability of the dimers found in experimentally determined crystalline structures. Considering that most of the future applications within the field would need a full understanding of the transport mechanism, we assess the reliability of the methods to be employed according to the nature of the problem. Finally, we perform a computationally guided molecular engineering of a new set of materials derived from tetracene (rubrene and highly twisted oligoacenes) which allows to robustly anticipate the reasons for their expected performance in organic-based electronic devices.     

First-principles investigation of the asymmetric contact effect on current-voltage characteristics of a molecular device

The properties of electronic transport in an electronic device composed of a spatially symmetric phenyldithiolate molecule sandwiched between two gold electrodes with asymmetric contact are investigated by the first-principles study. It is found that the I-V and G-V characteristics of a device show significant asymmetry and the magnitudes of current and conductance depend remarkably on the variation of molecule-metal distance at one of the two contacts. Namely, an asymmetric contact would lead to the weak rectifying effects on the current-voltage characteristics of a molecular device. We also calculate self-consistently other microscopic quantities such as the local density of states, the total density of states, and the distribution of charges in the asymmetric molecular models under the applied bias. The results show that the highest-occupied molecular orbital (HOMO) is responsible for the resonant tunneling and the shifting of the HOMO due to the charging of the device under the bias voltage is the intrinsic origin of asymmetric I(G)-V characteristics.     

Conformations and charge transport characteristics of biphenyldithiol self-assembled-monolayer molecular electronic devices: a multiscale computational study

We report a computational study of conformations and charge transport characteristics of biphenyldithiol (BPDT) monolayers in the (sqrt.3 x sqrt.3)R30 degrees packing ratio sandwiched between Au(111) electrodes. From force-field molecular-dynamics and annealing simulations of BPDT self-assembled monolayers (SAMs) with up to 100 molecules on a Au(111) substrate, we identify an energetically favorable herringbone-type SAM packing configuration and a less-stable parallel packing configuration. Both SAMs are described by the (2sqrt.3 x sqrt.3)R30 degrees unit cell including two molecules. With subsequent density-functional theory calculations of one unit cell of the (i) herringbone SAM with the molecular tilt angle theta approximately 15 degrees , (ii) herringbone SAM with theta approximately 30 degrees , and (iii) parallel SAM with theta approximately 30 degrees, we confirm that the herringbone packing configuration is more stable than the parallel one but find that the energy variation with respect to the molecule tilting within the herringbone packing is very small. Next, by capping these SAMs with the top Au(111) electrode, we prepare three molecular electronic device models and calculate their coherent charge transport properties within the matrix Green's function approach. Current-voltage (I-V) curves are then obtained via the Landauer-Buttiker formula. We find that at low-bias voltages (|V| < or = 0.2 V) the I-V characteristics of models (ii) and (iii) are similar and the current in model (i) is smaller than that in (ii) and (iii). On the other hand, at higher-bias voltages (|V| > or 0.5 V), the I-V characteristics of the three models show noticeable differences due to different phenyl band structures. We thus conclude that the BPDT SAM I-V characteristics in the low-bias voltage region are mainly determined by the -Au [corrected] interaction within the individual molecule-electrode contact, while both intramolecular conformation and intermolecular interaction can affect the BPDT SAM I-V characteristics in the high-bias voltage region.     

Molecular vibrations of [n]oligoacenes (n=2-5 and 10) and phonon dispersion relations of polyacene

As model compounds for nanosize carbon clusters, the phonon dispersion curves of polyacene are constructed based on density functional theory calculations for [n]oligoacenes (n=2-5, 10, and 15). Complete vibrational assignments are given for the observed Fourier-transform infrared and Raman spectra of [n]oligoacenes (n=2-5). Raman intensity distributions by the 1064-nm excitation are well reproduced by the polarizability-approximation calculations for naphthalene and anthracene, whereas several bands of naphthacene and pentacene at 1700-1100 cm(-1) are calculated to be enhanced by the resonance Raman effect. It is found from vibronic calculations that the coupled a(g) modes between the Kekulé deformation and joint CC stretching give rise to the Raman enhancements of the Franck-Condon type, and that the b(3g) mode corresponding to the graphite G mode is enhanced by vibronic coupling between the (1)L(a)((1)B(1u)) and (1)B(b)((1)B(2u)) states. The phonon dispersion curves of polyacene provide a uniform foundation for understanding molecular vibrations of the oligoacenes in terms of the phase difference. The mode correlated with the defect-sensitive D mode of the bulk carbon networks is also found for the present one-dimensional system.     

A Photosensitive Liquid Crystal Studied by 14N NMR, 2H NMR,and DFT Calculations

An azobenzene derivative, namely diheptylazobenzene, showing the nematic and smectic A liquid crystalline phases, was investigated by means of a combined approach based on NMR and DFT calculations. ¹⁴N NMR quadrupole‐ and chemical‐shift‐perturbed spectra were acquired in the whole mesophasic range, providing both experimental quadrupolar splittings and chemical shift anisotropy values. On the same mesogen, deuterium labelled at the α‐position of the hydrocarbon chain, ²H NMR quadrupole‐perturbed spectra were recorded. The analysis of these NMR data was performed with the help of ab initio calculations, in vacuo and by taking into account the effect of the anisotropic environment typical of liquid crystals, by using the IEF‐PCM model. The geometry optimizations of the azomesogen in the trans and cis configurations were performed by DFT calculations employing the combination of B3LYP functional with the 6‐311G(d) basis set. The analysis of experimental NMR data was performed by considering the trans configuration as the most populated one and the corresponding quadrupolar tensors and chemical shielding tensors were determined at the DFT level of theory. The main result of this work is the determination of a relatively high and temperature‐dependent molecular biaxiality of the trans state of this azomesogen.     

A Modiˉed Molecular Structural Mechanics Method for Analysis of Carbon Nanotubes

A modified molecular structural mechanics method, based on molecular mechanics and similar to the finite element method, was developed. The energy of a system was expressed by the force field functions of the molecular mechanics. Under the small deformation assumption and by the principle of minimum potential energy, the system function was established. The properties of tension and bending of single-walled carbon nanotubes were analyzed. The Young's modulus is about 0.36 TPa nm, which agrees perfectly with the results of previous analysis by other researchers. It is found, for the first time, that the Young's moduli, for Zigzag nanotubes, are different from each other when the system energy was expressed as the sum of two or three individual energy terms in molecular mechanics. Whereas, the Young's moduli were the same for the Armchair nanotubes. It is found, when simulating the bending, that the deflections are closer to the theoretical ones, of the classical elasticity, when the diameter of the carbon nanotube increases.     

Assessing ionic transport during apple juice electro-acidification: influence on system efficiency

The purpose of this work was to assess ionic transport, during electro-acidification of apple juice, and its influence on the acidification rate and energy usage. In order to fulfill this objective, advanced current–voltage characterization of monopolar (cationic and anionic) and bipolar membranes was carried out. Experiments were conducted using two electrodialysis configurations (bipolar-cationic membranes and bipolar-anionic membranes) with two different systems (KCl–KCl and KCl–juice). For the bipolar-anionic configuration, the system HCl–juice was also considered. The characteristic values of the transmembrane potential and of the estimated membrane resistance were correlated to the ionic transport through the membranes, and to the energy usage of the systems. Influence of the membrane boundary layers on the transmembrane potential was also investigated by working at different feed flow rates. It was shown that for the present operating conditions, the boundary layers do not affect the transmembrane potential. Furthermore, use of HCl with the bipolar-anionic configuration enables to benefit from the advantages of each configuration: a small anionic membrane resistance due to the transport of H⁺, as for the bipolar-cationic configuration, and a high acidification rate of the apple juice due the neutralization of OH⁻ by HCl. It was for this system that the acidification rate was the fastest.     

槲皮素与桑色素ESI-MS裂解行为的比较分析

运用量子化学方法辅助解析并比较槲皮素与桑色素在电喷雾离子阱质谱(ESI-MS)负离子模式下的裂解行为。依据密度泛函理论(Density functional theory,DFT),在B3LYP/6-31G(d)水平,对槲皮素与桑色素的分子空间构型进行优化,确定稳定的几何构型与去质子化位点,在RB3LYP/6-31+G(2d,2p)水平,计算相对碎裂电压下的二级质谱中碎片离子处于稳定状态时的能量,通过比较准分子离子稳定构型并结合基组重叠误差(Basis set superposition error,BSSE)校正后的键解离能(Bond dissociation energy,BDE),推导了质谱碎裂过程。结果显示:槲皮素的稳定构型为A,B,C环处于同一平面,桑色素上的2'-OH使得B环与AC环之间翻转一定角度,二面角D(1,2,1',6')为-134.662 4°。槲皮素与桑色素的质谱裂解过程主要通过C环跨环裂解产生,且具有多种开裂方式,开裂先后顺序为:1,2开裂、0,2开裂、1,3开裂、1,4开裂与0,4开裂,分别生成碎片离子1,2A-,0,2A-,1,3A-,1,4A-与0,4A-,并逐步进行后续裂解,而2'-OH的存在促进了桑色素的裂解。该研究为进一步揭示黄酮醇类化合物的质谱裂解规律提供了理论依据。