Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

The second law of thermodynamics and the global climate system : A review of the maximum entropy production principle

The long-term mean properties of the global climate system and those of turbulent fluid systems are reviewed from a thermodynamic viewpoint. Two general expressions are derived for a rate of entropy production due to thermal and viscous dissipation (turbulent dissipation) in a fluid system. It is shown with these expressions that maximum entropy production in the Earth s climate system suggested by Paltridge, as well as maximum transport properties of heat or momentum in a turbulent system suggested by Malkus and Busse, correspond to a state in which the rate of entropy production due to the turbulent dissipation is at a maximum. Entropy production due to absorption of solar radiation in the climate system is found to be irrelevant to the maximized properties associated with turbulence. The hypothesis of maximum entropy production also seems to be applicable to the planetary atmospheres of Mars and Titan and perhaps to mantle convection. Lorenz s conjecture on maximum generation of available potential energy is shown to be akin to this hypothesis with a few minor approximations. A possible mechanism by which turbulent fluid systems adjust themselves to the states of maximum entropy production is presented as a selffeedback mechanism for the generation of available potential energy. These results tend to support the hypothesis of maximum entropy production that underlies a wide variety of nonlinear fluid systems, including our planet as well as other planets and stars

Forecasting cosmic doomsday from CMB-LSS cross-correlations.

A broad class of dark energy models, which have been proposed in attempts at solving the cosmological constant problems, predict a late time variation of the equation of state with redshift. The variation occurs as a scalar field picks up speed on its way to negative values of the potential. The negative potential energy eventually turns the expansion into contraction and the local universe undergoes a big crunch. In this paper we show that cross-correlations of the cosmic microwave background anisotropy and matter distribution, in combination with other cosmological data, can be used to forecast the imminence of such cosmic doomsday.

Home Heating and Cooling, 2011

Every year, a typical family in the United States spends around half of its home energy budget on heating and cooling. In Iowa, that percentage can be higher, due to temperature extremes reached during the winter and summer months. Unfortunately, many of those dollars often are wasted, because conditioned air escapes through leaky ceilings, walls and foundations—or flows through inadequately insulated attics, exterior walls and basements. In addition, many heating systems and air conditioners aren’t properly maintained or are more than 10 years old and very inefficient, compared to models being sold today. As a result, it makes sense to analyze your home as a collection of systems that must work together in order to achieve peak energy savings. For example, you won’t get anywhere near the savings you’re expecting from a new furnace if your airhandling ducts are uninsulated and leak at every joint. The most energy-efficient central air-conditioning setup won’t perform to your expectations if your attic insulation is inadequate and can’t reduce solar heat gain to help keep your home cool. And planting the wrong types of trees or shrubs close to your home adversely can affect potential energy savings all year long.

Dimerization of polyacetylene treated as a spin-Peierls distortion of the Heisenberg Hamiltonian

Extracting a bond-length-dependent Heisenberg-like Hamiltonian from the potential-energy surfaces of the two lowest states of ethylene, it is possible to study the geometry of polyacetylene by minimization of the cohesive energy, using both variational-cluster and Rayleigh-Schrödinger perturbative expansions. The dimerization amplitude is satisfactorily reproduced. Optimizing the variational-cluster-expansion total energy with the equal-bond-length constraint, the barrier to reversal of alternation is obtained. The alternating-to-regular phase transition is treated from the Néel-state starting function and appears to be of second order.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

Energiansäästökohteiden kartoitus konepajateollisuuden yrityksessä

Tämän diplomityön tavoitteena on kartoittaa Vaahto Oy:n Hollolan tehtaan energiankulutus ja energiansäästökohteet. Ensin tutkittiin tehtaan energiankulutus ja energiankulutuksen jakautuminen. Tutkimuksessa käytettiin saatavilla olevia kulutustietoja, luettiin konekirjoja sekä haastateltiin tehtaan työntekijöitä. Lisäksi tehdashallien lämmityslaitteiden hyötysuhteet mitattiin. Tutkimuksen päätavoite oli selvittää, miksi tehtaan lämmitysenergiankulutus on kasvanut ja kannattaisiko rakennusten lämmittämiseen käyttää vaihtoehtoista lämmitysmuotoa öljylämmitykselle. Potentiaalisista energiansäästökohteista tehtiin investointilaskelmat ja toimenpide-ehdotukset. Kannattaviksi toimenpiteiksi tutkimuksessa todettiin: lämmityspolttoaineen vaihtaminen maakaasuun, nosto-ovien hankkiminen, paineilmaverkon huolto, paineilmakompressorin lämmöntalteenotto, tehdastilojen sisälämpötilan tarkastus ja työnjohtotilojen ilmanvaihdon käyntiaikojen muutos. Toimenpiteillä arvioidaan vuotuisten energiakustannusten pienenevän noin 34 000 euroa. Toimenpiteiden toteuttamisen arvioidaan maksavan 135 000 ¤, mistä lämmitysjärjestelmän vaihdon osuus on 100 000 ¤.

0 AND 02 ADSORPTION ON Ag/Cu[100]

