Small energy gap revealed in CrBr3 by scanning tunneling spectroscopy.

CrBr3 is a layered van der Waals material with magnetic ordering down to the 2D limit. For decades, based on optical measurements, it is believed that the energy gap of CrBr3 is in the range of 1.68-2.1 eV. However, controversial results have indicated that the band gap of CrBr3 is possibly smaller than that. An unambiguous determination of the energy gap is critical to the correct interpretations of the experimental results of CrBr3. Here, we present the scanning tunneling microscopy and spectroscopy (STM/S) results of CrBr3 thin and thick flakes exfoliated onto highly ordered pyrolytic graphite (HOPG) surfaces and density functional theory (DFT) calculations to reveal the small energy gap (peak-to-peak energy gap to be 0.57 ± 0.04 eV; or the onset signal energy gap to be 0.29 ± 0.05 eV from dI/dV spectra). Atomic resolution topography images show the defect-free crystal structure and the dI/dV spectra exhibit multiple peak features measured at 77 K. The conduction band - valence band peak pairs in the multi-peak dI/dV spectrum agree very well with all reported optical transitions. STM topography images of mono- and bi-layer CrBr3 flakes exhibit edge degradation due to short air exposure (∼15 min) during sample transfer. The unambiguously determined small energy gap settles the controversy and is the key in better understanding CrBr3 and similar materials.

High-temperature 2D ferromagnetism in conjugated microporous porphyrin-type polymers.

The need for magnetic 2D materials that are stable to the enviroment and have high Curie temperatures is very important for various electronic and spintronic applications. We have found that two-dimensional porphyrin-type aza-conjugated microporous polymer crystals are such a material (Fe-aza-CMPs). Fe-aza-CMPs are stable to CO, CO2, and O2 atmospheres and show unusual adsorption, electronic, and magnetic properties. Indeed, they are semiconductors with small energy band gaps ranging from 0.27 eV to 0.626 eV. CO, CO2, and O2 molecules can be attached in three different ways where single, double, or triple molecules are bound to iron atoms in Fe-aza-CMPs. For different attachment configurations we find that for CO and CO2 a uniform distribution of the molecules is most energetically favorable while for O2 molecules aggregation is most energetically preferable. The magnetic moments decrease from 4 to 2 to 0 for singly, doubly, triply occupied configurations for all gasses respectively. The most interesting magnetic properties are found for O2 molecules attached to the Fe-aza-CMP. For a single attachment configuration we find that an antiferromagnetic state is favorable. When two O2 molecules are attached, the calculations show the highest exchange integral with a value of J = 1071 µeV. This value has been verified by two independent methods where in the first method J is calculated by the energy difference between ferromagnetic and anitferromagnetic configurations. The second method is based on the frozen magnon approach where the magnon dispersion curve has been fitted by the Ising model. For the second method J has been estimated at J = 1100 µeV in excellent agreement with the first method.

Ferromagnetism in 2D organic iron hemoglobin crystals based on nitrogenated conjugated micropore materials.

In this work we study a low-cost two-dimensional ferromagnetic semiconductor with possible applications in biomedicine, solar cells, spintronics, and energy and hydrogen storage. From first principle calculations we describe the unique electronic, transport, optical, and magnetic properties of a π-conjugated micropore polymer (CMP) with three iron atoms placed in the middle of an isolated pore locally resembling heme complexes. This material exhibits strong Fe-localized dz2 bands. The bandgap is direct and equal to 0.28 eV. The valence band is doubly degenerate at the Γ-point and for larger k-wavevectors the HOMO band becomes flat with low contribution to charge mobility. The absorption coefficient is roughly isotropic. The conductivity is also isotropic with the nonzero contribution in the energy range 0.3-8 eV. The xy-component of the imaginary part of the dielectric tensor determines the magneto-optical Faraday and Kerr rotation. Nonvanishing rotation is observed in the interval of 0.5-5.0 eV. This material is found to be a ferromagnet of an Ising type with long-range exchange interactions with a very high magnetic moment per unit cell, m = 6 µB. The exchange integral is calculated by two independent methods: (a) from the energy difference between ferromagnetic and antiferromagnetic states and (b) from a magnon dispersion curve. In the former case Jnn = 27 µeV. In the latter case the magnon dispersion is fitted by the Ising model with the nearest and next-nearest neighbor spin interactions. From these estimations we find that Jnn = 19.5 µeV and Jnnn = -3 µeV. Despite the different nature of the calculations, the exchange integrals are only within 28% difference.

