Zeolites have been widely used in many different fields including catalysis, adsorption and separation, ion exchange, or gas storage. Conventional zeolites have three- dimensional (3D) structures with microporous channel system; typical pore sizes are well below 1 nanometer, therefore, diffusion limitation plays important role in many process and bulkier reactants (or products) cannot enter (or leave) the zeolite channel system. Two-dimensional (2D) zeolites prepared in last years can lift all diffusion limitation and they thus offer a very attractive alternative to conventional 3D zeolites. 2D zeolites attracted considerable attention on the experimental side; however, understanding of 2D zeolites based on computational investigation or on a combination of experimental and computational investigation is limited. A motivation for the computational work presented here is to improve our understanding of properties of 2D zeolites based on computational investigation. The originality of the research presented herein is in the strategy: we carried out systematic investigation of properties of corresponding 2D and 3D zeolites and we focus on the identification of similarities and differences. The most important zeolite properties, i.e., presence of Brønsted and Lewis acid sites, are investigated. A number of...

Nesta dissertação, estudamos cadeias unidimensionais antiferromagnéticas de spins 1/2 modeladas pelo Hamiltoniano de Heisenberg na presença de inomogeneidades causadas principalmente pela introdução de ligações substitucionais (defeitos nas ligações) e por efeitos de borda. Interessados então em determinar a energia do estado fundamental de sistemas com quaisquer distribuições das ligações, utilizamos o formalismo da Teoria do Funcional da Densidade (DFT) desenvolvido para o modelo de Heisenberg. O formalismo da DFT permite a estimativa da energia do estado fundamental de sistemas não-homogêneos conhecendo-se o sistema homogêneo. Construímos funcionais na aproximação da ligação local (LBA), proposta recentemente em analogia à já conhecida LSA (aproximação local para o spin). A obtenção dos funcionais se baseou no estudo do modelo de uma cadeia de spins em que as ligações são alternadas, isto é, a interação de troca se alterna em valor de sítio para sítio. Isso originou um funcional não-local na interação de troca da cadeia. Apesar disso, continuamos utilizando a nomenclatura LBA. Todos os resultados fornecidos pelos funcionais são comparados a dados provenientes de diagonalização numérica exata. / In this dissertation, we use the Heisenberg model to describe inhomogeneous antiferromagnetic spin 1/2 chains. The translational invariance is broken mainly due to the non-uniform distribution of bond interactions (defects) and the presence of boundaries. Interested in obtaining the ground-state energy of systems with any distribution of exchange couplings (Jij), we use the density-functional theory (DFT) formalism, developed for the Heisenberg model. The DFT formalism allows an estimate of the ground-state energy of inhomogeneous systems based on the homogeneous systems. We build functionals for the ground-state energy using a local bond approximation (LBA), recently proposed in analogy to the already known LSA (local spin approximation). To obtain the functionals we studied a model that describes an alternating chain, in which the exchange coupling alternates from site-to-site. This resulted in non-local functionals on the spin-spin exchange interaction. Nevertheless, we still call them LBA functionals. All the results from the functionals are compared with exact numerical data.

