CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments, we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by screened hybrid functional calculations.

Ionization levels of doped copper indium sulfide chalcopyrites

The electronic structure of modified chalcopyrite CuInS2 has been analyzed from first principles within the density functional theory. The host chalcopyrite has been modified by introducing atomic impurities M at substitutional sites in the lattice host with M = C, Si, Ge, Sn, Ti, V, Cr, Fe, Co, Ni, Rh, and Ir. Both substitutions M for In and M for Cu have been analyzed. The gap and ionization energies are obtained as a function of the M-S displacements. It is interesting for both spintronic and optoelectronic applications because it can provide significant information with respect to the pressure effect and the nonradiative recombination.

Double ion implantation and pulsed laser melting processes for third generation solar cells

In the framework of the third generation of photovoltaic devices, the intermediate band solar cell is one of the possible candidates to reach higher efficiencies with a lower processing cost. In this work, we introduce a novel processing method based on a double ion implantation and, subsequently, a pulsed laser melting (PLM) process to obtain thicker layers of Ti supersaturated Si. We perform ab initio theoretical calculations of Si impurified with Ti showing that Ti in Si is a good candidate to theoretically form an intermediate band material in the Ti supersaturated Si. From time-of-flight secondary ion mass spectroscopy measurements, we confirm that we have obtained a Ti implanted and PLM thicker layer of 135 nm. Transmission electron microscopy reveals a single crystalline structure whilst the electrical characterization confirms the transport properties of an intermediate band material/Si substrate junction. High subbandgap absorption has been measured, obtaining an approximate value of 104 cm−1 in the photons energy range from 1.1 to 0.6 eV.

Electronic and optical properties of substitutional V, Cr and Ir impurities in Cu2ZnSnS4

The substitution of cation atoms by V, Cr and It in the natural and synthetic quaternary Cu2ZnSnS4 semiconductor is analyzed using first-principles methods. In most of the substitutions, the electronic structure of these modified CZTS is characterized for intermediate bands with different occupation and position within of the energy band gap. A study of the symmetry and composition of these intermediate bands is carried out for all substitutions. These bands permit additional photon absorption and emission channels depending on their occupation. The optical properties are obtained and analyzed. The absorption coefficients are split into contributions from the different absorption channels and from the inter- and intra-atomic components. The sub bandgap transitions are significant in many cases because the anion states contribute to the valence, conduction and intermediates bands. These properties could therefore be used for novel optoelectronic devices.

Ionophores and receptors using cation-π interactions: Collarenes

Cation-π interactions are important forces in molecular recognition by biological receptors, enzyme catalysis, and crystal engineering. We have harnessed these interactions in designing molecular systems with circular arrangement of benzene units that are capable of acting as ionophores and models for biological receptors. [n]Collarenes are promising candidates with high selectivity for a specific cation, depending on n, because of their structural rigidity and well-defined cavity size. The interaction energies of [n]collarenes with cations have been evaluated by using ab initio calculations. The selectivity of these [n]collarenes in aqueous solution was revealed by using statistical perturbation theory in conjunction with Monte Carlo and molecular dynamics simulations. It has been observed that in [n]collarenes the ratio of the interaction energies of a cation with it and the cation with the basic building unit (benzene) can be correlated to its ion selectivity. We find that collarenes are excellent and efficient ionophores that bind cations through cation-π interactions. [6]Collarene is found to be a selective host for Li+ and Mg2+, [8]collarene for K+ and Sr2+, and [10]collarene for Cs+ and Ba2+. This finding indicates that [10]collarene and [8]collarene could be used for effective separation of highly radioactive isotopes, 137Cs and 90Sr, which are major constituents of nuclear wastes. More interestingly, collarenes of larger cavity size can be useful in capturing organic cations. [12]Collarene exhibits a pronounced affinity for tetramethylammonium cation and acetylcholine, which implies that it could serve as a model for acetylcholinestrase. Thus, collarenes can prove to be novel and effective ionophores/model-receptors capable of heralding a new direction in molecular recognition and host-guest chemistry.

Conformation of coenzyme pyrroloquinoline quinone and role of Ca2+ in the catalytic mechanism of quinoprotein methanol dehydrogenase

The ab initio structures of 2,7,9-tricarboxypyrroloquinoline quinone (PQQ), semiquinone (PQQH), and dihydroquinone (PQQH2) have been determined and compared with ab initio structures of the (PQQ)Ca2+, (PQQH)Ca2+, and (PQQH2)Ca2+ complexes as well as the x-ray structure of (PQQ)Ca2+ bound at the active site of the methanol dehydrogenase (MDH) of methyltropic bacteria. Plausible mechanisms for the MDH oxidation of methanol involving the (PQQ)Ca2+ complex are explored via ab initio computations and discussed. Considering the reaction of methanol with PQQ in the absence of Ca2+, nucleophilic addition of methanol to the PQQ C-5 carbonyl followed by a retro-ene elimination is deemed unlikely due to large energy barrier. A much more favorable disposition of the methanol C-5 adduct to provide formaldehyde involves proton ionization of the intermediate followed by elimination of methoxide concerted with hydride transfer to the oxygen of the C-4 carbonyl. Much the same transition state is reached if one searches for the transition state beginning with Asp-303–CO2−general-base removal of the methanol proton of the (PQQ)Ca2+O(H)CH3 complex concerted with hydride transfer to the oxygen at C-4. For such a mechanism the role of the Ca2+ moiety would be to (i) contribute to the formation of the ES complex (ii) provide a modest decrease in the pKa of methanol substrate,; and (iii) polarize the oxygen at C-5.

