Fully relativistic calculations on the potential energy surfaces of the lowest 23 states of molecular chlorine

The electronic structure and spectroscopic properties (R(e), omega(e), omega(e)x(e), beta(e), and T(e)) of the ground state and the 22 lowest excited states of chlorine molecule were studied within a four-component relativistic framework using the MOLFDIR program package. The potential energy curves of all possible 23 covalent states were calculated using relativistic complete open shell configuration interaction approach. In addition, four component multireference configuration interaction with single and double excitation calculations were performed in order to infer the effects due to dynamical correlation in vertical excitations. The calculated properties are in good agreement with the available experimental data.

Morse potential energy spectra through the variational method and supersymmetry

The Variational Method is applied within the context of Supersymmetric Quantum Mechanics to provide information about the energy and eigenfunction of the lowest levels of a Hamiltonian. The approach is illustrated by the case of the Morse potential applied to several diatomic molecules and the results are compared with stabilished results. (C) 2000 Elsevier Science B.V.

Unitarity and nonrelativistic potential energy in a higher-order Lorentz symmetry breaking electromagnetic model

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

The Dynamics of a Small Model Glass Former as Viewed from Its Potential Energy Landscape

Despite intensive research during the last decades, thetheoreticalunderstanding of supercooled liquids and the glasstransition is stillfar from being complete. Besides analytical investigations,theso-called energy-landscape approach has turned out to beveryfruitful. In the literature, many numerical studies havedemonstratedthat, at sufficiently low temperatures, all thermodynamicquantities can be predicted with the help of the propertiesof localminima in the potential-energy-landscape (PEL). The main purpose of this thesis is to strive for anunderstanding ofdynamics in terms of the potential energy landscape. Incontrast to the study of static quantities, this requirestheknowledge of barriers separating the minima.Up to now, it has been the general viewpoint that thermallyactivatedprocesses ('hopping') determine the dynamics only belowTc(the critical temperature of mode-coupling theory), in thesense that relaxation rates follow from local energybarriers.As we show here, this viewpoint should be revisedsince the temperature dependence of dynamics is governed byhoppingprocesses already below 1.5Tc.At the example of a binary mixture of Lennard-Jonesparticles (BMLJ),we establish a quantitative link from the diffusioncoefficient,D(T), to the PEL topology. This is achieved in three steps:First, we show that it is essential to consider wholesuperstructuresof many PEL minima, called metabasins, rather than singleminima. Thisis a consequence of strong correlations within groups of PELminima.Second, we show that D(T) is inversely proportional to theaverageresidence time in these metabasins. Third, the temperaturedependenceof the residence times is related to the depths of themetabasins, asgiven by the surrounding energy barriers. We further discuss that the study of small (but not toosmall) systemsis essential, in that one deals with a less complex energylandscapethan in large systems. In a detailed analysis of differentsystemsizes, we show that the small BMLJ system consideredthroughout thethesis is free of major finite-size-related artifacts.

AM1 and PM3 Calculations of the Potential Energy Surfaces for CH2OH Reactions with NO and NO2

The AM1 and PM3 molecular orbital methods have been utilized to investigate the reactions of CH20H with NO and NO2 PM3 and AM1 calculated heats of formation differ from experimental values by 8.6 and 18.8 kcal mol-', respectively. The dominant reaction of CH20H with NO is predicted to produce the adduct HOCH2N0, supporting the hypothesis of Pagsberg, Munk, Anastasi, and Simpson. Calculated activation energies for the NO2 system predict the formation of the adducts HOCH2N02 and HOCH20N0. In addition, the PM3 calculations predict that the abstraction reaction producing CH20 and HN02 is more likely than one producing CH20 and HONO from reactions of CH20H with NO2.

Exploration of the Potential Energy Surfaces, Prediction of Atmospheric Concentrations, and Prediction of Vibrational Spectra for the HO2···(H2O)n (n = 1−2) Hydrogen Bonded Complexes

The hydroperoxy radical (HO2) plays a critical role in Earth's atmospheric chemistry as a component of many important reactions. The self-reaction of hydroperoxy radicals in the gas phase is strongly affected by the presence of water vapor. In this work, we explore the potential energy surfaces of hydroperoxy radicals hydrogen bonded to one or two water molecules, and predict atmospheric concentrations and vibrational spectra of these complexes. We predict that when the HO2 concentration is on the order of 108molecules·cm-3 at 298 K, that the number of HO2···H2O complexes is on the order of 107molecules·cm-3 and the number of HO2···(H2O)2 complexes is on the order of 106molecules·cm-3. Using the computed abundance of HO2···H2O, we predict that, at 298 K, the bimolecular rate constant for HO2···H2O + HO2 is about 10 times that for HO2 + HO2.

An Appraisal of the Current Status and the Potential of Surface Water in Upper Anseba Catchment, Eritrea

The study deals with the status and potential of surface water resources in Upper Anseba, Central Highlands of Eritrea, one of the most densely populated regions in Eritrea, including small scale farming and the country's capital city. water demand is increasing rapidly for all uses. The area has no perennial water course and depends very largely on reservoirs for its water supply. The report finds that there are 74 reservoirs in the area, of which 49 are in Upper Anseba. Total reservoir capacity already corresponds to 70% of runoff. the capacity of some of the reservoirs already exceeds annual runoff of their catchment. Recommendations thus include the use of water saving technologies for irrigation; and above all, preparation of a regional master plan for development, including water allocation planning with a mid term perspective.

Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.

Universality of the self gravitational potential energy of any fundamental particle.

