Single molecular electronics is considered as a promising alternative to conventional electronics towards the ultimate miniaturization of electronic devices. The achievements in this area are reviewed in a number of works [1] -[3] . A considerable attention is devoted to design and research of molecular diodes. Since the fundamental work of Aviram and Ratner [4] , plenty of objects were investigated on the basis of organic [5] -[8] and hybrid [9] [10] compounds which demonstrate diode properties. However, the transfer mechanism of electrons in such molecules remains obscure. In the model of Aviram and Ratner [4] , charge transfer in molecules containing electron donor (D) and acceptor (A) moieties separated by an insulating bridge (A-b-D molecules) was carried out in the three stages. Under forward bias, electrons can tunnel from one electrode to an acceptor (A), then to donor, and to the second electrode. At the opposite bias, it would require a larger potential to bring the energy levels of the donor and acceptor into favorable alignment with the electrodes. In a number of works [11] , it is proved that this model is realized. In 2001, the model of the molecular diode based on an asymmetric molecule, where the electroactive part was connected to electrodes through the isolating chains was offered [12] . The asymmetry could be attained by alkyl chains of the different lengths [13] or by various types of links with electrodes, for example at the expense of a chemisorption on the one end of a molecule and physical adsorption on the other end [14] [15] . According to R.S. Williams and all [13] , only one electronic level of electroactive part of a molecule participates in electrons transfer. The single-level model recently was used in a number of works for determination of the position of the molecular energy level and coupling to the electrodes. These parameters can be determined from a form of current-voltage dependencies (J(V)-curves) [16] [17] . Recently [18] a new type of molecular diode (organo-inorganic hybrid systems), based on the organoimido derivatives of hexamolybdates, has been designed. The transport of electrons was investigated using density functional theory combined with the non-equilibrium Green’s function. High and robust rectification ratio was predicted for these compounds.

The first STM investigation of organo-inorganic hybrid compounds is presented in this work. We study the molecular diode, where the electroactive moiety consists of Keggin’s heteropolyanion adsorbed on highly oriented pyrolytic graphite (HOPG) and connected to the other electrode by an alkyl chain. The surface of HOPG is one of the most often used surfaces in the experiments with STM [19] . In this work, the applicability of single-level model is considered. The coupling between CTMVA and a surface of HOPG, and energy gap between HOMO and LUMO orbitals are determined. Similar structures were studied in the work [20] by means of X-ray photoelectron spectroscopy (XPS), scanning tunnel electron microscopy and electrochemical methods. These compounds form ordered structures on the surface of HOPG similar to those observed in [21] -[23] after adsorption of the Keggin-type heteropolyacids. As shown in [20] , heteropolyanion is adsorbed on the surface of graphite, and an alkyl chain located above it. The method of film formation is described below.

1 ml of water solution (50 mg, concentration 25 mM) containing the excess of Na₂H₃[PMo₁₀V₂O₄₀] was added to 5 ml of cetyltrimethylammonium bromide [(С₁₆H₃₃)N(Me₃)]Br in methylene chloride solution (9 mg, concentration 5 mM). After vigorous stirring for 1 hour the orange-and-yellow precipitate was subsequently transferred to the Eppendorf tube and centrifuged for 2 minutes at 6000 cycles per second. The 20 µl of the resulting solution was dropped on the surface of HOPG (3 × 5 mM). After evaporation the surface was twice washed with deionized water (20 µl each portion). Then the sample was placed in vacuum for further measurements. The molecules of CTMVA aggregated on a freshly cleaved hydrophobic HOPG surface therefore the surface was previously hydroxylated. For this purpose the sample was introduced in solution of three percent hydrogen peroxide and irradiated for two minutes by 365 nm UV light. After the pretreatment the surface became hydrophilic with molecules of CTMVA adsorbed on HOPG with [PMo₁₀V₂O₄₀] anions positioned on the surface and hydrophobic alkyl chains above them. In the course of further study the sample was examined by means of scanning tunnel microscopy that revealed a spot with the adsorbed molecules forming a monolayer structured film. The topographical image coincided with the one displayed in the study [20] . All spectroscopic results were obtained on this film.

Measurements were carried out on STM VT “Omicron” at the room temperature in the ultrahigh-vacuum installation (5 × 10¹⁰ mbar). Electrochemically etched W tip of STM was cleaned by prolonged argon ions bombardment and current impulses at the voltage of 10 V. During the measurements the tip was brought into the contact with alkyl tails of the adsorbed molecules. Formation of bonds between the tip and alkyl tails of the adsorbed molecule was determined by measuring the dependencies of current on distance between a tip and a surface. The conductance decays exponentially with increasing length. Decay constant (~0.9 Å⁻¹) was in accordance with the data obtained for saturated σ-bonded alkanes. It proves that only the part of alkyl chain between the tip and the surface participates in electron transport. J(V)-dependencies were measured with the currents two order higher than currents prior to the loss of contact of the tip with molecules.

Figure 1 demonstrate J(V) curves measured on the area 10 × 10 nm². The curves are well reproducible. Rectification ratio RR = J(V+)/J(V−) reaches maximum value RR @ 15 ÷ 18 at the voltage 1.5 V. Good reproducibility most likely is connected with interaction of several alkyl chains with the tip. The part of J(V) curve for negative bias is presented on Figure 2, where the features inherent to all other curves can be seen. The noticeable increase of current occurs at the zero point and at voltage approaching to minus 2 V. Rapid and smooth growth of