O efeito Hall quântico surge em gases de elétrons bidimensionais (2DEG) na presença de altos campos magnéticos B. O campo magnético quantiza o movimento planar dos elétrons em órbitas ciclotrônicas caracterizadas pelos níveis de Landau. Neste regime a resistividade transversal (ou Hall) ρxy em função de B exibe platôs em submúltiplos inteiros de e2/h, i.e., ρxy = ν-1 e2/h, sendo ν o fator de preenchimento dos níveis de Landau. Por sua vez, a resistividade longitudinal ρxx apresenta picos nas transições entre platôs de ρxy. Em primeira instância, ρxx é uma medida indireta da densidade de estados no nível de Fermi g(εF), e os picos dos mesmos indicam cruzamentos do nível de Fermi εF com niveis de Landau. Assim, o diagrama de densidade de elétrons n2D e B dos picos de ρxx ~ g(εF) fornece um mapa topológico da estrutura eletrônica do sistema. Em sistemas de duas subbandas, ρxx(n2D, B) exibe estruturas em forma de anel devido a cruzamentos de níveis de Landau de subbandas distintas [experimentos do grupo do Prof. Jiang (UCLA)]. Estes cruzamentos podem ainda levar a instabilidades ferromagnéticas. Investigamos estas instabilidades usando a teoria do funcional da densidade (DFT) para o cálculo da estrutura eletrônica, e o modelo de Ando (formalismo de Kubo) para o cálculo de ρxx e ρxy. Para temperaturas mais altas (340 mK) obtemos as estruturas em forma de anel em ρxx. Para temperaturas mais baixas (70 mK), observamos uma quebra dos anéis devido a transições de fase ferromagnéticas. Variando-se o ângulo θ de B com relação ao 2DEG observa-se o encolhimento do anel. Nossos resultados mostram que o ângulo de colapso total do anel depende de uma competição entre o termo de troca da interação de Coulomb (princípio de Pauli) e cruzamentos evitados devido ao ângulo θ finito. As transições de fase exibem ainda o fenômeno de histerese. Na região de instabilidade ferromagnética obtemos diferentes soluções variando B de forma crescente ou decrescente. Estas soluções possuem energias total diferentes, de forma que representam estados fundamental e excitado de muitos corpos. Esta observação, juntamente com resultados anteriores do grupo [Freire & Egues (2007)], representam as primeiras realizações teóricas da previsão da possibilidade de estados excitados como mínimos locais do funcional de energia do estado fundamental [Perdew & Levy (1985)]. O modelo aqui proposto fornece excelente acordo com os experimentos considerados. Adicionalmente, a observação sistemática e experimentalmente verificada dos estados excitados valida as previsões de Perdew & Levy. Aplicamos ainda estas mesmas ideias no cálculo da estrutura eletrônica e condutância de fios quânticos na presença de campos magnéticos, mostrando que cruzamentos de modos transversais também exibem instabilidades ferromagnéticas observadas em experimentos recentes [Dissertação de Mestrado de Filipe Sammarco, IFSC/USP], fortalecendo a validade do modelo apresentado nesta tese. / The quantum Hall effect arises in two dimensional electron gases (2DEG) under high magnetic fields B. The magnetic field quantizes the planar motion of the electrons into cyclotron orbits given by the Landau levels. In this regime the transversal (Hall) resistivity ρxy shows plateaus as a function of B at integer sub-multiples of e2/h, i.e., ρxy = ν-1 e2/h, where n is the filling factor of the Landau levels. The longitudinal resistivity ρxx shows peaks at the transition between the plateaus of ρxy. In principle, ρxx is an indirect measure of the density of states at the Fermi level g(εF), so that the peaks indicate when the Fermi level εF crosses a Landau level. Therefore, a density-B-field diagram n2D-B of the ρxx ~ g(εF) peaks shows a topological map of the electronic structure of the system. In two-subband systems, ρxx( n2D, B) shows ringlike structures due to crossings of spin-split Landau levels from distinct subbands [experiments from the group of Prof. Jiang (UCLA)] that could lead to ferromagnetic instabilities. We study these instabilities using the density functional theory (DFT) to calculate the electronic structure, and Ando\'s model (Kubo formalism) for ρxx and ρxy. At higher temperatures (340 mK) we also obtain the ringlike structures in ρxx. At lower temperatures (70 mK) we see broken rings due to quantum Hall ferromagnetic phase transitions. Tilting B by theta with respect to the 2DEG normal we find that the ring structure shrinks. Our results show that the angle of full collapse depends on a competition between the exchange term from the Coulomb interaction (Pauli principle) and the anticrossing of Landau levels due to the finite angle theta. Additionally, at the instabilities we observe hysteresis. Sweeping the B field up or down near these regions we obtain two different solutions with distinct total energies, corresponding to the ground state and an excited state of the many-body system. This result, together with previous results of our group [Freire & Egues (2007)], are the first realizations of the theoretical prediction of the possibility of excited states as local minima of the ground state energy functional [Perdew & Levy (1985)]. The model proposed here shows an excellent agreement with the experiments. Additionally, the systematic and experimentally verified observation of excited states corroborates the predictions of Perdew & Levy. Similar ideas as presented here when applied to the electronic structure and conductance of quantum wires with an in-plane magnetic field show ferromagnetic instabilities at crossings of the wire transverse modes [Master Thesis of Filipe Sammarco, IFSC/USP], also with excellent experimental agreement. This strengthen the range of validity of the model proposed in this Thesis.

