Tunable Mechanical and Thermal Properties of One-Dimensional

Sep 30, 2015 - Recently, carbyne chain, the one-dimensional sp-hybridized carbon allotrope in the form of either α-carbyne (polyyne), with alternatin...
2 downloads 11 Views 1MB Size
Subscriber access provided by NEW YORK MED COLL

Article

Tunable Mechanical and Thermal Properties of One-dimensional Carbyne Chain: Phase Transition and Microscopic Dynamics Xiangjun Liu, Gang Zhang, and Yong-Wei Zhang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b08026 • Publication Date (Web): 30 Sep 2015 Downloaded from http://pubs.acs.org on October 3, 2015

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 32

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

The Journal of Physical Chemistry

Tunable Mechanical and Thermal Properties of One-dimensional Carbyne Chain: Phase Transition and Microscopic Dynamics Xiangjun Liu†, Gang Zhang†*, Yong-Wei Zhang† †

Institute of High performance Computing, A*STAR, Singapore 138632 Abstract

Recently, carbyne chain, the one-dimensional sp-hybridized carbon allotrope in the form of either α-carbyne (polyyne) with alternating single and triple bonds or β-carbyne (cumulene) with repeating double bonds, has attracted more and more attention. However, the mechanical and thermal properties of individual phase, their phase transition dynamics and defect formation remain largely unknown. Our molecular dynamics simulations show that the critical temperature for the phase transition from cumulene to polyyne is 499 K, and the phase transition is ultrafast and completed within 150 fs. To achieve perfect polyyne, however, refined temperature control is needed so as to avoid defective bonds. The bending stiffness and Young’s modulus of cumulene are significantly higher than those of polyyne, while both of them are comparable to the hardest natural materials. The large difference in the stress-strain behavior between cumulene and polyyne provides a novel route for storing mechanical energy. Furthermore, the thermal conductivity of cumulene is found to be two times higher than that of polyyne, and the defective bonds can dramatically decrease the thermal conductivity of polyyne to only 13% of that of pristine polyyne. The significant changes in the mechanical and thermal properties between the two phases of carbyne chain present a great opportunity for its use as strain sensors, mechanical connectors and mechanical/thermal energy storage devices.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Introduction Carbon allotropes have attracted continuous interest due to their stability and various fascinating properties.1 Carbon atom is capable of forming various allotropes: three-dimensional sp3 diamond, two-dimensional sp2 graphene,2 quasi-one-dimensional carbon nanotube,3 or zerodimensional sp2 fullerenes.4 Each of them has notably different electronic and mechanical properties. Diamond is a wide band gap insulator and one of the hardest natural materials. On the contrary, graphene has a semimetal electronic structure with linear band dispersion and extraordinarily high electron mobility.5,6 One of the newest potential carbon allotropes is the monoatomistic chain of sp-hybridized carbon, and its bulk modification is known as carbyne. Although the existence of this exotic phase of bulk carbon is still debated until today,7 the isolated chains have been observed in several experiments, for example, in gas-phase deposition, epitaxial growth, electrochemical synthesis, or pulling the atomic chains from graphene or carbon nanotubes.8−19 The longest isolable carbon chains was produced in solution with a length of about 60 Å.17 Theoretical21−23 and experimental studies have revealed interesting properties for carbon chains. A sp carbon chain can have two possible types of structures: cumulene with repeating double bonds over the whole chain, ( = C = C =) n ; or polyyne with alternating single and triple bonds, ( −C ≡ C −) n . In the symmetric cumulene structure, it is characterized by σ-bonds along the axis (s-px orbitals), and two decoupled π-bonds per atom from the perpendicular py and pz orbitals. The π-electrons are uniformly distributed over the chain, which makes cumulene a metallic conductor. On the other hand, in the symmetry-breaking polyyne, the electrons are mostly localized at the triple bonds, thus polyyne is expected to be semiconducting. As a 1D conductor, cumulene can be destabilized by Peierls distortion, so that polyyne with alternating 2 ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

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

The Journal of Physical Chemistry

bond lengths becomes energetically favoured24−25 with a band gap at the edge of Brillouin zone to lower the total energy per atom by 2 meV. As the band gap increases with increasing strain,25−27 the Peierls distortion also increases in strained carbon chains compared with the freestanding ones, which has been observed experimentally in recent studies.28 These novel electronic properties and strain-induced metal-semiconductor transition are promising in the applications of electromechanical switching, strain sensor, and atomistic scale circuit26,29-31. To drive the progress of these applications, an in-depth theoretical understanding of the structural transition between cumulene and polyyne is necessary. Although polyyne is believed to be more energetically favorable, a recent theoretical study by Artyukhov et al. suggested that in unstrained carbyne chains, the energetic gain by the Peierls instability is so small that it can be overcompensated by zero-point vibrations.27 This finding is critically important since it predicts that metallic cumulene can exist under specified conditions. Thus the cumulene-polyyne transition temperature in mechanically relaxed chains is an important question that deserves further studies. Although many theoretical studies of carbon chains have been performed in the past decade, most studies focused on the strain-induced transition, while the temperature effect and dynamical information about the bond-rearranging process during phase transition is still very limited. Furthermore, compared with the large amount of studies of electronic properties, the thermal properties of 1D carbon chain remain elusive. Molecular dynamics (MD) simulation is a powerful tool to handle many-body problems at atomic level. Besides predicting the intrinsic thermal properties without the thermodynamic-limit assumption,32 the MD simulation is also very efficient for studying mechanical properties33 and structural changes34 of nanomaterials. In MD simulation, the Newton's equation of motion for each atom is numerically solved with the interatomic force derived from the system’s potential 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 32

