Subtleties in Quantum Mechanical Metastability

Total Page:16

File Type:pdf, Size:1020Kb

Subtleties in Quantum Mechanical Metastability View metadata, citation and similar papers atSubtleties core.ac.uk in Quantum Mechanical Metastability brought to you by CORE provided by CERN Document Server R. M. Cavalcantia) Departamento de F´ısica, Pontif´ıcia Universidade Cat´olica do Rio de Janeiro, Cx. Postal 38071, CEP 22452-970, Rio de Janeiro, RJ, Brazil C. A. A. de Carvalhob) Instituto de F´ısica, Universidade Federal do Rio de Janeiro, Cx. Postal 68528, CEP 21945-970, Rio de Janeiro, RJ, Brazil (April 28, 1997) We present a detailed discussion of some features of quantum mechanical metastability. We analyze the nature of decaying (quasistationary) states and the regime of validity of the exponencial law, as well as decays at finite temperature. We resort to very simple systems and elementary techniques to emphasize subtle aspects of the problem. I. INTRODUCTION method. Why this method works is the question we try to answer in this paper in a very elementary way.3 Thus, Complex-energy eigenfunctions made their d´ebut in in Section II, we show Gamow’s method in action for a Quantum Mechanics through the hands of Gamow, in very simple potential. Some of the results obtained there the theory of alpha-decay.1 Gamow imposed an “outgo- are used in Section III, where we study the time evolu- ing wave boundary condition” on the solutions of the tion of a wave packet initially confined in the potential Schr¨odinger equation for an alpha-particle trapped in the well defined in Section II. This is done with the help of 4 nucleus. Since there is only an outgoing flux of alpha- the propagator built with normalizable eigenfunctions, particles, the wave function ψ(r, t) must behave far from associated to real eigenenergies. As a bonus, we show the nucleus as (for simplicity, we consider an s-wave)2 that the exponential decay law is not valid for very small times or for very large times. This is the content of Sec- e−iEt+ikr tion IV, where the region of validity of the exponential ψ(r, t) ∼ (r →∞). (1) r decay law is roughly delimited. This boundary condition, together with the requirement Another topic we address in this paper is decay at fi- of finiteness of the wave function at the origin, gives rise nite temperature. This is done in Section V, where we to a quantization condition on the values of k (and, there- study, with the help of an exactly solvable toy model, fore, on the values of E = k2/2). It turns out that such how the decay of a metastable system is modified when values are complex: it is coupled to a heat bath. It is shown that, under suit- able conditions, the decay is exponential, with a decay kn = κn − iKn/2,En=n−iΓn/2. (2) rate Γ given by the thermal average of the Γn’s, It follows that −βn nΓne −Γn t+Knr Γ= . (4) 2 e −βn |ψ (r, t)| ∼ (r →∞). (3) n e n r2 P P Thus, if Γn > 0, the probability of finding the alpha- particle in the nucleus decays exponentially in time. The Although this result appears to be rather obvious, in fact lifetime of the nucleus would be given by τ =1/Γ ,and it is not: the decay of a metastable system is an intrinsi- n n cally non-equilibrium process and, so, there is no a priori the energy of the emitted alpha-particle by n. Although very natural, this interpretation suffers from reason for the decay rate to be given by (4). Finally, in some difficulties. How can the energy, which is an ob- Section VI, we discuss the concept of “free energy of a servable quantity, be complex? In other words, how can metastable phase,” a point where we think there is some the Hamiltonian, which is a Hermitean operator, have confusion in the literature. complex eigenvalues? Also, the eigenfunctions are not The results and ideas presented here are not really normalizable, since Γn positive implies Kn positive and, new, but discussions on these matters usually involve 2 therefore, according to (3), |ψn(r, t)| diverges exponen- the use of sophisticated mathematical techniques, such tially with r. as functional5,6 or complex analysis.7–11 For this reason, In spite of such problems (which, in fact, are closely we have tried to make the presentation as clear as possi- related), it is a fact of life that alpha-decay, as well as ble by resorting to very simple systems and elementary other types of decay, does obey an exponential decay techniques — in fact, techniques that can be found in law, with a rate close to that obtained using Gamow’s any standard Quantum Mechanics textbook.12 1 II. DECAYING STATES On the other hand, the derivative of the wave function has a discontinuity at x = a, which can be determined In order to exhibit Gamow’s method in action, we shall by integrating both sides of (7) from a − ε to a + ε,with + study the escape of a particle from the potential well ε→0 : given by: 2λ ϕ0(a)−ϕ0(a)= ϕ (a), (12) 2 1 a 2 V (x)= (λ/a) δ(x − a)forx>0, (5) +∞ for x<0. from which there follows another relation between A and n B: Motion in the region x<0 is forbidden because of the B e−ika cos ka infinite wall at the origin. The positive dimensionless = . (13) constant λ is a measure of the “opacity” of the barrier at A i − 2λ/ka x = a; in the limit λ , the barrier becomes impene- →∞ Combining (11) and (13), we obtain a quantization con- trable, and the energy levels inside the well are quantized. dition for k: If λ is finite, but very large, a particle is no more con- fined to the well, but it usually stays there for a long 2λ cotan ka = i − . (14) time before it escapes. If λ is not so large, the particle ka can easily tunnel through the barrier, and quickly escape The roots of Eq. (14) are complex; when λ 1, those from the potential well. Metastability, therefore, can only which are closest to the origin are given by8 be achieved if the barrier is very opaque, i.e., λ is very large. For this reason, we shall assume this to be the case 1 2nπλ nπ 2 k ≈ −i (nπ λ) (15) in what follows and, whenever possible, we shall retain n a 1+2λ 2λ only the first non-trivial term in a 1/λ expansion. ≡ κn − iKn/2. To find out how fast the particle escapes from the po- tential well, we must solve the Schr¨odinger equation The corresponding eigenenergies are 2 2 3 ∂ 1 ∂ λ 2 1 nπ (nπ) i ψ(x, t)=− ψ(x, t)+ δ(x−a)ψ(x, t). (6) En = kn/2 ≈ − i (16) ∂t 2∂x2 a 2 a (2λa)2 ≡ − iΓ /2. ψ(x, t)=exp(−iEt) ϕ(x) is a particular solution of this n n equation, provided ϕ(x) satisfies the time-independent (Note that Kn κn and Γn n; these results will be Schr¨odinger equation used later.) The imaginary part of En gives rise to an 2 exponential decay of |ψn(x, t)| , with lifetime equal to 1 d2 λ ϕ(x)+ δ(x a)ϕ(x)=Eϕ(x). (7) 2 − 2 − 2(λa) 2 dx a τ =1/Γ ≈ . (17) n n (nπ)3 Denoting the regions 0 <x<aand x>aby the in- dices 1 and 2, respectively, the corresponding wave func- Since the corresponding value of B/A is very small (∼ tions ϕj (x)(j=1,2) satisfy the free-particle Schr¨odinger n/λ), one may be tempted to say that the probability of equation: finding the particle outside the well is negligible in com- parison with the probability of finding the particle inside 2 1 d the well. Normalizing ψ in such a way that the latter − ϕ (x)=Eϕ (x). (8) n 2 dx2 j j equals one when t = 0, the probability of finding the particle inside the well at time t,ifitwereinthen-th Since the wall at the origin is impenetrable, ϕ (0) must 1 decaying state at t =0,wouldbe be zero; the solution of Eq. (8) which obeys this bound- a ary condition is 2 Pn(t)= |ψn(x, t)| dx =exp(−Γnt). (18) √ 0 ϕ1(x)=Asin kx (k = 2E ). (9) Z The trouble with this interpretation is that Im kn ≡ 1,6,13 To determine ϕ2 (x), we follow Gamow’s reasoning −Kn/2 < 0, and so ψn(x, t) diverges exponentially as and require ϕ2(x) to be an outgoing wave. Therefore, we x →∞, since, according to (10), select, from the admissible solutions of Eq. (8), 2 2 |ψn(x, t)| = |Bn| exp(−Γnt + Knx) (19) ikx ϕ2(x)=Be . (10) outside the well. Because of this “exponential catastro- phe”, the decaying states are nonnormalizible and, there- The wave function must be continuous at x = a,sothat fore, cannot be accepted as legitimate solutions of the ϕ(a)=ϕ (a), or 1 2 Schr¨odinger equation (although one can find in the liter- B ature the assertion that they are “rigorous” solutions of = e−ika sin ka. (11) 14 A the time-dependent Schr¨odinger equation ). 2 ∞ 2 III. TIME EVOLUTION OF A WAVE PACKET G(x, x0; t)= e−ik t/2 ϕ∗(k, x) ϕ(k, x0) dk. (27) Z0 We now return to Eq. (8) and write, for the solution The function ϕ(k, x) is the solution of Eq.
Recommended publications
  • Metastability of Anatase
    Metastability of anatase: size dependent and irreversible anatase-rutile phase SUBJECT AREAS: DENDRIMERS transition in atomic-level precise titania POLYMER CHEMISTRY Norifusa Satoh1, Toshio Nakashima2 & Kimihisa Yamamoto2 NANOSCIENCE AND TECHNOLOGY COORDINATION CHEMISTRY 1Photovoltaic Materials Unit, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki, 305-0047 Japan, 2Chemical Resources Laboratory, Tokyo Institute of Technology, 4259 Nagatsuta, Yokohama 226-8503 Japan. Received 12 March 2013 Since crystal phase dominantly affects the properties of nanocrystals, phase control is important for the applications. To demonstrate the size dependence in anatase-rutile phase transition of titania, we used Accepted quantum-size titania prepared from the restricted number of titanium ions within dendrimer templates for 21 May 2013 size precision purposes and optical wave guide spectroscopy for the detection. Contrary to some theoretical calculations, the observed irreversibility in the transition indicates the metastablity of anatase; Published thermodynamics cannot explain the formation of metastable states. Therefore, we take into account the 7 June 2013 kinetic control polymerization of TiO6 octahedral units to explain how the crystal phase of the crystal-nucleus-size titania is dependent on which coordination sites, cis-ortrans-, react in the TiO6 octahedra, suggesting possibilities for the synthetic phase control of nanocrystals. In short, the dendrimer Correspondence and templates give access to crystal nucleation chemistry. The paper will also contribute to the creation of artificial metastable nanostructures with atomic-level precision. requests for materials should be addressed to N.S. (SATOH. he size dependence of nanocrystals during a solid-solid phase transition has attracted much attention1,2, [email protected]) because phase control is a key step to improve functionalities in photovoltaics3,4, ferroelectricity5, magnetics6, 7 or K.Y.
    [Show full text]
  • Multiple Interfaces in Diffusional Phase Transitions in Binary Mesogen-Non-Mesogen Mixtures Undergoing Metastable Phase Separations
    Multiple interfaces in diffusional phase transitions in binary mesogen-non-mesogen mixtures undergoing metastable phase separations Ezequiel R. Soulé 1,2 , Cyrille Lavigne 2,3 , Linda Reven 3 and Alejandro D. Rey 2* 1. Institute of Materials Science and Technology (INTEMA), University of Mar del Plata and National Research Council (CONICET), J. B. Justo 4302, 7600 Mar del Plata, Argentina 2. Department of Chemical Engineering, McGill University, Montreal, Quebec H3A 2B2, Canada 3. Center for Self-Assembled Chemical Structure (CSACS), Chemistry Department, McGill University, 801 Sherbrooke St. West, Montreal, Quebec H3A 2K6, Canada * corresponding author Abstract Theory and simulations of simultaneous chemical demixing and phase ordering are performed for a mixed order parameter system with an isotropic-isotropic (I-I) phase separation that is metastable with respect to an isotropic-nematic (I-N) phase ordering transition. Under certain conditions, the disordered phase transforms into an ordered phase via the motion of a double front containing a metastable phase produced by I-I demixing, a thermodynamically driven mechanism not previously reported. Different kinetic regimes are found depending on the location of the initial conditions in the thermodynamic phase diagram and the ratio between diffusional and nematic phase ordering mobilities. For a diffusional process, depending if the temperature is above or below the critical co-dissolution point, an inflection point or a phase separation takes place in the depletion layer. This phase separation leads to the formation of a second interface where the separation of the two metastable isotropic phases grows monotonically with time. The observed deviations from the typical Fickian concentration profiles are associated with strong positive deviations of the mixture from ideality due to couplings between concentration and nematic ordering.
    [Show full text]
  • A Case Study of N-Doped (Tio2)N for Photocatalysis
    Metastability Triggered Reactivity in Clusters at Realistic Conditions: A Case Study of N-doped (TiO2)n for Photocatalysis Shikha Saini, Pooja Basera, Manish Kumar, Preeti Bhumla, Saswata Bhattacharya∗ Department of Physics, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India E-mail: [email protected][SB] Phone: +91-2659 1359. Fax: +91-2658 2037 Abstract Here we report a strategy, by taking a prototypical model system for photocatal- ysis (viz. N-doped (TiO2)n clusters), to accurately determine low energy metastable structures that can play a major role with enhanced catalytic reactivity. Computa- tional design of specific metastable photocatalyst with enhanced activity is never been easy due to plenty of isomers on potential energy surface. This requires fixing vari- ous parameters viz. (i) favorable formation energy, (ii) low fundamental gap, (iii) low excitation energy and (iv) high vertical electron affinity (VEA) and low vertical ion- arXiv:2006.07395v1 [cond-mat.mtrl-sci] 12 Jun 2020 ization potential (VIP). We validate here by integrating several first principles based methodologies that consideration of the global minimum structure alone can severely underestimate the activity. As a first step, we have used a suite of genetic algorithms [viz. searching clusters with conventional minimum total energy ((GA)E); searching EA IP clusters with specific property i.e. high VEA ((GA)P ), and low VIP ((GA)P )] to 1 model the N-doped (TiO2)n clusters. Following this, we have identified its free energy using ab initio thermodynamics to confirm that the metastable structures are not too far from the global minima.
    [Show full text]
  • Metastable–Solid Phase Diagrams Derived from Polymorphic
    Metastable–solid phase diagrams derived from polymorphic solidification kinetics Babak Sadigha,1 , Luis Zepeda-Ruiza, and Jonathan L. Belofa,1 aLawrence Livermore National Laboratory, Physical and Life Sciences Directorate, Livermore, CA 94550 Edited by Pablo G. Debenedetti, Princeton University, Princeton, NJ, and approved January 16, 2021 (received for review August 24, 2020) Nonequilibrium processes during solidification can lead to kinetic Through X-ray diffraction of the freely suspended droplets, the stabilization of metastable crystal phases. A general frame- dynamics of crystal nucleation during solidification have been work for predicting the solidification conditions that lead to investigated. As a result, solidification of diverse crystalline metastable-phase growth is developed and applied to a model phases such as fcc, bcc, icosahedral, and quasicrystalline has face-centered cubic (fcc) metal that undergoes phase transitions been observed, and emergence of metastable phases, often in to the body-centered cubic (bcc) as well as the hexagonal close- the bcc structure, in sufficiently undercooled liquids has been packed phases at high temperatures and pressures. Large-scale demonstrated. Through application of CNT, the undercooling molecular dynamics simulations of ultrarapid freezing show that necessary for metastable-phase growth has been rationalized. bcc nucleates and grows well outside of the region of its thermo- It is conjectured that phase selection takes place in the nucle- dynamic stability. An extensive study of crystal–liquid equilibria ation stage; the phase emerging is one with the smallest critical confirms that at any given pressure, there is a multitude of nucleation barrier. Consequently, models for solid–liquid inter- metastable solid phases that can coexist with the liquid phase.
    [Show full text]
  • The Metastability of an Electrochemically Controlled
    [28]W. L. Jorgensen, J. Tirado-Rives, J. Am. Chem. Soc. 1988, 110, 1657. guise of bistable [2]rotaxanes in which the ring component can [29]Gaussian 98 (Revision A.11.3), M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. be induced[5] to move relative to the dumbbell-shaped one by Scuseria, M. A. Robb, J. R. Cheeseman, V. G. Zakrzewski, J. A. Montgomery, R. E. Stratmann, J. C. Burant, S. Dapprich, J. M. Millam, A. D. Daniels, K. N. altering the redox characteristics of the molecules. Such Kudin, M. C. Strain, O. Farkas, J. Tomasi, V. Barone, M. Cossi, R. Cammi, B. precisely controllable nanoscale molecular machines and Mennucci, C. Pomelli, C. Adamo, S. Clifford, J. Ochterski, G. A. Petersson, switches have attracted a lot of attention[2, 3] because of their P. Y. Ayala, Q. Cui, K. Morokuma, D. K. Malick, A. D. Rabuck, K. Raghava- potential to meet the expectations of a visionary[6] and to act as chari, J. B. Foresman, J. Cioslowski, J. V. Ortiz, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. Gomperts, R. L. Martin, D. J. Fox, T. some of the smallest components for the engineering of Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, C. Gonzalez, M. nanoelectromechanical systems (NEMs) and the fabrication of Challacombe, P. M. W. Gill, B. G. Johnson, W. Chen, M. W. Wong, J. L. nanoelectronic devices.[7] Andres, M. Head-Gordon, E. S. Replogle, J. A. Pople, Gaussian, Inc., Although the redox-switching properties of numerous bista- Pittsburgh, PA, 2002.
    [Show full text]
  • The Disentangling of Hysteretic Spin Transition, Polymorphism and Metastability in Bistable Thin Films Formed by Sublimation of Bis(Scorpionate) Fe( II ) Molecules O
    The disentangling of hysteretic spin transition, polymorphism and metastability in bistable thin films formed by sublimation of bis(scorpionate) Fe( II ) molecules O. Iasco, Marie-Laure Boillot, A. Bellec, R. Guillot, E. Riviere, S. Mazerat, S. Nowak, D. Morineau, A. Brosseau, F. Miserque, et al. To cite this version: O. Iasco, Marie-Laure Boillot, A. Bellec, R. Guillot, E. Riviere, et al.. The disentangling of hys- teretic spin transition, polymorphism and metastability in bistable thin films formed by sublimation of bis(scorpionate) Fe( II ) molecules. Journal of Materials Chemistry C, Royal Society of Chemistry, 2017, 5 (42), pp.11067-11075. 10.1039/C7TC03276E. hal-01625120 HAL Id: hal-01625120 https://hal.archives-ouvertes.fr/hal-01625120 Submitted on 4 Dec 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. The disentangling of hysteretic spin transition, polymorphism and metastability in bistable thin films formed by sublimation of bis(scorpionate) Fe(II) molecules O. Iasco,a M.-L. Boillot,a* A. Bellec,b R. Guillot,a E. Rivière,a S. Mazerat,a S. Nowak,c D. Morineau,d A. Brosseau,e F. Miserque,F V. Repainb and T.
    [Show full text]
  • Arxiv:2101.05736V1 [Cond-Mat.Mtrl-Sci] 14 Jan 2021
    Metastable piezoelectric group IV monochalcogenide monolayers with a buckled honeycomb structure Shiva P. Poudel1, ∗ and Salvador Barraza-Lopez1, 2, y 1Department of Physics, University of Arkansas, Fayetteville, AR 72701, USA 2Institute for Nanoscience and Engineering, University of Arkansas, Fayetteville, Arkansas 72701, USA (Dated: January 15, 2021) Multiple two-dimensional materials are being na¨ıvely termed stable on the grounds of displaying phonon dispersions with no negative frequencies, and of not collapsing on molecular dynamics calcu- lations at fixed volume. But, if these phases do not possess the smallest possible structural energy, how does one understand and establish their actual meta-stability? To answer this question, twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure{belonging to sym- metry group P3m1{and an out-of-plane intrinsic electric polarization are shown to be metastable by three independendent methods. First, we uncover a coordination-preserving structural trans- formation from the low-buckled honeycomb structure onto the lower-energy Pnm21 (or Pmmn for PbS, PbSe, and PbTe) phase to estimate energy barriers EB that must be overcome during such structural transformation. Using the curvature of the local minima and EB as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (nevertheless, and with the exception of PbS and PbSe, these phases display es- cape times ranging from 700 years to multiple times the age of the universe, and can be considered \stable" for practical purposes only in that relative sense).
    [Show full text]
  • Study of the Phase Transitions in the Binary System NPG-TRIS for Thermal Energy Storage Applications
    materials Article Study of the Phase Transitions in the Binary System NPG-TRIS for Thermal Energy Storage Applications Sergio Santos-Moreno 1,2,3 , Stefania Doppiu 1,*, Gabriel A. Lopez 3 , Nevena Marinova 2, Ángel Serrano 1 , Elena Silveira 2 and Elena Palomo del Barrio 1,4 1 Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, 01510 Vitoria-Gasteiz, Spain; [email protected] (S.S.-M.); [email protected] (Á.S); [email protected] (E.P.d.B.) 2 TECNALIA, Basque Research and Technology Alliance (BRTA), Parque Tecnológico de San Sebastián, 20009 Donostia-San Sebastián, Spain; [email protected] (N.M.); [email protected] (E.S.) 3 Applied Physics II, University of the Basque Country UPV-EHU, 48940 Leioa, Spain; [email protected] 4 Ikerbasque, Basque Foundation for Science, 348013 Bilbao, Spain * Correspondence: [email protected] Received: 11 February 2020; Accepted: 3 March 2020; Published: 5 March 2020 Abstract: Neopentylglycol (NPG) and tris(hydroxymethyl)aminomethane (TRIS) are promising phase change materials (PCMs) for thermal energy storage (TES) applications. These molecules undergo reversible solid-solid phase transitions that are characterized by high enthalpy changes. This work investigates the NPG-TRIS binary system as a way to extend the use of both compounds in TES, looking for mixtures that cover transition temperatures different from those of pure compounds. The phase diagram of NPG-TRIS system has been established by thermal analysis. It reveals the existence of two eutectoids and one peritectic invariants, whose main properties as PCMs (transition temperature, enthalpy of phase transition, specific heat and density) have been determined.
    [Show full text]
  • Nuclear Spin Squeezing in Helium-3 by Continuous Quantum Nondemolition Measurement
    Nuclear spin squeezing in Helium-3 by continuous quantum nondemolition measurement Alan Serafin,1 Matteo Fadel,2 Philipp Treutlein,2 and Alice Sinatra1 1Laboratoire Kastler Brossel, ENS-Universit´ePSL, CNRS, Universit´ede la Sorbonne et Coll`egede France, 24 rue Lhomond, 75231 Paris, France 2Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland (Dated: December 15, 2020) We propose a technique to control the macroscopic collective nuclear spin of a Helium-3 vapor in the quantum regime using light. The scheme relies on metastability exchange collisions to mediate interactions between optically accessible metastable states and the ground-state nuclear spin, giving rise to an effective nuclear spin-light quantum nondemolition interaction of the Faraday form. Our technique enables measurement-based quantum control of nuclear spins, such as the preparation of spin-squeezed states. This, combined with the day-long coherence time of nuclear spin states in Helium-3, opens the possibility for a number of applications in quantum technology. Introduction. The nuclear spin of Helium-3 atoms in x a room-temperature gas is a very well isolated quantum B y z ϕ system featuring record-long coherence times of up to λ several days [1]. It is nowadays used in a variety of appli- κ 2 cations, such as magnetometry [2], gyroscopes for navi- 1083 nm gation [3], as target in particle physics experiments [1], and even in medicine for magnetic resonance imaging of the human respiratory system [4]. Moreover, Helium-3 gas cells are used for precision measurements in funda- mental physics, e.g. in the search for anomalous forces FIG.
    [Show full text]
  • A RELAXATION MODEL for LIQUID-VAPOR PHASE CHANGE with METASTABILITY Francois James, Hélène Mathis
    A RELAXATION MODEL FOR LIQUID-VAPOR PHASE CHANGE WITH METASTABILITY Francois James, Hélène Mathis To cite this version: Francois James, Hélène Mathis. A RELAXATION MODEL FOR LIQUID-VAPOR PHASE CHANGE WITH METASTABILITY. Communications in Mathematical Sciences, International Press, 2016, 74 (8), pp.2179-2214. 10.4310/CMS.2016.v14.n8.a4. hal-01178947v2 HAL Id: hal-01178947 https://hal.archives-ouvertes.fr/hal-01178947v2 Submitted on 14 Mar 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A RELAXATION MODEL FOR LIQUID-VAPOR PHASE CHANGE WITH METASTABILITY FRANC¸OIS JAMES AND HEL´ ENE` MATHIS Abstract. We propose a model that describes phase transition including metastable states present in the van der Waals Equation of State. From a convex optimization problem on the Helmoltz free energy of a mixture, we deduce a dynamical system that is able to depict the mass transfer between two phases, for which equilibrium states are either metastable states, stable states or a coexistent state. The dynamical system is then used as a relaxation source term in an isothermal 4×4 two-phase model. We use a Finite Volume scheme that treats the convective part and the source term in a fractional step way.
    [Show full text]
  • Energy Transfer from PO Excited States to Alkali Metal Atoms in the Phosphorus Chemiluminescence Flame [Metastable P0(4Jj,)/(P0.P0)* Excimerl AHSAN U
    Proc. Natl. Acad. Sci. USA Vol. 77, No. 12, pp. 6952-6955, December 1980 Chemistry Energy transfer from PO excited states to alkali metal atoms in the phosphorus chemiluminescence flame [metastable P0(4jj,)/(P0.P0)* excimerl AHSAN U. KHAN Institute of Molecular Biophysics and Department of Chemistry, Florida State University, Tallahassee, Florida 32306 Communicated by Michael Kasha, August 21, 1980 ABSTRACT Phosphorus chemiluminescence under ambient cence is a visible continuum upon which are superposed a conditions of a phosphorus oxidation flame is found to offer an number of sharp emission bands. In 1938 Rumpf (9) confirmed efficient electronic energy transferring system to alkali metal Ball regarding the ultraviolet atoms. The lowest resonance lines, 2P312,12-.S1/2, of potas- the observations of Ghosh and sium and sodium are excited by energy transfer when an argon emission from the reaction and attributed the visible emission stream at 800C carrying potassium or sodium atoms intersects also to PO. In 1957 Walsh (10) made a detailed investigation a phosphorus vapor stream, either at the flame or'in the post- of the ultraviolet bands originating from the A22+ and B2Z* flame region. The lowest electronically excited metastable 4fl; states of PO and suggested that an unknown 2z+ state was re- state of PO or the (P0OP)* excimer is considered to be the sponsible for the visible emission; In 1965 Cordes and Witschel probable energy donor. The (PO.;PO)* excimer resilts from the on the interaction of the "ll; stite of one PO molecule with the ground (11) pointed out that the sharp band systems superposed 211r state of another.
    [Show full text]
  • Clusters, Metastability, and Nucleation: Kinetics of First-Order Phase Transitions M
    Journal of Statistical Physics, VoL 18, No. 1, 1978 Clusters, Metastability, and Nucleation: Kinetics of First-Order Phase Transitions M. Kalos, 1 Joel L. Lebowitz, 2'3'40. Penrose, 5 and A. Sur 6 Received July 21, 1977 We describe and interpret computer simulations of the time evolution of a binary alloy on a cubic lattice, with nearest neighbor interactions favoring like pairs of atoms. Initially the atoms are arranged at random; the time evolution proceeds by random interchanges of nearest neighbor pairs, using probabilities compatible with the equilibrium Gibbs distribution at tem- perature T. For temperatures 0.59T~, 0.81To, and 0.89Tc, with density p of A atoms equal to that in the B-rich phase at coexistence, the density C~ of clusters of l A atoms approximately satisfies the following empirical formulas: C1 ~ w(1 - p)a and C~ ,.~ (1 - p)~Q,w ~ (2 ~< l ~< 10). Here w is a parameter and we define Q~ = ~K e-~K~, where the sum goes over all translationally nonequivalent /-particle clusters and E(K) is the energy of formation of the cluster K. For l > 10, Qz is not known exactly; so we use an extrapolation formula Q~ ~ AwZ~l -~ exp(-bP), where w~ is the value of w at coexistence. The same formula (with w > w~) also fits the observed values of C~ (for small values of l) at densities greater than the coexistence density (for T = 0.59T~): When the supersaturation is small, the simula- tions show apparently metastable states, a theoretical estimate of whose lifetime is compatible with the observations.
    [Show full text]