Properties of molecules in tunnel junctions

Molecular tunnel junctions involve studying the behaviour of a single molecule sandwiched between metal leads. When a molecule makes contact with electrodes, it becomes open to the environment which can heavily influence its properties, such as electronegativity and electron transport. While the most common computational approaches remain to be single particle approximations, in this thesis it is shown that a more explicit treatment of electron interactions can be required. By studying an open atomic chain junction, it is found that including electron correlations corrects the strong lead-molecule interaction seen by the ΔSCF approximation, and has an impact on junction I − V properties. The need for an accurate description of electronegativity is highlighted by studying a correlated model of hexatriene-di-thiol with a systematically varied correlation parameter and comparing the results to various electronic structure treatments. The results indicating an overestimation of the band gap and underestimation of charge transfer in the Hartree-Fock regime is equivalent to not treating electron-electron correlations. While in the opposite limit, over-compensating for electron-electron interaction leads to underestimated band gap and too high an electron current as seen in DFT/LDA treatment. It is emphasised in this thesis that correcting electronegativity is equivalent to maximising the overlap of the approximate density matrix to the exact reduced density matrix found at the exact many-body solution. In this work, the complex absorbing potential (CAP) formalism which allows for the inclusion metal electrodes into explicit wavefunction many-body formalisms is further developed. The CAP methodology is applied to study the electron state lifetimes and shifts as the junction is made open.

Role of Landau quantization on the neutron-drip transition in magnetar crusts

The role of a strong magnetic field on the neutron-drip transition in the crust of a magnetar is studied. The composition of the crust and the neutron-drip threshold are determined numerically for different magnetic field strengths using the experimental atomic mass measurements from the 2012 Atomic Mass Evaluation complemented with theoretical masses calculated from the Brussels-Montreal Hartree-Fock-Bogoliubov nuclear mass model HFB-24. The equilibrium nucleus at the neutron-drip point is found to be independent of the magnetic field strength. As demonstrated analytically, the neutron-drip density and pressure increase almost linearly with the magnetic field strength in the strongly quantizing regime for which electrons lie in the lowest Landau level. For weaker magnetic fields, the neutron-drip density exhibits typical quantum oscillations. In this case, the neutron-drip density can be either increased by about 14% or decreased by 25% depending on the magnetic field strength. These variations are shown to be almost universal, independently of the nuclear mass model employed. These results may have important implications for the physical interpretation of timing irregularities and quasiperiodic oscillations detected in soft gamma-ray repeaters and anomalous x-ray pulsars, as well as for the cooling of strongly magnetized neutron stars.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

Integer filling facfor phases and isospin in vertical diatomic artificial molecules

Integer filling factor phases of many-electron vertically coupled diatomic artificial quantum dot molecules are investigated for different values of the interdot coupling. The experimental results are analyzed within local-spin density functional theory for which we have determined a simple lateral confining potential law that can be scaled for the different coupling regimes, and Hartree-Fock theory. Maximum density droplets composed of electrons in both bonding and antibonding or just bonding states are revealed, and interesting isospin-flip physics appears for weak interdot coupling when the systematic depopulation of antibonding states leads to changes in isospin.

Electronic structure of a neutral oxygen vacancy in SrTiO3

The electronic structure of an isolated oxygen vacancy in SrTiO3 has been investigated with a variety of ab initio quantum mechanical approaches. In particular we compared pure density functional theory (DFT) approaches with the Hartree-Fock method, and with hybrid methods where the exchange term is treated in a mixed way. Both local cluster models and periodic calculations with large supercells containing up to 80 atoms have been performed. Both diamagnetic (singlet state) and paramagnetic (triplet state) solutions have been considered. We found that the formation of an O vacancy is accompanied by the transfer of two electrons to the 3d(z2) orbitals of the two Ti atoms along the Ti-Vac-Ti axis. The two electrons are spin coupled and the ground state is diamagnetic. New states associated with the defect center appear in the gap just below the conduction band edge. The formation energy computed with respect to an isolated oxygen atom in the triplet state is 9.4 eV.

Finite element method for the accurate solution of diatomic molecules

We present the Finite-Element-Method (FEM) in its application to quantum mechanical problems solving for diatomic molecules. Results for Hartree-Fock calculations of H_2 and Hartree-Fock-Slater calculations of molecules like N_2 and C0 have been obtained. The accuracy achieved with less then 5000 grid points for the total energies of these systems is 10_-8 a.u., which is demonstrated for N_2.

Crystal structure and theoretical calculations of Julocrotine, a natural product with antileishmanial activity

