Structure–Property of Lithium–Sulfur Nanoparticles via Molecular

Oct 9, 2018 - Lithium–sulfur (Li–S) batteries offer higher energy densities than most reported lithium-ion batteries. However, our understanding o...
0 downloads 0 Views 12MB Size
Subscriber access provided by UNIV OF LOUISIANA

Surfaces, Interfaces, and Applications

Structure-Property of Lithium-Sulfur Nanoparticles via Molecular Dynamics Simulation Ying Li, Nichols A. Romero, and Kah Chun Lau ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.8b09128 • Publication Date (Web): 09 Oct 2018 Downloaded from http://pubs.acs.org on October 11, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Structure-Property of Lithium-Sulfur Nanoparticles via Molecular Dynamics Simulation Ying Li,a,b,* Nichols A. Romero,a,b Kah Chun Lauc,* a

Computational Science Division, Argonne National Laboratory, IL 60439, USA Leadership Computing Facility, Argonne National Laboratory, IL 60439, USA c Department of Physics & Astronomy, California State University Northridge, CA 91330, USA * [email protected]; [email protected] b

Keywords: Molecular dynamics, Nanoparticles, Lithium sulfur materials, ReaxFF force field, Structural property Abstract: Lithium-sulfur (Li-S) batteries offer higher energy densities than most reported lithium-ion batteries. However, our understanding of Li-S battery is still largely unknown at the level of the nanoscale. The structural properties of Li-S materials were investigated via molecular dynamics (MD) simulations using the ReaxFF force field. everal Li-S nanoparticles with different Li:S composition ratios (2:1 and 2:8) and various structures are studied. Our MD simulations show that, among the four structures we constructed for Li2S8 nanoparticles, the core-shell structure is the most thermodynamically stable one during the charging (delithiation) process. In contrast to bulk crystal Li2S, we find the presence of mixed lithium sulfide and polysulfide species are common features for these Li-S (Li2S, Li2S8) nanoparticles. The complex distribution of these sulfide and polysulfide speciation are dictated by both stoichiometry and local atomic structures in the nanoparticle. These findings will provide insight into further development of functionalized lithium-sulfur cathodes. 1. Introduction As the demand for electrical energy continues to increase due to economic and population growth, the development of robust electrical energy storage technologies (such as batteries) is urgent to store sustainable energy resources (e.g., solar, wind) in a safe and environmentally friendly manner.1-3 To catch up to the efficiency of energy production with conventional fossil fuels technologies, a substantial enhancement of the energy density for rechargeable batteries is required.4 This requirement has attracted much attention to a new rechargeable battery design based on lithium-sulfur (Li-S) 5-6. The energy density of Li-S is several times higher than conventional Li-ion batteries.3,7 However, to date, the advances in Li-S battery technologies are hindered by low rechargeability, rapid capacity fading, polysulfide dissolution,8 insufficient control and lack of basic understandings of the electrode/electrolyte interface.5,9 Thus, fundamental and systematic studies of the structure and properties of Li-S batteries and related materials are critical for the development of a practical rechargeable Li-S battery. One of the serious challenges in Li-S batteries is the design of a robust cathode.10 Lithium sulfide (Li2S) is currently under intense study because it is the main insulating discharging product of Li-S batteries and an interesting candidate cathode material.11-13 When Li2S solid is employed as the cathode material, the mechanical vulnerability of the cathode due to volume expansion of sulfur particles caused by the lithiation in the discharging process can be reduced significantly in contrast to a pure sulfur cathode.3,14-16 The advantages of using Li2S as the cathode material are obvious, but there are still lots of fundamental problems waiting to be solved. One of the most common issue is that there are many higher order polysulfide (Li2Sn, n > 2) species found together with the Li2S species.17 In contrast to Li2S, the chemical properties of the polysulfide (Li2Sn, n > 2) are erratic due to the complicated molecular structures and unpredictable (electro)chemical reactions.18-20 Thus resulting in wide distributions of neutral or anionic Li2Sx, or radical LiSx species in electrolyte solutions.21-22 In attempt to improve the electronic conductivity and suppress polysulfide dissolution at cathode, methods of confining sulfur within the cathodes will continue to be an active research area.23-24 To fabricate composite cathodes for lithium–sulfur batteries (depicted in figure 1), ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 23