energy. In the present work, we perform molecular dynamics simulations using the firstprinciples based ReaxFF35,36 force field to address the questions raised above, aiming to not only reveal new insights into the fundamental properties of carbyne, but also provide useful design guidelines towards the functional use of this novel material. Computational Methods In this work, MD simulations are employed to study the mechanical and thermal properties of carbyne using the LAMMPS package.37 In all MD simulations performed here, the ReaxFF force field,35,36 an empirical bond-order-dependent potential, is used to describe the reactive, (single-, double-, and triple-) covalent bonding interactions, in which dissociation and reaction curves were obtained from fitting to quantum mechanics calculations. ReaxFF38 is based on a bond order/bond distance relationship, a concept introduced by Tersoff39 and first employed to carbon chemistry by Brenner40. Bond orders ( BOij' ), including contributions from σ, π, and π-π bonds, are calculated from the interatomic distances, using equation, '

σ

π

ππ

ij

ij

ij

ij

BO = BO + BO + BO

p p p     r ij  bo 2   r ij  bo 4   r ij  bo 6  .          = exp p + exp p + exp p  bo1  σ    bo 3  π    bo 5  ππ    r 0    r 0    r 0     

(1)

The instantaneous bond orders are updated every MD step, thus the bond-order-dependent ReaxFF provides accurate descriptions of bond breaking and bond formation during simulation. This allows us for the possibility to catch the instantaneous bond state in the simulation. In addition, ReaxFF also accounts for polarization effects by using a geometry-dependent charge calculation scheme. Thus, the bond-order based ReaxFF method is built to bridge the gap between Quantum mechanics method and conventional molecular dynamics method. ReaxFF partitions the overall system energy into contributions from various partial energy terms, include the terms of bond, lone pair, over-coordination, under-coordination, angle strain, 4 ACS Paragon Plus Environment

Page 5 of 32

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

The Journal of Physical Chemistry

angle conjugation, C2 correction, torsion, torsion conjugation, hydrogen bond to properly handle the nature of preferred configurations of atomic and molecular orbitals, and other terms of van der Waals and Coulomb interactions.38-40 The latter two non-bonded interactions are calculated between each atom pair, irrespective of connectivity, and are shielded to avoid excessive repulsion at short distances. This treatment of non-bonded interactions allows ReaxFF to describe ionic, covalent, and intermediate materials, thus, greatly improving its transferability. In this paper, the version of ReaxFF force field (ReaxFFCHO) and its parameters are adopted from Ref. 36, which is developed for carbon-rich hydrocarbon materials. The ReaxFFCHO expanded the ReaxFF reactive force field training set to include additional transition states and chemical reactivity of systems relevant to these reactions and optimized the force field parameters against a quantum mechanics-based training set. The ReaxFFCHO can accurately provide the chemical reaction process of dehydrogenation, e.g. from H3 C − C H2 − C H2 − C H3 to HC ≡ C − C ≡ CH .36 It has also been demonstrated that the ReaxFF potential provides a reasonable

account of the chemical and mechanical behavior of many carbon-based materials,41−44 including carbyne.44 The velocity Verlet algorithm is employed to integrate Newton’s equations of atom motion numerically with a time step of 0.1 fs. First, the system was equilibrated at a constant target temperature T for 100 ps using Nosé-Hoover temperature thermostat45 (NVT ensemble). After the constant temperature relaxation, we continued to relax the system with NVE (constant volume and no thermostat) ensemble for 100 ps. During this stage, the total energy and temperature of the system were monitored. We found that the total energy was conserved and the temperature of the entire system remained constant with fluctuations around T, which means that the system has reached equilibrium. It is worth pointing out that, in order to take the quantum effects of phonon occupation into account, the temperature calculated from the MD simulation 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 6 of 32

2 2 3 (TMD) has been corrected using the quantum correction scheme46,47 T MD = T2 ∫0T T dx xx , where e −1 TD D

TMD, TD (= 517 K for carbyne48) and T are the MD temperature, Debye temperature and corrected temperature, respectively. Results and Discussions Temperature Induced Phase Transition First, we explored the temperature induced phase transition of carbyne from cumulene to polyyne. In experiments, neutron-diffraction technique gives the most accessible information on the microscopic structures of a material. Usually, the bond-length distribution function calculated from MD simulations can provide comparable results to the neutron-diffraction data on atomistic structure.49 In our MD simulation, a cumulene chain is initially constructed and relaxed at 5 K. The temperature is then increased with a rate of approximately 0.04 K/fs. Figure 1a shows the bond-length distribution function of carbyne chain at 498 K (Similar distributions at temperature lower than 498K are not shown here). A single peak is clearly shown, with an equilibrium bond length of 1.37 Å. However, at a higher temperature such as 500 K, the bond-length distribution is quite different, as shown in Figure 1b (Similar distributions at temperature higher than 500K are not shown here). At this temperature, the distribution peak at 1.37 Å disappears, while two separated peaks appear at 1.21 and 1.58 Å, respectively. In DFT studies, the bond length alternation (BLA, defined as the difference between the long and short bonds) is usually used to describe the phase transition. An increase of BLA from cumulene (0.088Å) to polyyne (0.248Å) is predicted by HSE06 calculations.26 Our MD simulation results of single/double/triple bond lengths are in a reasonable agreement with those (1.54Å, 1.35Å, and 1.20Å for single, double, and triple bond, respectively) from first-principles calculations.50 6 ACS Paragon Plus Environment

Page 7 of 32

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

The Journal of Physical Chemistry

Clearly, the critical temperature of phase transition is between 498K and 500K. With further refined calculations, the critical temperature is found to be 499K, below which only carboncarbon double bonds are observed, while above which all double bonds change to alternating single and triple bonds. Using DFT calculations together with analytic modeling, Artyukhov et al.27 studied the strain-temperature “phase diagram” of the cumulene-polyyne transition and found in unstrained carbyne, the transition temperature is about 800K. The differences between our MD predicted transition temperature and the value from DFT could be due to the following reasons. First, in principle, each phonon mode is equally excited in classical MD simulations, which is different from the quantum statistical distribution. Furthermore, the periodic boundary condition is used along the axis of carbyne chain in our MD simulations. Note that adoption of different boundary conditions has been found to have certain influence on the atomic structure of carbyne chains in DFT calculations.24,51,52 Second, the ReaxFF potential could not exactly reproduce the force field predicted by first-principles calculations, while the strength of Peierls distortion is extremely sensitive to the force field.26 Third, in DFT calculations, the well-known underestimated band gap53,54 may result in a much weakened Peierls instability, which may lead to an overestimated transition temperature. Overall, the most important characteristic in the transition process is well revealed by our MD calculations. A significant advantage of MD study is its ability to provide a detailed description of dynamics process, including bond breaking and creating. To explore the dynamics picture of the phase transition process, we record the time-dependent bond-length distribution function at the critical temperature of 499 K. As shown in Figure 2a, initially, only double-bonds exist in the system. Local perturbations created by the vibrations of carbon atoms in the cumulene chain 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

induce the nucleation of a pair of single- and triple-bonds at one carbon atom. This leads to its adjacent double bonds to continue to convert to triple-/single-bonds, causing a spontaneous phase transition from cumulene to polyyne. This phase transition process is reflected in the decreased distribution peak of double-bonds and the increased distribution peaks of single- and triple-bonds. It is interesting to note that the phase transition process is ultrafast. For a carbyne chain with 50 atoms, it only needs 145 fs to complete the phase transition process, as shown in Figure 2a−2e. The ideal polyyne atomic structure with alternating single- and triple-bonds is predicted when the temperature is increased from 5K to 499K with a sufficiently slow heating rate. However, if the heating rate is high (e.g. 2.5K/fs) or if the cumulene chain is initialized under a temperature over 499K, several non-ideal bonds are observed as shown in Figure 2j. The phase transition process is shown in Figure 2f−2j. Unlike the representative bond-length distribution function in perfect polyyne, there exists one additional peak at 1.45 Å, which is shorter than that of the single-bond, but longer than the double-bond and triple-bond. In this work, for simplicity, this bond is termed as “defective bond”. Similar “soliton” defects have also been reported in polyethylene chains.55 This phenomenon indicates that, during the fabrication/formation process, temperature control is a critical factor to produce perfect polyyne. In previous DFT studies 26, 27, the bond length alternation is used to evaluate the strain induced phase change, which is defined as the difference between the long and short bonds. However, as the length of the defective bond is within the lengths of single- and triple-bond, this defective bond was ignored by previous DFT studies. Mechanical Properties of Carbyne

8 ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

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

The Journal of Physical Chemistry

Next, we investigated the bending rigidity of both cumulene and polyyne chains. As the focus here is the bending rigidity of carbyne chains, a curvature is imposed on the carbyne chain by forming a ring structure. The potential energy of the ring (Er) is obtained by energy minimization with an energy tolerance of 4.3×10-8 eV. To directly assess the bending energy of the ring, we calculated the potential energy of a linear chain with the same atom number (Ec), as shown in Figure 3a. The difference between these two energies is the bending energy (Eb = Er – Ec). A series of ring structures are constructed with carbon atoms N from 8 to 50, with associated ideal radii (R) of approximately 1.70 to 10.66 Å for cumulene rings and 1.82 to 11.07 Å for polyyne rings, respectively. The results of bending energy of cumulene and polyyne are shown in Figure 3b. These carbyne rings represented a set of structures with constant curvatures, ρ (= R-1). Combing the bending energy with the minimized perimeter L of each ring resulted in the normalized energy per length, EL (= Eb/L eV/Å), as shown in Figure 3c. It is found that the bending energy of cumulene chain is higher than that of polyyne chain with the same atom number N, but as ρ decreases (N increases), the difference gradually decreases and it disappears when ρ is below 0.09 Å-1. In continuum elastic theory, the elastic bending energy Eb of a curved element can be calculated as26,44

Eb =

1 L 2 ∫0 Dρ dsˆ, 2

(2)

where D is the bending stiffness and ρ is the curvature. Assuming a constant bending stiffness D, the normalized bending energy per length is:

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

EL =

E b = 1 Dρ 2 . L

Page 10 of 32

(3)

2

Figure 3d shows the normalized bending energy per length EL versus ρ2. When ρ2 is less than 0.1 Å-2, the calculated normalized bending energy depends linearly on ρ2, which is in good agreement with the relationship described in Eq. 3. In this curvature range, the fitted bending stiffness D are 8.5 eV·Å for cumulene chain, and 6.7 eV·Å for polyyne chain. Hence, cumulene has a 27% higher bending stiffness than polyyne, which results from the fact that the sp2 conjugated double-bond in cumulene forms a more rigid linear structure with a less amplitude of bending/torsional motion than the alternating single and triple bonds in polyyne. First-principles method has been used to calculate the bending stiffness of polyyne chain with DDFT =

1 ∂ 2 EL , where a = 2.565Å is the unit cell length, with DDFT=3.56 eV·Å.26 Using a ∂ρ 2

the same definition, the calculated bending stiffness in our work is 2.6 eV·Å , which is close to that presented in Ref. 26. Using MD simulation with ReaxFF potential, Nair et al. studied the length-dependent vibrational frequencies of carbyne chains.41 In their simulation, a carbyne chain with two ends fixed was excited by perturbing the central atoms and allowed to freely vibrate. The fundamental frequencies of vibration were obtained by analyzing the oscillations via a discrete Fourier transform. The highest frequency was found to be close to 6THz in the shortest carbyne chain of 5 Å. Then the chain was assumed to be an elastic beam, and the frequency f was related to the stiffness D by the relation 2πf = α 2 D m L3 , where m is the mass, α is a constant of 4.43 for the chain with a load at the center. As the vibration frequency f is dependent on the length of carbyne chain, the estimated bending stiffness is also length-dependent. For a chain with length of 8 Å, 10 ACS Paragon Plus Environment

Page 11 of 32

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

The Journal of Physical Chemistry

the bending stiffness was found to be about 8.5 eV·Å , which is close to our predicted value of 6.7 eV·Å for polyyne chain. To use carbyne as a structural connection in mechanical components, persistence length is an important parameter that dictates its usability. The bending stiffness can be evaluated using the concept of persistence length Lp in polymer physics: Lp = D/kBT, where kB is the Boltzmann constant and T is the temperature.56 The persistence length Lp for carbyne chain is found to be about 33 nm at 300 K, which is in the same order of magnitude as those of double strand DNA (45−50 nm)57 and graphene nanoribbons (10−100 nm).58 This result indicates that, for a carbyne chain shorter than Lp, it behaves like a flexible elastic rod; while for a carbyne chain longer than Lp, its properties can only be described statistically. Since the longest carbyne chain observed to date is only 6 nm, its mechanical properties can be well described by a flexible elastic rod model. Then we apply uniaxial tensile strain along the carbyne axis to determine the maximum strength/strain. Uniaxial tension strain is applied by displacing both ends of the chain in opposite directions, at a constant rate of 2.5×10-5 Å/fs. We also select two different strain rates: 2.5×10-4 Å/fs, and 2.5×10-6 Å/fs, and found the stress-strain behavior of the carbyne chain is independent of the strain rate in this range. To obtain a complete stress vs. strain curve, we increase the applied strain until failure occurs, resulting in the maximum strain and stress. In the calculations, stress is calculated using the virial stress formulation. The usual definition of stress faces difficulty for the mono-atomistic chain due to the ambiguity in defining its cross section area. Here, we adopt a diameter of 3.35 Å59 to calculate the cross sectional area of carbyne chain, similar to the approach in defining the cross section of graphene and carbon nanotubes.60,61 To alleviate this ambiguity, we also report the equivalent 1D stress (units of force), which does not need the cross-sectional area. 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 12 of 32

The stress-strain curve for polyyne chain at 500 K is shown in Figure 4a. It is found that the chain was broken when the strain reached about 8.0% with an ultimate stress of 67.1 GPa or a pulling force of 7.53 nN, which is in perfect match with the experimental value of ~7.916 nN reported in previous work 62. According to the linear elastic regime of the curves, the Young’s modulus γ can be calculated by: γ = (dσ/dԑ)|ԑ=0

(4)

where σ is the tensile stress, and ԑ is the extensional strain. Here the Young’s modulus is approximated by using a linear fit of the stress versus strain data up to 2.0% strain. The calculated Young’s modulus of polyyne is about 1345 GPa, which is higher than that of diamond (1220 GPa)63 and that of graphene (1000 GPa)26. The estimated breaking force for polyyne chain is 7.53 nN, or a specific strength of 5.7×107 N·m/kg, which again significantly outperforms graphene (4.7~5.5×107 N·m/kg), carbon nanotubes (4.3~5.0×107 N·m/kg), and diamond (2.5~6.5×107 N·m/kg).64−68 It should be noted that, to calculate the cross sectional area of carbyne chain, the diameter of the chain is adopted as 0.335 nm here, similar to the approach in defining the cross section of graphene and carbon nanotubes. Based on this cross sectional area, a breaking stress (67.1 GPa) is obtained. However, in a previous theoretical and experimental study of the mechanical properties of carbyne chain, the diameter of 0.2 nm was also adopted, which is the cutoff distance in Tersoff-Brenner bond order potential62. If we also adopt 0.2nm to calculate the cross sectional area, the obtained breaking stress is about 200 GPa, which is very close to the experimental value of 245 GPa69.

12 ACS Paragon Plus Environment

Page 13 of 32

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

The Journal of Physical Chemistry

Clearly, the perfect polyyne possesses an ultra-high Young’s modulus. However, as expected, defects in the polyyne chain are able to reduce the mechanical strength. For instance, the polyyne chain with 4% defective bonds (here, two defective bonds existed in the polyyne chain) exhibits the ultimate strain of about 7.3%, and the maximum stress of 66.3 GPa, which are lower than those of pristine polyyne at the same temperature. It is interesting to find that unlike the polyyne chain, the cumulene chain has a totally different mechanical response under strain at 66 K. As shown in Figure 4b, there are obviously three steps observed during the elongation process of cumulene. The first one is a linear elastic process, providing a Young’s modulus of about 1753 GPa, which is 30% higher than that of perfect polyyne. This elastic stage corresponds to the elastic deformations arising from chain uncoiling and stretching. Subsequently, the carbon-carbon double-bonds are being stretched, resulting in the phase transition from cumulene to polyyne, but the carbyne chain is not broken. It is known that the double bonds are made up of one σ bond and one π bond. Before transforming to singleand triple-bonds, π bond must be broken. Because carbon-carbon double-bond is energetically favorable compared to carbon-carbon single-bond, it is energetically costly to break the double bond, thus a large external work is required to overcome this energy barrier. Therefore, in the stress-strain curve, the stress increases sharply until the π bonds break and simultaneously rearrange themselves to form single- and triple-bonds. As the length of two carbon-carbon double-bonds is 2.74Å, while the length of a pair of single- and triple-bond is 2.79Å, there is a rapid decrease in stress during the transition from double-bonds to single- and triple-bonds. This is the origin of the sharp peak in the stress-strain curve. In the third step, since the chain has been transformed to polyyne with alternating single- and triple-bonds, upon further straining, the chain keeps its structure until it is broken at about 11%. The sharp increase in stress to resist the chain 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 14 of 32

deformation provides new physical insight into the bond breaking/rearrangement in monoatomic structures. The above results also indicate that cumulene and polyyne are super strong 1D materials along the axial direction. Thermal Conductivity of Carbyne Chains Finally, we explored the thermal transport in the carbyne chains using equilibrium molecular dynamics (EMD) simulations, which is based on the Green-Kubo formula (GKF) derived from the fluctuation-dissipation theorem and linear response theory. In GKF, the thermal conductivity is related to heat current autocorrelation function (HCACF),70

κ=

V

kBT

∫ J i (0 ) J i (t ) dt , ∞

2 0

(5)

where kB is the Boltzmann constant, V is the volume of the system, and J i (0) J i (t ) is the average heat flux autocorrelation function along the carbyne chain direction, respectively. The heat flux J is calculated by70 N

J = ∑ ε i vi + i

1 N 1 N ∑ (Fij ⋅ vi )rij + ∑ (Fijk ⋅ vi )(rij + rik ), ij ; i ≠ j 2 6 ijk ;i≠ j ≠k

(6)

where, ԑi and vi are the energy and velocity associated with atom i, respectively. Vector rij denotes the interatomic distance between two atoms, and Fij and Fijk denote the two-body and three-body force, respectively. The carbyne (both cumulene and polyyne) chains are constructed with 50 carbon atoms. Periodic boundary condition (PBC) is applied along the axis of the chain. Free boundary condition (FBC) is applied in directions perpendicular to the carbyne chain. Once the chain is relaxed and the temperature reaches the targeted value, the heat flux is collected at each step for 10 ns. After that,

14 ACS Paragon Plus Environment

Page 15 of 32

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

The Journal of Physical Chemistry

the thermal conductivity is calculated according to Eq. 5. It should be noted that, since PBC is adopted along the axis of the chain during the thermal transport simulation, the edge/end effect and length effect on the thermal transport are not considered in our results. The calculated thermal conductivity of cumulene is about 83 W/mK at 480 K, which is comparable to that of the widely used thermal and electric conductor, copper nanowire (75 W/mK for copper nanowire with 20 µm long and 180 nm wide71), but significantly higher than that of semiconductor silicon nanowires72,73. The thermal conductivity of the polyyne chain is about 42 W/mK at 500 K, which is about half of that of cumulene. For the polyyne with 4% defects, its thermal conductivity dramatically decreases to 5.5 W/mK, which is only about 13% of that of pristine polyyne. It should be noted that the thermal conductivity calculated from MD simulations only includes the phonon contribution and does not include the electronic contribution. Since cumulene is metallic, we estimate the electronic contribution to the thermal conductivity of cumulene using the Wiedemann-Franz law: κe = σLwT, where, Lw (= 2(kB/e)2 = 1.45×10-8 V2/K2) is the Lorenz number, T is 480 K, and σ is the electrical conductivity. According to experimental measurements28, σ = ~1.60×102 S/m. Thus, the electronic contribution to the thermal conductivity κe is about 1.11×10-3 W/mK, which is negligible compared to the phonon thermal conductivity. In order to understand the underlying mechanism of this remarkable thermal transport difference among cumulene, polyyne, and polyyne with 4% defects, we calculated the heat current autocorrelation function (HCACF), C JJ (t ) = J i (0) J i (t ) , of these three chains. It is known that the decay of HCACF in bulk material is exponential.74 Also, the initial fast decay is due to

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 32

the high frequency phonon modes while the slow decay arises from the low frequency phonon modes. From the HCACF of cumulene as shown in Figure 5a, it is seen that there is a rapid initial decay followed by a long time decay (15 ps) of the correlation function in cumulene. The HCACF of polyyne is similar to that of cumulene, but its decay time is shorter, which is about 10 ps. In general, high frequency (short wavelength) phonons have a limited contribution to thermal conductivity due to their low group velocity. Thus the large difference in relaxation time for the long wavelength phonons along the chain is responsible for the remarkable difference in thermal conductivity between cumulene and polyyne. It is also seen that the HCACF of polyyne with defective bonds is obviously lower than that of the pristine polyyne. Interestingly it shows a regular high frequency oscillation. We note that a similar oscillation was also observed in heterogeneous systems with different atomic masses.75,76 From Figure 5b, it is also seen that the decay time in the polyyne with defects is remarkably shorter than that of pristine cumulene and polyyne, which is clearly due to the presence of defective bonds. These defected atoms lead to increased phonon scattering, which introduces localized vibrational modes and thus reduces the thermal conductivity. Tunability of thermal conductivity, which can control or rectify heat fluxes, may have important applications in thermal devices.77,78 Recently, a graphene-based thermal modulator is predicted theoretically, in which the heat flux is controlled by applying a tunable external pressure.79 On the other hand, temperature-induced change in thermal conductivity is important for application in thermal memory. Perhaps a better choice for tuning thermal conductivity is through using phase change materials.80 In fact, phase change materials have been widely applied to store and release heat energy due to their high latent heat. Although the latent heat in existing phase change materials can be huge, the change in their thermal conductivity at phase 16 ACS Paragon Plus Environment

Page 17 of 32

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

The Journal of Physical Chemistry

change temperature is usually insignificant. For example, the thermal conductivity of PEG4000, a representative phase change material, is 0.308 W/mK in solid state and drops to 0.216 W/mK after solid-to-liquid phase transition, giving rise to the tunability (high/low ratio) of 143%.81,82 For the cumulene-polyyne phase transition, the thermal conductivity is reduced from 83 W/mK to 42 W/mK, which gives rise to the tunability of 200%. In addition, this type of phase transition can be induced by both temperature and tensile strain, thus carbyne chain is a promising candidate for the realization of mechanically and/or thermally tuned solid-state thermal memories. Conclusions In summary, to fundamentally understand the characteristics of carbyne chain and vigorously expedite its practical applications, here, we performed systematic molecular dynamics simulations to investigate the mechanical and thermal properties of carbyne with two different phase forms. Our results indicate that the critical temperature for carbyne to transform from cumulene to polyyne is 499 K. Both cumulene and polyyne have a high bending stiffness and Young’s modulus, which are comparable with those of the hardest natural materials. Furthermore, cumulene has a fairly high thermal conductivity, which is two times higher than that of polyyne. A small amount of defects in polyyne is able to significantly reduce its thermal conductivity. For example, a 4% of defective bonds can reduce the thermal conductivity dramatically to only 13% of that of pristine polyyne. These unusual mechanical and thermal properties suggest that carbyne is promising for applications in electromechanical switching, mechanical connections and storage of mechanical and thermal energy.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 18 of 32

AUTHOR INFORMATION Corresponding Author E-mail: [email protected]; Phone: 65-64191583 Notes The authors declare no competing financial interests. Acknowledgement: The authors gratefully acknowledge the financial support from the Agency for Science, Technology and Research (A*STAR), Singapore and the use of computing resources at the A*STAR Computational Resource Centre, Singapore.

18 ACS Paragon Plus Environment

Page 19 of 32

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

The Journal of Physical Chemistry

References (1) Hirsch, A. The era of carbon allotropes. Nat. Mater., 2010, 9, 868–871. (2) Geim, A. K.; Novoselov, K. S. The rise of graphene. Nature Materials 2007, 6, 183-191. (3) Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56-58. (4) Kroto, H. W.; Heath, J. R.; O’Brien, S. C.; Curl, R. F.; Smalley, R. E. C60: Buckminsterfullerene. Nature 1985, 318, 162–163. (5) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric field effect in atomically thin carbon films. Science 2004, 306, 666–669. (6) Zhang, Y.; Tan, Y.-W.; Stormer, H. L.; Kim, P. Experimental observation of the quantum Hall effect and Berry's phase in graphene. Nature 2005, 438, 201-204. (7) Baughman, R. Dangerously seeking linear carbon. Science 2006, 312, 1009–1010. (8) Heath, J. R.; Zhang, Q.; O’Brien, S. C.; Curl, R. F.; Kroto, H. W.; Smalley, R. E. The formation of long carbon chain molecules during laser vaporization of graphite. J. Am. Chem. Soc. 1987, 109, 359–363. (9) Gibtner, T., Hampel, F., Gisselbrecht, J. P., Hirsch, A. End-cap stabilized oligoynes: Model compounds for the linear sp carbon allotrope carbyne. Chem-Eur. J. 2002, 8, 408–432. (10)

Cataldo, F. A method for synthesizing polyynes in solution. Carbon 2005, 43, 2792–2800.

(11)

Eisler, S.; Slepkov, A. D.; Elliott, E.; Luu, T.; McDonald, R.; Hegmann, F. A.; Tykwinski,

R. R. Polyynes as a model for carbyne: synthesis, physical properties, and nonlinear optical response. J. Am. Chem. Soc. 2005, 127, 2666–2676. (12)

Jin, C.; Lan, H.; Peng, L.; Suenaga, K.; Iijima, S. Deriving carbon atomic chains from

graphene. Phys. Rev. Lett. 2009, 102, 205501.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(13)

Page 20 of 32

Tykwinski, R. R.; Chalifoux, W.; Eisler, S.; Lucotti, A.; Tommasini, M.; Fazzi, D.; Del

Zoppo, M.; Zerbi, G. Toward carbyne: Synthesis and stability of really long polyynes. Pure. Appl. Chem. 2010, 82, 891–904. (14)

Chalifoux, W. A.; Tykwinski, R. R. Synthesis of polyynes to model the sp-carbon

allotrope carbyne. Nat. Chem. 2010, 2, 967–971. (15)

Bӧrrnert, F.; Bӧrrnert, C.; Gorantla, S.; Liu, X.; Bachmatiuk, A.; Joswig, J.-O.; Wagner,

F. R.; Schäffel, F.; Warner, J. H.; Schönfelder, R.; Rellinghaus, B.; Gemming, T.; Thomas, J.; Knupfer, M.; Büchner, B.; Rümmeli, M. H. Single-wall-carbon-nanotube/single-carbonchain molecular junctions. Phys. Rev. B 2010, 81, 085439. (16)

Song, B.; Schneider, G. F.; Xu, Q.; Pandraud, G.; Dekker, C.; Zandbergen, H. Atomic-

scale electron-beam sculpting of near-defect-free graphene nanostructures. Nano Lett. 2011, 11, 2247–2250. (17)

Robertson, A. W.;Warner, J. H. Atomic resolution imaging of graphene by transmission

electron microscopy. Nanoscale 2013, 5, 4079–4093. (18)

Casillas, G.; Mayoral, A.; Liu, M.; Ponce, A.; Artyukhov, V. I.; Yakobson, B. I.; Jose-

Yacaman, M. New insights into the properties and interactions of carbon chains as revealed by HRTEM and DFT analysis. Carbon 2014, 66, 436–441. (19)

Kano, E.; Takeguchi, M.; Fujita, J.; Hashimoto, A. Direct observation of Pt-terminating

carbyne on graphene. Carbon 2014, 80, 382–386. (20)

Saito, R.; Dresselhaus, G.; Dresselhaus,M. S. Physical Properties of Carbon Nanotubes,

Imperial College Press, London, 1998. (21)

Yakobson, B. I.; Campbell, M. P.; Brabec, C. J.; Bernholc, J. High strain rate fracture and

C-chain unraveling in carbon nanotubes. Comput. Mater. Sci. 1997, 8, 341–348.

20 ACS Paragon Plus Environment

Page 21 of 32

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

The Journal of Physical Chemistry

(22)

Hobi, E.; Pontes, R. B.; Fazzio, A.; da Silva, A. J. R. Formation of atomic carbon chains

from graphene nanoribbons. Phys. Rev. B 2010, 81, 201406. (23)

Erdogan, E.; Popov, I.; Rocha, C. G.; Cuniberti, G.; Roche, S.; Seifert, G. Engineering

carbon chains from mechanically stretched graphene-based materials. Phys. Rev. B 2011, 83, 041401. (24)

Cahangirov, S.; Topsakal, M.; Ciraci, S. Long-range interactions in carbon atomic chains.

Phys.Rev. B 2010, 82, 195444. (25)

Cretu, O.; Botello-Mendez, A.; Janowska, I.; Pham-Huu, C.; Charlier, J.-C.; Banhart, F.

Electrical transport measured in atomic carbon chains. Nano Lett. 2013, 13, 3487–3493. (26)

Liu, M.; Artyukhov, V. I.; Lee, H.; Xu, F.; Yakobson, B. I. Carbyne from first principles:

chain of C atoms, a nanorod or a nanorope. ACS Nano, 2013, 7, 10075–10082. (27)

Artyukhov, V. I.; Liu, M.; Yakobson, B. I. Mechanically induced metal–insulator

transition in carbyne. Nano Lett. 2014, 14, 4224–4229. (28)

Torre, A. La; Botello-Mendez, A.; Baaziz, W.; Charlier, J.-C.; Banhart, F. Strain-induced

metal-semiconductor transition observed in atomic carbon chains. Nat. Commun. 2015, 6, 6636. (29)

Khoo, K. H.; Neaton, J.B.; Son, Y.W.; Cohen, M. L.; Louie, S. G. Negative differential

resistance in carbon atomic wire-carbon nanotube junctions. Nano Lett. 2008, 8, 2900–2905. (30)

Lin, Z. Z.; Ning, X. J. Controlling the electronic properties of monatomic carbon chains.

Epl-Europhys Lett. 2011, 95, 47012. (31)

Hu, Y. H. Bending effect of sp-hybridized carbon (carbyne) chains on their structures and

properties. J. Phys. Chem. C 2011, 115, 1843–1850.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(32)

Page 22 of 32

Zhang, G. Nanoscale Energy Transport and Harvesting: A Computational Study, Pan

Stanford Publishing, 2014. (33)

Li, N.-N.; Sha, Z.-D.; Pei, Q.-X.; Zhang, Y.-W. Hydrogenated grain boundaries control

the strength and ductility of polycrystalline graphene. J. Phys. Chem. C 2014, 118, 13769−13774. (34)

Sha, Z. D.; Quek, S. S.; Pei, Q. X.; Liu, Z. S.; Wang, T. J.; Shenoy, V. B.; Zhang, Y. W.

On the failure load and mechanism of polycrystalline graphene by nanoindentation. Sci. Rep. 2014, 4, 7437. (35)

Strachan, A.; Kober, E. M.; Van Duin, A. C. T.; Oxgaard, J.; Goddard, W. A. Thermal

decomposition of RDX from reactive molecular dynamics. J. Chem. Phys. 2005, 122, 054502. (36)

Chenoweth, K.; Van Duin, A. C. T.; Goddard, W. A. A ReaxFF reactive force field for

molecular dynamics simulations of hydrocarbon oxidation. J. Phys. Chem. A 2008, 112, 1040–1053. (37)

Plimpton, S. J. Fast parallel algorithms for short-range molecular dynamics. J. Comput.

Phys. 1995, 117, 1−19. (38)

Van Duin, A. C. T.; Dasgupta, S.; Lorant, F.; Goddard, W. A. ReaxFF: A reactive force

field for hydrocarbons. J. Phys. Chem. A 2001, 105, 9396–9409. (39)

Tersoff, J. New empirical approach for the structure and energy of covalent systems.

Phys. ReV. B 1988, 37, 6991. (40)

Brenner, D. W. Empirical potential for hydrocarbons for use in simulating the chemical

vapor deposition of diamond films. Phys. ReV. B 1990, 42, 9458.

22 ACS Paragon Plus Environment

Page 23 of 32

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

The Journal of Physical Chemistry

(41)

Roman, R. E.; Kwan, K.; and Cranford, S. W. Mechanical properties and defect

sensitivity of diamond nanothreads. Nano Lett. 2015, 15, 1585−1590. (42)

Nair, A. K.; Cranford, S. W.; Buehler, M. J. The minimal nanowire: Mechanical

properties of carbyne. Europhys. Lett. 2011, 95, 16002. (43)

Cranford, S. Thermal stability of idealized folded carbyne loops. Nanoscale Res. Lett.

2013, 8, 490. (44)

Kocsis, A. J.; Yedama, N. A. R.; Cranford, S. W. Confinement and controlling the

effective compressive stiffness of carbyne. Nanotech. 2014, 25, 335709. (45)

Nosé, S. A unified formulation of the constant temperature molecular-dynamics methods.

J. Chem. Phys. 1984, 81, 511–519. (46)

Hu, J.; Ruan, X.; and Chen, Y. P. Thermal conductivity and thermal rectification in

graphene nanoribbons: a molecular dynamics study. Nano Lett. 2009, 9, 2730–2735. (47)

Zhang, J.; Huang, X.; Yue, Y.; Wang, J.; Wang, X. Dynamic response of graphene to

thermal impulse. Phys. Rev. B 2011, 84, 235416. (48)

Lebedev, B. V., Chapter 5.2 in Heimann R. B.; Evsyukov S. E.; Kavan L. Carbyne and

carbynoid structures, Kluwer academic press, Dordrecht, pp. 317−331, 1999. (49)

Liu, X.; Zhang, G.; Pei, Q.; Zhang, Y.-W. Modulating the thermal conductivity of silicon

nanowires via surface amorphization. Sci. China Tech. Sci. 2014, 57, 699−705. (50)

Allcock, H. R. Introduction to materials chemistry. John Wiley & Sons, Inc., 2011.

(51)

Fan, X.; Liu, L.; Lin, J.; Shen, Z.; and Kuo, J.-L. Density functional theory study of finite

carbon chains. ACS Nano 2009, 3, 3788−3794. (52)

Timoshevskii, A.; Kotrechko, S.; and Matviychuk, Y. Atomic structure and mechanical

properties of carbyne. Phys. Rev. B 2015, 91, 245434.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(53)

Page 24 of 32

Pomerantz, A. E.; Han, J. H.; Musgrave, C. B. Calculating cumulene/poly-yne

isomerization energies. J. Phys. Chem. A 2004, 108, 4030–4035. (54)

Zhang, G.; Musgrave, C. B. Comparison of DFT methods for molecular orbital

eigenvalue calculations. J. Phys. Chem. A 2007, 111, 1554–1561. (55)

Roth, S.; Carroll, D. One-dimensional metals: Conjugated polymers, organics crystals,

carbon nanotubes. John Wiley & Sons Inc. 2004. (56)

Doi, M.; Edwards, S. F. The theory of polymer dynamics. Oxford: Clarendon Press, 1986.

(57)

Hagerman, P. J. Flexibility of DNA. Annu. Rev. Biophys. Biophys. Chem. 1988, 17, 265–

286. (58)

Bets, K. V.; Yakobson, B. I. Spontaneous twist and intrinsic instabilities of pristine

graphene nanoribbons. Nano Res. 2009, 2, 161–166. (59)

Cranford, S.; Sen, D.; Buehler, M. J. Meso-origami: folding multilayer graphene sheets.

Appl. Phys. Lett. 2009, 95, 123121. (60)

Zhang G.; Li, B. Anomalous vibrational energy diffusion in carbon nanotubes. J. Chem.

Phys. 2005,123, 014705. (61)

Liu, X.; Zhang, G.; Zhang, Y.-W. Thermal conduction across graphene cross-linkers. J.

Phys. Chem. C 2014, 118, 12541–12547. (62)

Mazilova, T. I.; Kotrechko, S.; Sadanov, E. V.; Ksenofontov, V. A.; Mikhailovskij, I. M.

High-field formation of linear carbon chains and atomic clusters. Int. J. Nanosci. 2010, 9, 151−157. (63)

Itzhaki, L. Altus, E.; Basch, H.; Hoz, S. Harder than diamond: Determining the cross-

sectional area and Young's modulus of molecular rods. Angewandte Chemie, 2005, 117, 7598–7601.

24 ACS Paragon Plus Environment

Page 25 of 32

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

The Journal of Physical Chemistry

(64)

Lee, C.; Wei, X.; Kysar, J. W.; Hone, J. Measurement of the elastic properties and

intrinsic strength of monolayer graphene. Science 2008, 321, 385–388. (65)

Liu, F.; Ming, P.; Li, J. Ab initio calculation of ideal strength and phonon instability of

graphene under tension. Phys. Rev. B 2007, 76, 064120. (66)

Ogata, S.; Shibutani, Y. Ideal tensile strength and band gap of single-walled carbon

nanotubes. Phys. Rev. B 2003, 68, 165409. (67)

Chang, C.-C.; Hsu, I.-K.; Aykol, M.; Hung, W.-H.; Chen, C.-C.; Cronin, S. B. A new

lower limit for the ultimate breaking strain of carbon nanotubes. ACS Nano 2010, 4, 5095– 5100. (68)

Telling, R. H.; Pickard, C. J.; Payne, M. C.; Field, J. E. Theoretical strength and cleavage

of diamond. Phys. Rev. Lett. 2000, 84, 5160–5163. (69)

Mikhailovskij, I. M.; Sadanov E. V.; Kotrechko, S. A.; Ksenofontov, V. A.;, Mazilova, T.

I. Measurement of the inherent strength of carbon atomic chains. Phys. Rev. B 2013, 87, 045410. (70)

Volz, S. G.; Chen, G. Molecular dynamics simulation of thermal conductivity of silicon

nanowires. Appl. Phys. Lett. 1999, 75, 2056–2058. (71)

Mehta, R.; Chugh, S.; Chen, Z. Enhanced electrical and thermal conduction in graphene-

encapsulated copper nanowires. Nano Lett. 2015, 15, 2024–2030. (72)

Hochbaum, Allon I.; Chen, R.; Delgado, R. D.; Liang, W.; Garnett, E. C.; Najarian, M.;

Majumdar, A.; Yang, P. Enhanced thermoelectric performance of rough silicon nanowires. Nature 2008, 451, 163−167. (73)

Boukai, A.; Bunimovich, I. Y.; Tahir-Kheli, J.; Yu, J. K.; Goddard, W. A.; and Heath, J.

R. Silicon nanowires as efficient thermoelectric materials. Nature 2008, 451, 168−171.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(74)

Page 26 of 32

Che, J.; Cagin, T.;, Deng, W.; Goddard, W. A. Thermal conductivity studies by

molecular dynamics simulations. J. Chem. Phys. 2000, 113, 6888–6900. (75)

Chen, J.; Zhang, G.; Li, B. Phonon coherent resonance and its effect on thermal transport

in core-shell nanowires. J. Chem. Phys. 2011, 135, 104508. (76)

Chen, J.; Zhang, G.; Li, B. Impacts of atomistic coating on thermal conductivity of

germanium nanowires. Nano Lett. 2012, 12, 2826–2832. (77)

Li, N.; Ren, J.; Wang, L.; Zhang, G.; Hänggi, P.; Li, B. Colloquium: Phononics:

manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys. 2012, 84, 1045– 1066. (78)

Tian, H.; Xie, D.; Yang, Y.; Ren, T.-L.; Zhang, G.; Wang, Y.-F.; Zhou, C.-J.; Peng, P.-G.;

Wang, L.-G.; Liu, L.-T. A novel solid-state thermal rectifier based on reduced graphene oxide. Sci. Rep. 2012, 2, 523. (79)

Liu, X.; Zhang, G.; Zhang, Y.-W. Graphene-based thermal modulators. Nano Res. 2015,

8, 2755−2762. (80)

Xie, R.; Bui, C. T.; Varghese, B.; Zhang, Q.; Sow, C. H.; Li, B.; Thong, T. L. An

electrically tuned solid-state thermal memory based on metal–insulator transition of singlecrystalline VO2 nanobeams. Adv. Funct. Mater. 2011, 21, 1602–1607. (81)

Zheng, R., Gao, J., Wang, J. & Chen, G. Reversible temperature regulation of electrical

and thermal conductivity using liquid-solid phase transitions. Nat. Commun. 2011, 2, 289– 295. (82)

Chen, R.; Cui, Y.; Tian, H.; Yao, R.; Liu, Z.; Shu, Y.; Li, C.; Yang, Y.; Ren, T.; Zhang,

G.; Zou, R. Controllable thermal rectification realized in binary phase change composites. Sci. Rep. 2015, 5, 8884.

26 ACS Paragon Plus Environment

Page 27 of 32

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

The Journal of Physical Chemistry

Figure 1. Bond-length distribution of carbyne chains with 50 atoms: (a) Cumulene and (b) Polyyne. Insets show the integral of bond number. It is clear that in polyyne, 50% bonds are single bonds and the rest 50% are triple bonds.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 28 of 32

Figure 2. Phase transition process of carbyne chains with 50 atoms at the critical temperature of 499 K. (a) – (e) Evolution of bond numbers from cumulene to pristine polyyne with a low heating rate of 0.04K/fs. (f) – (j) Evolution of bond numbers from cumulene to polyyne with defective bonds (22 single-bond, 24 triple-bond, and 2 defective bonds) with a heating rate of 2.5 K/fs.

28 ACS Paragon Plus Environment

Page 29 of 32

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

The Journal of Physical Chemistry

Figure 3. (a) A straight carbyne chain (lower panel) and a carbyne ring representing a system with a constant curvature (upper panel). (b) Bending energy of carbyne chain Eb as a function of number of atoms in the chain. (c) Bending energy EL as a function of the carbyne chain curvature ρ. (d) Bending energy EL as a function of ρ2.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 30 of 32

Figure 4. Stress-strain curves of carbyne chains. (a) Polyyne chain with or without defects at 500K. (b) Cumulene chain at 66K.

30 ACS Paragon Plus Environment

Page 31 of 32

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

The Journal of Physical Chemistry

Figure 5. The heat flux autocorrelation functions for carbyne chains. (a) Cumulene and perfect polyyne chains. (b) Polyyne chain with 4% defective bonds. The temperature of cumulene is 480K, and temperature of polyyne is at 500K.

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Table of contents graphic

32 ACS Paragon Plus Environment

Page 32 of 32