Room temperature d0 ferromagnetism in PbS films: nonuniform distribution of Pb vacancies.

Because of the importance of ferromagnetism at room temperature, we search for new materials that can exhibit a non-vanishing magnetic moment at room temperature and at the same time can be used in spintronics. The experimental results indicate that d0 ferromagnetism without any magnetic impurities takes place in PbS films made of close-packed lead sulfide nanoparticles of 30 nm. To explain the existence of the d0 ferromagnetism, we propose a model where various PbS bulk and surface configurations of Pb-vacancies are analyzed. The bulk configurations have a zero magnetic moment while the two surface configurations with Pb vacancies with the same non-vanishing magnetic moments and lowest ground state energies contribute to the total magnetization. Based on the experimental value of the saturation magnetization, 0.2 emu g-1, we have found that the calculated Pb vacancy concentration should be about 3.5%, which is close to typical experimental values. Besides being very important for applications, there is one feature of PbS d0 ferromagnetism that makes this material special for fundamental research: PbS ferromagnetism can exhibit topologically driven spatial magnetic moment distributions (e.g., magnetic skyrmions) due to large spin-orbit coupling.

Large enhancement in photocurrent by Mn doping in CdSe/ZTO quantum dot sensitized solar cells.

We find a large enhancement in the efficiency of CdSe quantum dot sensitized solar cells by doping with manganese. In the presence of Mn impurities in relatively small concentrations (2.3%) the photoelectric current increases by up to 190%. The average photocurrent enhancement is about 160%. This effect cannot be explained by a light absorption mechanism because the experimental and theoretical absorption spectra demonstrate that there is no change in the absorption coefficient in the presence of the Mn impurities. To explain such a large increase in the injection current we propose a tunneling mechanism of electron injection from the quantum dot LUMO state to the Zn2SnO4 (ZTO) semiconductor photoanode. The calculated enhancement is approximately equal to 150% which is very close to the experimental average value of 160%. The relative discrepancy between the calculated and experimentally measured ratios of the IPCE currents is only 6.25%. For other mechanisms (such as electron trapping, etc.) the remaining 6.25% cannot explain the large change in the experimental IPCE. Thus we have indirectly proved that electron tunneling is the major mechanism of photocurrent enhancement. This work proposes a new approach for a significant improvement in the efficiency of quantum dot sensitized solar cells.

Optical spectra of CdMnSe of nano-ferro- and antiferro-magnets.

We study optical transitions in CdSe quantum dots doped by Mn atoms. At low concentrations the transitions are spin-forbidden. Nevertheless, strong light absorption was experimentally found. To explain this effect we propose a new mechanism that includes two or more Mn atoms closely placed to each other containing the electrons with opposite spin projections. In this case the spin-flip is unnecessary. In addition we study absorption from quantum dots containing two Mn atoms with different multiplicities. We find that the strongest absorption from the gap is for an antiferromagnetic arrangement. The obtained results confirm the experimental concentration dependencies.

Spectroscopic and electronic structure properties of CdSe nanocrystals: spheres and cubes.

In this work we study the electronic structure of CdmSem quantum dots of various sizes and different shapes such as spheres and cubes using DFT, TDDFT, and CIS methods. This work requires a careful computational analysis where a proper exchange-correlation functional has to be chosen to fit the experimental optical gap. We find some differences in the optical and HOMO-LUMO gap values between spheres and cubes. In general, the gaps for the cubes have higher values than those for the spheres. We also calculate optical absorption spectra using the data for energy levels and oscillator strengths for optical transitions. We find that DFT yields some discrepancy in the density of states for the spheres and cubes. However, the density of states calculated by TDDFT and CIS provide better agreement. The results of the calculation can be useful for quantum dots synthesized in laser ablation experiments.

Tunneling current-controlled spin states in few-layer van der Waals magnets.