Julocrotine, N-(2,6-dioxo-l-phenethyl-piperidin-3-yl)-2-methyl-butyramide, is a potent antiproliferative agent against the promastigote and amastigote forms of Leishmania amazonensis (L.). In this work, the crystal structure of Julocrotine was solved by X-ray diffraction, and its geometrical parameters were compared with theoretical calculations at the B3LYP and HF level of theory. IR and NMR spectra also have been obtained and compared with theoretical calculations. IR absorptions calculated with the B3LYP level of theory employed together with the 6-311G+(d,p) basis set, are close to those observed experimentally. Theoretical NMR calculations show little deviation from experimental results. The results show that the theory is in accordance with the experimental data. (C) 2007 Wiley Periodicals, Inc.

Propriedades física de molécilas orgânicas e compostos do tipo perovskita CdSiO3 e CaPbO3

In the first part of this work our concern was to investigate the thermal effects in organic crystals using the theory of the polarons. To analyse such effect, we used the Fröhlich s Hamiltonian, that describes the dynamics of the polarons, using a treatment based on the quantum mechanics, to elucidate the electron-phonon interaction. Many are the forms to analyzing the polaronic phenomenon. However, the measure of the dielectric function can supply important information about the small polarons hopping process. Besides, the dielectric function measures the answer to an applied external electric field, and it is an important tool for the understanding of the many-body effects in the normal state of a polaronic system. We calculate the dielectric function and its dependence on temperature using the Hartree-Fock decoupling method. The dieletric function s dependence on the temperature is depicted by through a 3D graph. We also analyzed the so called Arrhenius resistivity, as a functionof the temperature, which is an important tool to characterize the conductivity of an organic molecule. In the second part we analyzed two perovskita type crystalline oxides, namely the cadmium silicate triclinic (CdSiO3) and the calcium plumbate orthorhombic (CaPbO3), respectively. These materials are normally denominated ABO3 and they have been especially investigated for displaying ferroelectric, piezoelectric, dielectrics, semiconductors and superconductors properties. We found our results through ab initio method within the functional density theory (DFT) in the GGA-PBE and LDA-CAPZ approximations. After the geometry optimization for the two structure using the in two approximations, we found the structure parameters and compared them with the experimental data. We still determined further the angles of connection for the two analyzed cases. Soon after the convergence of the energy, we determined their band structures, fundamental information to characterize the nature of the material, as well as their dielectrics functions, optical absorption, partial density of states and effective masses of electrons and holes

Estudo da série iso-eletrônica do átomo de hélio pelo método hiperesférico

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

GHOST POLES IN THE NUCLEON PROPAGATOR - VERTEX CORRECTIONS AND FORM-FACTORS

Vertex corrections are taken into account in the Schwinger-Dyson equation for the nucleon propagator in a relativistic field theory of fermions and mesons. The usual Hartree-Fock approximation for the nucleon propagator is known to produce the appearance of complex (ghost) poles which violate basic theorems of quantum field theory. In a theory with vector mesons there are vertex corrections that produce a strongly damped vertex function in the ultraviolet. One set of such corrections is known as the Sudakov form factor in quantum electrodynamics. When the Sudakov form factor generated by massive neutral vector mesons is included in the Hartree-Fock approximation to the Schwinger-Dyson equation for the nucleon propagator, the ghost poles disappear and consistency with basic requirements of quantum field theory is recovered.

Theoretical analysis on TiO(2)(110)/V surface

Theoretical analysis based on the Hartree-Fock method were performed in order to study the stoichiometric TiO(2) (110) surface and the vanadium substituted system. The Pople with polarization 3-21G* basis set level was used. The TiO(2) (110) surface was modeled using a (TiO(2))(15) cluster model. In order to take into account the finite size of the cluster, we have studied two different models: the point charge and the hydrogen saturated methodologies. The charge values used in the point charge calculations were optimized. The density of states, orbital self-consistend field (SCF) energies, and Mulliken charge values were analyzed. The method and model's dependence on the analyzed results are discussed. The theoretical results are compared with available experimental data. (C) 2001 John Wiley & Sons, Inc.

Ghost poles in the nucleon propagator in the linear sigma model approach and its role in pi N low-energy theorems

Complex mass poles, or ghost poles, are present in the Hartree-Fock solution of the Schwinger-Dyson equation for the nucleon propagator in renormalizable models with Yukawa-type meson-nucleon couplings, as shown many years ago by Brown, Puff and Wilets (BPW), These ghosts violate basic theorems of quantum field theory and their origin is related to the ultraviolet behavior of the model interactions, Recently, Krein et.al, proved that the ghosts disappear when vertex corrections are included in a self-consistent way, softening the interaction sufficiently in the ultraviolet region. In previous studies of pi N scattering using ''dressed'' nucleon propagator and bare vertices, did by Nutt and Wilets in the 70's (NW), it was found that if these poles are explicitly included, the value of the isospin-even amplitude A((+)) is satisfied within 20% at threshold. The absence of a theoretical explanation for the ghosts and the lack of chiral symmetry in these previous studies led us to re-investigate the subject using the approach of the linear sigma-model and study the interplay of low-energy theorems for pi N scattering and ghost poles. For bare interaction vertices we find that ghosts are present in this model as well and that the A((+)) value is badly described, As a first approach to remove these complex poles, we dress the vertices with phenomenological form factors and a reasonable agreement with experiment is achieved, In order to fix the two cutoff parameters, we use the A((+)) value for the chiral limit (m(pi) --> 0) and the experimental value of the isoscalar scattering length, Finally, we test our model by calculating the phase shifts for the S waves and we find a good agreement at threshold. (C) 1997 Elsevier B.V. B.V.