Femtosecond observation of benzyne intermediates in a molecular beam: Bergman rearrangement in the isolated molecule

In this communication, we report our femtosecond real-time observation of the dynamics for the three didehydrobenzene molecules (p-, m-, and o-benzyne) generated from 1,4-, 1,3-, and 1,2-dibromobenzene, respectively, in a molecular beam, by using femtosecond time-resolved mass spectrometry. The time required for the first and the second C-Br bond breakage is less than 100 fs; the benzyne molecules are produced within 100 fs and then decay with a lifetime of 400 ps or more. Density functional theory and high-level ab initio calculations are also reported herein to elucidate the energetics along the reaction path. We discuss the dynamics and possible reaction mechanisms for the disappearance of benzyne intermediates. Our effort focuses on the isolated molecule dynamics of the three isomers on the femtosecond time scale.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Análise Conformacional e estudo das interações eletrônicas de algumas α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas

A presente Dissertação relata a síntese e o estudo conformacional das α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas p-X-Φ-C(O)CH[SeΦ][P(O)(OEt2] (X=OMe 1, Me 2, H 3, F 4, Cl 5, Br 6 e NO2 7) através da banda de estiramento da carbonila no infravermelho, em solventes de polaridade crescente apoiado por cálculos ab initio HF/6-3IG**. A comparação entre a freqüência e a intensidade relativa dos componentes do dubleto, para os derivados 6 e 7, e do singleto para os derivados 1-5, no solvente apolar tetracloreto de carbono, e dos componentes do dubleto, nos solventes de polaridade crescente (clorofórmio, diclorometano e acetonitrila), para os derivados 1-7, com os dados do cálculo ab initio de 3 (composto de referência), indicou que ambas as conformações estáveis (g1 e g2) apresentam a ligação C-Se na geometria anti-clinal (gauche) em relação à carbonila (C=O), enquanto que a ligação C-P assume uma geometria sin-periplanar (cis) em relação à carbonila. A análise dos contatos interatômicos de átomos relevante em comparação com a soma de seus raios de van der Waals, indicou que ambas as conformações g1 e g2 são fortemente estabilizadas pelo sinergismo das interações orbitalares e eletrostáticas π*(CO) / nSe e Oδ-[CO].....Pδ+[PO]. Analogamente, as interações mais fracas Oδ-[OR]..... Cδ+[CO], 0-Hδ+[SeΦ]....Oδ-[PO] e 0-Hδ+[ΦC(O)]....Oδ-[CO] estabilizam as conformações g1 e g2, aproximadamente na mesma extensão. No entanto, somente a conformação g1 é estabilizada pela interação eletrostática (ligação de hidrogênio) Hδ+[α-CH].....Oδ-[OR], enquanto que sómente a conformação g2 é desestabilizada pelo Efeito de Campo Repulsivo entre os dipolos Cδ+=.Oδ- e Pδ+-ORδ- Assim sendo, pode-se concluir que no dubleto de VCO no IV, o componente de maior freqüência e de menor intensidade corresponde à conformação menos estável g2 (do cálculo) enquanto que o componente de menor freqüência e mais intenso corresponde à conformação mais estável g1 (do cálculo). Estes dados estão de pleno acordo com os deslocamentos de freqüência mais negativos da carbonila (ΔVCO) do confôrmero mais estável g1 em relação ao menos estável g2.

Estudo da cristalinidade de filmes finos de nitreto de índio e simulado pelo pacote de programas Wien2k

Neste trabalho, foi utilizado o método de deposição assistida por feixe de íons (IBAD na sigla em inglês) para produção de filmes finos de nitreto de índio em substratos de silício (111) e Safira-C. Variando as condições de deposição e utlilizando a técnica de difração de raios-X, investigou-se com o intuito de obter os parâmetros que resultam em filmes finos com melhor grau de cristalinidade. Os filmes produzidos a 380C apresentaram alta cristalinidade, superior àqueles a 250C. Temperaturas muito superiores a 380C não ocasionam a formação de filme cristalino de InN, como foi observado ao utilizar a temperatura de 480C; o mesmo se observa ao utilizar temperatura ambiente. Na temperatura considerada adequada ,de 380C, obteve-se que a utilização de Ra, ou seja, a razão de fluxo de partículas entre o nitrogênio e índio, em torno de 2,3 permite obter um melhor grau de cristalinização, o qual decresce conforme se diverge desse valor. A comparação entre difratogramas de amostras produzidas com e sem a evaporação prévia de titânio, o qual é possível observar um deslocamento dos picos do InN, indicam que o efeito Gettering permite a redução de impurezas no filme, principalmente de oxigênio. Utilizou-se a técnica de Retroespalhamento de Rutherford para obtenção da composição dos elementos e o perfil de profundidade. Notou-se uma forte mistura dos elementos do substrato de silício e safira com o nitreto de índio mesmo próximos a superfície. A presença indesejável de impurezas, principalmente o oxigênio, durante a deposição de filmes finos é praticamente inevitável. Desta forma, cálculos ab initio baseados na Teoria do Funcional da Densidade (DFT) foram realizados para investigar defeitos isolados e complexos de oxigênio no nitreto de índio e a sua influência nas propriedades óticas. Considerou-se diferentes concentrações de oxigênio (x=2,76, 8,32, 11,11 e 22,22%) aplicando-se o método PBEsolGGA e TB-mBJ para o tratamento da energia e potencial de troca e correlação. Obteve-se que é energeticamente favorável o oxigênio existir principalmente como defeito carregado e isolado. Os resultados utilizando a aproximação de TB-mBJ indicam um estreitamento do bandgap conforme a concentração de oxigênio aumenta. Entretanto, a alta contribuição do efeito de Moss-Burstein resulta num efetivo alargamento do band gap, gerando valores de band gap ótico maiores que no do bulk de nitreto de índio.

Propriedades magnéticas de trímeros de FexCo1-x depositados em Pt(111)

Quais propriedades magnéticas são modificadas quando se agrupam átomos de Fe/Co para formar estruturas quasi-2D, se comparadas aos nanofios (quasi-1D) de FexCo1-x? E como estas propriedades reagem com a variação da proporção de Fe/Co nos aglomerados? A fim de responder a estas questões, trímeros de FexCo1-x depositados em Pt(111) são investigados utilizando o método de primeiros princípios Real Space-Linear Muffin-Tin Orbital-Atomic Sphere Approximation (RS-LMTO-ASA) no âmbito da Teoria do Funcional da Densidade (DFT). Diferentes configurações de trímeros triangulares são consideradas, variando-se as posições e a concentração dos átomos de Fe/Co. Neste trabalho, demonstra-se a ocorrência de uma tendência não-linear estritamente decrescente dos momentos orbitais médios como função da concentração de Fe, distinta do encontrado tanto para os nanofios de FexCo1-x (dependência linear) quanto para a monocamada correspondente (dependência não-linear). Os resultados obtidos mostram ainda que os momentos orbitais variam com o ambiente local e com a direção de magnetização, especialmente quando associados aos átomos de Co, em concordância com publicações anteriores. A mudança de dimensionalidade quasi-1D (nanofios) para quasi-2D (trímeros compactos) não afeta o comportamento dos momentos de spin, que permanecem descritos por uma função linear com respeito à proporção de Fe/Co. Ambos o formato e a concentração de Fe nos sistemas apresentam um papel importante nos valores de energia de anisotropia magnética. Em adição, observou-se que o subtrato de Pt opera ativamente na definição das propriedades magnéticas dos aglomerados. Embora todas as configurações lineares e compactas dos aglomerados de FexCo1-x sejam estáveis e exibam interações fortemente ferromagnéticas entre os primeiros vizinhos, nem todas revelaram o ordenamento colinear como estado fundamental, apresentando uma interação de Dzyaloshinskii-Moriya não-desprezível induzida pelo acoplamento spin-órbita. Estes casos específicos são: o trímero triangular de Co puro e o trímero linear (nanofio) de Fe puro, para o qual foi verificado o acoplamento do tipo Ruderman-Kittel-Kasuya-Yosida entre os átomos de Fe constituintes. Os resultados obtidos contribuem para o entendimento de quais mecanismos definem o magnetismo nos trímeros de FexCo1-x/Pt(111), e discutem as questões presentes atualmente na literatura no contexto destes sistemas.

Dinâmica molecular de sistemas iônicos poliatômicos: modelos polarizável e não-polarizável

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito.

Prediction of hidden multiferroic order in graphene zigzag ribbons

An electronic phase with coexisting magnetic and ferroelectric order is predicted for graphene ribbons with zigzag edges. The electronic structure of the system is described with a mean-field Hubbard model that yields results very similar to those of density functional calculations. Without further approximations, the mean-field theory is recasted in terms of a BCS wave function for electron-hole pairs in the edge bands. The BCS coherence present in each spin channel is related to spin-resolved electric polarization. Although the total electric polarization vanishes, due to an internal phase locking of the BCS state, strong magnetoelectric effects are expected in this system. The formulation naturally accounts for the two gaps in the quasiparticle spectrun, Δ0 and Δ1, and relates them to the intraband and interband self-energies.