-
Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
Authors:
Penghua Ying,
Wenjiang Zhou,
Lucas Svensson,
Erik Fransson,
Fredrik Eriksson,
Ke Xu,
Ting Liang,
Bai Song,
Shunda Chen,
Paul Erhart,
Zheyong Fan
Abstract:
Path-integral molecular dynamics (PIMD) simulations are crucial for accurately capturing nuclear quantum effects in materials. However, their computational intensity and reliance on multiple software packages often limit their applicability at large scales. Here, we present an integration of PIMD methods, including thermostatted ring-polymer molecular dynamics (TRPMD), into the open-source GPUMD p…
▽ More
Path-integral molecular dynamics (PIMD) simulations are crucial for accurately capturing nuclear quantum effects in materials. However, their computational intensity and reliance on multiple software packages often limit their applicability at large scales. Here, we present an integration of PIMD methods, including thermostatted ring-polymer molecular dynamics (TRPMD), into the open-source GPUMD package, combined with highly accurate and efficient machine-learned neuroevolution potential (NEP) models. This approach achieves almost the accuracy of first-principles calculations with the computational efficiency of empirical potentials, enabling large-scale atomistic simulations that incorporate nuclear quantum effects. We demonstrate the efficacy of the combined NEP-PIMD approach by examining various thermal properties of diverse materials, including lithium hydride (LiH), three porous metal-organic frameworks (MOFs), and elemental aluminum. For LiH, our NEP-PIMD simulations successfully capture the isotope effect, reproducing the experimentally observed dependence of the lattice parameter on the reduced mass. For MOFs, our results reveal that achieving good agreement with experimental data requires consideration of both nuclear quantum effects and dispersive interactions. For aluminum, the TRPMD method effectively captures thermal expansion and phonon properties, aligning well with quantum mechanical predictions. This efficient NEP-PIMD approach opens new avenues for exploring complex material properties influenced by nuclear quantum effects, with potential applications across a broad range of materials.
△ Less
Submitted 6 September, 2024;
originally announced September 2024.
-
Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
Authors:
Haikuan Dong,
Yongbo Shi,
Penghua Ying,
Ke Xu,
Ting Liang,
Yanzhou Wang,
Zezhu Zeng,
Xin Wu,
Wenjiang Zhou,
Shiyun Xiong,
Shunda Chen,
Zheyong Fan
Abstract:
Molecular dynamics (MD) simulations play an important role in understanding and engineering heat transport properties of complex materials. An essential requirement for reliably predicting heat transport properties is the use of accurate and efficient interatomic potentials. Recently, machine-learned potentials (MLPs) have shown great promise in providing the required accuracy for a broad range of…
▽ More
Molecular dynamics (MD) simulations play an important role in understanding and engineering heat transport properties of complex materials. An essential requirement for reliably predicting heat transport properties is the use of accurate and efficient interatomic potentials. Recently, machine-learned potentials (MLPs) have shown great promise in providing the required accuracy for a broad range of materials. In this mini review and tutorial, we delve into the fundamentals of heat transport, explore pertinent MD simulation methods, and survey the applications of MLPs in MD simulations of heat transport. Furthermore, we provide a step-by-step tutorial on developing MLPs for highly efficient and predictive heat transport simulations, utilizing the neuroevolution potentials (NEPs) as implemented in the GPUMD package. Our aim with this mini review and tutorial is to empower researchers with valuable insights into cutting-edge methodologies that can significantly enhance the accuracy and efficiency of MD simulations for heat transport studies.
△ Less
Submitted 24 April, 2024; v1 submitted 29 January, 2024;
originally announced January 2024.
-
Correcting force error-induced underestimation of lattice thermal conductivity in machine learning molecular dynamics
Authors:
Xiguang Wu,
Wenjiang Zhou,
Haikuang Dong,
Penghua Ying,
Yanzhou Wang,
Bai Song,
Zheyong Fan,
Shiyun Xiong
Abstract:
Machine learned potentials (MLPs) have been widely employed in molecular dynamics (MD) simulations to study thermal transport. However, literature results indicate that MLPs generally underestimate the lattice thermal conductivity (LTC) of typical solids. Here, we quantitatively analyze this underestimation in the context of the neuroevolution potential (NEP), which is a representative MLP that ba…
▽ More
Machine learned potentials (MLPs) have been widely employed in molecular dynamics (MD) simulations to study thermal transport. However, literature results indicate that MLPs generally underestimate the lattice thermal conductivity (LTC) of typical solids. Here, we quantitatively analyze this underestimation in the context of the neuroevolution potential (NEP), which is a representative MLP that balances efficiency and accuracy. Taking crystalline silicon, GaAs, graphene, and PbTe as examples, we reveal that the fitting errors in the machine-learned forces against the reference ones are responsible for the underestimated LTC as they constitute external perturbations to the interatomic forces. Since the force errors of a NEP model and the random forces in the Langevin thermostat both follow a Gaussian distribution, we propose an approach to correcting the LTC by intentionally introducing different levels of force noises via the Langevin thermostat and then extrapolating to the limit of zero force error. Excellent agreement with experiments is obtained by using this correction for all the prototypical materials over a wide range of temperatures. Based on spectral analyses, we find that the LTC underestimation mainly arises from increased phonon scatterings in the low-frequency region caused by the random force errors.
△ Less
Submitted 26 May, 2024; v1 submitted 21 January, 2024;
originally announced January 2024.
-
General-purpose machine-learned potential for 16 elemental metals and their alloys
Authors:
Keke Song,
Rui Zhao,
Jiahui Liu,
Yanzhou Wang,
Eric Lindgren,
Yong Wang,
Shunda Chen,
Ke Xu,
Ting Liang,
Penghua Ying,
Nan Xu,
Zhiqiang Zhao,
Jiuyang Shi,
Junjie Wang,
Shuang Lyu,
Zezhu Zeng,
Shirong Liang,
Haikuan Dong,
Ligang Sun,
Yue Chen,
Zhuhua Zhang,
Wanlin Guo,
Ping Qian,
Jian Sun,
Paul Erhart
, et al. (3 additional authors not shown)
Abstract:
Machine-learned potentials (MLPs) have exhibited remarkable accuracy, yet the lack of general-purpose MLPs for a broad spectrum of elements and their alloys limits their applicability. Here, we present a feasible approach for constructing a unified general-purpose MLP for numerous elements, demonstrated through a model (UNEP-v1) for 16 elemental metals and their alloys. To achieve a complete repre…
▽ More
Machine-learned potentials (MLPs) have exhibited remarkable accuracy, yet the lack of general-purpose MLPs for a broad spectrum of elements and their alloys limits their applicability. Here, we present a feasible approach for constructing a unified general-purpose MLP for numerous elements, demonstrated through a model (UNEP-v1) for 16 elemental metals and their alloys. To achieve a complete representation of the chemical space, we show, via principal component analysis and diverse test datasets, that employing one-component and two-component systems suffices. Our unified UNEP-v1 model exhibits superior performance across various physical properties compared to a widely used embedded-atom method potential, while maintaining remarkable efficiency. We demonstrate our approach's effectiveness through reproducing experimentally observed chemical order and stable phases, and large-scale simulations of plasticity and primary radiation damage in MoTaVW alloys. This work represents a significant leap towards a unified general-purpose MLP encompassing the periodic table, with profound implications for materials science.
△ Less
Submitted 12 June, 2024; v1 submitted 8 November, 2023;
originally announced November 2023.
-
Dissimilar thermal transport properties in $κ$-Ga$_2$O$_3$ and $β$-Ga$_2$O$_3$ revealed by machine-learning homogeneous nonequilibrium molecular dynamics simulations
Authors:
Xiaonan Wang,
Jinfeng Yang,
Penghua Ying,
Zheyong Fan,
Jin Zhang,
Huarui Sun
Abstract:
The lattice thermal conductivity (LTC) of Ga$_2$O$_3$ is an important property due to the challenge in the thermal management of high-power devices. We develop machine-learned neuroevolution potentials for single-crystalline $β$-Ga$_2$O$_3$ and $κ$-Ga$_2$O$_3$, and apply them to perform homogeneous nonequilibrium molecular dynamics simulations to predict their LTCs. The LTC of $β$-Ga$_2$O$_3$ was…
▽ More
The lattice thermal conductivity (LTC) of Ga$_2$O$_3$ is an important property due to the challenge in the thermal management of high-power devices. We develop machine-learned neuroevolution potentials for single-crystalline $β$-Ga$_2$O$_3$ and $κ$-Ga$_2$O$_3$, and apply them to perform homogeneous nonequilibrium molecular dynamics simulations to predict their LTCs. The LTC of $β$-Ga$_2$O$_3$ was determined to be 10.3 $\pm$ 0.2 W/(m K), 19.9 $\pm$ 0.2 W/(m K), and 12.6 $\pm$ 0.2 W/(m K) along [100], [010], and [001], respectively, aligning with previous experimental measurements. For the first time, we predict the LTC of $κ$-Ga$_2$O$_3$ along [100], [010], and [001] to be 4.5 $\pm$ 0.0 W/(m K), 3.9 $\pm$ 0.0 W/(m K), and 4.0 $\pm$ 0.1 W/(m K), respectively, showing a nearly isotropic thermal transport property. The reduced LTC of $κ$-Ga$_2$O$_3$ versus $β$-Ga$_2$O$_3$ stems from its restricted low-frequency phonons up to 5 THz. Furthermore, we find that the $β$ phase exhibits a typical temperature dependence slightly stronger than $\sim T^{-1}$, whereas the $κ$ phase shows a weaker temperature dependence, ranging from $\sim T^{-0.5}$ to $\sim T^{-0.7}$.
△ Less
Submitted 2 November, 2023;
originally announced November 2023.
-
Combining linear-scaling quantum transport and machine-learning molecular dynamics to study thermal and electronic transports in complex materials
Authors:
Zheyong Fan,
Yang Xiao,
Yanzhou Wang,
Penghua Ying,
Shunda Chen,
Haikuan Dong
Abstract:
We propose an efficient approach for simultaneous prediction of thermal and electronic transport properties in complex materials. Firstly, a highly efficient machine-learned neuroevolution potential is trained using reference data from quantum-mechanical density-functional theory calculations. This trained potential is then applied in large-scale molecular dynamics simulations, enabling the genera…
▽ More
We propose an efficient approach for simultaneous prediction of thermal and electronic transport properties in complex materials. Firstly, a highly efficient machine-learned neuroevolution potential is trained using reference data from quantum-mechanical density-functional theory calculations. This trained potential is then applied in large-scale molecular dynamics simulations, enabling the generation of realistic structures and accurate characterization of thermal transport properties. In addition, molecular dynamics simulations of atoms and linear-scaling quantum transport calculations of electrons are coupled to account for the electron-phonon scattering and other disorders that affect the charge carriers governing the electronic transport properties. We demonstrate the usefulness of this unified approach by studying thermoelectric transport properties of a graphene antidot lattice.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
Mechanisms of temperature-dependent thermal transport in amorphous silica from machine-learning molecular dynamics
Authors:
Ting Liang,
Penghua Ying,
Ke Xu,
Zhenqiang Ye,
Chao Ling,
Zheyong Fan,
Jianbin Xu
Abstract:
Amorphous silica (a-SiO$_2$) is a foundational disordered material for which the thermal transport properties are important for various applications. To accurately model the interatomic interactions in classical molecular dynamics (MD) simulations of thermal transport in a-SiO$_2$, we herein develop an accurate yet highly efficient machine-learned potential model that allowed us to generate a-SiO…
▽ More
Amorphous silica (a-SiO$_2$) is a foundational disordered material for which the thermal transport properties are important for various applications. To accurately model the interatomic interactions in classical molecular dynamics (MD) simulations of thermal transport in a-SiO$_2$, we herein develop an accurate yet highly efficient machine-learned potential model that allowed us to generate a-SiO$_2$ samples closely resembling experimentally produced ones. Using the homogeneous nonequilibrium MD method and a proper quantum-statistical correction to the classical MD results, quantitative agreement with experiments is achieved for the thermal conductivities of bulk and 190 nm-thick a-SiO$_2$ films over a wide range of temperatures. To interrogate the thermal vibrations at different temperatures, we calculated the current correlation functions corresponding to the transverse acoustic (TA) and longitudinal acoustic (LA) collective vibrations. The results reveal that below the Ioffe-Regel crossover frequency, phonons as well-defined excitations, remain applicable in a-SiO$_2$ and play a predominant role at low temperatures, resulting in a temperature-dependent increase in thermal conductivity. In the high-temperature region, more phonons are excited, accompanied by a more intense liquid-like diffusion event. We attribute the temperature-independent thermal conductivity in the high-temperature range of a-SiO$_2$ to the collaborative involvement of excited phonon scattering and liquid-like diffusion in heat conduction. These findings provide physical insights into the thermal transport of a-SiO$_2$ and are expected to be applied to a vast range of amorphous materials.
△ Less
Submitted 1 November, 2023; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Combining the D3 dispersion correction with the neuroevolution machine-learned potential
Authors:
Penghua Ying,
Zheyong Fan
Abstract:
Machine-learned potentials (MLPs) have become a popular approach of modelling interatomic interactions in atomistic simulations, but to keep the computational cost under control, a relatively short cutoff must be imposed, which put serious restrictions on the capability of the MLPs for modelling relatively long-ranged dispersion interactions. In this paper, we propose to combine the neuroevolution…
▽ More
Machine-learned potentials (MLPs) have become a popular approach of modelling interatomic interactions in atomistic simulations, but to keep the computational cost under control, a relatively short cutoff must be imposed, which put serious restrictions on the capability of the MLPs for modelling relatively long-ranged dispersion interactions. In this paper, we propose to combine the neuroevolution potential (NEP) with the popular D3 correction to achieve a unified NEP-D3 model that can simultaneously model relatively short-ranged bonded interactions and relatively long-ranged dispersion interactions. We show the improved descriptions of the binding and sliding energies in bilayer graphene can be obtained by the NEP-D3 approach compared to the pure NEP approach. We implement the D3 part into the GPUMD package such that it can be used out of the box for many exchange-correlation functionals. As a realistic application, we show that dispersion interactions result in approximately a 10% reduction in thermal conductivity for three typical metal-organic frameworks.
△ Less
Submitted 8 October, 2023;
originally announced October 2023.
-
Sub-micrometer phonon mean free paths in metal-organic frameworks revealed by machine-learning molecular dynamics simulations
Authors:
Penghua Ying,
Ting Liang,
Ke Xu,
Jin Zhang,
Jianbin Xu,
Zheng Zhong,
Zheyong Fan
Abstract:
Metal-organic frameworks (MOFs) are a family of materials that have high porosity and structural tunability and hold great potential in various applications, many of which requiring a proper understanding of the thermal transport properties. Molecular dynamics (MD) simulations play an important role in characterizing the thermal transport properties of various materials. However, due to the comple…
▽ More
Metal-organic frameworks (MOFs) are a family of materials that have high porosity and structural tunability and hold great potential in various applications, many of which requiring a proper understanding of the thermal transport properties. Molecular dynamics (MD) simulations play an important role in characterizing the thermal transport properties of various materials. However, due to the complexity of the structures, it is difficult to construct accurate empirical interatomic potentials for reliable MD simulations of MOFs. To this end, we develop a set of accurate yet highly efficient machine-learned potentials for three typical MOFs, including MOF-5, HKUST-1, and ZIF-8, using the neuroevolution potential approach as implemented in the GPUMD package, and perform extensive MD simulations to study thermal transport in the three MOFs. Although the lattice thermal conductivity (LTC) values of the three MOFs are all predicted to be smaller than 1 $\rm{W/(m\ K)}$ at room temperature, the phonon mean free paths (MFPs) are found to reach the sub-micrometer scale in the low-frequency region. As a consequence, the apparent LTC only converges to the diffusive limit for micrometer single crystals, which means that the LTC is heavily reduced in nanocrystalline MOFs. The sub-micrometer phonon MFPs are also found to be correlated with a moderate temperature dependence of LTC between those in typical crystalline and amorphous materials. Both the large phonon MFPs and the moderate temperature dependence of LTC fundamentally change our understanding of thermal transport in MOFs.
△ Less
Submitted 3 June, 2023;
originally announced June 2023.
-
Anisotropic and high thermal conductivity in monolayer quasi-hexagonal fullerene: A comparative study against bulk phase fullerene
Authors:
Haikuan Dong,
Chenyang Cao,
Penghua Ying,
Zheyong Fan,
Ping Qian,
Yanjing Su
Abstract:
Recently a novel two-dimensional (2D) C$_{60}$ based crystal called quasi-hexagonal-phase fullerene (QHPF) has been fabricated and demonstrated to be a promising candidate for 2D electronic devices [Hou et al. Nature 606, 507-510 (2022)]. We construct an accurate and transferable machine-learned potential to study heat transport and related properties of this material, with a comparison to the fac…
▽ More
Recently a novel two-dimensional (2D) C$_{60}$ based crystal called quasi-hexagonal-phase fullerene (QHPF) has been fabricated and demonstrated to be a promising candidate for 2D electronic devices [Hou et al. Nature 606, 507-510 (2022)]. We construct an accurate and transferable machine-learned potential to study heat transport and related properties of this material, with a comparison to the face-centered-cubic bulk-phase fullerene (BPF). Using the homogeneous nonequilibrium molecular dynamics and the related spectral decomposition methods, we show that the thermal conductivity in QHPF is anisotropic, which is 137(7) W/mK at 300 K in the direction parallel to the cycloaddition bonds and 102(3) W/mK in the perpendicular in-plane direction. By contrast, the thermal conductivity in BPF is isotropic and is only 0.45(5) W/mK. We show that the inter-molecular covalent bonding in QHPF plays a crucial role in enhancing the thermal conductivity in QHPF as compared to that in BPF. The heat transport properties as characterized in this work will be useful for the application of QHPF as novel 2D electronic devices.
△ Less
Submitted 9 February, 2023; v1 submitted 8 August, 2022;
originally announced August 2022.
-
GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
Authors:
Zheyong Fan,
Yanzhou Wang,
Penghua Ying,
Keke Song,
Junjie Wang,
Yong Wang,
Zezhu Zeng,
Ke Xu,
Eric Lindgren,
J. Magnus Rahm,
Alexander J. Gabourie,
Jiahui Liu,
Haikuan Dong,
Jianyang Wu,
Yue Chen,
Zheng Zhong,
Jian Sun,
Paul Erhart,
Yanjing Su,
Tapio Ala-Nissila
Abstract:
We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in [Fan et al., Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package GPUMD. We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev…
▽ More
We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in [Fan et al., Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package GPUMD. We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev basis functions and by extending the angular descriptor with some four-body and five-body contributions as in the atomic cluster expansion approach. We also detail our efficient implementation of the NEP approach in graphics processing units as well as our workflow for the construction of NEP models, and we demonstrate their application in large-scale atomistic simulations. By comparing to state-of-the-art MLPs, we show that the NEP approach not only achieves above-average accuracy but also is far more computationally efficient. These results demonstrate that the GPUMD package is a promising tool for solving challenging problems requiring highly accurate, large-scale atomistic simulations. To enable the construction of MLPs using a minimal training set, we propose an active-learning scheme based on the latent space of a pre-trained NEP model. Finally, we introduce three separate Python packages, GPYUMD, CALORINE, and PYNEP, which enable the integration of GPUMD into Python workflows.
△ Less
Submitted 29 June, 2022; v1 submitted 20 May, 2022;
originally announced May 2022.
-
Investigation of Deformation and Fracture Mechanisms in Two-dimensional Gallium Telluride Multilayers Using Nanoindentation
Authors:
Yan Zhou,
Shi Zhou,
Penghua Ying,
Qinghua Zhao,
Yong Xie,
Mingming Gong,
Pisu Jiang,
Hui Cai,
Bin Chen,
Sefaattin Tongay,
Wanqi Jie,
Jin Zhang,
Tao Wang,
Dong Liu,
Martin Kuball
Abstract:
Two-dimensional (2D) materials possess great potential for flexible devices, ascribing to their outstanding electrical, optical, and mechanical properties. However, their mechanical deformation property and fracture mechanism, which are inescapable in many applications like flexible optoelectronics, are still unclear or not thoroughly investigated due methodology limitations. In light of this, suc…
▽ More
Two-dimensional (2D) materials possess great potential for flexible devices, ascribing to their outstanding electrical, optical, and mechanical properties. However, their mechanical deformation property and fracture mechanism, which are inescapable in many applications like flexible optoelectronics, are still unclear or not thoroughly investigated due methodology limitations. In light of this, such mechanical properties and mechanisms are explored on example of gallium telluride (GaTe), a promising optoelectronic candidate with an ultrahigh photo-responsibility and a high plasticity within 2D family. Considering the driving force insufficient in atomic force microscopy (AFM)-based nanoindentation method, here the mechanical properties of both substrate-supported and suspended GaTe multilayers were systematically investigated through full-scale Berkovich-tip nanoindentation, micro-Raman spectroscopy, AFM, and scanning electron microscopy. An unusual concurrence of multiple pop-in and load-drop events in loading curve was observed. By further correlating to molecular dynamics calculations, this concurrence was unveiled originating from the interlayer sliding mediated layers-by-layers fracture mechanism within GaTe multilayers. The van der Waals force between GaTe multilayers and substrates was revealed much stronger than that between GaTe interlayers, resulting in the easy sliding and fracture of multilayers within GaTe. This work provides new insights into the deformation and fracture mechanisms of GaTe and other similar 2D multilayers in flexible applications.
△ Less
Submitted 23 April, 2022;
originally announced April 2022.
-
$s$-wave superconductivity in the noncentrosymmetric W$_3$Al$_2$C superconductor: An NMR study
Authors:
D. Tay,
T. Shang,
Y. P. Qi,
T. P. Ying,
H. Hosono,
H-R. Ott,
T. Shiroka
Abstract:
We report on a microscopic study of the noncentrosymmetric superconductor W$_3$Al$_2$C (with $T_c$ = 7.6 K), mostly by means of $^{27}$Al- and $^{13}$C nuclear magnetic resonance (NMR). Since in this material the density of states at the Fermi level is dominated by the tungsten's 5$d$ orbitals, we expect a sizeable spin-orbit coupling (SOC) effect. The normal-state electronic properties of W$_3$Al…
▽ More
We report on a microscopic study of the noncentrosymmetric superconductor W$_3$Al$_2$C (with $T_c$ = 7.6 K), mostly by means of $^{27}$Al- and $^{13}$C nuclear magnetic resonance (NMR). Since in this material the density of states at the Fermi level is dominated by the tungsten's 5$d$ orbitals, we expect a sizeable spin-orbit coupling (SOC) effect. The normal-state electronic properties of W$_3$Al$_2$C resemble those of a standard metal, but with a Korringa product $1/(T_{1}T)$ significantly smaller than that of metallic Al, reflecting the marginal role played by $s$-electrons. In the superconducting state, we observe a reduction of the Knight shift and an exponential decrease of the NMR relaxation rate $1/T_1$, typical of $s$-wave superconductivity. This is further supported by the observation of a small but distinct coherence peak just below $T_c$ in the $^{13}$C NMR relaxation-rate, in agreement with the fully-gapped superconducting state inferred from the electronic specific-heat data well below $T_c$. The above features are compared to those of members of the same family, in particular, Mo$_3$Al$_2$C, often claimed to exhibit unconventional superconductivity. We discuss why, despite the enhanced SOC, W$_3$Al$_2$C does not show spin-triplet features in its superconducting state and consider the broader consequences of our results for noncentrosymmetric superconductors in general.
△ Less
Submitted 6 February, 2022;
originally announced February 2022.
-
Narrow-gap Semiconducting Superhard Amorphous Carbon with Superior Toughness
Authors:
Shuangshuang Zhang,
Yingju Wu,
Kun Luo,
Bing Liu,
Yu Shu,
Yang Zhang,
Lei Sun,
Yufei Gao,
Mengdong Ma,
Zihe Li,
Baozhong Li,
Pan Ying,
Zhisheng Zhao,
Wentao Hu,
Vicente Benavides,
Olga P. Chernogorova,
Alexander V. Soldatov,
Julong He,
Dongli Yu,
Bo Xu,
Yongjun Tian
Abstract:
New carbon forms exhibiting extraordinary physico-chemical properties can be generated from nanostructured precursors under extreme pressure. Nevertheless, synthesis of such fascinating materials is often not well understood that results, as is the case of C60 precursor, in irreproducibility of the results and impeding further progress in the materials design. Here the semiconducting amorphous car…
▽ More
New carbon forms exhibiting extraordinary physico-chemical properties can be generated from nanostructured precursors under extreme pressure. Nevertheless, synthesis of such fascinating materials is often not well understood that results, as is the case of C60 precursor, in irreproducibility of the results and impeding further progress in the materials design. Here the semiconducting amorphous carbon having bandgaps of 0.1-0.3 eV and the advantages of isotropic superhardness and superior toughness over single-crystal diamond and inorganic glasses are produced from transformation of fullerene at high pressure and moderate temperatures. A systematic investigation of the structure and bonding evolution was carried out by using rich arsenal of complimentary characterization methods, which helps to build a model of the transformation that can be used in further high p,T synthesis of novel nanocarbon systems for advanced applications. The produced amorphous carbon materials have the potential of demanding optoelectronic applications that diamond and graphene cannot achieve
△ Less
Submitted 15 June, 2021;
originally announced June 2021.
-
Superconducting gap symmetry of the noncentrosymmetric superconductor W3Al2C
Authors:
R. Gupta,
T. P. Ying,
Y. P. Qi,
H. Hosono,
R. Khasanov
Abstract:
A detailed zero-field and transverse-field muon spin relaxation/rotation ($μ$SR) experiemnts have been carried out on the recently discovered non-centrosymmetric superconductor W$_3$Al$_2$C to speculate about its superconducting ground state. Bulk nature of superconductivity below 7.6 K is confirmed through magnetization measurements. No change in the $μ$SR spectra collected above and below $T_c$…
▽ More
A detailed zero-field and transverse-field muon spin relaxation/rotation ($μ$SR) experiemnts have been carried out on the recently discovered non-centrosymmetric superconductor W$_3$Al$_2$C to speculate about its superconducting ground state. Bulk nature of superconductivity below 7.6 K is confirmed through magnetization measurements. No change in the $μ$SR spectra collected above and below $T_c$ is visible, ruling out the possibility of spontaneous magnetic field below $T_c$. This confirms that time-reversal symmetry is preserved for W$_3$Al$_2$C upon entering in the superconducting ground state. Temperature dependent superfluid density [$ρ_s(T)$], which directly reflects the superconducting gap symmetry is obtained by the analysis of spectra obtained from the transverse-field $μ$SR experiments. Despite a non-centrosymmetric structure, W$_3$Al$_2$C adopts a fully gaped spin-singlet superconducting ground state with a zero temperature value of gap $Δ_0$ = 1.158(8) meV with gap-to-$T_c$ ratio 2$Δ_0/k_BT_c\approx$3.54, classifying this material as a weakly-coupled superconductors.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Nodeless superconductivity in the SnAs-based van der Waals type superconductor NaSn2As2
Authors:
E. J. Cheng,
J. M. Ni,
F. Q. Meng,
T. P. Ying,
B. L. Pan,
Y. Y. Huang,
Darren Peets,
Q. H. Zhang,
S. Y. Li
Abstract:
We grew the single crystals of the SnAs-based van der Waals (vdW)-type superconductor NaSn$_2$As$_2$ and systematically measured its resistivity, specific heat, and ultralow-temperature thermal conductivity. The superconducting transition temperature $T_c$ = 1.60 K of our single crystal is 0.3 K higher than that previously reported. A weak but intrinsic anomaly situated at 193 K is observed in bot…
▽ More
We grew the single crystals of the SnAs-based van der Waals (vdW)-type superconductor NaSn$_2$As$_2$ and systematically measured its resistivity, specific heat, and ultralow-temperature thermal conductivity. The superconducting transition temperature $T_c$ = 1.60 K of our single crystal is 0.3 K higher than that previously reported. A weak but intrinsic anomaly situated at 193 K is observed in both resistivity and specific heat, which likely arises from a charge-density-wave (CDW) instability. Ultralow-temperature thermal conductivity measurements reveal a fully-gapped superconducting state with a negligible residual linear term in zero magnetic field, and the field dependence of $κ_0 / T$ further suggests NaSn$_2$As$_2$ is an $s$-wave superconductor.
△ Less
Submitted 4 June, 2018;
originally announced June 2018.
-
Discrete superconducting phases in FeSe-derived superconductors
Authors:
T. P. Ying,
M. X. Wang,
Z. Y. Zhao,
Z. Z. Zhang,
X. Y. Jia,
Y. C. Li,
B. Lei,
Q. Li,
Y. Yu,
E. J. Cheng,
Z. H. An,
Y. Zhang,
W. Yang,
X. H. Chen,
S. Y. Li
Abstract:
A general feature of unconventional superconductors is the existence of a superconducting dome in the phase diagram as a function of carrier concentration. For the simplest iron-based superconductor FeSe (with transition temperature Tc ~ 8 K), its Tc can be greatly enhanced by doping electrons via many routes, even up to 65 K in monolayer FeSe/SiTiO3. However, a clear phase diagram with carrier co…
▽ More
A general feature of unconventional superconductors is the existence of a superconducting dome in the phase diagram as a function of carrier concentration. For the simplest iron-based superconductor FeSe (with transition temperature Tc ~ 8 K), its Tc can be greatly enhanced by doping electrons via many routes, even up to 65 K in monolayer FeSe/SiTiO3. However, a clear phase diagram with carrier concentration for FeSe-derived superconductors is still lacking. Here, we report the observation of a series of discrete superconducting phases in FeSe thin flakes by continuously tuning carrier concentration through the intercalation of Li and Na ions with a solid ionic gating technique. Such discrete superconducting phases are robust against the substitution of Se by 20% S, but are vulnerable to the substitution of Fe by 2% Cu, highlighting the importance of the iron site being intact. A complete superconducting phase diagram for FeSe-derivatives is given, which is distinct from other unconventional superconductors.
△ Less
Submitted 5 February, 2018;
originally announced February 2018.
-
Electronic structure of FeS
Authors:
J. Miao,
X. H. Niu,
D. F. Xu,
Q. Yao,
Q. Y. Chen,
T. P. Ying,
S. Y. Li,
Y. F. Fang,
J. C. Zhang,
S. Ideta,
K. Tanaka,
B. P. Xie,
D. L. Feng,
Fei Chen
Abstract:
Here we report the electronic structure of FeS, a recently identified iron-based superconductor. Our high-resolution angle-resolved photoemission spectroscopy studies show two hole-like ($α$ and $β$) and two electron-like ($η$ and $δ$) Fermi pockets around the Brillouin zone center and corner, respectively, all of which exhibit moderate dispersion along $k_z$. However, a third hole-like band (…
▽ More
Here we report the electronic structure of FeS, a recently identified iron-based superconductor. Our high-resolution angle-resolved photoemission spectroscopy studies show two hole-like ($α$ and $β$) and two electron-like ($η$ and $δ$) Fermi pockets around the Brillouin zone center and corner, respectively, all of which exhibit moderate dispersion along $k_z$. However, a third hole-like band ($γ$) is not observed, which is expected around the zone center from band calculations and is common in iron-based superconductors. Since this band has the highest renormalization factor and is known to be the most vulnerable to defects, its absence in our data is likely due to defect scattering --- and yet superconductivity can exist without coherent quasiparticles in the $γ$ band. This may help resolve the current controversy on the superconducting gap structure of FeS. Moreover, by comparing the $β$ bandwidths of various iron chalcogenides, including FeS, FeSe$_{1-x}$S$_x$, FeSe, and FeSe$_{1-x}$ Te$_x$, we find that the $β$ bandwidth of FeS is the broadest. However, the band renormalization factor of FeS is still quite large, when compared with the band calculations, which indicates sizable electron correlations. This explains why the unconventional superconductivity can persist over such a broad range of isovalent substitution in FeSe$_{1-x}$Te$_{x}$ and FeSe$_{1-x}$S$_{x}$.
△ Less
Submitted 25 March, 2017;
originally announced March 2017.
-
Fully gapped superconducting state in Au2Pb: a natural candidate for topological superconductor
Authors:
Y. J. Yu,
Y. Xu,
Y. Xing,
J. Zhang,
T. P. Ying,
X. C. Hong,
M. X. Wang,
X. Zhang,
S. Jia,
J. Wang,
S. Y. Li
Abstract:
We measured the ultra-low-temperature specific heat and thermal conductivity of Au$_2$Pb single crystal, a possible three-dimensional Dirac semimetal with a superconducting transition temperature $T_c \approx$ 1.05 K. The electronic specific heat can be fitted by a two-band s-wave model, which gives the gap amplitudes $Δ_1$(0)/$k_BT_c$ = 1.38 and $Δ_2$(0)/$k_BT_c$ = 5.25. From the thermal conducti…
▽ More
We measured the ultra-low-temperature specific heat and thermal conductivity of Au$_2$Pb single crystal, a possible three-dimensional Dirac semimetal with a superconducting transition temperature $T_c \approx$ 1.05 K. The electronic specific heat can be fitted by a two-band s-wave model, which gives the gap amplitudes $Δ_1$(0)/$k_BT_c$ = 1.38 and $Δ_2$(0)/$k_BT_c$ = 5.25. From the thermal conductivity measurements, a negligible residual linear term $κ_0/T$ in zero field and a slow field dependence of $κ_0/T$ at low field are obtained. These results suggest that Au$_2$Pb has a fully gapped superconducting state in the bulk, which is a necessary condition for topological superconductor if Au$_2$Pb is indeed one.
△ Less
Submitted 20 September, 2016;
originally announced September 2016.
-
Observation of two distinct superconducting domes under pressure in tetragonal FeS
Authors:
J. Zhang,
F. L. Liu,
T. P. Ying,
N. N. Li,
Y. Xu,
L. P. He,
X. C. Hong,
Y. J. Yu,
M. X. Wang,
J. Shen,
W. G. Yang,
S. Y. Li
Abstract:
As the simplest iron-based superconductor, FeSe forms a tetragonal structure with transition temperature Tc ~ 8 K. With assistance of pressure, or other techniques, Tc can be greatly enhanced, even to above liquid nitrogen temperature. The newly discovered superconducting tetragonal FeS (Tc ~ 4.5 K), a sulfide counterpart of FeSe, promotes us on its high pressure investigation. The transport and s…
▽ More
As the simplest iron-based superconductor, FeSe forms a tetragonal structure with transition temperature Tc ~ 8 K. With assistance of pressure, or other techniques, Tc can be greatly enhanced, even to above liquid nitrogen temperature. The newly discovered superconducting tetragonal FeS (Tc ~ 4.5 K), a sulfide counterpart of FeSe, promotes us on its high pressure investigation. The transport and structure evolution of FeS with pressure have been studied. A rapid suppression of Tc and vanishing of superconductivity at 4.0 GPa are observed, followed by a second superconducting dome with a 30% enhancement in maximum Tc. An onsite tetragonal to hexagonal phase transition occurs around 7.0 GPa, followed by a broad pressure range of phase coexistence. The residual deformed tetragonal phase is considered as the source of second superconducting dome. The observation of two superconducting domes in iron-based superconductors poses great challenges for understanding their pairing mechanism.
△ Less
Submitted 7 June, 2016; v1 submitted 18 April, 2016;
originally announced April 2016.
-
Nodal superconductivity in FeS: Evidence from quasiparticle heat transport
Authors:
T. P. Ying,
X. F. Lai,
X. C. Hong,
Y. Xu,
L. P. He,
J. Zhang,
M. X. Wang,
Y. J. Yu,
F. Q. Huang,
S. Y. Li
Abstract:
We report low-temperature heat transport measurements on superconducting iron sulfide FeS with $T_c \approx$ 5 K, which has the same crystal structure and similar electronic band structure to the superconducting iron selenide FeSe. In zero magnetic field, a significant residual linear term $κ_0/T$ is observed. At low field, $κ_0/T$ increases rapidly with the increase of field. These results provid…
▽ More
We report low-temperature heat transport measurements on superconducting iron sulfide FeS with $T_c \approx$ 5 K, which has the same crystal structure and similar electronic band structure to the superconducting iron selenide FeSe. In zero magnetic field, a significant residual linear term $κ_0/T$ is observed. At low field, $κ_0/T$ increases rapidly with the increase of field. These results provide strong evidence for nodal superconducting gap in FeS. We compare it with the sister compound FeSe, and other iron-based superconductors with nodal gap.
△ Less
Submitted 13 June, 2016; v1 submitted 24 November, 2015;
originally announced November 2015.
-
Observation of superconductivity at 30 K~46 K in AxFe2Se2 (A = Li, Na, Ba, Sr, Ca, Yb, and Eu)
Authors:
T. P. Ying,
X. L. Chen,
G. Wang,
S. F. Jin,
T. T. Zhou,
X. F. Lai,
H. Zhang,
W. Y. Wang
Abstract:
New iron selenide superconductors by intercalating smaller-sized alkali metals (Li, Na) and alkaline earths using high-temperature routes have been pursued ever since the discovery of superconductivity at about 30 K in KFe2Se2, but all have failed so far. Here we demonstrate that a series of superconductors with enhanced Tc=30~46 K can be obtained by intercalating metals, Li, Na, Ba, Sr, Ca, Yb, a…
▽ More
New iron selenide superconductors by intercalating smaller-sized alkali metals (Li, Na) and alkaline earths using high-temperature routes have been pursued ever since the discovery of superconductivity at about 30 K in KFe2Se2, but all have failed so far. Here we demonstrate that a series of superconductors with enhanced Tc=30~46 K can be obtained by intercalating metals, Li, Na, Ba, Sr, Ca, Yb, and Eu in between FeSe layers by the ammonothermal method at room temperature. Analysis on their powder X-ray diffraction patterns reveals that all the main phases can be indexed based on body-centered tetragonal lattices with a~3.755-3.831 Å while c~15.99-20.54 Å. Resistivities show the corresponding sharp transitions at 45 K and 39 K for NaFe2Se2 and Ba0.8Fe2Se2, respectively, confirming their bulk superconductivity. These findings provide a new starting point for studying the properties of these superconductors and an effective synthetic route for the exploration of new superconductors as well.
△ Less
Submitted 29 May, 2012; v1 submitted 20 February, 2012;
originally announced February 2012.
-
Finite Size Scaling, Fisher Zeroes and N=4 Super Yang-Mills
Authors:
P. R. Crompton,
W. Janke,
Z. X. Xu,
H. P. Ying
Abstract:
We investigate critical slowing down in the local updating continuous-time Quantum Monte Carlo method by relating the finite size scaling of Fisher Zeroes to the dynamically generated gap, through the scaling of their respective critical exponents. As we comment, the nonlinear sigma model representation derived through the hamiltonian of our lattice spin model can also be used to give a effectiv…
▽ More
We investigate critical slowing down in the local updating continuous-time Quantum Monte Carlo method by relating the finite size scaling of Fisher Zeroes to the dynamically generated gap, through the scaling of their respective critical exponents. As we comment, the nonlinear sigma model representation derived through the hamiltonian of our lattice spin model can also be used to give a effective treatment of planar anomalous dimensions in N=4 SYM. We present scaling arguments from our FSS analysis to discuss quantum corrections and recent 2-loop results, and further comment on the prospects of extending this approach for calculating higher twist parton distributions.
△ Less
Submitted 25 August, 2004;
originally announced August 2004.
-
Critical dynamics and universality of the random-bond Potts ferromagnet with tri-distributed quenched disorders
Authors:
H. P. Ying,
B. J. Bian,
D. R. Ji,
H. J. Luo,
L. Schuelke
Abstract:
Critical behavior in short-time dynamics is investigated by a Monte Carlo study for the random-bond Potts ferromagnet with a trinary distribution of quenched disorders on two-dimensional triangular lattices. The dynamic scaling is verified and applied to estimate critical exponents $θ$, $z$ and $β/ν$ for several realizations of the trinary distribution. Our critical scaling analysis strongly ind…
▽ More
Critical behavior in short-time dynamics is investigated by a Monte Carlo study for the random-bond Potts ferromagnet with a trinary distribution of quenched disorders on two-dimensional triangular lattices. The dynamic scaling is verified and applied to estimate critical exponents $θ$, $z$ and $β/ν$ for several realizations of the trinary distribution. Our critical scaling analysis strongly indicates that the bond randomness influences the critical universality.
△ Less
Submitted 17 March, 2001; v1 submitted 16 March, 2001;
originally announced March 2001.
-
Short-time dynamics of random-bond Potts ferromagnet with continuous self-dual quenched disorders
Authors:
Z. Q. Pan,
H. P. Ying,
D. W. Gu
Abstract:
We present Monte Carlo simulation results of random-bond Potts ferromagnet with the Olson-Young self-dual distribution of quenched disorders in two-dimensions. By exploring the short-time scaling dynamics, we find universal power-law critical behavior of the magnetization and Binder cumulant at the critical point, and thus obtain estimates of the dynamic exponent $z$ and magnetic exponent $η$, a…
▽ More
We present Monte Carlo simulation results of random-bond Potts ferromagnet with the Olson-Young self-dual distribution of quenched disorders in two-dimensions. By exploring the short-time scaling dynamics, we find universal power-law critical behavior of the magnetization and Binder cumulant at the critical point, and thus obtain estimates of the dynamic exponent $z$ and magnetic exponent $η$, as well as the exponent $θ$. Our special attention is paid to the dynamic process for the $q$=8 Potts model.
△ Less
Submitted 6 March, 2001;
originally announced March 2001.
-
Corrections to Scaling for the Two-dimensional Dynamic XY Model
Authors:
H. P. Ying,
B. Zheng,
Y. Yu,
S. Trimper
Abstract:
With large-scale Monte Carlo simulations, we confirm that for the two-dimensional XY model, there is a logarithmic correction to scaling in the dynamic relaxation starting from a completely disordered state, while only an inverse power law correction in the case of starting from an ordered state. The dynamic exponent $z$ is $z=2.04(1)$.
With large-scale Monte Carlo simulations, we confirm that for the two-dimensional XY model, there is a logarithmic correction to scaling in the dynamic relaxation starting from a completely disordered state, while only an inverse power law correction in the case of starting from an ordered state. The dynamic exponent $z$ is $z=2.04(1)$.
△ Less
Submitted 5 March, 2001;
originally announced March 2001.
-
Short-time critical dynamics and universality on a two-dimensional Triangular Lattice
Authors:
H. P. Ying,
L. Wang,
J. B Zhang,
M. Jiang,
J. Hu
Abstract:
Critical scaling and universality in short-time dynamics for spin models on a two-dimensional triangular lattice are investigated by using Monte Carlo simulation. Emphasis is placed on the dynamic evolution from fully ordered initialstates to show that universal scaling exists already in the short-time regime in form of power-law behavior of the magnetization and Binder cumulant. The results mea…
▽ More
Critical scaling and universality in short-time dynamics for spin models on a two-dimensional triangular lattice are investigated by using Monte Carlo simulation. Emphasis is placed on the dynamic evolution from fully ordered initialstates to show that universal scaling exists already in the short-time regime in form of power-law behavior of the magnetization and Binder cumulant. The results measured for the dynamic and static critical exponents, $θ$, $z$, $β$ and $ν$, confirm explicitly that the Potts models on the triangular lattice and square lattice belong to the same universality class. Our critical scaling analysis strongly suggests that the simulation for the dynamic relaxation can be used to determine numerically the universality.
△ Less
Submitted 27 February, 2001;
originally announced February 2001.
-
Dynamic Monte Carlo Study of the Two-Dimensional Quantum XY Model
Authors:
H. P. Ying,
H. J. Luo,
L. Schuelke,
B. Zheng
Abstract:
We present a dynamic Monte Carlo study of the Kosterlitz-Thouless phase transition for the spin-1/2 quantum XY model in two dimensions. The short-time dynamic scaling behaviour is found and the dynamical exponent $θ$, $z$ and the static exponent $η$ are determined at the transition temperature.
We present a dynamic Monte Carlo study of the Kosterlitz-Thouless phase transition for the spin-1/2 quantum XY model in two dimensions. The short-time dynamic scaling behaviour is found and the dynamical exponent $θ$, $z$ and the static exponent $η$ are determined at the transition temperature.
△ Less
Submitted 9 July, 1998;
originally announced July 1998.