Effective control of magnetic phases in two-dimensional magnets would constitute crucial progress in spintronics, holding great potential for future computing technologies. Here, we report a new approach of leveraging tunneling current as a tool for controlling spin states in CrI3. We reveal that a tunneling current can deterministically switch between spin-parallel and spin-antiparallel states in few-layer CrI3, depending on the polarity and amplitude of the current. We propose a mechanism involving nonequilibrium spin accumulation in the graphene electrodes in contact with the CrI3 layers. We further demonstrate tunneling current-tunable stochastic switching between multiple spin states of the CrI3 tunnel devices, which goes beyond conventional bi-stable stochastic magnetic tunnel junctions and has not been documented in two-dimensional magnets. Our findings not only address the existing knowledge gap concerning the influence of tunneling currents in controlling the magnetism in two-dimensional magnets, but also unlock possibilities for energy-efficient probabilistic and neuromorphic computing.

Internal relaxation in dye sensitized solar cells based on Zn2SnO4 nanostructures.

In this work we study the effect of internal relaxation in a (Bu(4)N)(2)Ru(dcbpyH)(2)(NCS)(2) (N719) dye molecule in a dye sensitized solar cell. Experimentally measured light intensity dependencies of short circuit current and open circuit voltage for two different types of photoanodes, ZTO (Zn(2)SnO(4)) nanorods and nanoparticles, are explained in the framework of the proposed microscopic theory. This theory is based on a density matrix equation with a Markovian relaxation term. The computational results are in favor of the fast relaxation inside the unoccupied and occupied bands rather than slow interband electron-hole recombination. The difference in experimental dependencies for ZTO nanorods and nanoparticles is explained by the difference in the electron transfer matrix elements, and therefore, the electron transfer injection constants for the different morphologies of the photoanodes.

The theory of transport in helical spin-structure crystals.

We study helical structures in spin-spiral single crystals. In the continuum approach for the helicity potential energy the simple electronic band splits into two non-parabolic bands. For low exchange integrals, the lower band is described by a surface with a saddle shape in the direction of the helicity axis. Using the Boltzmann equation with the relaxation due to acoustic phonons, we discover the dependence of the current on the angle between the electric field and helicity axis leading to the both parallel and perpendicular to the electric field components in the electroconductivity. The latter can be interpreted as a planar Hall effect. In addition, we find that the transition rates depend on an electron spin allowing the transition between the bands. The electric conductivities exhibit nonlinear behaviors with respect to chemical potentialµ. We explain this effect as the interference of the band anisotropy, spin conservation, and interband transitions. The proposed theory with the spherical model in the effective mass approximation for conduction electrons can elucidate nonlinear dependencies that can be identified in experiments. We find the excellent agreement between the theoretical and experimental data for parallel resistivity depending on temperature at the phase transition from helical to ferromagnetic state in aMnPsingle crystal. In addition, we predict that the perpendicular resistivity abruptly drops to zero in the ferromagnetic phase.

The adsorption energy and diffusion of a pentacene molecule on a gold surface.

The nature of the chemical bonding of a pentacene molecule to a gold surface is studied. The calculations are carried out using two very different methodologies, the ab inito gaussian molecular orbital method and a numerical atomic orbital method, developed from the well tested SIESTA approach. Using the GAUSSIAN 09 package, we employ both local density B3LYP, and long-range correlated functionals CAM-B3LYP, ωB97, and ωB97X. For comparison, we also calculate the adsorption energy using the ATOMISTIX TOOLKIT with the revised PBE functional. Within computational and experimental errors we find that the best description of the binding energies can be obtained from GAUSSIAN calculations using long-range ωB97 and ωB97X exchange functionals. Thus the nature of chemical bonding of a pentacene to gold is a van der Waals type. To understand the large variation in the geometries computed by different methods, we calculate energy profiles in both X- and Y-directions. The energy barriers appear to be very small and comparable with the value of room temperature. Thus a pentacene molecule moves on a gold surface with almost no friction at room temperatures. An estimation of the work function is often obtained from a simple electrostatic approach. We test this estimation and find that this approach cannot be used because it significantly underestimates the work function. This investigation gives insights into the structure and bonding of pentacene to a gold surface and provides ideas for the improvement of methodologies for computing the properties of van der Waals adsorbates.

Spin filtering and spin separation in 2D materials by topological spin Hall effect.