Simulações de propriedades de estrutura eletrônica possuem fundamental importância para a física do estado sólido e química quântica. A teoria do funcional da densidade (DFT) é atualmente o método de estrutura eletrônica mais empregado, desde escalas atômicas e nanoscópicas até aglomerados biomoleculares. A acurácia da DFT depende essencialmente de aproximações para os efeitos de troca e correlação, para as quais existem vínculos a serem satisfeitos como forma de controlar sua construção. Esse é um tópico de grande importância, pois a construção de melhores funcionais é necessária para uma descrição cada vez mais precisa dos efeitos de muitos corpos na DFT. No presente trabalho, investigamos o comportamento da energia de troca e correlação e o desenvolvimento de funcionais aproximados sob a ótica de um vínculo universal de sistemas de interação Coulombiana, o limite inferior de Lieb-Oxford. Primeiramente apresentamos evidências de que em diversas classes de sistemas a energia de troca e correlação é distante do limite de Lieb-Oxford. A redução do limite foi implementada nos funcionais Perdew-Burke-Erzenhof (PBE), porém a forma com que o vínculo é implementado apenas aumentou a energia de troca. Propusemos em seguida que o limite de Lieb-Oxford não fosse utilizado apenas para determinar o valor de um parâmetro, como em PBE, mas que fosse ponto-de-partida de uma nova forma família de funcionais, do tipo hiper-GGA. Exploramos uma construção não-empírica, com implementação pós-autoconsistente. A particular forma proposta se beneficiou da redução do limite Lieb-Oxford, obtendo resultados satisfatórios para as energias de correlação. / Electronic-structure calculations play a fundamental role in solid-state physics and quantum chemistry. Density-functional theory (DFT) is today the most-widely used electronic-structure method, from atomic and nanoscopic scales to biomolecular aggregates. The accuracy of DFT depends essentially on approximations to the exchange and correlation energy, which are controlled by exact constraints. This is a very important issue, since the improvement of functionals is the key to a better description of many-body effects. In the present work, we investigate the exchange-correlation energy and approximate functionals from the viewpoint of an universal constraint on interacting Coulomb systems: the Lieb-Oxford lower bound. Initially we present evidence that for several classes of systems (atoms, ions, molecules and solids), the actual exchange-correlation energies are far from the Lieb-Oxford lower bound. A tighter form of this bound was conjectured; implemented in the Perdew-Burke-Erzenhof (PBE) functionals, and tested for atoms, molecules and solids. Finally, we propose to use the Lieb-Oxford bound not just to fix the value of a parameter as in PBE, but as a starting point for a new family of hyper-GGA functionals. For these, we explored a non-empirical construction, investigating its performance for atoms and small molecules post-selfconsistently. The particular HGGA proposed benefited from the tightening of the Lieb-Oxford bound and exhibited satisfactory correlation energies.

