Melting Line of the Lennard-Jones System, Infinite Size, and Full Potential

Total Page:16

File Type:pdf, Size:1020Kb

Melting Line of the Lennard-Jones System, Infinite Size, and Full Potential THE JOURNAL OF CHEMICAL PHYSICS 127, 104504 ͑2007͒ Melting line of the Lennard-Jones system, infinite size, and full potential ͒ ͒ Ethan A. Mastnya and Juan J. de Pablob Chemical and Biological Engineering Department, University of Wisconsin-Madison, 1415 Engineering Drive, Madison, Wisconsin 53706–1691, USA ͑Received 9 February 2007; accepted 5 June 2007; published online 13 September 2007͒ Literature estimates of the melting curve of the Lennard-Jones system vary by as much as 10%. The origin of such discrepancies remains unclear. We present precise values for the Lennard-Jones melting temperature, and we examine possible sources of systematic errors in the prediction of melting points, including finite-size and interaction-cutoff effects. A hypothetical thermodynamic integration path is used to find the relative free energies of the solid and liquid phases, for various system sizes, at constant cutoff radius. The solid-liquid relative free energy and melting temperature scale linearly as the inverse of the number of particles, and it is shown that finite-size effects can account for deviations in the melting temperature ͑from the infinite-size limit͒ of up to 5%. An extended-ensemble density-of-states method is used to determine free energy changes for each phase as a continuous function of the cutoff radius. The resulting melting temperature predictions exhibit an oscillatory behavior as the cutoff radius is increased. Deviations in the melting temperature ͑from the full potential limit͒ arising from a finite cutoff radius are shown to be of comparable magnitude as those resulting from finite-size effects. This method is used to identify melting temperatures at five different pressures, for the infinite-size and full potential Lennard-Jones system. We use our simulation results as references to connect the Lennard-Jones solid equation of state of van der Hoef with the Lennard-Jones fluid equation of state of Johnson. Once the references are applied the two equations of state are used to identify a melting curve. An empirical equation that fits this melting curve is provided. We also report a reduced triple point temperature ͓ ͔ Ttr=0.694. © 2007 American Institute of Physics. DOI: 10.1063/1.2753149 I. INTRODUCTION proaches unity. Thermodynamic properties are subsequently corrected for the truncated portion of the interaction potential The Lennard-Jones system is arguably one of the sim- by assuming a value of unity for the radial distribution func- plest models capable of reproducing the complete thermody- tion at distances greater than rc. For crystalline solids, how- namic behavior of classical fluids. Over the last two decades, ever, the radial distribution does not approach unity in an a consensus has emerged over the precise coordinates of its asymptotic manner for short to intermediate length scales. It vapor-liquid coexistence ͑or binodal͒ curve and its critical is therefore of interest to determine the nature and magnitude point. In contrast, available estimates of the melting curve of the systematic errors that arise in simulations of solid- vary considerably. Simulations of solid-liquid equilibria are liquid equilibria as a result of finite-size and cutoff-radius considerably more challenging than those of vapor-liquid effects. equilibria. It is unclear whether discrepancies between differ- Theoretical values for the melting curve of the Lennard- ent results are due to methodological challenges or to finite- Jones system are not available. A definitive standard for size or finite-cutoff radius effects. The aim of this article is to comparison of the results of molecular simulations does not evaluate the magnitude of such effects and to present results exist. One must therefore rely on somewhat subjective as- for the melting curve of Lennard-Jonesium in the thermody- sessments of existing calculations to determine the accuracy namic limit ͑N→ϱ, where N is the number of particles in the of a particular value of the melting point. Published esti- system͒. mates of melting points of the Lennard-Jones system vary by Molecular simulations of simple or complex fluids are as much as 10%. The work presented in this article is moti- generally conducted on systems consisting of a few hundred vated by a desire to establish a precise standard for the melt- or thousand particles. Such simulations typically adopt a ing temperature of Lennard-Jonesium. We examine in con- truncated interaction potential, whose value is assumed to siderable detail some of the more common sources of vanish beyond a specified cutoff radius. For the Lennard- systematic errors in the prediction of melting points, namely, ␴ Jones system, that cutoff radius is often set to rc =2.5 , finite-size and cutoff radius effects, in the hope that the where ␴ is the size or length scale associated with the analysis provided in this work will permit more objective ␴ Lennard-Jones potential energy function. At rc =2.5 , the ra- consideration of melting point estimates obtained by means dial distribution function g͑r͒ of a Lennard-Jones liquid ap- of molecular simulation. The melting points presented in this work are generated ͒ 1,2 a Electronic mail: [email protected] by thermodynamic integration along a hypothetical path, ͒ b Electronic mail: [email protected] that integration yields the relative free energies of the solid 0021-9606/2007/127͑10͒/104504/8/$23.00127, 104504-1 © 2007 American Institute of Physics Reuse of AIP Publishing content is subject to the terms: https://publishing.aip.org/authors/rights-and-permissions. Downloaded to IP: 128.112.35.141 On: Fri, 12 Feb 2016 14:58:30 104504-2 E. A. Mastny and J. J. de Pablo J. Chem. Phys. 127, 104504 ͑2007͒ and liquid phases, of various system sizes, at constant cutoff simulating 108-, 256-, and 500-particle systems. He also radius. The solid-liquid relative free energy and melting tem- considered two cutoff schemes, namely, a standard long- perature are shown to scale as the inverse size of the system, range correction beyond half the box length and a perfect- 1/N. It is also shown that finite-size effects from this inte- lattice correction for the solid phase beyond half the box gration path can account for melting temperature errors of up length. He noted that the coexistence pressures did not ex- to 5%. An extended-ensemble density-of-states method is hibit any apparent scaling relationship for either type of cor- then used to determine the free energy change in each phase rection, and he was therefore unable to extrapolate his results as a continuous function of the cutoff radius. Our results to infinite system size. indicate that the melting temperature exhibits oscillations In a previous study,10 we used a density-of-states Monte whose magnitude decreases as the cutoff radius is increased. Carlo sampling method with order parameters of potential Melting temperature prediction errors due to a finite cutoff energy and volume. In this work, neither finite-size nor cut- radius can also amount to 5% of the value corresponding to off effects were investigated ͑500 particles and 2.5␴ were the full Lennard-Jones potential. used͒. The results of this work are also sensitive to the We use the relative free energies determined from the method used for matching the free energy functions in the thermodynamic integrations and corrected for large cutoff energy and volume space between the equilibrium solid and radius as references to connect the Lennard-Jones solid equa- liquid regions. tion of state of van der Hoef3 with the Lennard-Jones fluid More recently, McNeil-Watson and Wilding11 also used equation of state of Johnson.4 Upon applying the references the phase-switch Monte Carlo method described by Err- the two equations of state were used to identify a melting ington to predict melting points. They then reweighted their curve. data to obtain a portion of the melting line. Their results were quantitatively similar to those of Errington. They explained that the reason they did not see a 1/N scaling of the melting II. LITERATURE REVIEW temperature was due to the fact that the choice of a “fluctu- In 1969 Hansen and Verlet reported the first melting ating” cutoff at one-half the box length introduces a coupling point estimates of the Lennard-Jones system.5 They evalu- between the cutoff and system size. ated the free energy of the solid and liquid phases indepen- The results of Errington’s and McNeil-Watson and Wild- dently, and calculated coexistence properties at four different ing do not scale to infinite-size systems, and it is unclear temperatures. The thermodynamic integration path used to whether a cutoff radius that changed with system size had a calculate the solid free energy, however, includes a phase significant effect. The 1/N scaling analysis presented by transition. these authors compared different systems because the cutoff Agrawal and Kofke6,7 carried out a subsequent study of was different. One would not expect different systems to the Lennard-Jones melting curve. They used a Gibbs-Duhem have any predictable scaling behavior. van der Hoef published an equation of state for the integration to track the melting point as the potential of in- 3 teraction was first changed from purely repulsive ͑hard Lennard-Jones solid. In his study almost 900 NVT simula- sphere͒ to the repulsive part of the Lennard-Jones potential tions in the temperature range of 0.1–2.0 and density range energy function, and then to the complete Lennard-Jones en- of 0.94–1.2 were performed. He used 2048 particles and a ␴ ergy function. The starting point for these calculations was a cutoff radius of 6.0 . The measured internal energies and single melting point corresponding to a system of hard pressures from these simulations were fitted to an empirical spheres.
Recommended publications
  • LATENT HEAT of FUSION INTRODUCTION When a Solid Has Reached Its Melting Point, Additional Heating Melts the Solid Without a Temperature Change
    LATENT HEAT OF FUSION INTRODUCTION When a solid has reached its melting point, additional heating melts the solid without a temperature change. The temperature will remain constant at the melting point until ALL of the solid has melted. The amount of heat needed to melt the solid depends upon both the amount and type of matter that is being melted. Therefore: Q = ML Eq. 1 f where Q is the amount of heat absorbed by the solid, M is the mass of the solid that was melted and Lf is the latent heat of fusion for the type of material that was melted, which is measured in J/kg, NOTE: to fuse means to melt. In this experiment the latent heat of fusion of water will be determined by using the method of mixtures ΣQ=0, or QGained +QLost =0. Ice will be added to a calorimeter containing warm water. The heat energy lost by the water and calorimeter does two things: 1. It melts the ice; 2. It warms the water formed by the melting ice from zero to the final equilibrium temperature of the mixture. Heat gained + Heat lost = 0 Heat needed to melt ice + Heat needed to warm the melt water + Heat lost by warn water = 0 MiceLf + MiceCw (Tf - 0) + MwCw (Tf – Tw) = 0 Eq. 2. * Note: The mass of the melted water is the same as the mass of the ice. where M = mass of warm water initially in calorimeter w Mice = mass of ice and water from the melted ice Cw = specific heat of water Lf = latent heat of fusion of water Tw = initial temperature water Tf = equilibrium temperature of mixture APPARATUS: Calorimeter, thermometer, balance, ice, water OBJECTIVE: To determine the latent heat of fusion of water PROCEDURE: 1.
    [Show full text]
  • Determination of the Identity of an Unknown Liquid Group # My Name the Date My Period Partner #1 Name Partner #2 Name
    Determination of the Identity of an unknown liquid Group # My Name The date My period Partner #1 name Partner #2 name Purpose: The purpose of this lab is to determine the identity of an unknown liquid by measuring its density, melting point, boiling point, and solubility in both water and alcohol, and then comparing the results to the values for known substances. Procedure: 1) Density determination Obtain a 10mL sample of the unknown liquid using a graduated cylinder Determine the mass of the 10mL sample Save the sample for further use 2) Melting point determination Set up an ice bath using a 600mL beaker Obtain a ~5mL sample of the unknown liquid in a clean dry test tube Place a thermometer in the test tube with the sample Place the test tube in the ice water bath Watch for signs of crystallization, noting the temperature of the sample when it occurs Save the sample for further use 3) Boiling point determination Set up a hot water bath using a 250mL beaker Begin heating the water in the beaker Obtain a ~10mL sample of the unknown in a clean, dry test tube Add a boiling stone to the test tube with the unknown Open the computer interface software, using a graph and digit display Place the temperature sensor in the test tube so it is in the unknown liquid Record the temperature of the sample in the test tube using the computer interface Watch for signs of boiling, noting the temperature of the unknown Dispose of the sample in the assigned waste container 4) Solubility determination Obtain two small (~1mL) samples of the unknown in two small test tubes Add an equal amount of deionized into one of the samples Add an equal amount of ethanol into the other Mix both samples thoroughly Compare the samples for solubility Dispose of the samples in the assigned waste container Observations: The unknown is a clear, colorless liquid.
    [Show full text]
  • Single Crystal Growth of Mgb2 by Using Low Melting Point Alloy Flux
    Single Crystal Growth of Magnesium Diboride by Using Low Melting Point Alloy Flux Wei Du*, Xugang Xian, Xiangjin Zhao, Li Liu, Zhuo Wang, Rengen Xu School of Environment and Material Engineering, Yantai University, Yantai 264005, P. R. China Abstract Magnesium diboride crystals were grown with low melting alloy flux. The sample with superconductivity at 38.3 K was prepared in the zinc-magnesium flux and another sample with superconductivity at 34.11 K was prepared in the cadmium-magnesium flux. MgB4 was found in both samples, and MgB4’s magnetization curve being above zero line had not superconductivity. Magnesium diboride prepared by these two methods is all flaky and irregular shape and low quality of single crystal. However, these two metals should be the optional flux for magnesium diboride crystal growth. The study of single crystal growth is very helpful for the future applications. Keywords: Crystal morphology; Low melting point alloy flux; Magnesium diboride; Superconducting materials Corresponding author. Tel: +86-13235351100 Fax: +86-535-6706038 Email-address: [email protected] (W. Du) 1. Introduction Since the discovery of superconductivity in Magnesium diboride (MgB2) with a transition temperature (Tc) of 39K [1], a great concern has been given to this material structure, just as high-temperature copper oxide can be doped by elements to improve its superconducting temperature. Up to now, it is failure to improve Tc of doped MgB2 through C [2-3], Al [4], Li [5], Ir [6], Ag [7], Ti [8], Ti, Zr and Hf [9], and Be [10] on the contrary, the Tc is reduced or even disappeared.
    [Show full text]
  • Fast Synthesis of Fe1.1Se1-Xtex Superconductors in a Self-Heating
    Fast synthesis of Fe1.1Se1-xTe x superconductors in a self-heating and furnace-free way Guang-Hua Liu1,*, Tian-Long Xia2,*, Xuan-Yi Yuan2, Jiang-Tao Li1,*, Ke-Xin Chen3, Lai-Feng Li1 1. Key Laboratory of Cryogenics, Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China 2.Department of Physics, Beijing Key Laboratory of Opto-electronic Functional Materials & Micro-nano Devices, Renmin University of China, Beijing 100872, P. R. China 3. State Key Laboratory of New Ceramics & Fine Processing, School of Materials Science and Engineering, Tsinghua University, Beijing 100084, P. R. China *Corresponding authors: Guang-Hua Liu, E-mail: [email protected] Tian-Long Xia, E-mail: [email protected] Jiang-Tao Li, E-mail: [email protected] Abstract A fast and furnace-free method of combustion synthesis is employed for the first time to synthesize iron-based superconductors. Using this method, Fe1.1Se1-xTe x (0≤x≤1) samples can be prepared from self-sustained reactions of element powders in only tens of seconds. The obtained Fe1.1Se1-xTe x samples show clear zero resistivity and corresponding magnetic susceptibility drop at around 10-14 K. The Fe1.1Se0.33Te 0.67 sample shows the highest onset Tc of about 14 K, and its upper critical field is estimated to be approximately 54 T. Compared with conventional solid state reaction for preparing polycrystalline FeSe samples, combustion synthesis exhibits much-reduced time and energy consumption, but offers comparable superconducting properties. It is expected that the combustion synthesis method is available for preparing plenty of iron-based superconductors, and in this direction further related work is in progress.
    [Show full text]
  • Problem Set #10 Assigned November 8, 2013 – Due Friday, November 15, 2013 Please Show All Work for Credit
    Problem Set #10 Assigned November 8, 2013 – Due Friday, November 15, 2013 Please show all work for credit To Hand in 1. 1 2. –1 A least squares fit of ln P versus 1/T gives the result Hvaporization = 25.28 kJ mol . 3. Assuming constant pressure and temperature, and that the surface area of the protein is reduced by 25% due to the hydrophobic interaction: 2 G 0.25 N A 4 r Convert to per mole, ↓determine size per molecule 4 r 3 V M N 0.73mL/ g 60000g / mol (6.02 1023) 3 2 2 A r 2.52 109 m 2 23 9 2 G 0.25 N A 4 r 0.0720N / m 0.25 6.02 10 / mol (4 ) (2.52 10 m) 865kJ / mol We think this is a reasonable approach, but the value seems high 2 4. The vapor pressure of an unknown solid is approximately given by ln(P/Torr) = 22.413 – 2035(K/T), and the vapor pressure of the liquid phase of the same substance is approximately given by ln(P/Torr) = 18.352 – 1736(K/T). a. Calculate Hvaporization and Hsublimation. b. Calculate Hfusion. c. Calculate the triple point temperature and pressure. a) Calculate Hvaporization and Hsublimation. From Equation (8.16) dPln H sublimation dT RT 2 dln P d ln P dT d ln P H T 2 sublimation 11 dT dT R dd TT For this specific case H sublimation 2035 H 16.92 103 J mol –1 R sublimation Following the same proedure as above, H vaporization 1736 H 14.43 103 J mol –1 R vaporization b.
    [Show full text]
  • Phase Diagrams a Phase Diagram Is Used to Show the Relationship Between Temperature, Pressure and State of Matter
    Phase Diagrams A phase diagram is used to show the relationship between temperature, pressure and state of matter. Before moving ahead, let us review some vocabulary and particle diagrams. States of Matter Solid: rigid, has definite volume and definite shape Liquid: flows, has definite volume, but takes the shape of the container Gas: flows, no definite volume or shape, shape and volume are determined by container Plasma: atoms are separated into nuclei (neutrons and protons) and electrons, no definite volume or shape Changes of States of Matter Freezing start as a liquid, end as a solid, slowing particle motion, forming more intermolecular bonds Melting start as a solid, end as a liquid, increasing particle motion, break some intermolecular bonds Condensation start as a gas, end as a liquid, decreasing particle motion, form intermolecular bonds Evaporation/Boiling/Vaporization start as a liquid, end as a gas, increasing particle motion, break intermolecular bonds Sublimation Starts as a solid, ends as a gas, increases particle speed, breaks intermolecular bonds Deposition Starts as a gas, ends as a solid, decreases particle speed, forms intermolecular bonds http://phet.colorado.edu/en/simulation/states- of-matter The flat sections on the graph are the points where a phase change is occurring. Both states of matter are present at the same time. In the flat sections, heat is being removed by the formation of intermolecular bonds. The flat points are phase changes. The heat added to the system are being used to break intermolecular bonds. PHASE DIAGRAMS Phase diagrams are used to show when a specific substance will change its state of matter (alignment of particles and distance between particles).
    [Show full text]
  • Chapter 10: Melting Point
    Chapter 10: Melting Point The melting point of a compound is the temperature or the range of temperature at which it melts. For an organic compound, the melting point is well defined and used both to identify and to assess the purity of a solid compound. 10.1 Compound Identification An organic compound’s melting point is one of several physical properties by which it is identified. A physical property is a property that is intrinsic to a compound when it is pure, such as its color, boiling point, density, refractive index, optical rotation, and spectra (IR, NMR, UV-VIS, and MS). A chemist must measure several physical properties of a compound to determine its identity. Since melting points are relatively easy and inexpensive to determine, they are handy identification tools to the organic chemist. The graph in Figure 10-1 illustrates the ideal melting behavior of a solid compound. At a temperature below the melting point, only solid is present. As heat is applied, the temperature of the solid initially rises and the intermolecular vibrations in the crystal increase. When the melting point is reached, additional heat input goes into separating molecules from the crystal rather than raising the temperature of the solid: the temperature remains constant with heat input at the melting point. When all of the solid has been changed to liquid, heat input raises the temperature of the liquid and molecular motion within the liquid increases. Figure 10-1: Melting behavior of solid compounds. Because the temperature of the solid remains constant with heat input at the transition between solid and liquid phases, the observer sees a sharp melting point.
    [Show full text]
  • The Melting Point Experimental Procedure
    Melting Point1 The temperature at which a solid melts is known as the melting point (MP) of that substance. The melting point is a physical property of a solid and can be used to help identify a substance. In practice, a solid usually melts over a range of temperatures rather than at one specific temperature. For this reason, it is more useful to speak of a melting point range. Although the term "melting point" is usually used, what is meant by this term is "melting point range." If the compound melts over a very narrow range, it can usually be assumed that the compound is relatively pure. Conversely, compounds that melt over a wide range are assumed to be relatively impure. Besides melting over a wide range, impure solids also melt at a temperature lower than that of the pure compound. For our purposes, a range greater than 2 °C is considered to be wide. For example, if an unknown solid melts at 102 – 106 °C, the 4 °C range suggests that the sample is impure. If the unknown is one of four possible compounds which melt at 102, 104, 106, and 108 °C, it is most likely that which melts at 108 °C. To summarize, an impure solid melts over a wide range and at a temperature lower than that of the pure solid. It should be noted that “insoluble” impurities such as bits of filter paper or dust have no effect on the MP of a substance. To impact the MP the impurity must be soluble in the solid.
    [Show full text]
  • 21/04/2021 1.Write the Different Units of Temperature. Commo
    Class: IX Topic: Matter in our surrounding Subject: Chemistry Date -21/04/2021 1.Write the different units of temperature. Common unit: Degree Celsius (oC) SI unit: Kelvin (K) Relation between common unit and SI unit of temperature: 0 oC = 273K 2. How can the different states of matter be changed? Physical states of matter can be inter converted into each other by following two ways: 1. By changing the temperature 2. By changing the pressure 1. Effect of Change of Temperature: a.Solid to liquid: • On increasing the temperature of solids, the kinetic energy of the particles increases which overcomes the forces of attraction between the particles thereby solid gets converted to a liquid. b. Liquid to gas: On heating a liquid like water, the kinetic energy of its particles increases as high as in a gas, thus causing the liquid to change to a gas. c. Gas to liquid: On cooling a gas like steam (or water vapour), the kinetic energy of its particles is lowered down, causing them to move slowly and bringing them closer, forming a liquid. d. Liquid to solid: When a liquid is cooled down by lowering its temperature, its particles lose the kinetic energy and come to a stationary position, start attracting each other more strongly causing the liquid to turn to soilid. 2. Effect of change of pressure Gas to liquid: Gases can be liquefied by applying pressure and reducing the temperature. When a high pressure is applied to a gas, it gets compressed and if the temperature is lowered, the gas is liquefied.
    [Show full text]
  • (PCM) Candidates for Thermal Energy Storage (TES) Applications Judith C
    High-Temperature Phase Change Materials (PCM) Candidates for Thermal Energy Storage (TES) Applications Judith C. Gomez NREL is a national laboratory of the U.S. Department of Energy, Office of Energy Efficiency & Renewable Energy, operated by the Alliance for Sustainable Energy, LLC. Milestone Report NREL/TP-5500-51446 September 2011 Contract No. DE-AC36-08GO28308 High-Temperature Phase Change Materials (PCM) Candidates for Thermal Energy Storage (TES) Applications Judith C. Gomez Prepared under Task No. CP09.2201 NREL is a national laboratory of the U.S. Department of Energy, Office of Energy Efficiency & Renewable Energy, operated by the Alliance for Sustainable Energy, LLC. National Renewable Energy Laboratory Milestone Report 1617 Cole Boulevard NREL/TP-5500-51446 Golden, Colorado 80401 September 2011 303-275-3000 • www.nrel.gov Contract No. DE-AC36-08GO28308 NOTICE This report was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or any agency thereof.
    [Show full text]
  • States, Boiling Point, Melting Point, and Solubility
    States, Boiling Point, Melting Point, and Solubility Solid Liquid Gas Defining States of Matter ● States of matter are defined by whether they hold SHAPE and VOLUME Element (Au) Compound (NaCl) Mixture (Milk, Salt, etc) ALL KEEP THE SAME SHAPE AND VOLUME = Solids Particle View of a Solid ● Particles in a solid are PACKED CLOSELY together and they are in a FIXED POSITION. Particles vibrate in place Liquids ● Liquids – has definite VOLUME but no defined SHAPE 100 ml Particle View of a Liquid ● Packed CLOSELY (like a solid), but move FREELY around each other (must stay in contact). Gases ● Gases - do NOT have definite SHAPE or VOLUME. Bromine gas fills up the entire volume of the container. Particle view of a Gas ● Particles can MOVE FREELY and will either fill up or squeeze into available space. Task ● Draw a diagram of ● A) Gas particles ● B) Liquid particles ● C) Solid particles Changes in States of Matter ● Thermal Energy – heat energy. ● More thermal energy = More particle movement Changing States Increase Thermal Energy (Heat up) point Solid Liquid Gas Melting Melting Boiling Point Decrease Thermal Energy (Cool off) Melting point ● Melting - change from solid to liquid ● Melting point - SPECIFIC temperature when melting occurs. ● Each pure substance has a SPECIFIC melting point. ● Examples: ● M.P. of Water = 0°C (32°F) ● M.P. of Nitrogen = -209.9 °C (-345.81998 °F) ● M.P. of Silver = 961.93 °C (1763.474 °F) ● M.P. of Carbon = 3500.0 °C (6332.0 °F) Melting Point ● Particles of a solid vibrate so fast that they break free from their fixed positions.
    [Show full text]
  • (2020) Origin of Superconductivity at Nickel-Bismuth Interfaces
    PHYSICAL REVIEW RESEARCH 2, 013270 (2020) Origin of superconductivity at nickel-bismuth interfaces Matthew Vaughan ,1 Nathan Satchell ,1 Mannan Ali,1 Christian J. Kinane,2 Gavin B. G. Stenning,2 Sean Langridge ,2 and Gavin Burnell 1,* 1School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom 2ISIS Neutron and Muon Source, Rutherford Appleton Laboratory, Harwell, Oxon OX11 0QX, United Kingdom (Received 22 November 2019; accepted 4 February 2020; published 6 March 2020) Unconventional superconductivity has been suggested to be present at the interface between bismuth and nickel in thin-film bilayers. In this work, we study the structural, magnetic, and superconducting properties of sputter deposited Bi/Ni bilayers. As-grown, our films do not display a superconducting transition; however, when stored at room temperature, after about 14 days our bilayers develop a superconducting transition up to 3.8 K. To systematically study the effect of low temperature annealing on our bilayers, we perform structural characterization with x-ray diffraction and polarized neutron reflectometry, along with magnetometry and low-temperature electrical transport measurements on samples annealed at 70◦ C. We show that the onset of superconductivity in our samples is coincident with the formation of ordered NiBi3 intermetallic alloy, a known s-wave superconductor. We calculate that the annealing process has an activation energy of (0.86 ± 0.06) eV. As a consequence, gentle heating of the bilayers will cause formation of the superconducting NiBi3 at the Ni/Bi interface, which poses a challenge to studying any distinct properties of Bi/Ni bilayers without degrading that interface.
    [Show full text]