Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Measurement of Bottom Versus Charm as a Function of Transverse Momentum with Electron-Hadron Correlations in p plus p Collisions at s=200 GeV

The momentum distribution of electrons from semileptonic decays of charm and bottom quarks for midrapidity |y|< 0.35 in p+p collisions at s=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range 2 < p(T)< 7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D -> e(+/-)K(-/+)X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in p(T). A fixed-order-plus-next-to-leading-log perturbative quantum chromodynamics calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is sigma(bb)=3.2(-1.1)(+1.2)(stat)(-1.3)(+1.4)(syst)mu b.

Semiclassical scattering in two dimensions

The semiclassical limit of quantum mechanical scattering in two dimensions is developed and the Wentzel-Kramers-Brillouin and eikonal results for two-dimensional scattering is derived. No backward or forward glory scattering is present in two dimensions. Other phenomena, such as rainbows and orbiting, do occur. (C) 2008 American Association of Physics Teachers.

Optical third-harmonic spectroscopy of the magnetic semiconductor EuTe

EuTe possesses the centrosymmetric crystal structure m3m of rocksalt type in which the second-harmonic generation is forbidden in electric dipole approximation but the third-harmonic generation (THG) is allowed. We studied the THG spectra of this material and observed several resonances in the vicinity of the band gap at 2.2-2.5 eV and at higher energies up to 4 eV, which are related to four-photon THG processes. The observed resonances are assigned to specific combinations of electronic transitions between the ground 4f(7) state at the top of the valence band and excited 4f(6)5d(1) states of Eu(2+) ions, which form the lowest energy conduction band. Temperature, magnetic field, and rotational anisotropy studies allowed us to distinguish crystallographic and magnetic-field-induced contributions to the THG. A strong modification of THG intensity for the 2.4 eV band and suppression of the THG for the 3.15 eV band was observed in applied magnetic field. Two main features of the THG spectra were assigned to 5d(t(2g)) and 5d(e(g)) subbands at 2.4 eV and 3.15 eV, respectively. A microscopic quantum-mechanical model of the THG response was developed and its conclusions are in qualitative agreement with the experimental results.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

The role of oxygen vacancy in the photoluminescence property at room temperature of the CaTiO(3)

In this paper, electron paramagnetic resonance, photoluminescence (PL) emission, and quantum mechanical calculations were used to observe and understand the structural order-disorder of CaTiO(3), paying special attention to the role of oxygen vacancy. The PL phenomenon at room temperature of CaTiO(3) is directly influenced by the presence of oxygen vacancies that yield structural order-disorder. These oxygen vacancies bonded at Ti and/or Ca induce new electronic states inside the band gap. Ordered and disordered CaTiO(3) was obtained by the polymeric precursor method. (C) 2009 American Institute of Physics. [DOI: 10.1063/1.3190524]

NMR quadrupolar system described as Bose-Einstein-condensate-like system

This paper presents a description of nuclear magnetic resonance (NMR) of quadrupolar systems using the Holstein-Primakoff (HP) formalism and its analogy with a Bose-Einstein condensate (BEC) system. Two nuclear spin systems constituted of quadrupolar nuclei I=3/2 ((23)Na) and I=7/2 ((133)Cs) in lyotropic liquid crystals were used for experimental demonstrations. Specifically, we derived the conditions necessary for accomplishing the analogy, executed the proper experiments, and compared with quantum mechanical prediction for a Bose system. The NMR description in the HP representation could be applied in the future as a workbench for BEC-like systems, where the statistical properties may be obtained using the intermediate statistic, first established by Gentile. The description can be applied for any quadrupolar systems, including new developed solid-state NMR GaAS nanodevices.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Entangled States and Collective Nonclassical Effects in Two-Atom Systems

We propose a review of recent developments on entanglement and nonclassical effects in collective two-atom systems and present a uniform physical picture of the many predicted phenomena. The collective effects have brought into sharp focus some of the most basic features of quantum theory, such as nonclassical states of light and entangled states of multiatom systems. The entangled states are linear superpositions of the internal states of the system which cannot be separated into product states of the individual atoms. This property is recognized as entirely quantum-mechanical effect and have played a crucial role in many discussions of the nature of quantum measurements and, in particular, in the developments of quantum communications. Much of the fundamental interest in entangled states is connected with its practical application ranging from quantum computation, information processing, cryptography, and interferometry to atomic spectroscopy.

Group theory and quasiprobability integrals of Wigner functions

The integral of the Wigner function of a quantum-mechanical system over a region or its boundary in the classical phase plane, is called a quasiprobability integral. Unlike a true probability integral, its value may lie outside the interval [0, 1]. It is characterized by a corresponding selfadjoint operator, to be called a region or contour operator as appropriate, which is determined by the characteristic function of that region or contour. The spectral problem is studied for commuting families of region and contour operators associated with concentric discs and circles of given radius a. Their respective eigenvalues are determined as functions of a, in terms of the Gauss-Laguerre polynomials. These polynomials provide a basis of vectors in a Hilbert space carrying the positive discrete series representation of the algebra su(1, 1) approximate to so(2, 1). The explicit relation between the spectra of operators associated with discs and circles with proportional radii, is given in terms of the discrete variable Meixner polynomials.

Photochemical Reactions of fac-[Mn(CO)(3)(phen)imidazole](+): Evidence for Long-Lived Radical Species Intermediates