O uso de nanomateriais em aplicações como catálise e medicina, despertou nos últimos anos o interesse no estudo das propriedades de nanoclusters. O estudo das propriedades termodinâmicas desses sistemas é essencial, pois mudanças estruturais originadas de mudanças de fase podem alterar propriedades como atividade catalítica, momento magnético e propriedades óticas. A dinâmica molecular vêm sendo utilizada para o estudo computacional das propriedades termodinâmicas de diversos nanomateriais, enquanto o uso de métodos de Monte Carlo (MC), nesse contexto, tem se restringido ao estudo de nanoclusters de Lennard-Jones (LJ). Para avaliar a viabilidade do uso de métodos de MC no estudo de propriedades de sistemas reais, uma implementação do método Parallel Tempering Monte Carlo (PTMC) utilizando algoritmos do estado da arte para realização de trocas, determinação de temperaturas e ajuste de deslocamentos foi construída. Através de testes, é mostrado que alguns dos algoritmos implementados podem não ser adequados ao estudo do problema em questão. A implementação foi validada com o estudo das propriedades termodinâmicas de nanoclusters de LJ com 38, 55 e 147 átomos, que possuem resultados conhecidos na literatura. Além disso, resultados para propriedades do nanocluster LJ98 são apresentados, e devido à características estruturais desse sistema, é observada uma transição sólido-sólido entre as estruturas tetraédricas e icosaédricas em temperatura abaixo da de fusão. A hipótese do uso do algoritmo PTMC para o estudo de propriedades de materiais reais, foi testada nas nanoligas (PtCo)55 e (PtNi)55, descritas pelo potencial de Gupta. Através da comparação das estruturas de mais baixa energia com resultados de teoria do funcional da densidade (DFT, do inglês), é mostrado que o uso do potencial de Gupta pode ser justificado, dados os baixos desvios no comprimento de ligação (menores que 2.4%) e a semelhança de outras características estruturais. Os resultados obtidos indicam que o método PTMC é capaz de identificar as mudanças de fase das nanoligas estudadas. Essas mudanças são ilustradas e analisadas com o uso de um algoritmo para a comparação da similaridade de estruturas, com o qual foi possível analisar a fusão dos nanoclusters Co55, Ni55, Pt30Co25 e Pt40Ni15 (obtida em temperaturas entre 900 e 1100 K); além da fusão, a 727 K, e transição sólido-sólido, a 300 K, para a Pt55. Com as estruturas mais frequentes, encontradas através da análise de similaridade, e com a realização de cálculos DFT, foi possível observar um deslocamento do centro da banda d em direção ao HOMO (Highest Occupied Molecular Orbital) causado pelo aumento da temperatura. Esse deslocamento, segundo o modelo da banda d válido para superfícies, pode indicar uma maior reatividade dos nanoclusters nesses casos. / The use of nanomaterials in applications such as catalysis and medicine, aroused in the last years interest in studying properties of nanoclusters. The study of thermodynamic properties of these systems is essential, since structural changes originated from phase changes can alter properties such as catalytic activity, magnetic moment and optical properties. Molecular dynamics have been used for the computational study of thermodynamic properties of various nanomaterials, while the use of Monte Carlo methods (MC), in this context, has been restricted to the study of Lennard-Jones (LJ) nanoclusters. To evaluate the feasibility of using MC methods to study properties of real systems, an implementation of the Parallel Tempering Monte Carlo (PTMC) method using state of the art algorithms to perform exchanges, determine the temperature set and adjust the maximum displacement, was built. Through testing, it is shown that some of the implemented algorithms may not be suitable for the study of the problem in question. The implementation was validated by studying the thermodynamic properties of LJ nanoclusters with 38, 55 and 147 atoms, which have results known in the literature. In addition, results for the properties of the LJ98 nanocluster are reported, and due to the structural features of this system, a solid-solid transition between the tetrahedral and icosahedral structures in a temperature below melting is observed. The possibility of using the PTMC algorithm in the study of properties of real materials, is tested in the (PtCo)55 and (PtNi)55 nanoalloys, described by the Gupta potential. By comparing the lowest energy structures with density functional theory (DFT) results, it is shown that the use of the Gupta potential can be justified, given the small deviation in the bond lenght (less than 2.4%) and the similarity of other structural features. The results indicate that the PTMC method is able to identify the phase changes in the studied nanoalloys. These changes are illustrated and analyzed with the use of an algorithm for comparing the structure similarity, which made possible the analysis of the melting of the Co55, Ni55, Pt30Co25 e Pt40Ni15 nanoclusters (obtained at temperatures between 900 e 1100 K); and the melting at 727 K, and solid-solid transition at 300 K, for Pt55. With the most frequent structures, obtained by the similarity analysis, and through DFT calculations, it was possible to observe a shift in the d band center to the HOMO (Highest Occupied Molecular Orbital) caused by the temperature increase. This shift, following the d band model valid for surfaces, may indicate a higher reactivity of the nanocluster in these cases.