Investigação da transição isolante-metal do CDM via algoritmos computacionais e estudo de nano-dispositivos orgânicos

Apresentamos neste trabalho um estudo teórico sobre polímeros orgânicos conjugados. É conhecido que estes sistemas, em geral semicondutores ou isolantes, sob dopagem química podem vir a adquirir propriedades elétricas de material condutor. E ainda, sob ação de campo elétrico, pequenos oligômeros podem apresentar comportamento equivalente ao de dispositivos usuais, mas com inúmeras vantagens como, por exemplo, tamanho extremamente reduzido (alguns nanômetros). Dessa forma no primeiro capítulo faremos uma breve introdução sobre polímeros orgânicos conjugados mostrando alguns resultados experimentais obtidos para o polímero 4-dicianometileno-4,4-ciclopenta [2,1-b: 3,4b’] ditiofeno – CDM, que é o objeto central de estudo desta dissertação. O capítulo 2 trata dos métodos quânticos utilizados. Citaremos a Teoria de Hartre-Fock (HF) e suas derivações semi-empíricas. A técnica de Interação de configuração (CI) e a Teoria do Funcional da Densidade (DFT) também serão tratadas neste capítulo. O capítulo 3 é dedicado a descrever as características de alguns dispositivos usuais como diodos e transistores. Aqui o fundamental é entender a composição, o funcionamento e principalmente, como se comportam suas curvas características corrente versus voltagem (IxV). Citaremos ainda alguns dispositivos eletrônicos extremamente pequenos. No capítulo 4 começa nossos resultados e discussões referentes a análise da transição isolante-metal em CDM sob ação de dopagem. Primeiramente a nível semiempírico, obtivemos a caracterização ótica de oligômeros de CDM neutro e na presença de defeitos conformacionais do tipo bipólarons negativo e positivo. Partindo de geometrias otimizadas via métodos AM1 e PM3 obtivemos o espectro de absorção para sistemas com e sem carga. A nível Hartree-Fock calculamos a Estrutura de Bandas e a Densidade de Estados (DOS) para o PCDM no estado neutro e dopado. O cálculo da DOS e da Dispersão foram realizados através de programas computacionais desenvolvidos aqui no Grupo de Física de Materiais da Amazônia (GFMA). Apresentamos ainda neste capítulo o espectro de absorção teórico para oligômeros de CDM com diversas configurações com geometrias totalmente otimizadas pelo DFT. No capítulo 5 temos os resultados relativos à análise de nanodispositivos baseados em tetrâmeros de CDM com e sem carga. As curvas do deslocamento de carga versus voltagem apresentam características de curvas de dispositivos usuais. Analisamos também o espectro de absorção teórico dos nanodispositivos para valores de tensão nula e em pontos de saturação de corrente nas regiões direta e reversa.

Transporte eletrônico e quiralidade molecular: um estudo de dispositivos orgânicos em sistemas de dois terminais

No presente trabalho, investigamos o transporte eletrônico molecular em dois compostos orgânicos, o Ponceau SS (PSS) e o Oligo-(para)fenileno-vinileno (PPV) através de cálculos ab initio e função de Green de não equilíbrio (FGNE). Estes métodos demonstraram equivalência para a descrição destes dispositivos moleculares. Fizemos cálculos quânticos para o Hamiltoniano derivado de Hartree-Fock (HF) e obtivemos as propriedades de corrente-voltagem (I-V) para as duas estruturas moleculares. Com o método FGNE conseguimos modelar o transporte através de um sistema de multiníveis eletrônicos obtendo a corrente descrevendo as regiões de ressonância e a assimetria do sistema. Como resposta o PSS demonstrou assimetria para polarizações direta e reversa e a ressonância é alcançada mostrando que o dispositivo opere como um transistor molecular bi-direcional. Para o PPV investigamos também as propriedades geométricas através da conexão entre transporte eletrônico e o grau de quiralidade molecular que foi calculado usando o índice quiral que depende apenas das posições atômicas. Obtivemos que moléculas quirais e propriedades estruturais podem induzir uma assimetria no transporte eletrônico, resultando num processo de retificação. Também obtivemos que a resposta elétrica (I-V) e momento de dipolo elétrico são proporcionais ao grau de quiralidade molecular. Estes resultados sugerem que o transporte eletrônico neste sistema pode ser explorado na avaliação do seu grau de quiralidade.