tunnel current with increase of voltage can be observed at a positive bias on a sample.

Molecules under investigation belong to the asymmetric molecules which consist of electroactive moiety and an alkyl chain. Rectification properties of such molecules are actively discussed in a literature. HOMO and LUMO orbitals of electroactive part were used for an explanation of rectification, and it was shown that the location of the electroactive part of molecule within the gap between electrodes is very important [8] . Kornilovich et al. [13] and Baranger et al. [24] considered single-level model where transport of electrons was carried out through molecular orbitals energetically located near Fermi levels of adjacent electrodes and asymmetrically coupled to them. These orbitals could be LUMO [15] or HOMO [24] orbitals. In these works the main attention was paid to definition of parameters that provide the maximum values of rectification ratio. According to the work [13] RR had the maximum values for molecules where the electroactive part was separated from electrodes by the alkyl chains containing two atoms of carbon on one side and 10 atoms on the other side. The authors assumed that the potential drop across the electroactive part of molecular rectifier is insignificant. In the work [24] authors excluded a short chain from model to achieve the maximum current and obtain RR @ 10. They showed that the potential drop across the electroactive moiety is insignificant only for the resonant voltages where molecular level is in alignment with Fermi energy of a lead. Stadler et al. [25] made calculations for various types of molecular diodes, including those proposed by Aviram-Ratner, and concluded that RR can’t exceed 20. However, in the latest experiments [26] it was shown that RR reaches 100. The measurements were conducted on self-assembled monolayers of alkanethiolates with head group (SC11Fc). The authors made the conclusion that only one molecular orbital (the HOMO of the Fc) participated in electron transport. Such large rectification ratio is not predicted by theoretical models.

Thus, in these works the main attention was paid to a qualitative explanation of the mechanism of electron transport, but not to attempts of determination of parameters of model by comparison with experimental data. The single-level model allows to restore parameters of molecular orbitals participating in electron transfer (such as, relative position E₀ of the energy of molecular orbital with respect to the Fermi energy of metal and coupling Г_L,R of molecular states to left and right electrodes). This approach developed in works [16] [17] is based on application of Landauer formula for current through metal-molecule-metal junction:

where is the transmission function that depends on the energy E and voltage V, is the Fermi function, m is the lead chemical potential and kT is the thermal energy. Transmission function in the single-level model is given by the Breit-Wigner formula:

In the case of asymmetric coupling the position of energy level is a function of the applied voltage. It was assumed [16] that the dependence on Г_L,_R can be expressed as:

The main assumptions for the validity of this model are coherent tunneling and electron transport through a single molecular resonance in the entire voltage range. The validity of the first assumption seems reasonable for short molecules, the second assumption needs justification. The above-mentioned works [13] [24] [26] justify the applicability of single-level model for an explanation of rectification properties of asymmetric molecules. But our attempts to fit experimental J(V) curves to single-level model described by formulas (1)-(3) in the voltage range of ±2 V were unsuccessful. The reasons of inapplicability of single-level model can be connected with an electronic structure of electroactive moiety of molecules. In the case under consideration the electroactive part consists of Keggin’s anion which electronic structure is considered in work [27] . It is known that it consists from a group of free molecular orbitals mostly d_xy orbitals of metal atoms and group of filled orbitals relating to atoms of oxygen. The size of band gap is approximately 2 eV (for [PMo₁₂] − moeities LUMOHOMO gap = 2.03 eV). It is need to take into account that molecules of CTMVA are chemisorbed on hydroxylated hydrophilic graphite surface on one side and connected with the tip through physically adsorbed alkyl chain. Therefore, the electronic structure of Keggin’s anion is closely tied with Fermi level of HOPG. Consequently the Fermi level of the tip can be in alignment with a number of closely spaced LUMO levels at the positive bias. On the contrary electrons tunnel from surface to the tip through the nearest LUMO level of heteropolyanion located slightly above Fermi level of graphite at negative surface bias. For voltages comparable with the band gap, the contribution of HOMO orbital into electron transfer has to be taken into account. Therefore we carried out the fitting using formulas (1)-(3) for negative bias considering only frontier HOMO and LUMO orbitals. The model curve for only LUMO orbital (Figure 2(a)) does not describe the noticeable increase of a current for voltages approaching to minus 2 V. The model curve that takes into account the contribution of HOMO and LUMO orbitals presented on Figure 2(b). Experimental data are well described by this curve. The parameters describing coupling of CTMVA molecules with one electrode have higher values in comparison with parameters of coupling with another electrode.

The histogram of values of the HOMO-LUMO gap obtained from fitting of 54 J(V)-curves, presented on Figure 3. The measured value of the band gap is in a good agreement with DFT calculations for Keggin’s heteropolyacides [26] .

A new type of inorganic-organic hybrid system based on cetyltrimethylammonium derivatives of 10-molibdo-2- vanadophosphate-anions [PMo₁₀V₂O₄₀]⁻⁵ reveals the good rectification properties. But rectification ratio turned to be lower than RR predicted for organoimido derivatives of hexamolybdates [18] . This may be due to the replacement of a conjugate organic part to an organic part with alkyl chain.

This work was supported by the Russian Foundation for Basic Research (project №11-03-00494-a).