When doing computations on the negatively (positively) charged NV-center in diamond, the common procedure is to add (subtract) an electron from the system. However, when using periodic boundary conditions, this addition/subtraction of an electron from the supercell would result in a divergent electrostatic energy. So an artificial background jellium charge of opposite charge that compensate the electronic charge to make the supercell neutral is needed. This introduces further problems that needs corrections. And this method is especially problematic for slab supercells, as the compensating background charge leads to a dipole, which diverges as the vacuum between the slab images increases. An alternative, recently proposed way of charging the NV-center is to introduce electron donors/acceptors in the form of nitrogen/boron atoms (at substitutional sites in the diamond lattice). In this way, we keep the supercell/slab neutral, and avoid correction schemes. In this work we verify that the addition of a substitutional nitrogen atom indeed has the same effect on the NV-center as the more traditional method of adding an extra electron to the system. Further, we investigate the effects of 1. Adding two substitutional nitrogen atoms to the system (3 nitrogen atoms in total, neutral supercell), 2. Adding a substitutional nitrogen atom and an electron to the system (2 nitrogen atom in total, negatively charged supercell), 3. Adding two electrons to the system (1 nitrogen atom, doubly negatively charged supercell). Additionally, we investigate the addition of acceptor dopants (boron) in order to analyze the effect on the electronic structure of the NV-center and diamond.

No description available.

Methoden der Dichtefunktionaltheorie (DFT) sind beliebte quantenmechanische (QM) Verfahren. Einige Größen, wie z.B. Dispersionsenergien oder Reaktionsbarrieren, werden mit Standardfunktionalen jedoch schlecht beschrieben. Solche Probleme sind für kleine Systeme mit Elektronenkorrelationsmethoden lösbar. In dieser Arbeit wird ein Hybridansatz vorgestellt, der Vorzüge von DFT und Elektronenkorrelationsmethoden vereint. Dazu erfolgen DFT-Rechnungen für das Gesamtsystem und ein darin eingebettetes Modell des aktiven Zentrums. Dieses wird außerdem mit Moller–Plesset Störungstheorie zweiter Ordnung (MP2) zur Korrektur entsprechender DFT-Werte beschrieben. Die Zuverlässigkeit so erzielter Ergebnisse wird durch Extrapolation auf den Grenzfall einer vollständigen Orbitalbasis und eines unverkleinerten Modells erhöht. MP2-Korrekturen werden für eine Reihe unterschiedlich großer Modelle ermittelt. Daran wird ein C6-Paarpotential angepasst, das auch für das (periodische) Gesamtsystem die Ableitung von Werten in MP2-Qualität erlaubt. Zur Anwendung der Hybridmethode in der Katalyseforschung werden Protonensprünge im Zeolith Chabasit und Protonierungsreaktionen von Isobuten im Zeolith Ferrierit untersucht. Die MP2/DFT-Ergebnisse bestätigen, dass Aktivierungsbarrieren mit einfachen Dichtefunktionalen unterschätzt werden. Geschwindigkeitskonstanten sind nach Korrektur um ein bis zwei Größenordnungen kleiner. DFT-Adsorptions- und -Chemisorptionsenergien für Kohlenwasserstoffe in Zeolithen sind nicht verlässlich. Unterschiedlich große MP2-Korrekturen für die untersuchten Produkte liegen in einem Bereich von –29 bis –70 kJ/mol. Gegenüber dem wasserstoffgebundenen tert-Butylcarbeniumion profitieren Oberflächenalkoxide energetisch am meisten von der Erfassung der Dispersion. Die berechnete MP2-Adsorptionswärme des Isobutens im Zeolith Ferrierit (–74+–10 kJ/mol) entspricht experimentellen Anhaltspunkten und zeigt die Zuverlässigkeit des in dieser Arbeit vorgestellten Hybridansatzes. / Density functional theory (DFT) belongs to the most popular computational approaches in quantum mechanics (QM). For certain terms, e.g., dispersion energies or reaction barriers, widely used functionals do not yield reliable results. For small systems these problems can be overcome by using electron correlation methods. In this work a hybrid approach is presented combining advantages of both DFT and electron correlation methods. DFT calculations are performed for the full system and an embedded model representing the active site. The embedded model is also described by second-order Moller–Plesset perturbation theory (MP2) to obtain corrections to corresponding DFT values. The reliability of the results obtained with this hybrid method is further improved by extrapolation to the limiting case of a complete orbital basis set and an unreduced model system. MP2 corrections are obtained for a series of models of increasing size to fit a C6 pair potential. The fitted potential is applied to the full (periodic) system yielding results of MP2 quality. The hybrid method is applied in the field of catalysis research to investigate proton jump reactions in the zeolite chabazite and protonation reactions of isobutene in the zeolite ferrierite. The MP2/DFT results obtained confirm that with purely gradient corrected density functionals reaction barriers are clearly underestimated. Rate constants are corrected by one to two orders of magnitude. DFT reaction energies for adsorption and chemisorption of hydrocarbons in zeolites are not reliable. MP2 corrections obtained range between –29 and –70 kJ/mol and are not the same for different products. Compared to the hydrogen-bonded tert-butyl carbenium ion the surface alkoxides benefit most from treating dispersion explicitly. The calculated MP2 heat of adsorption of isobutene in ferrierite (–74+–10 kJ/mol) corresponds to experimental clues and underlines the reliability of the hybrid MP2/DFT approach presented in this work.