fundamental studies of the structural properties, speciation of lithium sulfide, polysulfide species, and the nanoscale of Li-S system are needed. To investigate the basic structural property of these lithium sulfide and polysulfide species at the nanoscale, it is important to study these systems (i.e., lithium sulfide and polysulfide species) in their nanoparticle morphology by focusing solely on nanoparticles (Fig. 1) using atomistic simulation. This in turn will provide us detailed molecular species characterization and spatial distributions of lithium sulfide and polysulfide species that generally present in Li-S battery. Thus, the study here is an important baseline study for future investigations of functionalized lithium-sulfur cathodes. 2. Method and Model We perform MD simulations using LAMMPS25 with ReaxFF potential,26 where the interaction of LiLi, Li-S and S-S atoms are based on the interatomic force field parameters proposed by van Duin et. al.27 The reactive force-field (ReaxFF) interatomic potential is a powerful computational tool for exploring material properties. Approaching from the classical side, ReaxFF casts the empirical interatomic potential within a bond-order formalism, thus implicitly describing chemical bonding without expensive quantum mechanics calculations. This allows ReaxFF to accurately model both covalent and electrostatic interactions for a diverse range of materials. The ReaxFF force field is an empirical interatomic potential that is based on classical principles require reasonable computational resources, which enables simulations to better describe dynamic processes over long timeframes and on large scales. Energy contributions to the ReaxFF potential are summarized by the following: 𝐸"#"$%& = 𝐸()*+ + 𝐸)-%. + 𝐸/*01% + 𝐸$)." + 𝐸-+2 + 𝐸3)41)&( + 𝐸56%78987 Ebond is a continuous function of interatomic distance and describes the energy associated with forming bonds between atoms. Eangle and Etors are the energies associated with three-body valence angle strain and four-body torsional angle strain. Eover is an energy penalty preventing the over coordination of atoms, which is based on atomic valence rules (e.g., a stiff energy penalty is applied if a carbon atom forms more than four bonds). ECoulomb and EvdW are electrostatic and dispersive contributions calculated between all atoms, regardless of connectivity and bond-order. ESpecific represents system specific terms that are not generally included, unless required to capture properties particular to the system of interest, such as lone-pair, conjugation, hydrogen bonding, and C2 corrections. Full functional forms can be found in the Supplementary Information of the 2008-C/H/O publication.28 Different from previous work on the effect of lithiation on Li-S material bulk mechanical properties (e.g., ultimate strength, yield strength, and Young’s modulus) reported by van Duin et. al. 27, the focus of current work is the large-scale Li-S nanoparticles structural property that presents as the discharging and charging products in Li-S battery. Studies on the structures are conducted on the crystalline and amorphous Li2S nanoparticles, as well as various stoichiometric Li2S8 nanoparticles. The ring and fragment analysis are performed by analyzing various LixSy species and other intermediates. In the ring analysis, a group of end-to-end bonded atoms is counted as a ring. In the fragment analysis, a group of non-end-to-end bonded atoms is counted as a fragment. 2.1. Crystalline and Amorphous Li2S Nanoparticles First, the system studied are crystalline and amorphous Li2S nanoparticles with diameter D = 10.0 nm consisting of 32,832 atoms. To construct the crystalline Li2S nanoparticle, the initial configuration of the nanoparticle is obtained by cutting bulk crystal Li2S. To guarantee the stoichiometry of the LixSy nanoparticle (i.e., x:y=2:1), the crystalline Li2S nanoparticle is maintained by removing the outer-most redundant atoms, which are defined as any Li atom has a coordination number less than 3 or any S atom has a coordination number less than 5, after cutting the bulk crystal Li2S structure into a sphere. The amorphous (Li2S)10944 nanoparticle is constructed by placing the crystalline (Li2S)10944 nanoparticle in a ACS Paragon Plus Environment

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

vacuum environment with periodic boundary conditions in all three dimensions in the MD simulation, which initially relaxed using a conjugate gradient method, then subjected to NPT simulations for 1 ns at T = 1400 K such that all three dimensions of the external stress tensors equal zero. Then the subsequent thermal equilibration process is followed by NVT simulations for an additional 4 ns at T = 1400 K. To obtain a stable structure of amorphous (Li2S)10944 nanoparticle, the system is then quenched to T = 300 K for additional 2 ns to ensure the proper amorphization at T = 300 K with NVT conditions. Throughout the simulation, the velocity-Verlet time stepping method is adopted and the integration time step is Dt = 1 fs. In the thermal equilibration process (i.e., NPT and NVT), the fluctuation (to the mean value) in total energy and temperature is smaller than 0.005 eV/atom and 0.50 K, respectively. Throughout the system preparation, the temperature and pressure are controlled by a Nosé-Hoover thermostat and barostat 29-30 with damping constants of 100 fs and 1000 fs, respectively. 2.2. Stoichiometric Li2S8 Nanoparticles The fundamental structure and properties of Li2S8 in the condensed phase and its agglomeration, aggregation, and precipitation at confined nanoscale in Li-S battery remain elusive31-33. In this work, we assume a model nanoparticle with a Li2S8 stoichiometry obtained by maintaining the number of sulfur atoms the same as the Li2S nanoparticle discharging product, i.e., (Li2S)10944 amorphous nanoparticle that subsequently leads to stoichiometric Li2S8 (i.e., (Li2S8)1368) nanoparticle due to the delithiation process during charging process. To construct a model system for Li2S8 aggregates in the nanoparticle morphology, a hybrid Monte Carlo/molecular dynamics (MC/MD) approach is used. Based on this approach, four different strategies are adopted to construct the Li2S8 nanoparticles: (1) extracting lithium atoms randomly from the outmost layer of the (Li2S)10944 amorphous nanoparticle to form (Li2S8)1368 nanoparticle, which has the distinct core-shell structure, where the outside shell is sulfur-rich, and the inner core is still Li2Salike; (2) extracting lithium atoms from the (Li2S)10944 amorphous nanoparticle according to a spindle distribution such that the lithium concentration is higher in the core region than the outer shell region to form an indistinguishable core-shell boundary throughout the structure of nanoparticle (where the spindle distribution is constructed by randomly extracting lithium atoms within the unit volume in the spherical coordinate of the (Li2S)10944 amorphous nanoparticle using Monte Carlo method); (3) maintaining sulfur atoms uniformly throughout the (Li2S)10944 amorphous nanoparticle to form a nearly homogeneous (Li2S8)1368 nanoparticle, (4) randomly packing the Li2S8 molecules obtained from density functional theory 34-36 to form an assembled (Li2S8)1368 nanoparticle. From all these four approaches, the stoichiometry of these LixSy nanoparticles is constrained to be x:y=2:8. We refer to these four (Li2S8)1368 nanoparticles with the following names: core-shell, spindle distribution, homogenous and assembled Li2S8 structures. The initial configurations of these nanoparticles are placed in vacuum environments with periodic boundary conditions in all three dimensions, subsequently relaxed using the conjugate gradient method and then subjected to NPT simulations for 1.32 ns at T = 300 K, such that all three dimensions of the external stress tensor are set to zero. The thermal equilibration is then followed by NVT simulations for an additional 2 ns at T = 300 K. For these simulation, the velocity-Verlet time stepping method is adopted and the integration time step is Dt = 1 fs. During the thermal equilibration process (i.e., NPT and NVT), the fluctuation (to the mean value) in total energy and temperature is smaller than 0.005 eV/atom and 0.50 K, respectively. The temperature and pressure are controlled by a Nosé-Hoover thermostat and barostat29-30 with damping constants of 100 fs and 1000 fs, respectively. 3. Results and Discussion The focus of this study is on the structural property of Li-S nanoparticle apart from the idealized stoichiometric lithium sulfide (Li2S) crystal with accurate atomistic descriptions. Reactive force field (ReaxFF) is applied in our large-scale molecular dynamics (MD) simulations to probe the structural properties of Li-S nanoparticles during the battery charging process. The path of improving the performance ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

of Li-S batteries begins with an in-depth structural property understanding of the Li-S system, which can be investigated either by advanced experimental characterization techniques or the state-of-the-art simulation methods. In an operating Li-S battery, systematic studies on the evolution of all possible LixSy speciation during discharging and charging process are not trivial. For Li2S solid embedded in conducting carbon electrodes (e.g., Fig. 1), the charging process happens when Li2S is delithiated, which leads to the presence of lithium-poor or sulfur-rich lithium polysulfide species (e.g., Li2Sn, n > 2). For these lithium polysulfide Li2Sn molecules and their intermediaries, the energetics and thermodynamic stability can be determined by quantum chemistry calculations at the molecular level.37 However, the condensed phase of those species remain elusive, especially the detailed structural property for nano-aggregates of these lithium polysulfide LixSy species. To study the structural property of LixSy nano-aggregates at the nanoscale, we focus our studies on these Li-S model systems: crystalline and amorphous Li2S nanoparticles, and various stoichiometric Li2S8 nanoparticles as described in Section 2.1 and 2.2. Figure 2 depicts the crystalline and amorphous configurations of the Li2S nanoparticles, along with the computed radial distribution function g(r) for different atomic pairs. From Fig. 2(c) and 2(d), the similarity between the crystalline and amorphous configuration of the Li2S nanoparticle can only be found within the first coordination shell for Li-Li, Li-S, and S-S bonds. The long-range crystallinity disappears beyond ~ 6.0 Å for the amorphous Li2S nanoparticle. Figure 3 shows the thermally equilibrated configurations of the various stoichiometric Li2S8 model nanoparticles (core-shell, spindle distribution, homogeneous and assembled), along with the computed radial density distribution of each element (i.e., Li and S) and radial pair distribution functions g(r) of the nanoparticles in the Li2S8 model systems. 3.1. Nanoparticle Energetic Analysis To understand how the structures of these nanoparticles correlated with thermodynamic stability, we compare the potential energy of amorphous and crystalline Li2S nanoparticles and the various Li2S8 nanoparticles along with the bulk 𝛼-sulfur. Heat of formation of the lithiated sulfur configuration as a function of lithium content with respect to 𝛼-sulfur and bcc-Li bulk were calculated using the relation ∆𝐸 = 𝐸 20 atoms) species are relatively small. 4. Conclusion We carried out MD simulations of the thermodynamic stability, structural analysis of mixed lithium sulfide and polysulfide speciation and its complex distribution in Li-S (stoichiometric Li2S and Li2S8) nanoparticles. These complex nanostructures are related to discharging and charging products in Li-S batteries. Our simulations show that the local atomic structure in a Li-S nanoparticle may prescribe its stability despite the similarity in stoichiometric compositions. In contrast to end discharge products, i.e., Li2S nanoparticles, our thermodynamic stability analysis suggests that the core-shell Li2S8 nanoparticles are most likely present during the charging process due to delithiation of Li2S nanoparticles. Based on our structural (ring and fragment) analysis, a broad distribution of sulfur, lithium sulfide and lithium polysulfide (i.e., Sn, LixS, and LixSy(y>1)) species are found in both stoichiometric Li2S and Li2S8 nanoparticles. A wide distribution in the population of sulfur Sn rings, and lithium polysulfide LixSy rings for the Li2S8 nanoparticles suggests that there is a complex energy landscape of competing intermediates within a narrow voltage window. Overall, in contrast to bulk crystal Li2S, the presence of mixed lithium sulfide and lithium polysulfide species and its complex distribution are found to be common features for these Li-S (Li2S and Li2S8) nanoparticles. To validate this prediction, characterization by operando x-ray absorption spectroscopy might be needed. The complex distribution of these sulfide and polysulfide speciation (i.e., sulfur (Sn), lithium sulfides (LixS) and lithium polysulfides (LixSy(y>1)) rings and fragments) are dictated by both stoichiometry and local atomic structures in the nanoparticle. The complicated Li-S speciation distribution and local atomic structures within these nanoparticles cannot be accessed via standard XRD spectra analysis. We believe the current systematic study will provide useful structuralproperty insights to the further development of functionalized cathode in Li-S battery. Acknowledgment This research used resources of the Argonne Leadership Computing Facility, which is a U.S. Department of Energy Office of Science User Facility operated under Contract DE-AC02-06CH11357. Y. Li acknowledges support from the Margaret Butler Postdoctoral Fellowship at Argonne National

ACS Paragon Plus Environment

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Laboratory. K. C. Lau acknowledges the support of California State University Northridge Faculty start-up fund.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

Reference (1) Goodenough, J. B. Electrochemical Energy Storage in a Sustainable Modern Society. Energ Environ Sci 2014, 7 (1), 14-18, DOI: 10.1039/c3ee42613k. (2) Luo, X.; Wang, J. H.; Dooner, M.; Clarke, J. Overview of Current Development in Electrical Energy Storage Technologies And The Application Potential in Power System Operation. Appl Energ 2015, 137, 511-536, DOI: 10.1016/j.apenergy.2014.09.081. (3) Manthiram, A.; Fu, Y. Z.; Chung, S. H.; Zu, C. X.; Su, Y. S. Rechargeable Lithium-Sulfur Batteries. Chem Rev 2014, 114 (23), 11751-11787, DOI: 10.1021/cr500062v. (4) Chu, S.; Cui, Y.; Liu, N. The Path Towards Sustainable Energy. Nat Mater 2017, 16 (1), 16-22, DOI: 10.1038/nmat4834. (5) Manthiram, A.; Fu, Y. Z.; Su, Y. S. Challenges and Prospects of Lithium-Sulfur Batteries. Accounts Chem Res 2013, 46 (5), 1125-1134, DOI: 10.1021/ar300179v. (6) Ji, X. L.; Nazar, L. F. Advances in Li-S Batteries. J Mater Chem 2010, 20 (44), 9821-9826, DOI: 10.1039/b925751a. (7) Bruce, P. G.; Freunberger, S. A.; Hardwick, L. J.; Tarascon, J. M. Li-O-2 and Li-S Batteries With High Energy Storage. Nat Mater 2012, 11 (1), 19-29, DOI: 10.1038/nmat3191. (8) Mikhaylik, Y. V.; Akridge, J. R. Polysulfide Shuttle Study In the Li/S Battery System. J Electrochem Soc 2004, 151 (11), A1969-A1976, DOI: 10.1149/1.1806394. (9) Li, C.; Wang, Z. B.; Wang, Q.; Gu, D. M. Recent Advances in Cathode Materials for Li-S Battery: Structure and Performance. Rare Metals 2017, 36 (5), 365-380, DOI: 10.1007/s12598-017-0900-2. (10) Son, Y.; Lee, J. S.; Son, Y.; Jang, J. H.; Cho, J. Recent Advances in Lithium Sulfide Cathode Materials and Their Use in Lithium Sulfur Batteries. Adv Energy Mater 2015, 5 (16), DOI: 10.1002/aenm.201500110. (11) Cai, K. P.; Song, M. K.; Cairns, E. J.; Zhang, Y. G. Nanostructured Li2S-C Composites as Cathode Material for High-Energy Lithium/Sulfur Batteries. Nano Lett 2012, 12 (12), 6474-6479, DOI: 10.1021/nl303965a. (12) Seh, Z. W.; Yu, J. H.; Li, W. Y.; Hsu, P. C.; Wang, H. T.; Sun, Y. M.; Yao, H. B.; Zhang, Q. F.; Cui, Y. Two-Dimensional Layered Transition Metal Disulphides for Effective Encapsulation of High-Capacity Lithium Sulphide Cathodes. Nat Commun 2014, 5, DOI: 10.1038/Ncomms6017. (13) Zu, C. X.; Klein, M.; Manthiram, A. Activated Li2S as a High-Performance Cathode for Rechargeable Lithium-Sulfur Batteries. J Phys Chem Lett 2014, 5 (22), 3986-3991, DOI: 10.1021/jz5021108. (14) Yin, Y. X.; Xin, S.; Guo, Y. G.; Wan, L. J. Lithium-Sulfur Batteries: Electrochemistry, Materials, and Prospects. Angew Chem Int Edit 2013, 52 (50), 13186-13200, DOI: 10.1002/anie.201304762. (15) Song, M. K.; Zhang, Y. G.; Cairns, E. J. A Long-Life, High-Rate Lithium/Sulfur Cell: A Multifaceted Approach to Enhancing Cell Performance. Nano Lett 2013, 13 (12), 5891-5899, DOI: 10.1021/nl402793z. (16) Pope, M. A.; Aksay, I. A. Structural Design of Cathodes for Li-S Batteries. Adv Energy Mater 2015, 5 (16), DOI: 10.1002/aenm.201500124. (17) Klein, M. J.; Veith, G. M.; Manthiram, A. Chemistry of Sputter-Deposited Lithium Sulfide Films. J Am Chem Soc 2017, 139 (31), 10669-10676, DOI: 10.1021/jacs.7b03379. (18) Yamin, H.; Gorenshtein, A.; Penciner, J.; Sternberg, Y.; Peled, E. Lithium Sulfur Battery - Oxidation Reduction-Mechanisms of Polysulfides in Thf Solutions. J Electrochem Soc 1988, 135 (5), 1045-1048, DOI: 10.1149/1.2095868. (19) Barchasz, C.; Molton, F.; Duboc, C.; Lepretre, J. C.; Patoux, S.; Alloin, F. Lithium/Sulfur Cell Discharge Mechanism: An Original Approach for Intermediate Species Identification. Anal Chem 2012, 84 (9), 3973-3980, DOI: 10.1021/ac2032244. (20) Patel, M. U. M.; Demir-Cakan, R.; Morcrette, M.; Tarascon, J. M.; Gaberscek, M.; Dominko, R. Li-S Battery Analyzed by UV/Vis in Operando Mode. Chemsuschem 2013, 6 (7), 1177-1181, DOI: 10.1002/cssc.201300142. ACS Paragon Plus Environment

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(21) Cuisinier, M.; Cabelguen, P. E.; Adams, B. D.; Garsuch, A.; Balasubramanian, M.; Nazar, L. F. Unique Behaviour of Nonsolvents For Polysulphides in Lithium-Sulphur Batteries. Energ Environ Sci 2014, 7 (8), 2697-2705, DOI: 10.1039/c4ee00372a. (22) Wujcik, K. H.; Pascal, T. A.; Pemmaraju, C. D.; Devaux, D.; Stolte, W. C.; Balsara, N. P.; Prendergast, D. Characterization of Polysulfide Radicals Present in an Ether-Based Electrolyte of a Lithium-Sulfur Battery During Initial Discharge Using In Situ X-Ray Absorption Spectroscopy Experiments and First-Principles Calculations. Adv Energy Mater 2015, 5 (16), DOI: 10.1002/aenm.201500285. (23) Cuisinier, M.; Cabelguen, P. E.; Evers, S.; He, G.; Kolbeck, M.; Garsuch, A.; Bolin, T.; Balasubramanian, M.; Nazar, L. F. Sulfur Speciation in Li-S Batteries Determined by Operando X-ray Absorption Spectroscopy. J Phys Chem Lett 2013, 4 (19), 3227-3232, DOI: 10.1021/jz401763d. (24) Pascal, T. A.; Villaluenga, I.; Wujcik, K. H.; Devaux, D.; Jiang, X.; Wang, D. R.; Balsara, N.; Prendergast, D. Liquid Sulfur Impregnation of Microporous Carbon Accelerated by Nanoscale Interfacial Effects. Nano Lett 2017, 17 (4), 2517-2523, DOI: 10.1021/acs.nanolett.7b00249. (25) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular-Dynamics. J Comput Phys 1995, 117 (1), 1-19, DOI: 10.1006/jcph.1995.1039. (26) Senftle, T. P.; Hong, S.; Islam, M. M.; Kylasa, S. B.; Zheng, Y. X.; Shin, Y. K.; Junkermeier, C.; Engel-Herbert, R.; Janik, M. J.; Aktulga, H. M.; Verstraelen, T.; Grama, A.; van Duin, A. C. T. The ReaxFF Reactive Force-Feld: Development, Applications and Future Directions. Npj Comput Mater 2016, 2, DOI: 10.1038/npjcompumats.2015.11. (27) Islam, M. M.; Ostadhossein, A.; Borodin, O.; Yeates, A. T.; Tipton, W. W.; Hennig, R. G.; Kumar, N.; van Duin, A. C. T. ReaxFF Molecular Dynamics Simulations on Lithiated Sulfur Cathode Materials. Phys Chem Chem Phys 2015, 17 (5), 3383-3393, DOI: 10.1039/c4cp04532g. (28) Chenoweth, K.; van Duin, A. C. T.; Goddard, W. A. ReaxFF Reactive Force Field for Molecular Dynamics Simulations of Hydrocarbon Oxidation. J Phys Chem A 2008, 112 (5), 1040-1053, DOI: 10.1021/jp709896w. (29) Nose, S. A Unified Formulation of the Constant Temperature Molecular-Dynamics Methods. J Chem Phys 1984, 81 (1), 511-519, DOI: 10.1063/1.447334. (30) Hoover, W. G. Canonical Dynamics - Equilibrium Phase-Space Distributions. Phys Rev A 1985, 31 (3), 1695-1697, DOI: 10.1103/PhysRevA.31.1695. (31) Visco, S. J.; Goncharenko, N.; Nimon, V.; Petrov, A.; Nimon, Y. S.; De Jonghe, L. C.; Katz, B. D.; Loginova, V. Lithium Sulfur Batteries and Electrolytes and Sulfur Cathodes Thereof. 9,660,265; Other: 14/657,695, 2017. (32) Tien Q. Duong, A. B. M. R. P. M.; Vehicle Technologies Office, E. E. a. R. E. Progress Report Advanced Battery Materials Research (BMR) Program. 2017.09. (33) Wang, M. C.; Yu, J. G.; Lin, S. C. Lithiation-Assisted Strengthening Effect and Reactive Flow in Bulk and Nanoconfined Sulfur Cathodes of Lithium-Sulfur Batteries. J Phys Chem C 2017, 121 (31), 17029-17037, DOI: 10.1021/acs.jpcc.7b05446. (34) Capelle, K. A Bird's-eye View of Density-Functional Theory. Braz J Phys 2006, 36 (4a), 1318-1343, DOI: 10.1590/S0103-97332006000700035. (35) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys Rev B 1964, 136 (3b), B864-+, DOI: 10.1103/PhysRev.136.B864. (36) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys Rev 1965, 140 (4a), 1133-&, DOI: 10.1103/PhysRev.140.A1133. (37) Assary, R. S.; Curtiss, L. A.; Moore, J. S. Toward a Molecular Understanding of Energetics in Li-S Batteries Using Nonaqueous Electrolytes: A High-Level Quantum Chemical Study. J Phys Chem C 2014, 118 (22), 11545-11558, DOI: 10.1021/jp5015466. (38) Li, X. M.; Wolden, C. A.; Ban, C. M.; Yang, Y. G. Facile Synthesis of Lithium Sulfide Nanocrystals for Use in Advanced Rechargeable Batteries. Acs Appl Mater Inter 2015, 7 (51), 28444-28451, DOI: 10.1021/acsami.5b09367. ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(39) Lee, S. K.; Lee, Y. J.; Sun, Y. K. Nanostructured Lithium Sulfide Materials for Lithium-Sulfur Batteries. J Power Sources 2016, 323, 174-188, DOI: 10.1016/j.jpowsour.2016.05.037. (40) van Duin, A. C. T.; Zeiri, Y.; Dubnikova, F.; Kosloff, R.; Goddard, W. A. Atomistic-Scale Simulations of The Initial Chemical Events In The Thermal Initiation of Triacetonetriperoxide. J Am Chem Soc 2005, 127 (31), 11053-11062, DOI: 10.1021/ja052067y.

ACS Paragon Plus Environment

Page 12 of 23

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Graphical Table of Contents

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. A schematic description of lithium sulfide/polysulfide nanoparticles that present in Li-S batteries cathode with carbon interfaces inspired by Nazar et. al. Reprinted in part with permission from reference 23 and reference 24.


ACS Paragon Plus Environment

Page 14 of 23

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces



Figure 2. Schematic visualizations of (a) crystalline and (b) amorphous configuration of a (Li2S)10944 nanoparticle with 10 nm in diameter. Sulfur and lithium atoms are colored yellow and pink, respectively. Radial pair distribution function g(r) for (c) crystalline and (d) amorphous configurations. gLi-Li(r), gLi-S(r) and gS-S(r) is pink, orange and yellow, respectively.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23



Figure 3. Cross-section of (Li2S8)1368 nanoparticles in (a) core-shell (b) spindle distribution (c) homogeneous and (d) assembled structures. Sulfur and lithium atoms are colored yellow and pink, respectively. Radial distribution of density d(r) for (e) core-shell (f) spindle distribution (g) homogeneous and (h) assembled structures of (Li2S8)1368 nanoparticles. d Li(r), dS(r) and dtot(r) is pink, blue and cyan, respectively. Radial pair distribution function g(r) for (i) core-shell (j) spindle distribution (k) homogeneous and (l) assembled structures of (Li2S8)1638 nanoparticles. g Li-Li(r), gLi-S(r) and gS-S(r) is pink, orange and yellow, respectively.

ACS Paragon Plus Environment

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces



Figure 4. Energy profile of various model Li-S nanoparticles in the unit of electronic voltage per atom as a function of simulation time.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. The simulated XRD spectra of various Li2S and Li2 S8 nanoparticles (plotted in lines) in comparison with the known bulk crystal Li2S (plotted in dots).

ACS Paragon Plus Environment

Page 18 of 23



Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 6. The fraction of atoms in the ring (red bars) and in the fragment (blue bars) for various Li-S nanoparticles.

ACS Paragon Plus Environment



ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23



Figure 7. The number of various lithium polysulfide rings (red bars) and positions of rings (blue bars) for the (a) crystalline Li2S and (b) amorphous Li2S nanoparticles. The position of rings is defined as the averaged distance of the center of the ring to the center of the nanoparticle (with standard deviation represented as the error bar). Various lithium polysulfide rings in (c) crystalline and (d) amorphous Li2S nanoparticle on real space, with ring sizes color coded as the color bar shown. The circles outline the shape of the nanoparticles.

ACS Paragon Plus Environment

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces



Figure 8. The number of lithium polysulfide and sulfur rings (red bars) and position of rings (blue bars) for the (a) core-shell (b) spindle distribution (c) homogeneous and (d) assembled Li2S8 nanoparticles. Various lithium polysulfide and sulfur rings for the (e) core-shell (f) spindle distribution (g) homogeneous and (h) assembled Li2S8 nanoparticles in real space, with stoichiometry ratio of rings color coded as shown in the color bar. The circles outline the shape of the crosscut nanoparticles.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 23



Figure 9. The number of sulfides (blue bars) and polysulfide fragments (green bars) for the (a) crystal and (b) amorphous Li2S nanoparticles. Various sulfide and polysulfide fragments for a (c) crystal and (d) amorphous Li2S nanoparticles in real space, with the fragment sizes color coded as shown in the color bar. The circles outline the shape of the crosscut nanoparticles.

ACS Paragon Plus Environment

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces



Figure 10. The number of various sulfur fragments (red bars), sulfide fragments (blue bars) and lithium polysulfide fragments (green bars) in (a) core-shell (b) spindle distribution (c) homogeneous and (d) assembled Li2S8 nanoparticles. Various fragments distribution in the crosscut (e) core-shell (f) spindle distribution (g) homogeneous and (h) assembled Li2S8 nanoparticle in real space, with the fragment size color coded as the color bar shown. The circles outline the shape of the cross-cut nanoparticles.

ACS Paragon Plus Environment