CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Density profiles in the raise and peel model with and without a wall; physics and combinatorics

We consider the raise and peel model of a one-dimensional fluctuating interface in the presence of an attractive wall. The model can also describe a pair annihilation process in disordered unquenched media with a source at one end of the system. For the stationary states, several density profiles are studied using Monte Carlo simulations. We point out a deep connection between some profiles seen in the presence of the wall and in its absence. Our results are discussed in the context of conformal invariance ( c = 0 theory). We discover some unexpected values for the critical exponents, which are obtained using combinatorial methods. We have solved known ( Pascal`s hexagon) and new (split-hexagon) bilinear recurrence relations. The solutions of these equations are interesting in their own right since they give information on certain classes of alternating sign matrices.

The reaction of methyl radical with nitrogen atom on the triplet potential energy surface: A CCSD(T)/CBS characterization

CCSD(T)/cc-pVnZ (n = D, T, Q) calculations followed by extrapolations to the CBS limit are used to characterize stationary states of species participating in the N((4)S) + CH(3) (2A ``) reaction on the triplet PES. A mechanistic model is investigated and reaction rates are computed for every step and the overall reaction. Our best CBS estimate (1.93 x 10(10) cm(3) molecule(1) s(1)) for the overall rate constant leading to the formation of H(2)CN + H compares well with the experimental values (8.5 x 10 (11) and 1.3 x 10(10) cm(3) molecule(1) s(1)), thus reducing significantly the discrepancy of a previous theoretical result (9.1 x 10(12) cm(3) molecule(1) s(1)). (C) 2008 Elsevier B.V. All rights reserved.

Tunneling time through a barrier using the local value of a time operator

A time for a quantum particle to traverse a barrier is obtained for stationary states by setting the local value of a time operator equal to a constant. This time operator, called the tempus operator because it is distinct from the time of evolution, is defined as the operator canonically conjugate to the energy operator. The local value of the tempus operator gives a complex time for a particle to traverse a barrier. The method is applied to a particle with a semiclassical wave function, which gives, in the classical limit, the correct classical traversal time. It is also applied to a quantum particle tunneling through a rectangular barrier. The resulting complex tunneling time is compared with complex tunneling times from other methods.

Homoclinic bifurcations in reversible Hamiltonian systems

We study the existence of homoclic solutions for reversible Hamiltonian systems taking the family of differential equations u(iv) + au - u +f(u, b) = 0 as a model, where fis an analytic function and a, b real parameters. These equations are important in several physical situations such as solitons and in the existence of finite energy stationary states of partial differential equations, but no assumptions of any kind of discrete symmetry is made and the analysis here developed can be extended to others Hamiltonian systems and successfully employed in situations where standard methods fail. We reduce the problem of computing these orbits to that of finding the intersection of the unstable manifold with a suitable set and then apply it to concrete situations. We also plot the homoclinic values configuration in parameters space, giving a picture of the structural distribution and a geometrical view of homoclinic bifurcations. (c) 2005 Published by Elsevier B.V.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Uma análise da dinâmica do pêndulo eletromecânico utilizando a teoria de pertubações

Pós-graduação em Matemática - IBILCE

Dinâmica de vórtices em filmes finos supercondutores de superfície variável

Pós-graduação em Ciência e Tecnologia de Materiais - FC

A Laplace transform approach to the quantum harmonic oscillator

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Exact correlation functions in particle-reaction models with immobile particles

Exact results on particle densities as well as correlators in two models of immobile particles, containing either a single species or else two distinct species, are derived. The models evolve following a descent dynamics through pair annihilation where each particle interacts once at most throughout its entire history. The resulting large number of stationary states leads to a non-vanishing configurational entropy. Our results are established for arbitrary initial conditions and are derived via a generating function method. The single-species model is the dual of the 1D zero-temperature kinetic Ising model with Kimball-Deker-Haake dynamics. In this way, both in finite and semi-infinite chains and also the Bethe lattice can be analysed. The relationship with the random sequential adsorption of dimers and weakly tapped granular materials is discussed.

Irreversibility by competition on a Glauber-Ising model by means of the entropy production.

An out of equilibrium Ising model subjected to an irreversible dynamics is analyzed by means of a stochastic dynamics, on a effort that aims to understand the observed critical behavior as consequence of the intrinsic microscopic characteristics. The study focus on the kinetic phase transitions that take place by assuming a lattice model with inversion symmetry and under the influence of two competing Glauber dynamics, intended to describe the stationary states using the entropy production, which characterize the system behavior and clarifies its reversibility conditions. Thus, it is considered a square lattice formed by two sublattices interconnected, each one of which is in contact with a heat bath at different temperature from the other. Analytical and numerical treatments are faced, using mean-field approximations and Monte Carlo simulations. For the one dimensional model exact results for the entropy production were obtained, though in this case the phase transition that takes place in the two dimensional counterpart is not observed, fact which is in accordance with the behavior shared by lattice models presenting inversion symmetry. Results found for the stationary state show a critical behavior of the same class as the equilibrium Ising model with a phase transition of the second order, which is evidenced by a divergence with an exponent µ ¼ 0:003 of the entropy production derivative.

Beeinflussung der Grenzflächenspannung zwischen den Schmelzen "unverträglicher" Homopolymere durch Copolymere

Im Rahmen dieser Arbeit wurde am System Polyethylenoxid / Polypropylenoxid (PEO / PPO) der Einfluß von Copolymeren auf die Grenzflächenspannung Sigma von Homopolymerblends untersucht. Als Additive dienten Triblockcopolymere EO-block-PO-block-EO bzw. PO-block-EO-block-PO, Diblockcopolymere S-block-EO sowie statistische Copolymere EO-ran-PO. Die Additive wurden so ausgewählt, daß sich Paare von Additiven jeweils in genau einer Eigenschaft (Zusammensetzung, Kettenlänge, Blockanordnung) unterscheiden, in allen anderen Parametern jedoch vergleichbar sind. Die Grenzflächenspannung wurde experimentell mit Hilfe der Pendant-Drop-Methode in Abhängigkeit von der Temperatur ermittelt, wobei das Polymer mit der höheren Dichte, PEO, die Tropfenphase und PPO die Matrixphase bildet. Das Additiv wurde bei Messung der Grenzflächenspannung der ternären Systeme in unterschiedlichen Konzentrationen entweder einer oder beiden Homopolymerphasen zugegeben. Die Konzentrationsabhängigkeit von Sigma lässt sich sowohl mit dem Modell von Tang und Huang als auch mit einem Langmuir-analogen Ansatz gut beschreiben.Um den Zusammenhang zwischen sigma und dem Phasenverhalten zu untersuchen, wurden für einige der ternären Systeme Trübungskurven bei 100°C aufgenommen. Der Vergleich zwischen den Phasendiagrammen und den korrespondierenden Werten von sigma weist darauf hin, dass ein Additiv sigma gerade dann wirksam reduziert, wenn es einem Homopolymer zugefügt wird, mit dem es nur begrenzt verträglich ist, da dann die Triebkraft zur Anlagerung an der Grenzfläche besonders ausgeprägt ist. Das bereits bekannte Phänomen, wonach der Wert der Grenzflächenspannung davon abhängig sein kann, in welcher der Phasen das Additiv zu Beginn der Messung vorliegt, wurde ausführlich untersucht. Es wird angenommen, dass das System nicht in jedem Fall das thermodynamische Gleichgewicht erlangt und der beobachtete Effekt auf das Erreichen stationärer Zustände zurückzuführen ist. Dieses Verhalten kann mit einem Modell beschrieben werden, in welches das Viskositätsverhältnis der Homopolymere sowie der Verteilungskoeffizient des Copolymers zwischen den Homopolymerphasen eingehen. Aus Löslichkeitsparametern wurde der binäre Wechselwirkungsparameter Chi PEO/PPO = 0.18 abgeschätzt und mit diesem die theoretischen Werte für sigma zwischen PEO und PPO nach den Modellen von Roe bzw. Helfand und Tagami berechnet. Der Vergleich mit den experimentellen Daten des binären Systems zeigt, dass beide Ansätze sigma-Werte liefern, die in der Größenordnung der experimentellen Daten liegen, hierbei erweist sich der Ansatz von Roe als besonders geeignet. Die Temperaturabhängigkeit der Grenzflächenspannung wird jedoch durch beide Ansätze unzutreffend wiedergegeben. Mit dem Modell von Helfand und Tagami wurden eine Grenzflächendicke von 7.9 Å und das Dichteprofil der Grenzfläche berechnet. Für die Copolymere EO92PO56EO92 und S9EO22 (die Indices geben die Zahl der Monomereinheiten an) können die Grenzflächenüberschusskonzentrationen, die kritische Mizellenkonzentration sowie der einem Additivmolekül an der Grenzschicht zur Verfügung stehende Platz bestimmt werden.Der Vergleich unterschiedlicher Copolymere hinsichtlich ihrer Fähigkeit, sigma wirkungsvoll herabzusetzen, zeigt, dass im Fall von Triblockcopolymeren die Anordnung der Blöcke gegenüber der Zusammensetzung eine untergeordnete Rolle spielt. Mit zunehmender Kettenlänge nimmt die Effektivität als Compatibilizer sowohl bei Blockcopolymeren als auch bei statistischen Copolymeren zu.

Aplicação e avaliação de desempenho de um sistema de otimização em tempo real em uma unidade de produção de propeno.

Com o objetivo de aumentar o lucro de plantas químicas, a Otimização em Tempo Real (RTO) é uma ferramenta que busca determinar as condições ótimas operacionais do processo em estado estacionário, respeitando as restrições operacionais estabelecidas. Neste trabalho foi realizada a implementação prática de um ciclo RTO em um processo de destilação por recompressão de vapor (VRD), propileno-propano, da Refinaria de Paulínia (Petrobras S.A.), a partir de dados históricos da planta. Foram consideradas as principais etapas de um ciclo clássico de RTO: identificação de estado estacionário, reconciliação de dados, estimação de parâmetros e otimização econômica. Essa unidade foi modelada, simulada e otimizada em EMSO (Environment for Modeling, Simulation and Optimization), um simulador de processos orientado a equações desenvolvido no Brasil. Foram analisados e comparados dois métodos de identificação de estado estacionário, um baseado no teste estatístico F e outro baseado em wavelets. Ambos os métodos tiveram resultados semelhantes e mostraram-se capazes de identificar os estados estacionários de forma satisfatória, embora seja necessário o ajuste de parâmetros na sua implementação. Foram identificados alguns pontos estacionários para serem submetidos ao ciclo RTO e foi possível verificar a importância de partir de um estado estacionário para a continuidade do ciclo, já que essa é uma premissa do método. A partir dos pontos analisados, os resultados deste estudo mostram que o RTO é capaz de aumentar o ganho econômico entre 2,5-24%, dependendo das condições iniciais consideradas, o que pode representar ganhos de até 18 milhões de dólares por ano. Além disso, para essa unidade, verificou-se que o compressor é um equipamento limitante no aumento de ganho econômico do processo.

Dynamics of Boolean networks: an exact solution

The dynamics of Boolean networks (BN) with quenched disorder and thermal noise is studied via the generating functional method. A general formulation, suitable for BN with any distribution of Boolean functions, is developed. It provides exact solutions and insight into the evolution of order parameters and properties of the stationary states, which are inaccessible via existing methodology. We identify cases where the commonly used annealed approximation is valid and others where it breaks down. Broader links between BN and general Boolean formulas are highlighted.