The needs of high speed performance electronic devices for various applications require novel materials and new physical phenomena. For these purposes we propose to study new physical effects based on electron scattering on magnetic skyrmions and vortices distributed in a 2D ferromagnetic material. We show that the topological spin Hall effect can be efficiently employed for the filtering, switching, and separation of spin currents. For some values of the parameters (conduction electron concentrations, and skyrmion/vortex sizes) it is possible to separate Hall currents for different electron spin projections as it is like for different carrier charges (electrons and holes) in the normal Hall effect. The calculations are performed using the Boltzmann kinetic equation for the nonequilibrium distribution function and the Lippmann-Schwinger equation for the transition matrix in the whole range of the adiabaticity parameter. The spin filtering due to the skyrmion/vortex scattering can be several orders of magnitude more efficient in the narrow range of the electron concentrations than that of the ordinary ferromagnetic spin polarization in spintronics.

Weak d 0 ferromagnetism: Zn vacancy condensation in ZnS nanocrystals.

We provide the explanation of the large discrepancy of three orders of magnitude between the experimentally measured and theoretically calculated magnetic moments in ZnS nanocrystals. We assume that the condensation of Zn vacancies into a single droplet takes place. The energy calculations reveal that the droplet phase is more favorable than the uniformly distributed vacancy configuration. The other assumption made is that a small magnetic moment could arise at the interface between the ZnS crystal and vacancy cluster. The calculations however dismiss this hypothesis because the magnetization of the layered system also vanishes. Thus we suggest that the experimentally low magnetization values could be explained from one of the two following pictures: (a) there are two phases where the vacancy cluster with the zero magnetic moment coexists along with the other phase, in which there are uniformly distributed Zn vacancies with low concentrations or (b) there is only a single vacancy phase-a vacancy droplet being in the metastable state with a weak nonvanishing magnetic moment.

Coexistence of Two Electronic Nano-Phases on a CH3NH3PbI3-xClx Surface Observed in STM Measurements.

Scanning tunneling microscopy is utilized to investigate the local density of states of a CH3NH3PbI3-xClx perovskite in cross-sectional geometry. Two electronic phases, 10-20 nm in size, with different electronic properties inside the CH3NH3PbI3-xClx perovskite layer are observed by the dI/dV mapping and point spectra. A power law dependence of the dI/dV point spectra is revealed. In addition, the distinct electronic phases are found to have preferential orientations close to the normal direction of the film surface. Density functional theory calculations indicate that the observed electronic phases are associated with local deviation of I/Cl ratio, rather than different orientations of the electric dipole moments in the ferroelectric phases. By comparing the calculated results with experimental data we conclude that phase A (lower contrast in dI/dV mapping at -2.0 V bias) contains a lower I/Cl ratio than that in phase B (higher contrast in dI/dV).

Ab initio electron propagators in molecules with strong electron-phonon interaction: II. Electron Green's function.

Ab initio electron propagator methods are developed to study electronic properties of molecular systems with strong electron-electron and electron-phonon interactions. For the calculation of electron Green's functions we apply a canonical small polaron transformation that intrinsically contains strong electron-phonon effects. In the transformed Hamiltonian, the energy levels for the noninteracting particles are shifted down by the relaxation (solvation) energies. The Coulomb integrals are also renormalized by the electron-phonon interaction. For certain values of the electron-phonon coupling constants, the renormalized Coulomb integrals can be negative which implies the attraction between two electrons. Within the small polaron transformation we develop a diagrammatic technique for the calculation of electron Green's function in which the electron-phonon interaction is already included into the multiple phonon correlation functions. Since the decoupling of the phonon correlation functions is impossible, and therefore, a Wick's theorem for such correlation functions is invalid, there is no Dyson equation for the electron Green's function. To find the electron Green's function, we use different approximations. One of them is a link-cluster approximation that includes diagonal transitions for the renormalized zeroth Green's function. In the linked-cluster approach the Dyson equation is derived in the most general case, where the self-energy operator is an arbitrary functional (not only in the Hartree-Fock approximation). It is shown that even a Hartree-Fock electron (hole) is not a particle any longer. It is a quasiparticle with a finite lifetime that depends on energy of particle and hole states in different ways. As a consequence of this, a standard description of a Hartree-Fock approximation in terms of wave functions becomes inappropriate in this problem. To challenge the linked-cluster approximation we develop a different approach: a sequential propagation approximation where scattering events occur only for sequential transitions. A self-consistent Hartree-Fock equation for a four-index Green's function matrix is derived. In conclusion, the proposed schemes can be considered for future method developments for quantum chemical calculations for large molecules with strong nonadiabatic effects, e-e correlated electron transfer reactions, and electron transport in molecular transport junctions.