The electronic absorption spectrum of fac[Mn(CO)(3)(phen)imH](+), fac-1 in CH(2)Cl(2) is characterized by a strong absorption band at 378 nm (epsilon(max) = 3200 mol(-1) L cm(-1)). On the basis of quantum mechanical calculations, the visible absorption band has been assigned to ligand-to-ligand charge-transfer (LLCT, im -> phen) and metal-to-ligand charge-transfer (MLCT, Mn -> phen) charge transfer transition. When fac-1 in CH(2)Cl(2) is irradiated with 350 nm continuous light, the absorption features are gradually shifted to represent those of the meridional complex mer-[Mn(CO)(3)(phen)imH](+), mer-1 (lambda(max) = 556 nm). The net photoreaction under these conditions is a photoisomerization, although, the presence of the long-lived radical species was also detected by (1)H NMR and FTIR spectroscopy. 355 nm continuous photolysis of fac-1 in CH(3)CN solution also gives the long-lived intermediate which is readly trapped by metylviologen (MV(2+)) giving rise to the formation of the one-electron reduced methyl viologen (MV(center dot+)). The UV-vis spectra monitored during the slow (45 min) thermal back reaction exhibited isosbestic conversion at 426 nm. On the basis of spectroscopic techniques and quantum mechanical calculations, the role of the radicals produced is analyzed.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

Utilização de Cabos de alta Temperatura em linhas Aéreas em Média Tensão

A presente dissertação teve como objetivo fazer uma análise da viabilidade técnica da utilização dos condutores de alta temperatura nas linhas aéreas de MT, identificar vantagens, analisar inconvenientes, e estabelecer um comparativo a custos médios com as soluções convencionais. Foi efetuado o estudo de um caso real da EDP Distribuição que consistia na necessidade do aumento da capacidade de transporte de energia da linha aérea a 15 kV Espinho-Sanguedo. Neste foi ponderada a solução onde se poderia efetuar passagem de linha simples para linha dupla em alumínio-aço (AA) 160 mm2 ou a solução alternativa e inovadora de substituição dos condutores existentes por condutores de alta temperatura ACCC 182 mm2. Para isso foram efetuados cálculos e também criada uma ferramenta de apoio à decisão, para validação dos mesmos, com o intuito de mais tarde poder ser aplicada nas linhas aéreas em Média Tensão em todo o país e, sempre que necessário, se possa fazer um estudo de ponderação técnica de forma sistemática e estruturada. Neste trabalho estão identificadas as vantagens, foram relatados os inconvenientes, e estabeleceu-se um comparativo a custos médios da utilização de condutores de alta temperatura com as soluções convencionais. Antes de poder ser realizado um estudo do caso concreto da Linha aérea Espinho-Sanguedo foi necessário um aprofundamento do estado da arte no que diz respeito à comparação entre o cabo de alta temperatura ACCC e o cabo convencional ACSR, sendo este o mais utilizado nas linhas aéreas em MT. Os cabos de alta temperatura trouxeram inovações neste tema de transporte de energia, e como tal surgiu a necessidade de um estudo mais aprofundado da sua constituição, destacando o seu núcleo formado pelo compósito de fibra de carbono e fibra de vidro. Foi também analisado vantagens e desvantagens do cabo de alta temperatura e até mesmo situações onde a sua aplicação poderá ser vantajosa, de modo a tirar proveito das suas caraterísticas em que se destacam altas temperaturas de funcionamento e flechas reduzidas. Para elaborar um projeto de uma linha aérea em média tensão é necessário considerar a legislação em vigor, os aspetos ambientais e económicos, respeitando e garantindo as premissas do cálculo elétrico e mecânico. Economicamente este tipo de cabo (ACCC) é mais dispendioso do que os convencionais, no entanto o estudo realizado permitiu perceber que a sua implementação técnica é vantajosa em linhas aéreas de elevada capacidade de transporte de energia, sobretudo nos casos onde serão necessárias instalar linhas duplas ou linhas simples de seções elevadas. Devido às suas caraterísticas mecânicas, estes cabos permitem melhorar as linhas na sua dimensão, podendo diminuir o número de apoios a instalar, podendo diminuir a robustez dos apoios e permitir maior facilidade na montagem. Estas vantagens traduzem-se em menores impactos ambientais e permitem sobretudo reduzir os constrangimentos com os proprietários dos terrenos onde os apoios são implantados.

Electrically tunable resonant scattering in fluorinated bilayer graphene

Adatom-decorated graphene offers a promising new path towards spintronics in the ultrathin limit. We combine experiment and theory to investigate the electronic properties of dilutely fluorinated bilayer graphene, where the fluorine adatoms covalently bond to the top graphene layer. We show that fluorine adatoms give rise to resonant impurity states near the charge neutrality point of the bilayer, leading to strong scattering of charge carriers and hopping conduction inside a field-induced band gap. Remarkably, the application of an electric field across the layers is shown to tune the resonant scattering amplitude from fluorine adatoms by nearly twofold. The experimental observations are well explained by a theoretical analysis combining Boltzmann transport equations and fully quantum-mechanical methods. This paradigm can be generalized to many bilayer graphene-adatom materials, and we envision that the realization of electrically tunable resonance may be a key advantage in graphene-based spintronic devices.