Using the relation proposed by Weinberg in 1972, combining quantum and cosmological parameters, we prove that the self gravitational potential energy of any fundamental particle is a quantum, with physical properties independent of the mass of the particle. It is a universal quantum of gravitational energy, and its physical properties depend only on the cosmological scale factor R and the physical constants ℏ and c. We propose a modification of the Weinberg’s relation, keeping the same numerical value, but substituting the cosmological parameter H/c by 1/R.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Theoretical study of two states reactivity of methane activation on iron atom and iron dimer

Density functional theory (DFT) calculations have been carried out to explore the catalytic activation of C–H bonds in methane by the iron atom, Fe, and the iron dimer, Fe2. For methane activation on an Fe atom, the calculations suggest that the activation of the first C–H bond is mediated via the triplet excited-state potential energy surface (PES), with initial excitation of Fe to the triplet state being necessary for the reaction to be energetically feasible. Compared with the breaking of the first C–H bond, the cleavage of the second C–H bond is predicted to involve a significantly higher barrier, which could explain experimental observations of the HFeCH3 complex rather than CH2FeH2 in the activation of methane by an Fe atom. For methane activation on an iron dimer, the cleavage of the first C–H bond is quite facile with a barrier only 11.2, 15.8 and 8.4 kcal/mol on the septet state energy surface at the B3LYP/6-311+G(2df,2dp), BPW91/6-311+G(2df,2dp) and M06/B3LYP level, respectively. Cleavage of the second C–H bond from HFe2CH3 involves a barrier calculated respectively as 18.0, 10.7 and 12.4 kcal/mol at the three levels. The results suggest that the elimination of hydrogen from the dihydrogen complex is a rate-determining step. Overall, our results indicate that the iron dimer Fe2 has a stronger catalytic effect on the activation of methane than the iron atom.

Structure, magnetism and colour in simple bis­(phosphine)nickel(II) dihalide complexes : an experimental and theoretical investigation

The complex [1,2-bis­(di-tert-butyl­phosphan­yl)ethane-[kappa]2P,P']di­iodido­nickel(II), [NiI2(C18H40P2] or (dtbpe-[kappa]2P)NiI2, [dtbpe is 1,2-bis­(di-tert-butyl­phosphan­yl)ethane], is bright blue-green in the solid state and in solution, but, contrary to the structure predicted for a blue or green nickel(II) bis­(phos­phine) complex, it is found to be close to square planar in the solid state. The solution structure is deduced to be similar, because the optical spectra measured in solution and in the solid state contain similar absorptions. In solution at room temperature, no 31P{1H} NMR resonance is observed, but the very small solid-state magnetic moment at temperatures down to 4 K indicates that the weak paramagnetism of this nickel(II) complex can be ascribed to temperature independent paramagnetism, and that the complex has no unpaired electrons. The red [1,2-bis­(di-tert-butyl­phosphan­yl)ethane-[kappa]2P,P']di­chlorido­nickel(II), [NiCl2(C18H40P2] or (dtbpe-[kappa]2P)NiCl2, is very close to square planar and very weakly paramagnetic in the solid state and in solution, while the maroon [1,2-bis­(di-tert-butyl­phosphan­yl)ethane-[kappa]2P,P']di­bromido­nickel(II), [NiBr2(C18H40P2] or (dtbpe-[kappa]2P)NiBr2, is isostructural with the diiodide in the solid state, and displays paramagnetism inter­mediate between that of the dichloride and the diiodide in the solid state and in solution. Density functional calculations demonstrate that distortion from an ideal square plane for these complexes occurs on a flat potential energy surface. The calculations reproduce the observed structures and colours, and explain the trends observed for these and similar complexes. Although theoretical investigation identified magnetic-dipole-allowed excitations that are characteristic for temperature-independent paramagnetism (TIP), theory predicts the mol­ecules to be diamagnetic.

Reactions of simple and peptidic alpha-carboxylate radical anions with dioxygen in the gas phase

alpha-Carboxylate radical anions are potential reactive intermediates in the free radical oxidation of biological molecules (e. g., fatty acids, peptides and proteins). We have synthesised well-defined alpha-carboxylate radical anions in the gas phase by UV laser photolysis of halogenated precursors in an ion-trap mass spectrometer. Reactions of isolated acetate ((center dot)CH(2)CO(2)) and 1-carboxylatobutyl (CH(3)CH(2)CH(2)(center dot)CHCO(2)(-)) radical anions with dioxygen yield carbonate (CO(3)(center dot-)) radical anions and this chemistry is shown to be a hallmark of oxidation in simple and alkyl-substituted cross-conjugated species. Previous solution phase studies have shown that C(alpha)-radicals in peptides, formed from free radical damage, combine with dioxygen to form peroxyl radicals that subsequently decompose into imine and keto acid products. Here, we demonstrate that a novel alternative pathway exists for two alpha-carboxylate C(alpha)-radical anions: the acetylglycinate radical anion (CH(3)C(O)NH(center dot)CHCO(2)(-)) and the model peptide radical anion, YGGFG(center dot-). Reaction of these radical anions with dioxygen results in concerted loss of carbon dioxide and hydroxyl radical. The reaction of the acetylglycinate radical anion with dioxygen reveals a two-stage process involving a slow, followed by a fast kinetic regime. Computational modelling suggests the reversible formation of the C(alpha) peroxyl radical facilitates proton transfer from the amide to the carboxylate group, a process reminiscent of, but distinctive from, classical proton-transfer catalysis. Interestingly, inclusion of this isomerization step in the RRKM/ME modelling of a G3SX level potential energy surface enables recapitulation of the experimentally observed two-stage kinetics.