Partículas nanométricas de óxidos têm despertado um grande interesse devido a ampla gama de aplicações, por exemplo em catálise, óptica, sensores de gases, semicondutores, entre outros. Por isso, há vários estudos para uma grande variedade de óxidos com composição MO2, em função do tamanho e terminação da superfície. Também existe um grande interesse no estudo de óxidos mistos devido a possibilidade de combinar dois ou mais óxidos em escala nanométrica, com objetivo de controlar as propriedades físicas e químicas em função da composição. Neste trabalho, utilizou-se os metais Ce e Zr, uma vez que os óxidos desses elementos têm um grande potencial de aplicação na indústria como em células combustíveis, reações catalíticas entre outros. Os óxidos mistos Cem-nZrnO2mpossuem a capacidade de armazenar ou liberar oxigênio ao longo de reações químicas, e portanto, são largamente empregados em catalisadores de três-vias na indústria automobilística. Neste estudo foi empregada a teoria do funcional da densidade (DFT) com o funcional de troca-correlação proposto por Perdew, Burke e Ernzerhof (PBE). Investigou-se as propriedades estruturais, energéticas e eletrônicas dos clusters puros e então, selecionou-se quatro estruturas, as quais foram aplicadas na construção dos clusters de óxidos mistos de céria-zircônia. Por meio dos estudos de eletronegatividade de Mulliken, verificou-se que os clusters de zircônia apresentam uma tendência maior em atrair elétrons, quando comparados aos clusters de céria. Observou-se também que as estruturas de mais baixa energia de céria e zircônia apresentam diferentes números de coordenação (CN), no caso 4 e 3 NNN (vizinhos mais próximos), respectivamente. Isso se deve ao fato de tais estruturas não serem esféricas e também aos diferentes tamanhos atômicos de Ce e Zr. Além disso, notou-se que há uma maior compactação dos clusters de menor energia se comparados aos de maior energia, tanto para o céria quanto para zircônia, obteve-se os seguintes valores de diâmetro 1 nm e 2 nm para essa variação da energia relativa. Com esses resultados tem-se que há uma propensão dos clusters mistos, manterem propriedades semelhantes dos cluster puros, devido a similaridade entre o comportamento das estruturas de céria e da zircônia, os quais foram os agentes geradores dos clusters mistos. Observou-se que as propriedades estruturais dos cluster puros como CN, dav, σ e raio, mantiveram-se com a mesmas tendências nos clusters mistos. Para a obtenção do cluster mais estável analisou-se a energia de formação dos sistemas, sendo que a configuração mais favorável de ser formada é a Ce0,27Zr0,73O2, a qual considerou-se a estrutura pGMC representativa das configurações estudadas. Levando-se em consideração os resultados obtidos, na análise de formação dos clusters mistos delimitou-se a região de maior probabilidade de formação de clusters mistos estáveis entre 25% a 60%, de cério na composição. / Oxide nanomaterials have aroused great interest due to their application in several fields, such as in catalysis, optics, gas sensors, semiconductors, among others. Therefore, there are several studies for a wide variety of oxides with composition, and MO2, depending on the size and surface termination. There is also great interest in the study of mixed oxides due to the possibility of combining two or more oxides on a nanometric scale, in order to control their physical and chemical properties. In this work the metals Ce and Zr were used, since the oxides of these elements have a promising application in the industry as in fuel cells, catalytic reactions among others. The mixed oxides Cem-nZrnO2m present the ability to store or release oxygen during chemical reactions and therefore are widely used in three-way catalysts in the automotive industry. In this study, we used density functional theory (DFT) with PBE for the exchange-correlation functional. To investigate the structural, electronic and energetic properties of clusters, four selected structures were applied to build up the mixed ceria-zirconia oxides. Through the Mulliken electronegativity studies, it was found that zirconia clusters have a greater tendency to attract electrons, than the ceria clusters. Also observed that the lowest energy structures of ceria and zirconia have different coordination numbers, 4 and 3 NNN (Number of Nearest Neighbours), respectively. This is due to the fact that such structures are not spherical and have different atomic sizes. In addition, it was noted that there is a higher compaction of the clusters of lower energy compared to higher energy, for both ceria and zirconia, it was obtained the following values of diameter 1 nm and 2 nm for this variation of relative energy change. These results show that there is a trend in which mixed clusters have similar properties of the pure cluster, due to the similarity between the behavior of the structures of ceria and zirconia, which were the generating agents of the mixed clusters. It was observed that the structural properties of the pure clusters like CN, dav, σ and cluster radii, remained with the same trends in the mixed clusters. To obtain the most stable clusters we analyzed the formation energy of the systems, the most favorable configuration to be obtained is Ce0,27Zr0,73O2, which was considered the representative pGMC structure of the studied configurations. Taking into consideration the results, obtained in the analysis of the formation of the mixed clusters, delimited the higher probability of formation of mixed stable clusters between 25% to 60%, of cerium in composition.