Ab initio electron propagators in molecules with strong electron-phonon interaction. I. Phonon averages.

Ab initio electron propagators in molecular systems with strong electron-electron and electron-phonon interactions are considered to study molecular electronic properties. This research is important in electron transfer reactions where the electron transition is not considered any longer as a single electron transfer process or in temperature dependences of current-voltage characteristics in molecular wires or aggregates. To calculate electron Green's functions, the authors apply a small polaron canonical transformation that intrinsically contains strong electron-phonon effects. According to this transformation, the excitation energies of the noninteracting Hamiltonian are shifted down by the relaxation (solvation) energy for each state. The electron-electron interaction is also renormalized by the electron-phonon coupling. For some values of the electron-phonon coupling constants, the renormalized Coulomb integrals can be negative resulting in the attraction between two electrons. Within this transformation, they develop a diagrammatic expansion for electron Green's function in which the electron-phonon interaction is included into the multiple phonon correlation functions. The multiple phonon correlation functions are exactly found. It is pointed out that Wick's theorem for such correlation functions is invalid. Consequently, there is no Dyson equation for electron Green's functions. The proposed approach can be considered for future method developments for quantum chemical calculations that include strong nonadiabatic (non-Born-Oppenheimer) effects.

Ab initio electron propagator calculations in molecular transport junctions: predictions of negative differential resistance.

In this work we study current-voltage characteristics in transport molecular junctions with a 1,4-benzene dithiol molecule as a bridge by using different ab initio electron propagator methods such as OVGF and P3 which are both programs in a Gaussian software package. The current-voltage characteristics are calculated for different values of Fermi energy in various basis sets such as 6-311++G(p,d) and cc-pVDZ and are compared with the experimental data. A good agreement is found in almost the entire voltage range. In addition, the results of our calculations indicate that the accuracy of ab initio electron propagator methods is in the range of 0.2-0.3 eV. Since the computational methods are truly ab initio, implying no adjustable parameters, functions, or functionals, the theoretical predictions can be improved only by changing the model of a transport device. The current-voltage characteristics predict peaks, i.e., negative differential resistances, for the various values of Fermi energy. As shown, the origin of the negative differential resistances is related to the voltage dependences of overlap integrals for the active terminal orbitals, expansion coefficients of partial atomic wavefunctions in Dyson orbitals, and the voltage dependences of Dyson poles (ionization potentials). We find that two peak behavior in the current-voltage characteristics can be explained by the anharmonicity of potential energy surfaces. As a result of our studies, we predict that negative differential resistances can be experimentally found by changing a position of Fermi level, i.e., by using different metal electrodes coated by a gold atomic monolayer.

Surface Green's function calculations: a nonrecursive scheme with an infinite number of principal layers.

A novel computational method for a surface Green's function matrix is introduced for the calculation of electrical current in molecular wires. The proposed nonrecursive approach includes an infinite number of principal layers and yields the second-order matrix equation for the transformed Green's function matrix. The solution is found by the direct diagonalization of the auxiliary matrix without any iteration process. As soon as complex roots of the auxiliary matrix (approximately GS) are calculated, the gaps and the bands in the surface electronic structure are found. It is shown that the solution of a second-order matrix equation determines the spectral density matrix, that is, the density of states for noninteracting electrons. Single and double principal layer models are studied both analytically and numerically. The energy interval for nonvanishing spectral matrices is determined. This method is applicable to matrices of any rank.

Electron tunneling dynamics in anharmonic bath.

Tunneling transition probability for a particle interacting with an anharmonic bath is found in a time-dependent Hartree approximation. The general expression is presented in terms of medium Keldysh functions that are assumed to be known. Furthermore, the transition probability is calculated in the noninteracting-blip approximation where the rate constant does not exhibit an activation dependence at high temperatures. The reorganization energy E(r) and the renormalized reaction heat epsilon are expressed in terms of the correlation matrix for a solvent and internal modes in both quantum and classical regimes. It is shown that E(r) and epsilon are temperature dependent.