In this work the adsorption mechanisms of atomic and molecular oxygen on Cu(100) surface are studied using ab initio simulation methods. Through the atomistic scale under-standing of the elementary oxidation processes we can further understand the large-scale oxidation. Copper is a material widely used in industry which makes it an interesting subject, and also understanding the oxidation of copper helps us understand the oxidation mechanism of other metals. First we have a look on some theory on surface alloys in general and behaviour of Ag on Cu(100) surface. After that the physical background there is behind the methods of density functional calculations are discussed, and some methods, namely potential energy surfaces and molecular dynamics, are introduced. Then there is a brief look on the numerical details used in the calculations, and after that, the results of the simulations are exhibited.

Oxygen adsorption on Cu(510)

Tässä työssä on tutkittu kuparin (510)-askelpinnan reaktiivisuutta käyttäen apuna kvanttimekaanisia ab initio laskentamenetelmiä. Tutkimus on toteutettu laskemalla happiatomin adsorptioenergia ja tilatiheys erilaisissa potentiaalisissa adsorptiopaikoissa pinnalla. Myös happimolekyylin adsorptiota ja hajoamista ontarkasteltu laskemalla pintaa lähestyvälle molekyylille potentiaalienergiapintoja. Energiapintojen tuloksia on myös täydennetty kvanttimekaanisilla molekyylidynamiikkalaskuilla. Metallisia askelpintoja pidetään yleisesti sileitä pintoja reaktiivisempina happea kohtaan, johtuen askeleen reunan pienentävästä vaikutuksesta molekyylin hajoamisen tiellä olevaan energiavaliin. On kuitenkin olemassa myös tuloksia, jotka osoittavat hapen tarttumisprosessin olevan hallitseva juuri terassialueella, askeleen reunan sijasta. Tässä työssä on todettu hapen adsorboituvan Cu(510)-pinnalla tehokkaimmin juuri terassilla olevaan hollow-paikkaan. Myös adsorptioenergiat ovat tällä pinnalla pienempiä kuin sileällä (100)-pinnalla. Potentiaalienergiapintojen perusteella Cu(510)-pinnan todetaan myös olevan vähemmän reaktiivinen kuin askelpintojen yleisesti odotetaan olevan, vaikka askeleen reunan todetaankin pienentävän happiatominhajoamisen esteenä olevaa energiavallia.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

Local thermodynamics and the generalized Gibbs-Duhem equation in systems with long-range interactions

The local thermodynamics of a system with long-range interactions in d dimensions is studied using the mean-field approximation. Long-range interactions are introduced through pair interaction potentials that decay as a power law in the interparticle distance. We compute the local entropy, Helmholtz free energy, and grand potential per particle in the microcanonical, canonical, and grand canonical ensembles, respectively. From the local entropy per particle we obtain the local equation of state of the system by using the condition of local thermodynamic equilibrium. This local equation of state has the form of the ideal gas equation of state, but with the density depending on the potential characterizing long-range interactions. By volume integration of the relation between the different thermodynamic potentials at the local level, we find the corresponding equation satisfied by the potentials at the global level. It is shown that the potential energy enters as a thermodynamic variable that modifies the global thermodynamic potentials. As a result, we find a generalized Gibbs-Duhem equation that relates the potential energy to the temperature, pressure, and chemical potential. For the marginal case where the power of the decaying interaction potential is equal to the dimension of the space, the usual Gibbs-Duhem equation is recovered. As examples of the application of this equation, we consider spatially uniform interaction potentials and the self-gravitating gas. We also point out a close relationship with the thermodynamics of small systems.

Íon Be H+ em campo elétrico uniforme intenso

The influence of a strong external electric field in chemical bonding is not extensively studied. In this work, the behavior of the potential energy curves of the BeH+ ion in an external electric field is investigated by means of variational calculations using molecular orbitals. For simplicity we consider the case where the external field is parallel to the molecular axis, in the direction of the dipole moment, which raises the polarization of the ion. The effects on dissociation are investigated.

Espectro infravermelho e análise conformacional do composto 3-fenil-2-oxo-1,2,3-oxatiazolidina

The infrared (IR) spectra of the four distict conformers located on the multidimensional potential energy surface (PES) for the 3-phenyl-1,2,3-oxathiazolidine 2-oxide compound have been calculated using the semiempirical quantum-mechanical method PM3. The band spectra are reported and compared directly with the experimental spectrum. The IR intensities are shown to be much more sensitive to conformational changes than the vibrational frequencies and so, the theoretical analysis of the IR spectrum can be used as a tool for helping in the elucidation of the structure of heterocyclic compounds.

Elementos da dinâmica química ao nível da teoria variacional do estado de transição com correções interpoladas

In this article are presented some fundamental elements of the conventional and of the variational transition state theories which are needed to carried out calculations of semi-classical chemical dynamics. Some important bottlenecks in building reliable potential energy surfaces using electronic structure calculations are also discussed. It is put emphasis on the methodology of the variational transition state theory with interpolated corrections (VTST-IC), and its application in the calculations of the rate constants and of the kinetic isotope effect (KIE) of CH4 + Cl « CH3 + HCl reaction.