Neste trabalho estudamos as propriedades estruturais, energéticas e eletrônicas das nanopartículas puras de prata (Ag) e de ouro (Au) e estruturas do tipo core-shell com número total de átomos variando de 147 à 923, no formato cubo-octaédrico. Estudamos também a adsorção da molécula de metanotiol (SCH4) sobre os sítios de coordenação dessas nanoestruturas, analisando, entre outros aspectos, os efeitos da interação de van de Waals. Para tanto, foram feitos cálculos teóricos de primeiros princípios dentro da Teoria do Funcional da Densidade (DFT) usando a Aproximação do Gradiente Generalizado (GGA) e Pseudopotenciais Ultrassuaves (USPP). Concluímos que as maiores nanopartículas puras e core-shell apresentam uma superfície mais esférica e suas energias de formação tendem às energias das superfícies [001] e [111] e dos bulks de Ag e de Au. Uma única camada de shell de ouro ou de prata na core-shell já determina praticamente o comportamento energético e as propriedades da nanopartícula. A inclusão da interação de van de Waals nos cálculos uniformiza, de certa forma, os padrões de deslocamento atômico das superfícies das nanopartículas e o comportamento energético das core-shell, sem entretanto alterar o perfil das densidades de estado. A adsorção da molécula de metanotiol nas nanopartículas puras de Ag e de Au e suas core-shell foi analisada verificando-se que ela praticamente não perturba os estados eletrônicos das nanopartículas e que sua estrutura molecular é preservada. Nas nanopartículas maiores verifica-se um único padrão de adsorção independente do número de camadas de shell nas estruturas core-shell. / In this work we study the structural, energetic and electronic properties of the pure nanoparticles of silver (Ag) and gold (Au) and their core-shell with total number of atoms ranging from 147 to 923 in cube-octahedral shape. We also investigated the adsorption of the methanethiol molecule (SCH4) in the coordination sites of these nanoparticles, analyzing, among other things, the influence of dispersion(van der Waals) interactions. Our simulations are performed using first principles theoretical calculations within of the Density Functional Theory (DFT) framework, described in terms of the Generalized Gradient Approximation (GGA), and by using Ultra-Soft Pseudo-potentials (USPP). We conclude that the largest pure nanoparticles and core-shell have a more spherical surface and their formation energies tend to formation energies of bare surfaces [001] and [111] and bulks of Ag and Au. A single layer of gold or silver shell already determines the properties and energetic behavior of the nanoparticles. The inclusion of van der Waals\' dispersion interaction in the calculations makes uniform, in certain way, the atomic displacement patterns surfaces of the nanoparticles and energetic behavior of core-shell, without change the form of the density of states. The adsorption of methanethiol molecule on the surface of the Ag and Au pure nanoparticles and their core-shell was analyzed and we verified that it almost does not disturb the electronic states of the nanoparticles and their molecular structure is preserved. In the largest nanoparticles we checking only one pattern of adsorption independent of the number of layers shell in the core-shell.