Showing 1–50 of 93 results for author: Tanaka, I

Superconducting phase diagram of Sb and Se substitution in CeOBiS2 single crystals

Authors: Tatsuya Suzuki, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sb- and Se-substituted CeOBiS2 single crystals have been successfully grown using CsCl/KCl flux. Sb and Se substitution dependence of the superconductivity on CeOBiS2 was investigated and the superconducting phase diagram at above 0.3 K was described. The non-linear boundary between superconductivity and non-superconductivity was revealed. Sb- and Se-substituted CeOBiS2 single crystals have been successfully grown using CsCl/KCl flux. Sb and Se substitution dependence of the superconductivity on CeOBiS2 was investigated and the superconducting phase diagram at above 0.3 K was described. The non-linear boundary between superconductivity and non-superconductivity was revealed. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 16 pages, 2 figures

Journal ref: Journal of the Physical Society of Japan 93, 054708 (2024)

Imaging quantum interference in a monolayer Kitaev quantum spin liquid candidate

Authors: Y. Kohsaka, S. Akutagawa, S. Omachi, Y. Iwamichi, T. Ono, I. Tanaka, S. Tateishi, H. Murayama, S. Suetsugu, K. Hashimoto, T. Shibauchi, M. O. Takahashi, M. G. Yamada, S. Nikolaev, T. Mizushima, S. Fujimoto, T. Terashima, T. Asaba, Y. Kasahara, Y. Matsuda

Abstract: Single atomic defects are prominent windows to look into host quantum states because collective responses from the host states emerge as localized states around the defects. Friedel oscillations and Kondo clouds in Fermi liquids are quintessential examples. However, the situation is quite different for quantum spin liquid (QSL), an exotic state of matter with fractionalized quasiparticles and topo… ▽ More Single atomic defects are prominent windows to look into host quantum states because collective responses from the host states emerge as localized states around the defects. Friedel oscillations and Kondo clouds in Fermi liquids are quintessential examples. However, the situation is quite different for quantum spin liquid (QSL), an exotic state of matter with fractionalized quasiparticles and topological order arising from a profound impact of quantum entanglement. Elucidating the underlying local electronic property has been challenging due to the charge neutrality of fractionalized quasiparticles and the insulating nature of QSLs. Here, using spectroscopic-imaging scanning tunneling microscopy, we report atomically resolved images of monolayer $αあるふぁ$-RuCl$_3$, the most promising Kitaev QSL candidate, on metallic substrates. We find quantum interference in the insulator manifesting as incommensurate and decaying spatial oscillations of the local density of states around defects with a characteristic bias dependence. The oscillation differs from any known spatial structures in its nature and does not exist in other Mott insulators, implying it is an exotic oscillation involved with excitations unique to $αあるふぁ$-RuCl$_3$. Numerical simulations can reproduce the observed oscillation by assuming that itinerant Majorana fermions of Kitaev QSL are scattered across the Majorana Fermi surface. The oscillation provides a new approach to exploring Kitaev QSLs through the local response against defects like Friedel oscillations in metals. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Predictive power of polynomial machine learning potentials for liquid states in 22 elemental systems

Authors: Hayato Wakai, Atsuto Seko, Hirosato Izuta, Takayuki Nishiyama, Isao Tanaka

Abstract: The polynomial machine learning potentials (MLPs) described by polynomial rotational invariants have been systematically developed for various systems and used in diverse applications in crystalline states. In this study, we systematically investigate the predictive power of the polynomial MLPs for liquid structural properties in 22 elemental systems with diverse chemical bonding properties, inclu… ▽ More The polynomial machine learning potentials (MLPs) described by polynomial rotational invariants have been systematically developed for various systems and used in diverse applications in crystalline states. In this study, we systematically investigate the predictive power of the polynomial MLPs for liquid structural properties in 22 elemental systems with diverse chemical bonding properties, including those showing anomalous melting behavior, such as Si, Ge, and Bi. We compare liquid structural properties obtained from molecular dynamics simulations using the density functional theory (DFT) calculation, the polynomial MLPs, and other interatomic potentials in the literature. The current results demonstrate that the polynomial MLPs consistently exhibit high predictive power for liquid structural properties with the same accuracy as that of typical DFT calculations. △ Less

Submitted 11 June, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: REVTeX 4-2, 30 pages with 35 figures

Growth and Characterization of Superconducting Bulk Crystal [(SnSe)$_{1+δでるた}$]$_m$(NbSe$_2$) Misfit Layer Compounds

Authors: Ryufa Shu, Masanori Nagao, Chiaya Yamamoto, Keisuke Arimoto, Junji Yamanaka, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: [(SnSe)$_{1+δでるた}$]$_m$(NbSe$_2$) ($m$ = 1-6, 8, and 12) highly orientated crystals 1-2 mm in size and well-defined c-planes were successfully grown using CsCl/KCl flux, including the first growth of crystals with $m = 12$. The stacked layers along the $c$ axis in the obtained crystals were directly observed by transmission electron microscopy as m alternating layers of SnSe and single layers of NbSe… ▽ More [(SnSe)$_{1+δでるた}$]$_m$(NbSe$_2$) ($m$ = 1-6, 8, and 12) highly orientated crystals 1-2 mm in size and well-defined c-planes were successfully grown using CsCl/KCl flux, including the first growth of crystals with $m = 12$. The stacked layers along the $c$ axis in the obtained crystals were directly observed by transmission electron microscopy as m alternating layers of SnSe and single layers of NbSe$_2$. The superconducting transition temperature of the obtained [(SnSe)$_{1+δでるた}$]$_m$(NbSe$_2$) crystals decreased with an increase in the number of SnSe layers per unit cell. As the superconducting anisotropy parameters increase, a significant increase is observed between $m = 4$ and 5. This indicates that the superconducting dimensionality becomes more two-dimensional with an increasing $m$. △ Less

Submitted 14 January, 2024; originally announced January 2024.

Report number: Journal of Alloys and Compounds, vol.978 (2024) 173486

Journal ref: Journal of Alloys and Compounds, vol.978 (2024) 173486

Suppression of superconductivity by Sb-substitution into CeOBiS2 single crystals

Authors: Tatsuya Suzuki, Yuji Hanada, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sb-substituted CeOBiS2 single crystals with 0.2-1.0 mm size have been successfully grown using CsCl/KCl flux. Sb substitution strongly suppressed the superconductivity in CeOBiS2, and the substitution of more than approximately 4 atomic percent in Bi-site disappeared the superconductivity at a measurement range of above 0.36 K. Furthermore, the drastic change of c-axis lattice parameter was observ… ▽ More Sb-substituted CeOBiS2 single crystals with 0.2-1.0 mm size have been successfully grown using CsCl/KCl flux. Sb substitution strongly suppressed the superconductivity in CeOBiS2, and the substitution of more than approximately 4 atomic percent in Bi-site disappeared the superconductivity at a measurement range of above 0.36 K. Furthermore, the drastic change of c-axis lattice parameter was observed in a neighborhood to Sb substitution amount, which maybe suggests the structural transition in CeOBiS2. The Ce valence state of the obtained single crystals was evaluated. However, the established correlations between Ce valence state and superconductivity could not be revealed. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Journal ref: Solid State Communications, vol.378 (2024) 115414

Algorithm for spin symmetry operation search

Authors: Kohei Shinohara, Atsushi Togo, Hikaru Watanabe, Takuya Nomoto, Isao Tanaka, Ryotaro Arita

Abstract: A spin space group provides a suitable way to fully exploit the symmetry of a spin arrangement with a negligible spin-orbit coupling. There has been a growing interest in applying spin symmetry analysis with the spin space group in the field of magnetism. However, there is no established algorithm to search for spin symmetry operations of the spin space group. This paper presents an exhaustive alg… ▽ More A spin space group provides a suitable way to fully exploit the symmetry of a spin arrangement with a negligible spin-orbit coupling. There has been a growing interest in applying spin symmetry analysis with the spin space group in the field of magnetism. However, there is no established algorithm to search for spin symmetry operations of the spin space group. This paper presents an exhaustive algorithm for determining spin symmetry operations of commensurate spin arrangements. The present algorithm searches for spin symmetry operations from the symmetry operations of a corresponding nonmagnetic crystal structure and determines their spin-rotation parts by solving a Procrustes problem. An implementation is distributed under a permissive free software license in spinspg v0.1.1: https://github.com/spglib/spinspg. △ Less

Submitted 21 October, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Growth of self-integrated atomic quantum wires and junctions of a Mott semiconductor

Authors: Tomoya Asaba, Lang Peng, Takahiro Ono, Satoru Akutagawa, Ibuki Tanaka, Hinako Murayama, Shota Suetsugu, Aleksandar Razpopov, Yuichi Kasahara, Takahito Terashima, Yuhki Kohsaka, Takasada Shibauchi, Masatoshi Ichikawa, Roser Valentí, Shin-ichi Sasa, Yuji Matsuda

Abstract: Continued advances in quantum technologies rely on producing nanometer-scale wires. Although several state-of-the-art nanolithographic technologies and bottom-up synthesis processes have been used to engineer such wires, critical challenges remain in growing uniform atomic-scale crystalline wires and constructing their network structures. Here we discover a simple method to fabricate atomic-scale… ▽ More Continued advances in quantum technologies rely on producing nanometer-scale wires. Although several state-of-the-art nanolithographic technologies and bottom-up synthesis processes have been used to engineer such wires, critical challenges remain in growing uniform atomic-scale crystalline wires and constructing their network structures. Here we discover a simple method to fabricate atomic-scale wires with various arrangements, including stripes, X-, Y-junctions, and nanorings. Single-crystalline atomic-scale wires of a Mott insulator, whose band gap is comparable to those of wide-gap semiconductors, are spontaneously grown on graphite substrates \DEL{and epitaxial monolayer graphene on SiC }by pulsed-laser deposition. These wires are one-unit-cell-thick and have an exact width of two- and four-unit-cells (1.4 and 2.8\,nm) and lengths up to a few $μみゅーm$. We show that the non-equilibrium reaction-diffusion processes may play an essential role in atomic pattern formation. Our findings offer a new perspective on the non-equilibrium self-organization phenomena on an atomic scale, paving a unique way for the quantum architecture of nano-network. △ Less

Submitted 22 May, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Journal ref: Sci. Adv.9,eabq5561(2023)

Effects of equivalent composition on superconducting properties of high-entropy REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, Gd) single crystals

Authors: Yuma Fujita, Masanori Nagao, Akira Miura, Daisuke Urushihara, Yoshikazu Mizuguchi, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Superconductors are influenced by high-entropy alloys (HEAs); these have been investigated in various functional materials. REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, and Gd in different combinations) single crystals with HEAs at the RE-site were successfully grown using the flux method. The obtained crystals were plate-shaped (1 mm$^2$) with a well-developed c-plane. Ce was present in both trivalent (C… ▽ More Superconductors are influenced by high-entropy alloys (HEAs); these have been investigated in various functional materials. REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, and Gd in different combinations) single crystals with HEAs at the RE-site were successfully grown using the flux method. The obtained crystals were plate-shaped (1 mm$^2$) with a well-developed c-plane. Ce was present in both trivalent (Ce$^{3+}$) and tetravalent (Ce$^{4+}$) electronic configurations; the concentration of Ce$^{4+}$ at the RE-site was approximately 10 at% in all single crystals. The single crystals showed superconducting transition temperature with zero resistivity within 1.2-4.2 K. The superconducting transition temperature, superconducting anisotropy, electronic specific heat coefficient, and Debye temperature of the crystals were not correlated with the mixed entropy at the RE-site. Except for the electronic specific heat coefficient, the variation of these parameters as a function of mixed entropy showed different trends for equivalent and non-equivalent RE element compositions. Thus, the configuration of RE elements influences the superconducting properties of REOBiS$_2$ single crystals, alluding to a method of modulating transition temperatures. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Journal ref: Physica C, vol.608 (2023) 1354254

Growth of REBa2Cu3Ox single-crystal whiskers utilizing the concept of high-entropy alloys

Authors: Yasunori Suzuki, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Single-crystal whiskers of REBa2Cu3Ox (RE-123, RE: rare earth element) cuprate superconductors were successfully grown using Te-doped precursors utilizing the concept of high-entropy-alloys (HEAs) at the RE site. The obtained whiskers were 0.5-1.0 mm length with approximately 10 um thickness. The width of the flat surface, which corresponds to the ab-plane, was 10-50 um. The RE site of the grown w… ▽ More Single-crystal whiskers of REBa2Cu3Ox (RE-123, RE: rare earth element) cuprate superconductors were successfully grown using Te-doped precursors utilizing the concept of high-entropy-alloys (HEAs) at the RE site. The obtained whiskers were 0.5-1.0 mm length with approximately 10 um thickness. The width of the flat surface, which corresponds to the ab-plane, was 10-50 um. The RE site of the grown whiskers could be easily substituted by large ionic radius RE elements (Gd, Dy, and Ho). Substitution with small ionic radius RE elements (Tm, Yb, and Lu) resulted in the opposite trend. The superconducting transition temperature and superconducting anisotropy of the grown RE-123 whiskers were 89-93 K and 6-10, respectively. The mixed entropy at the RE site in the RE-123 whiskers did not affect the superconducting transition temperature and superconducting anisotropy. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 27 pages, 8 figures, 2 tables

Journal ref: Jpn. J. Appl. Phys. 62 (2023) 033001

Implementation strategies in phonopy and phono3py

Authors: Atsushi Togo, Laurent Chaput, Terumasa Tadano, Isao Tanaka

Abstract: Scientific simulation codes are public property sustained by the community. Modern technology allows anyone to join scientific software projects, from anywhere, remotely via the internet. The phonopy and phono3py codes are widely used open source phonon calculation codes. This review describes a collection of computational methods and techniques as implemented in these codes and shows their impl… ▽ More Scientific simulation codes are public property sustained by the community. Modern technology allows anyone to join scientific software projects, from anywhere, remotely via the internet. The phonopy and phono3py codes are widely used open source phonon calculation codes. This review describes a collection of computational methods and techniques as implemented in these codes and shows their implementation strategies as a whole, aiming to be useful for the community. Some of the techniques presented here are not limited to phonon calculations and may therefore be useful in other area of condensed matter physics. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Algorithms for magnetic symmetry operation search and identification of magnetic space group from magnetic crystal structure

Authors: Kohei Shinohara, Atsushi Togo, Isao Tanaka

Abstract: A crystal symmetry search is crucial for computational crystallography and materials science. Although algorithms and implementations for the crystal symmetry search have been developed, their extension to magnetic space groups (MSGs) remains limited. In this paper, algorithms for determining magnetic symmetry operations of magnetic crystal structures, identifying magnetic space-group types of giv… ▽ More A crystal symmetry search is crucial for computational crystallography and materials science. Although algorithms and implementations for the crystal symmetry search have been developed, their extension to magnetic space groups (MSGs) remains limited. In this paper, algorithms for determining magnetic symmetry operations of magnetic crystal structures, identifying magnetic space-group types of given MSGs, searching for transformations to a BNS setting, and symmetrizing the magnetic crystal structures using the MSGs are presented. The determination of magnetic symmetry operations is numerically stable and is implemented with minimal modifications from the existing crystal symmetry search. Magnetic space-group types and transformations to the BNS setting are identified by a two-step approach combining space-group type identification and the use of affine normalizers. Point coordinates and magnetic moments of the magnetic crystal structures are symmetrized by projection operators for the MSGs. An implementation is distributed with a permissive free software license in spglib v2.0.2: https://github.com/spglib/spglib. △ Less

Submitted 7 June, 2023; v1 submitted 27 November, 2022; originally announced November 2022.

A Monomer Mott Insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$ as a Potential Nodal Line System

Authors: N. Yoneyama, M. K. Nuryadin, T. Tsumuraya, S. Iguchi, T. Takei, N. Kumada, M. Nagao, I. Tanaka, T. Sasaki

Abstract: We report the band structure calculations and the experimental results of resistivity and magnetic susceptibility in a spin-1/2 (BEDT-TTF)$^{\bullet +}$ monomer Mott insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$. The band calculations indicate a Dirac semimetal state with nodal lines at the Fermi level. The resistivity and the magnetic susceptibility as functions of temperature are well interpreted in ter… ▽ More We report the band structure calculations and the experimental results of resistivity and magnetic susceptibility in a spin-1/2 (BEDT-TTF)$^{\bullet +}$ monomer Mott insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$. The band calculations indicate a Dirac semimetal state with nodal lines at the Fermi level. The resistivity and the magnetic susceptibility as functions of temperature are well interpreted in terms of the monomer Mott insulating state instead of the expected semimetal state probably owing to strong electron correlation. In addition, we find that an Arrhenius-type steep reduction of the paramagnetic susceptibility appears below approximately 25 K, which indicates a spin-singlet ground state. △ Less

Submitted 28 September, 2022; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 4 pages, 5 figures. J. Phys. Soc. Jpn (2022) in press

Journal ref: J. Phys. Soc. Jpn. vol. 91 113704 (2022)

Investigation of superconductivity in Ce-doped (La,Pr)OBiS2 single crystals

Authors: Masanori Nagao, Yuji Hanada, Akira Miura, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Single crystals of Ce-doped (La,Pr)OBiS2 superconductors, multinary rare-earth elements substituted ROBiS2, were successfully grown. The grown crystals typically had a size of 1-2 mm and a plate-like shape with a well-developed c-plane. The c-axis lattice constants of the obtained (La,Ce,Pr)OBiS2 single crystals were approximately 13.6-13.7 A, and the superconducting transition temperature was 1.2… ▽ More Single crystals of Ce-doped (La,Pr)OBiS2 superconductors, multinary rare-earth elements substituted ROBiS2, were successfully grown. The grown crystals typically had a size of 1-2 mm and a plate-like shape with a well-developed c-plane. The c-axis lattice constants of the obtained (La,Ce,Pr)OBiS2 single crystals were approximately 13.6-13.7 A, and the superconducting transition temperature was 1.23-2.18 K. Valence fluctuations of Ce and Pr were detected through X-ray absorption spectroscopy analysis. In contrast to (Ce,Pr)OBiS2 and (La,Ce)OBiS2, the superconducting transition temperature of (La,Ce,Pr)OBiS2 increased with increasing concentrations of the tetravalent state at the R-site. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Journal ref: Materials, vol.15 (2022) 2977

LO-mode phonon of KCl and NaCl at 300 K by inelastic X ray scattering measurements and first principles calculations

Authors: Atsushi Togo, Hiroyuki Hayashi, Terumasa Tadano, Satoshi Tsutsui, Isao Tanaka

Abstract: Longitudinal-optical (LO) mode phonon branches of KCl and NaCl were measured using inelastic X-ray scattering (IXS) at 300 K and calculated by the first-principles phonon calculation with the stochastic self-consistent harmonic approximation. Spectral shapes of the IXS measurements and calculated spectral functions agreed well. We analyzed the calculated spectral functions that provide higher reso… ▽ More Longitudinal-optical (LO) mode phonon branches of KCl and NaCl were measured using inelastic X-ray scattering (IXS) at 300 K and calculated by the first-principles phonon calculation with the stochastic self-consistent harmonic approximation. Spectral shapes of the IXS measurements and calculated spectral functions agreed well. We analyzed the calculated spectral functions that provide higher resolutions of the spectra than the IXS measurements. Due to strong anharmonicity, the spectral functions of these phonon branches have several peaks and the LO modes along $Γがんま$--L paths are disconnected. △ Less

Submitted 22 April, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Intrinsic suppression of topological thermal Hall effect in an exactly solvable quantum magnet

Authors: S. Suetsugu, T. Yokoi, K. Totsuka, T. Ono, I. Tanaka, S. Kasahara, Y. Kasahara, Z. Chengchao, H. Kageyama, Y. Matsuda

Abstract: In contrast to electron (fermion) systems, topological phases of charge neutral bosons have been poorly understood despite recent extensive research on insulating magnets. The most important unresolved issue is how the inevitable inter-bosonic interactions influence the topological properties. It has been proposed that the quantum magnet SrCu$_2$(BO$_3$)$_2$ with an exact ground state serves as an… ▽ More In contrast to electron (fermion) systems, topological phases of charge neutral bosons have been poorly understood despite recent extensive research on insulating magnets. The most important unresolved issue is how the inevitable inter-bosonic interactions influence the topological properties. It has been proposed that the quantum magnet SrCu$_2$(BO$_3$)$_2$ with an exact ground state serves as an ideal platform for this investigation, as the system is expected to be a magnetic analogue of a Chern insulator with electrons replaced by bosonic magnetic excitations (triplons). Here, in order to examine topologically protected triplon chiral edge modes in SrCu$_2$(BO$_3$)$_2$, we measured and calculated the thermal Hall conductivity $κかっぱ_{xy}$. Our calculations show that the sign of $κかっぱ_{xy}$ is negative, which is opposite to the previous calculations, and its magnitude is 2$πぱい$ times smaller. No discernible $κかっぱ_{xy}$ was observed, which is at most 20-30% of our calculations if present. This implies that even relatively weak inter-particle interactions seriously influence the topological transport properties at finite temperatures. These demonstrate that, in contrast to fermionic cases, the picture of non-interacting topological quasi-particles cannot be naively applied to bosonic systems, calling special attention to the interpretation of the topological bosonic excitations reported for various insulating magnets. △ Less

Submitted 11 January, 2022; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 11 pages, 11 figures

Journal ref: Phys. Rev. B 105, 024415 (2022)

Finding well-optimized special quasirandom structures with decision diagram

Authors: Kohei Shinohara, Atsuto Seko, Takashi Horiyama, Isao Tanaka

Abstract: The advanced data structure of the zero-suppressed binary decision diagram (ZDD) enables us to efficiently enumerate nonequivalent substitutional structures. Not only can the ZDD store a vast number of structures in a compressed manner, but also can a set of structures satisfying given constraints be extracted from the ZDD efficiently. Here, we present a ZDD-based efficient algorithm for exhaustiv… ▽ More The advanced data structure of the zero-suppressed binary decision diagram (ZDD) enables us to efficiently enumerate nonequivalent substitutional structures. Not only can the ZDD store a vast number of structures in a compressed manner, but also can a set of structures satisfying given constraints be extracted from the ZDD efficiently. Here, we present a ZDD-based efficient algorithm for exhaustively searching for special quasirandom structures (SQSs) that mimic the perfectly random substitutional structure. We demonstrate that the current approach can extract only a tiny number of SQSs from a ZDD composed of many substitutional structures (>$10^{12}$). As a result, we find SQSs that are optimized better than those proposed in the literature. A series of SQSs should be helpful for estimating the properties of substitutional solid solutions. Furthermore, the present ZDD-based algorithm should be useful for applying the ZDD to the other structure enumeration problems. △ Less

Submitted 29 November, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Fluorine solubility and superconducting properties of Sm(O,F)BiS$_2$ single crystals

Authors: Koki Kinami, Yuji Hanada, Masanori Nagao, Akira Miura, Shigeto Hirai, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sm(O,F)BiS2 superconducting single crystals with various nominal F contents have been grown using KI-KCl flux. The solid solution limit of F at the O-site in Sm(O,F)BiS2 single crystals was approximately less than 20 at%. F concentrations of Sm(O,F)BiS2 single crystals were hardly controlled by nominal F contents. It suggests the existence of intrinsic stable phase. Superconductivity was shown at… ▽ More Sm(O,F)BiS2 superconducting single crystals with various nominal F contents have been grown using KI-KCl flux. The solid solution limit of F at the O-site in Sm(O,F)BiS2 single crystals was approximately less than 20 at%. F concentrations of Sm(O,F)BiS2 single crystals were hardly controlled by nominal F contents. It suggests the existence of intrinsic stable phase. Superconductivity was shown at around 5 K. Superconducting transition temperature with zero resistivity were increased with increasing nominal F contents within less than 70 at% in O-site, despite similar F concentration and c-axis lattice constants. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Journal ref: Journal of Alloys and Compounds vol.883 (2021) 160812

Nematicity Liquid in a Trimerized-Kagome Antiferromagnet

Authors: Ichiro Tanaka, Hirokazu Tsunetsugu

Abstract: We theoretically study low-temperature properties of the antiferromagnetic Heisenberg model with half-integer spin $S$ on the Kagome lattice with large trimerization. We have derived a low-energy effective model for general $S$ in terms of spin and nematicity operators in triangular units, and studied their low-temperature correlations for S=3/2 by classical Monte Carlo simulations. The previous s… ▽ More We theoretically study low-temperature properties of the antiferromagnetic Heisenberg model with half-integer spin $S$ on the Kagome lattice with large trimerization. We have derived a low-energy effective model for general $S$ in terms of spin and nematicity operators in triangular units, and studied their low-temperature correlations for S=3/2 by classical Monte Carlo simulations. The previous study for the S=1/2 case [M. Ferrero et al., Phys. Rev. B 68, 214431 (2003)] reported a spin liquid state at low temperatures driven by a glassy behavior of isolated dimers and trimers of nematicities. The results for S=3/2 show that nematicity dimers and trimers are connected by weak links to form a defective planar network. At very low temperatures, nematicities show glassy slow dynamics, and cluster spins continue to fluctuate in a nonperiodically frozen background of nematicities. The characteristic time of glassy dynamics scales with temperature following a power law. △ Less

Submitted 20 May, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 5 pages, 6 figures

NaInX2 (X = S, Se) layered materials for energy harvesting applications: First-principles insights into optoelectronic and thermoelectric properties

Authors: M. M. Hossain, M. A. Hossain, S. A. Moon, M. A. Ali, M. M. Uddin, S. H. Naqib, A. K. M. A. Islam M. Nagao S. Watauchi, I. Tanaka

Abstract: In the present study, the structural, electronic, optical and thermoelectric properties of two isostructural chalcogenide materials, NaInS2 and NaInSe2 with hexagonal symmetry (R-3m) have been studied using the first principles method. A very good agreement has been found between our results with the available experimental and theoretical ones. The studied materials are semiconducting in nature as… ▽ More In the present study, the structural, electronic, optical and thermoelectric properties of two isostructural chalcogenide materials, NaInS2 and NaInSe2 with hexagonal symmetry (R-3m) have been studied using the first principles method. A very good agreement has been found between our results with the available experimental and theoretical ones. The studied materials are semiconducting in nature as confirmed from the electronic band structure and optical properties.The strong hybridizations among s orbitals of Na, In and Se atoms push the bottom of the conduction band downward resulting in a narrower band gap of NaInSe2 compared to that of NaInS2 compound. Different optical (dielectric function, photoconductivity, absorption coefficient, reflectivity, refractive index and loss function) and thermoelectric (Seebeck coefficient, electrical conductivity, power factor and thermal conductivity) properties of NaInX2 (X = S, Se) have been studied in detail for the first time. It is found that all these properties are significantly anisotropic due to the strongly layered structure of NaInX2 (X = S, Se). Strong optical absorption with sharp peaks is found in the far visible to mid ultraviolet (UV) regions while the reflectivity is low in the UV region for both the compounds. Such features indicate feasibility of applications in optoelectronic sector.The calculated thermoelectric power factors at 1000 K for NaInS2 and NaInSe2 along a-axis are found to be 151.5 micro Watt /cmK2 and 154 micro Watt/cmK2, respectively and the corresponding ZT values are ~0.70. The obtained thermal conductivity along a-axis for both compounds is high (~22 W/mK).This suggests that the reduction of such high thermal conductivity is important to achieve higher ZT values of the NaInX2(X = S, Se) compounds. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 23 pages, 8 figures, 4 Tables

Influence of Se doping in recently synthesized NaInS2-xSex solid solutions for potential thermo-mechanical applications studied via first-principles method

Authors: M. M. Hossain, M. A. Ali, M. M. Uddin, M. A. Hossain, M. Rasadujjaman, S. H. Naqib, M. Nagao, S. Watauchi, I. Tanaka

Abstract: In the present work, the structural and hitherto unexplored thermal and mechanical properties of NaInS2-xSex (x = 0, 0.5, 1.0, 1.5 and 2.0) compounds have been studied using the density functional theory. Besides, the elastic anisotropy indices and hardnesses of NaInS2-xSex have been investigated as Se content is varied. The mechanical stability of all the compounds under study has been confirmed.… ▽ More In the present work, the structural and hitherto unexplored thermal and mechanical properties of NaInS2-xSex (x = 0, 0.5, 1.0, 1.5 and 2.0) compounds have been studied using the density functional theory. Besides, the elastic anisotropy indices and hardnesses of NaInS2-xSex have been investigated as Se content is varied. The mechanical stability of all the compounds under study has been confirmed. The ratio of shear to bulk modulus (G/B) is low suggesting that the NaInS2-xSex (x = 0.5 and 1.5) compounds exhibit damage tolerant (ductility) properties while rest of the compositions are brittle in nature. The predicted hardness (H) values are also influenced with the Se content in the following order: H (NaInSSe) > H (NaInS2) > H (NaInSe2) > H (NaInS1.5Se0.5) > H (NaInS0.5Se1.5). All the anisotropic indices under study indicate that NaInS2-xSex compounds are anisotropic in nature. The Mulliken bond population analysis suggests that the degree of covalency of In-S/Se bonds decreases when S is substituted by Se. The origin of low Debye temperature (ΘしーたD) and low minimum thermal conductivity (Kmin) have been successfully explained by considering the mean atomic weight (M/n) and average bond strength of the compounds. Temperature dependence of heat capacities (Cv, Cp) and linear thermal expansion coefficient (αあるふぁ ) are also estimated using the quasi-harmonic Debye model and discussed. The low values of Kmin, ΘしーたD and αあるふぁ and damage tolerant behavior clearly indicate that the NaInS2-xSex (x = 0.5 and 1.5) compounds can be used as promising thermal barrier coating materials for high temperature applications. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 23 pages, 11 figures, 7 Tables

Growth and anisotropy evaluation of NbBiCh$_3$ (Ch = S, Se) misfit-layered superconducting single crystals

Authors: Masanori Nagao, Akira Miura, Yoichi Horibe, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: NbBiCh$_3$ (Ch = S, Se) misfit-layered superconducting single crystals were successfully grown using a CsCl/KCl flux for the first time. The obtained crystals had a well-developed habit parallel to the c-plane with a typical width of 1-2 mm and thickness of 10-40 um. The superconducting transition temperatures with zero resistivity of NbBiS$_3$ single crystals obtained from the nominal composition… ▽ More NbBiCh$_3$ (Ch = S, Se) misfit-layered superconducting single crystals were successfully grown using a CsCl/KCl flux for the first time. The obtained crystals had a well-developed habit parallel to the c-plane with a typical width of 1-2 mm and thickness of 10-40 um. The superconducting transition temperatures with zero resistivity of NbBiS$_3$ single crystals obtained from the nominal composition of Nb0.9Bi1.2S3 was 0.31 K, and that value of the NbBiSe$_3$ single crystals grown from the stoichiometry composition (NbBiSe$_3$) was 2.3 K. Sharp decreases in electric resistivity and magnetic susceptibility at approximately 3 K suggested a possible superconducting transition temperature of NbBiSe$_3$. The normal-state anisotropy values of grown NbBiS$_3$ and NbBiSe$_3$ single crystals were 2.2-2.4 and 1.5-1.6, respectively. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: 4 figures + 3 figures (Supplementary Materials)

Journal ref: Solid State Communications 321 (2020) 114051

Application of machine learning potentials to predict grain boundary properties in fcc elemental metals

Authors: Takayuki Nishiyama, Atsuto Seko, Isao Tanaka

Abstract: Accurate interatomic potentials are in high demand for large-scale atomistic simulations of materials that are prohibitively expensive by density functional theory (DFT) calculation. In this study, we apply machine learning potentials in a recently constructed repository to the prediction of the grain boundary energy in face-centered-cubic elemental metals, i.e., Ag, Al, Au, Cu, Pd, and Pt. The sy… ▽ More Accurate interatomic potentials are in high demand for large-scale atomistic simulations of materials that are prohibitively expensive by density functional theory (DFT) calculation. In this study, we apply machine learning potentials in a recently constructed repository to the prediction of the grain boundary energy in face-centered-cubic elemental metals, i.e., Ag, Al, Au, Cu, Pd, and Pt. The systematic application of machine learning potentials shows that they enable us to predict grain boundary structures and their energies accurately. The grain boundary energies predicted by the MLPs are in agreement with those calculated by DFT, although no grain boundary structures were included in training datasets of the present MLPs. △ Less

Submitted 13 August, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. Materials 4, 123607 (2020)

Phonon structure of titanium under shear deformation along $\{10\bar{1}2\}$ twinning mode

Authors: Atsushi Togo, Yuta Inoue, Isao Tanaka

Abstract: We investigated phonon behavior of hexagonal close packed titanium under homogeneous shear deformation corresponding to the $\{10\bar{1}2\}$ twinning mode using first-principles calculation and phonon calculation. By this deformation, we found that a phonon mode located at a point on Brillouin zone boundary is drastically soften increasing the shear and finally it triggers a spontaneous structural… ▽ More We investigated phonon behavior of hexagonal close packed titanium under homogeneous shear deformation corresponding to the $\{10\bar{1}2\}$ twinning mode using first-principles calculation and phonon calculation. By this deformation, we found that a phonon mode located at a point on Brillouin zone boundary is drastically soften increasing the shear and finally it triggers a spontaneous structural transition by breaking the crystal symmetry toward twin from parent. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Journal ref: Phys. Rev. B 102, 024106 (2020)

Derivative structure enumeration using binary decision diagram

Authors: Kohei Shinohara, Atsuto Seko, Takashi Horiyama, Masakazu Ishihata, Junya Honda, Isao Tanaka

Abstract: A derivative structure is a nonequivalent substitutional atomic configuration derived from a given primitive cell. The enumeration of derivative structures plays an essential role in searching for the ground states in multicomponent systems. However, it is computationally hard to enumerate derivative structures if the number of derivative structures of a target system becomes huge. In the present… ▽ More A derivative structure is a nonequivalent substitutional atomic configuration derived from a given primitive cell. The enumeration of derivative structures plays an essential role in searching for the ground states in multicomponent systems. However, it is computationally hard to enumerate derivative structures if the number of derivative structures of a target system becomes huge. In the present study, we introduce a novel compact data structure of the zero-suppressed binary decision diagram (ZDD) to enumerate derivative structures much more efficiently. The present study shows its simple applications to the enumeration of structures derived from the face-centered cubic and hexagonal close-packed lattices in binary, ternary, and quaternary systems. The present ZDD-based procedure should significantly contribute not only to various computational approaches based on the derivative structures but also to a wide range of combinatorial issues in physics and materials science. △ Less

Submitted 13 July, 2020; v1 submitted 28 February, 2020; originally announced February 2020.

Growth of superconducting Sm(O,F)BiS2 single crystals

Authors: Koki Kinami, Yuji Hanada, Masanori Nagao, Akira Miura, Yosuke Goto, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Superconducting Sm(O,F)BiS2 single crystals were successfully grown using KI-KCl flux. The obtained single crystals were flat-shaped and their size was 500um. The structure and composition of the obtained single crystals were obtained using X-ray diffraction and X-ray spectroscopy. The lattice constants of the obtained Sm(O,F)BiS2 single crystals were evaluated to be a = 3.9672(4) A and c = 13.417… ▽ More Superconducting Sm(O,F)BiS2 single crystals were successfully grown using KI-KCl flux. The obtained single crystals were flat-shaped and their size was 500um. The structure and composition of the obtained single crystals were obtained using X-ray diffraction and X-ray spectroscopy. The lattice constants of the obtained Sm(O,F)BiS2 single crystals were evaluated to be a = 3.9672(4) A and c = 13.417(6) A, while their F/(O + F) molar ratio was approximately 0.17. The superconducting transition temperature of the obtained Sm(O,F)BiS2 single crystals was 4.8 K. This was different from the reported Sm(O,F)BiS2 with no superconducting transition above 2 K. The superconducting anisotropies of the upper critical field and effective mass model of the obtained Sm(O,F)BiS2 single crystals were estimated to be approximately 26 and 10, respectively. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: 22 pages, 7 figures

Journal ref: Cryst.Growth Des. vol.19, pp.6136-6140(2019)

Growth and characterization of (La,Ce)OBiS2 single crystals

Authors: Yuji Hanada, Masanori Nagao, Akira Miura, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: (La,Ce)OBiS2 single crystals are successfully grown using CsCl flux. The obtained single crystals have a plate-like shape with a size of 0.5-2.0 mm and a well-developed ab-plane. The thickness of the crystals increases with a decrease in the Ce content. The (La,Ce)OBiS2 single crystals exhibit superconducting properties. Moreover, the superconducting transition temperature increases with an increa… ▽ More (La,Ce)OBiS2 single crystals are successfully grown using CsCl flux. The obtained single crystals have a plate-like shape with a size of 0.5-2.0 mm and a well-developed ab-plane. The thickness of the crystals increases with a decrease in the Ce content. The (La,Ce)OBiS2 single crystals exhibit superconducting properties. Moreover, the superconducting transition temperature increases with an increase in the Ce content. The obtained (La,Ce)OBiS2 single crystals exhibit a tetragonal structure and the valence state of Ce in the crystals is determined to be mixed with Ce3+ and Ce4+ based on X-ray absorption fine structure spectroscopy. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 12 pages, 6 figures

Journal ref: Japanese Journal of Applied Physics 58, 063001 (2019)

Vibrational Effects in X-ray Absorption Spectra of 2D Layered Materials

Authors: Weine Olovsson, Teruyasu Mizoguchi, Martin Magnuson, Stefan Kontur, Olle Hellman, Isao Tanaka, Claudia Draxl

Abstract: With the examples of the C $K$-edge in graphite and the B $K$-edge in hexagonal BN, we demonstrate the impact of vibrational coupling and lattice distortions on the X-ray absorption near-edge structure (XANES) in 2D layered materials. Theoretical XANES spectra are obtained by solving the Bethe-Salpeter equation of many-body perturbation theory, including excitonic effects through the correlated mo… ▽ More With the examples of the C $K$-edge in graphite and the B $K$-edge in hexagonal BN, we demonstrate the impact of vibrational coupling and lattice distortions on the X-ray absorption near-edge structure (XANES) in 2D layered materials. Theoretical XANES spectra are obtained by solving the Bethe-Salpeter equation of many-body perturbation theory, including excitonic effects through the correlated motion of core-hole and excited electron. We show that accounting for zero-point motion is important for the interpretation and understanding of the measured X-ray absorption fine structure in both materials, in particular for describing the $σしぐま^*$-peak structure. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: & pages, 4 figures

Journal ref: J. Phys. Chem. C 2019, 123, 15, 9688-9692

Growth and transport properties under high pressure of PrOBiS2 single crystals

Authors: Masanori Nagao, Akira Miura, Ryo Matsumoto, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Kiyoharu Tadanaga, Isao Tanaka

Abstract: PrOBiS2 single crystals were successfully grown using a KCl flux. The obtained crystals had a plate-like shape with a typical size of 0.5-1.0 mm and well-developed ab-plane. The Pr valence of the grown crystals, as determined by an X-ray absorption fine structure spectroscopy analysis, was almost trivalent. The PrOBiS2 single crystals did not exhibit superconductivity down to 0.25 K. The transport… ▽ More PrOBiS2 single crystals were successfully grown using a KCl flux. The obtained crystals had a plate-like shape with a typical size of 0.5-1.0 mm and well-developed ab-plane. The Pr valence of the grown crystals, as determined by an X-ray absorption fine structure spectroscopy analysis, was almost trivalent. The PrOBiS2 single crystals did not exhibit superconductivity down to 0.25 K. The transport properties under high pressures up to 50 GPa ensured semiconducting behavior for 2-300 K using a diamond anvil cell. △ Less

Submitted 20 April, 2019; originally announced April 2019.

Comments: 19 pages, 5 figures

Journal ref: Solid State Communications 296 (2019) 17-20

Group-theoretical high-order rotational invariants for structural representations: Application to linearized machine learning interatomic potential

Authors: Atsuto Seko, Atsushi Togo, Isao Tanaka

Abstract: Many rotational invariants for crystal structure representations have been used to describe the structure-property relationship by machine learning. The machine learning interatomic potential (MLIP) is one of the applications of rotational invariants, which provides the relationship between the energy and the crystal structure. Since the MLIP requires the highest accuracy among machine learning es… ▽ More Many rotational invariants for crystal structure representations have been used to describe the structure-property relationship by machine learning. The machine learning interatomic potential (MLIP) is one of the applications of rotational invariants, which provides the relationship between the energy and the crystal structure. Since the MLIP requires the highest accuracy among machine learning estimations of the structure-property relationship, the enumeration of rotational invariants is useful for constructing MLIPs with the desired accuracy. In this study, we introduce high-order linearly independent rotational invariants up to the sixth order based on spherical harmonics and apply them to linearized MLIPs for elemental aluminum. A set of rotational invariants is derived by the general process of reducing the Kronecker products of irreducible representations (Irreps) for the SO(3) group using a group-theoretical projector method. A high predictive power for a wide range of structures is accomplished by using high-order invariants with low-order invariants equivalent to pair and angular structural features. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Comments: 15 pages, 4 figures

Journal ref: Phys. Rev. B 99, 214108 (2019)

Growth and physical properties of Ce(O,F)Sb(S,Se)2 single crystals with site-selected chalcogen atoms

Authors: Masanori Nagao, Masashi Tanaka, Akira Miura, Miho Kitamura, Koji Horiba, Satoshi Watauchi, Yoshihiko Takano, Hiroshi Kumigashira, Isao Tanaka

Abstract: Ce(O,F)Sb(S,Se)2 single crystals were successfully grown using a CsCl/KCl flux method. The obtained crystals have a plate-like shape with the typical size of 1-2 mm and well-developed ab-plane, which enables X-ray single crystal structural analysis. The Ce(O,F)Sb(S,Se)2 crystallizes in a monoclinic space group, P21/m, with lattice parameters of a = 4.121(7) A, b = 4.109(7) A, c = 13.233(15) A, bet… ▽ More Ce(O,F)Sb(S,Se)2 single crystals were successfully grown using a CsCl/KCl flux method. The obtained crystals have a plate-like shape with the typical size of 1-2 mm and well-developed ab-plane, which enables X-ray single crystal structural analysis. The Ce(O,F)Sb(S,Se)2 crystallizes in a monoclinic space group, P21/m, with lattice parameters of a = 4.121(7) A, b = 4.109(7) A, c = 13.233(15) A, beta = 97.94(7) deg. It is composed of alternate stacking of Ce-(O,F) and Sb-SSe layers, and the Sb-SSe layer includes selective occupation of Se atoms in its in-plane site. The valence state of Ce is estimated to be Ce3+ by X-ray absorption fine spectroscopy analysis. The single crystals show an insulating behavior, and a magnetic ordering around 6 K. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 27 pages, 7 figures, 3 Tables

Journal ref: Solid State Communications vol.289(2019)38-42

Pressure-induced Insulator to Metal Transition of Mixed Valence Compound Ce(O,F)SbS$_{2}$

Authors: Ryo Matsumoto, Masanori Nagao, Masayuki Ochi, Hiromi Tanaka, Hiroshi Hara, Shintaro Adachi, Kazuki Nakamura, Ryo Murakami, Sayaka Yamamoto, Tetsuo Irifune, Hiroyuki Takeya, Isao Tanaka, Kazuhiko Kuroki, Yoshihiko Takano

Abstract: Transport properties of Ce$_{0.85}$F0.15SbS$_{2}$ and undoped CeOSbS$_{2}$ under high pressure were investigated experimentally and theoretically. Electrical resistivity measurements of the Ce$_{0.85}$F0.15SbS$_{2}$ single crystals were performed under various high pressures using a diamond anvil cell with boron-doped diamond electrodes. The samples showed the insulator to metal transition by appl… ▽ More Transport properties of Ce$_{0.85}$F0.15SbS$_{2}$ and undoped CeOSbS$_{2}$ under high pressure were investigated experimentally and theoretically. Electrical resistivity measurements of the Ce$_{0.85}$F0.15SbS$_{2}$ single crystals were performed under various high pressures using a diamond anvil cell with boron-doped diamond electrodes. The samples showed the insulator to metal transition by applying high pressure up to 30-40 GPa. On the other hand, the undoped CeOSbS$_{2}$ showed almost same transport property with the F-doped sample under high pressure. The valence state analysis using X-ray photoelectron spectroscopy revealed a simple valence state of Ce3+ in Ce$_{0.85}$F0.15SbS$_{2}$ and mixed valence state between Ce3+ and Ce4+ in undoped CeOSbS$_{2}$. The valence fluctuation in Ce carried out the comparable transport nature in the both samples. A band calculation suggests that the undoped CeOSbS$_{2}$ could be metallic under high pressure of 30 GPa in accordance with the experimental results. A superior thermoelectric property of power factor in CeOSbS$_{2}$ was estimated under high pressure around 20 GPa in comparison with that of ambient pressure. △ Less

Submitted 2 November, 2018; originally announced November 2018.

Journal ref: J.Appl. Phys.125,075102, (2019)

Bulk superconductivity in La2O2M4S6-type layered oxychalcogenide La2O2Bi3Ag0.6Sn0.4S5.7Se0.3

Authors: Rajveer Jha, Yosuke Goto, Tatsuma D. Matsuda, Yuji Aoki, Masanori Nagao, Isao Tanaka, Yoshikazu Mizuguchi

Abstract: Recently, we reported the observation of superconductivity at ~0.5 K in a La2O2M4S6-type (M: metal) layered oxychalcogenide La2O2Bi3AgS6, which is a layered compound related to the BiS2-based superconductor system but possesses a thicker Bi3AgS6-type conducting layer. In this study, we have developed the La2O2Bi3AgS6-type materials by element substitutions to increase the transition temperature (T… ▽ More Recently, we reported the observation of superconductivity at ~0.5 K in a La2O2M4S6-type (M: metal) layered oxychalcogenide La2O2Bi3AgS6, which is a layered compound related to the BiS2-based superconductor system but possesses a thicker Bi3AgS6-type conducting layer. In this study, we have developed the La2O2Bi3AgS6-type materials by element substitutions to increase the transition temperature (Tc) and to induce bulk nature of superconductivity. A resistivity anomaly observed at 180 K in La2O2Bi3AgS6 was systematically suppressed by Sn substitution for the Ag site. By the Sn substitution, Tc increased, and the shielding volume fraction estimated from magnetization measurements also increased. The highest Tc (= 2.3 K) and the highest shielding volume fraction (~20%) was observed for La2O2Bi3Ag0.6Sn0.4S6. The superconducting properties were further improved by Se substitutions for the S site. By the combinational substitutions of Sn and Se, bulk-superconducting phase of La2O2Bi3Ag0.6Sn0.4S5.7Se0.3 with a Tc of 3.0 K (Tconset = 3.6 K) was obtained. △ Less

Submitted 15 April, 2019; v1 submitted 19 October, 2018; originally announced October 2018.

Comments: 18 Pages, 8 figures. The title has been changed from "Sn substitution effect on the superconductivity of layered oxychalcogenide La2O2Bi3Ag1-xSnxS6". Two authors have been added

Journal ref: Scientific Reports 9, 13346 (2019)

$\texttt{Spglib}$: a software library for crystal symmetry search

Authors: Atsushi Togo, Kohei Shinohara, Isao Tanaka

Abstract: A computer algorithm to search symmetries of crystal structures as implemented in the \texttt{spglib} code is described. An iterative algorithm is employed to robustly identify space group types tolerating a certain amount of distortion in the crystal structures. The source code is distributed under the 3-Clause BSD License, a permissive open-source software license. This paper focuses on the algo… ▽ More A computer algorithm to search symmetries of crystal structures as implemented in the \texttt{spglib} code is described. An iterative algorithm is employed to robustly identify space group types tolerating a certain amount of distortion in the crystal structures. The source code is distributed under the 3-Clause BSD License, a permissive open-source software license. This paper focuses on the algorithm for identifying the space group symmetry of the crystal structures. △ Less

Submitted 13 March, 2024; v1 submitted 5 August, 2018; originally announced August 2018.

Lattice thermal conductivities of two SiO$_2$ polymorphs by first-principles calculation and phonon Boltzmann transport equation

Authors: Keiyu Mizokami, Atsushi Togo, Isao Tanaka

Abstract: Lattice thermal conductivities of two SiO$_2$ polymorphs, i.e., $αあるふぁ$-quartz (low) and $αあるふぁ$-cristobalite (low), were studied using first-principles anharmonic phonon calculation and linearized phonon Boltzmann transport equation. Although $αあるふぁ$-quartz and $αあるふぁ$-cristobalite have similar phonon densities of states, phonon frequency dependencies of phonon group velocities and lifetimes are dissimilar, whic… ▽ More Lattice thermal conductivities of two SiO$_2$ polymorphs, i.e., $αあるふぁ$-quartz (low) and $αあるふぁ$-cristobalite (low), were studied using first-principles anharmonic phonon calculation and linearized phonon Boltzmann transport equation. Although $αあるふぁ$-quartz and $αあるふぁ$-cristobalite have similar phonon densities of states, phonon frequency dependencies of phonon group velocities and lifetimes are dissimilar, which results in largely different anisotropies of the lattice thermal conductivities. For $αあるふぁ$-quartz and $αあるふぁ$-cristobalite, distributions of the phonon lifetimes effective to determine the lattice thermal conductivities are well described by energy and momentum conservations of three phonon scatterings weighted by phonon occupation numbers and one parameter that represents the phonon-phonon interaction strengths. △ Less

Submitted 9 April, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Journal ref: Phys. Rev. B 97, 224306 (2018)

Coexistence of superconductivity and charge-density wave in the quasi-one-dimensional material HfTe3

Authors: Saleem J. Denholme, Akinori Yukawa, Kohei Tsumura, Masanori Nagao, Ryuji Tamura, Satoshi Watauchi, Isao Tanaka, Hideaki Takayanagi, Nobuaki Miyakawa

Abstract: We present the first experimental evidence for metallicity, superconductivity (SC) and the co-existence of charge density waves (CDW) in the quasi-one-dimensional material HfTe3. The existence of such phenomena is a typical characteristic of the transition metal chalcogenides however, without the application of hydrostatic pressure/chemical doping, it is rare for a material to exhibit the co-exist… ▽ More We present the first experimental evidence for metallicity, superconductivity (SC) and the co-existence of charge density waves (CDW) in the quasi-one-dimensional material HfTe3. The existence of such phenomena is a typical characteristic of the transition metal chalcogenides however, without the application of hydrostatic pressure/chemical doping, it is rare for a material to exhibit the co-existence of both states. Materials such as HfTe3 can therefore provide us with a unique insight into the relationship between these multiple ordered states. By improving on the original synthesis conditions, we have successfully synthesised single phase HfTe3 and confirmed the resultant structure by performing Rietveld refinement. Using low temperature resistivity measurements, we provide the first experimental evidence of SC at ~1.4 K as well as a resistive anomaly indicative of a CDW formation at ~82 K. By the application of hydrostatic-pressure, the resistivity anomaly shifts to higher temperature. The results show that HfTe3 is a promising new material to help study the relationship between SC and CDW. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 9 pages, 4 figures, includes supplementary information

Journal ref: Scientific Reports 7, Article number: 45217 (2017)

Compositional descriptor-based recommender system accelerating the materials discovery

Authors: Atsuto Seko, Hiroyuki Hayashi, Isao Tanaka

Abstract: Structures and properties of many inorganic compounds have been collected historically. However, it only covers a very small portion of possible inorganic crystals, which implies the presence of numerous currently unknown compounds. A powerful machine-learning strategy is mandatory to discover new inorganic compounds from all chemical combinations. Herein we propose a descriptor-based recommender-… ▽ More Structures and properties of many inorganic compounds have been collected historically. However, it only covers a very small portion of possible inorganic crystals, which implies the presence of numerous currently unknown compounds. A powerful machine-learning strategy is mandatory to discover new inorganic compounds from all chemical combinations. Herein we propose a descriptor-based recommender-system approach to estimate the relevance of chemical compositions where stable crystals can be formed [i.e., chemically relevant compositions (CRCs)]. As well as data-driven compositional similarity used in the literature, the use of compositional descriptors as a prior knowledge can accelerate the discovery of new compounds. We validate our recommender systems in two ways. Firstly, one database is used to construct a model, while another is used for the validation. Secondly, we estimate the phase stability for compounds at expected CRCs using density functional theory calculations. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 8 pages, 7 figures

Linearized machine-learning interatomic potentials for non-magnetic elemental metals: Limitation of pairwise descriptors and trend of predictive power

Authors: Akira Takahashi, Atsuto Seko, Isao Tanaka

Abstract: Machine-learning interatomic potential (MLIP) has been of growing interest as a useful method to describe the energetics of systems of interest. In the present study, we examine the accuracy of linearized pairwise MLIPs and angular-dependent MLIPs for 31 elemental metals. Using all of the optimal MLIPs for 31 elemental metals, we show the robustness of the linearized frameworks, the general trend… ▽ More Machine-learning interatomic potential (MLIP) has been of growing interest as a useful method to describe the energetics of systems of interest. In the present study, we examine the accuracy of linearized pairwise MLIPs and angular-dependent MLIPs for 31 elemental metals. Using all of the optimal MLIPs for 31 elemental metals, we show the robustness of the linearized frameworks, the general trend of the predictive power of MLIPs and the limitation of pairwise MLIPs. As a result, we obtain accurate MLIPs for all 31 elements using the same linearized framework. This indicates that the use of numerous descriptors is the most important practical feature for constructing MLIPs with high accuracy. An accurate MLIP can be constructed using only pairwise descriptors for most non-transition metals, whereas it is very important to consider angular-dependent descriptors when expressing interatomic interactions of transition metals. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 14 figures

Matrix- and tensor-based recommender systems for the discovery of currently unknown inorganic compounds

Authors: Atsuto Seko, Hiroyuki Hayashi, Hisashi Kashima, Isao Tanaka

Abstract: Chemically relevant compositions (CRCs) and atomic arrangements of inorganic compounds have been collected as inorganic crystal structure databases. Machine-learning is a unique approach to search for currently unknown CRCs from vast candidates. Herein we propose matrix- and tensor-based recommender system approaches to predict currently unknown CRCs from database entries of CRCs. Firstly, the per… ▽ More Chemically relevant compositions (CRCs) and atomic arrangements of inorganic compounds have been collected as inorganic crystal structure databases. Machine-learning is a unique approach to search for currently unknown CRCs from vast candidates. Herein we propose matrix- and tensor-based recommender system approaches to predict currently unknown CRCs from database entries of CRCs. Firstly, the performance of the recommender system approaches to discover currently unknown CRCs is examined. A Tucker decomposition recommender system shows the best discovery rate of CRCs as the majority of the top 100 recommended ternary and quaternary compositions correspond to CRCs. Secondly, systematic density functional theory (DFT) calculations are performed to investigate the phase stability of the recommended currently unknown CRCs. The phase stability of the 27 compositions reveals that 23 currently unknown compounds are stable. These results indicate that the recommender system has great potential to accelerate the discovery of new compounds. △ Less

Submitted 20 September, 2017; originally announced October 2017.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. Materials 2, 013805 (2018)

Descriptors for Machine Learning of Materials Data

Authors: Atsuto Seko, Atsushi Togo, Isao Tanaka

Abstract: Descriptors, which are representations of compounds, play an essential role in machine learning of materials data. Although many representations of elements and structures of compounds are known, these representations are difficult to use as descriptors in their unchanged forms. This chapter shows how compounds in a dataset can be represented as descriptors and applied to machine-learning models f… ▽ More Descriptors, which are representations of compounds, play an essential role in machine learning of materials data. Although many representations of elements and structures of compounds are known, these representations are difficult to use as descriptors in their unchanged forms. This chapter shows how compounds in a dataset can be represented as descriptors and applied to machine-learning models for materials datasets. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 12 pages, 10 figures

DynaPhoPy: A code for extracting phonon quasiparticles from molecular dynamics simulations

Authors: Abel Carreras, Atsushi Togo, Isao Tanaka

Abstract: We have developed a computational code, DynaPhoPy, that allow us to extract the microscopic anharmonic phonon properties from molecular dynamics (MD) simulations using the normal-mode-decomposition technique as presented by Sun et al. [T. Sun, D. Zhang, R. Wentzcovitch, 2014]. Using this code we calculated the quasiparticle phonon frequencies and linewidths of crystalline silicon at different temp… ▽ More We have developed a computational code, DynaPhoPy, that allow us to extract the microscopic anharmonic phonon properties from molecular dynamics (MD) simulations using the normal-mode-decomposition technique as presented by Sun et al. [T. Sun, D. Zhang, R. Wentzcovitch, 2014]. Using this code we calculated the quasiparticle phonon frequencies and linewidths of crystalline silicon at different temperatures using both of first-principles and the Tersoff empirical potential approaches. In this work we show the dependence of these properties on the temperature using both approaches and compare them with reported experimental data obtained by Raman spectroscopy [M. Balkanski, R. Wallis, E. Haro, 1983 and R. Tsu, J. G. Hernandez, 1982]. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: 14 pages, 12 figures

Conceptual and practical bases for the high accuracy of machine learning interatomic potential

Authors: Akira Takahashi, Atsuto Seko, Isao Tanaka

Abstract: Machine learning interatomic potentials (MLIPs) based on a large dataset obtained by density functional theory (DFT) calculation have been developed recently. This study gives both conceptual and practical bases for the high accuracy of MLIPs, although MLIPs have been considered to be simply an accurate black-box description of atomic energy. We also construct the most accurate MLIP of the element… ▽ More Machine learning interatomic potentials (MLIPs) based on a large dataset obtained by density functional theory (DFT) calculation have been developed recently. This study gives both conceptual and practical bases for the high accuracy of MLIPs, although MLIPs have been considered to be simply an accurate black-box description of atomic energy. We also construct the most accurate MLIP of the elemental Ti ever reported using a linearized MLIP framework and many angular-dependent descriptors, which also corresponds to a generalization of the modified embedded atom method (MEAM) potential. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Journal ref: Phys. Rev. Materials 1, 063801 (2017)

Superconductivity and its Enhancement under High Pressure in "F-free" Single Crystals of CeOBiS$_2$

Authors: Masashi Tanaka, Masanori Nagao, Ryo Matsumoto, Noriyuki Kataoka, Ikuo Ueta, Hiromi Tanaka, Satoshi Watauchi, Isao Tanaka, Yoshihiko Takano

Abstract: "F-free" CeOBiS$_2$ single crystals have successfully grown, thoroughly eliminating a concern about F-contamination by using a high-purity CsCl flux. The obtained crystals have a plate-like shape with a size of 1.0-3.0 mm in the well-developed plane. Single crystal X-ray structural analysis clearly revealed that the CeOBiS$_2$ crystallizes with a space group $P4/nmm$ (with lattice parameters of… ▽ More "F-free" CeOBiS$_2$ single crystals have successfully grown, thoroughly eliminating a concern about F-contamination by using a high-purity CsCl flux. The obtained crystals have a plate-like shape with a size of 1.0-3.0 mm in the well-developed plane. Single crystal X-ray structural analysis clearly revealed that the CeOBiS$_2$ crystallizes with a space group $P4/nmm$ (with lattice parameters of $a$ = 4.0189(6) Å, $c$ = 13.573(2) Å). The bond valence sum estimation and X-ray photoelectron spectroscopy analysis showed that the chemical state of Ce was in the mixed valence of Ce$^{3+}$ and Ce$^{4+}$. The single crystals show superconductivity with zero resistivity at $\sim$1.3 K. This is the first conclusive evidence of superconductivity driven by Ce valence fluctuation in surely non-doped CeOBiS$_2$. The superconducting transition temperature was enhanced up to $\sim$3.8 K by applying hydrostatic pressure. △ Less

Submitted 12 June, 2017; originally announced June 2017.

Comments: 27 pages, 8 figures, J. Alloys and Compd. (2017)

Journal ref: J. Alloys and Compd. 722, 467-473 (2017)

Growth and superconducting properties of Cd-doped La(O,F)BiS2 single crystals

Authors: Masanori Nagao, Masashi Tanaka, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Cd-doped La(O,F)BiS2 single crystals were grown using a CsCl/KCl flux. The grown crystals have a plate-like shape with ~1-2 mm square size in a well-developed ab-plane. The Cd doping in the crystals was successfully characterized by single crystal X-ray diffraction and electron probe microanalysis techniques. The superconductivity of La(O,F)BiS2 was gradually suppressed with Cd doping. The superco… ▽ More Cd-doped La(O,F)BiS2 single crystals were grown using a CsCl/KCl flux. The grown crystals have a plate-like shape with ~1-2 mm square size in a well-developed ab-plane. The Cd doping in the crystals was successfully characterized by single crystal X-ray diffraction and electron probe microanalysis techniques. The superconductivity of La(O,F)BiS2 was gradually suppressed with Cd doping. The superconducting transition temperature with zero resistivity of La(O0.54F0.46)(Bi0.92Cd0.08)S2 was 2.3 K. The Cd doping does not change the superconducting anisotropy so much, albeit the considerable suppression of Tc. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: In Press, Solid State Communications 2017

Journal ref: Solid State Communications 261 (2017) 32-36

Direct evidence of hidden local spin polarization in a centrosymmetric superconductor LaO$_{0.55}$F$_{0.45}$BiS$_2$

Authors: Shi-Long Wu, Kazuki Sumida, Koji Miyamoto, Kazuaki Taguchi, Tomoki Yoshikawa, Akio Kimura, Yoshifumi Ueda, Masashi Arita, Masanori Nagao, Satoshi Watauchi, Isao Tanaka, Taichi Okuda

Abstract: Conventional Rashba spin polarization is caused by the combination of strong spin-orbit interaction (SOI) and spatial inversion asymmetry. However, Rashba- and Dresselhaus-type spin-split states are predicted in LaOBiS$_2$ system by recent theory even though the crystal structure is centrosymmetric, which stem from the local inversion asymmetry of active BiS$_2$ layer. By performing high-resolutio… ▽ More Conventional Rashba spin polarization is caused by the combination of strong spin-orbit interaction (SOI) and spatial inversion asymmetry. However, Rashba- and Dresselhaus-type spin-split states are predicted in LaOBiS$_2$ system by recent theory even though the crystal structure is centrosymmetric, which stem from the local inversion asymmetry of active BiS$_2$ layer. By performing high-resolution spin- and angle-resolved photoemission spectroscopy, we have investigated the electronic band structure and spin texture of superconductor LaO$_{0.55}$F$_{0.45}$BiS$_2$. Our studies present direct spectroscopic evidence for the local spin polarization in the vicinity of X point of both valence band and conduction band. Especially the coexistence of Rashba-like and Dresselhaus-like spin textures has been observed in the conduction band for the first time. The finding is of key importance for fabrication of proposed dual-gated spin-field effect transistor (SFET). Moreover, the spin-split band leads to a spin-momentum locking Fermi surface from which novel superconductivity emerges. Our demonstration not only expands the scope of spintronic materials but also enhances the understanding of SOI related superconductivity. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: 12 pages, 6 figures

Thermoelectric phase diagram of the SrTiO3-SrNbO3 solid solution system

Authors: Yuqiao Zhang, Bin Feng, Hiroyuki Hayashi, Tetsuya Tohei, Isao Tanaka, Yuichi Ikuhara, Hiromichi Ohta

Abstract: Thermoelectric energy conversion - the exploitation of the Seebeck effect to convert waste heat into electricity - has attracted an increasing amount of research attention for energy harvesting technology. Niobium-doped strontium titanate (SrTi1-xNbxO3) is one of the most promising thermoelectric material candidates, particularly as it poses a much lesser environmental risk in comparison to materi… ▽ More Thermoelectric energy conversion - the exploitation of the Seebeck effect to convert waste heat into electricity - has attracted an increasing amount of research attention for energy harvesting technology. Niobium-doped strontium titanate (SrTi1-xNbxO3) is one of the most promising thermoelectric material candidates, particularly as it poses a much lesser environmental risk in comparison to materials based on heavy metal elements. Two-dimensional electron confinement, e.g. through the formation of superlattices or two-dimensional electron gases, is recognized as an effective strategy to improve the thermoelectric performance of SrTi1-xNbxO3. Although electron confinement is closely related to the electronic structure, the fundamental electronic phase behavior of the SrTi1-xNbxO3 solid solution system has yet to be comprehensively investigated. Here, we present a thermoelectric phase diagram for the SrTi1-xNbxO3 (0.05 =< x =< 1) solid solution system, which we derived from the characterization of epitaxial films. We observed two thermoelectric phase boundaries in the system, which originate from the step-like decrease in carrier effective mass at x ~ 0.3, and from a local minimum in carrier relaxation time at x ~ 0.5. The origins of these phase boundaries are considered to be related to isovalent/heterovalent B-site substitution: parabolic Ti 3d orbitals dominate electron conduction for compositions with x < 0.3, whereas the Nb 4d orbital dominates when x > 0.3. At x ~ 0.5, a tetragonal distortion of the lattice, in which the B-site is composed of Ti4+ and Nb4+ ions, leads to the formation of tail-like impurity bands, which maximizes the electron scattering. These results provide a foundation for further research into improving the thermoelectric performance of SrTi1-xNbxO3. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: 20 pages, 6 figures

Journal ref: J. Appl. Phys. 121, 185102-1-7(2017)

Unconventional superconductivity in the BiS$_2$-based layered superconductor NdO$_{0.71}$F$_{0.29}$BiS$_2$

Authors: Yuichi Ota, Kozo Okazaki, Haruyoshi Q. Yamamoto, Takashi Yamamoto, Shuntaro Watanabe, Chuangtian Chen, Masanori Nagao, Satoshi Watauchi, Isao Tanaka, Yoshihiko Takano, Shik Shin

Abstract: We investigate the superconducting-gap anisotropy in one of the recently discovered BiS$_2$-based superconductors, NdO$_{0.71}$F$_{0.29}$BiS$_2$ ($T_c$ $\sim$ 5 K), using laser-based angle-resolved photoemission spectroscopy. Whereas the previously discovered high-$T_c$ superconductors such as copper oxides and iron-based superconductors, which are believed to have unconventional superconducting m… ▽ More We investigate the superconducting-gap anisotropy in one of the recently discovered BiS$_2$-based superconductors, NdO$_{0.71}$F$_{0.29}$BiS$_2$ ($T_c$ $\sim$ 5 K), using laser-based angle-resolved photoemission spectroscopy. Whereas the previously discovered high-$T_c$ superconductors such as copper oxides and iron-based superconductors, which are believed to have unconventional superconducting mechanisms, have $3d$ electrons in their conduction bands, the conduction band of BiS$_2$-based superconductors mainly consists of Bi 6$p$ electrons, and hence the conventional superconducting mechanism might be expected. Contrary to this expectation, we observe a strongly anisotropic superconducting gap. This result strongly suggests that the pairing mechanism for NdO$_{0.71}$F$_{0.29}$BiS$_2$ is unconventional one and we attribute the observed anisotropy to competitive or cooperative multiple paring interactions. △ Less

Submitted 24 March, 2017; originally announced March 2017.

Journal ref: Phys. Rev. Lett. 118, 167002 (2017)

First-principles high-throughput screening of shape-memory alloys based on energetic, dynamical, and structural properties

Authors: Joohwi Lee, Yuji Ikeda, Isao Tanaka

Abstract: First-principles-based materials screening is systematically performed to discover new combinations of chemical elements possibly making shape-memory alloys (SMAs). The B2, D03, and L21 crystal structures are considered as the parent phases, and the 2H and 6M structures are considered as the martensitic phases. 3,384 binary and 3,243 ternary alloys (6,627 in total) with stoichiometric composition… ▽ More First-principles-based materials screening is systematically performed to discover new combinations of chemical elements possibly making shape-memory alloys (SMAs). The B2, D03, and L21 crystal structures are considered as the parent phases, and the 2H and 6M structures are considered as the martensitic phases. 3,384 binary and 3,243 ternary alloys (6,627 in total) with stoichiometric composition ratios are investigated by the materials screening in terms of energetic and dynamical stabilities of the martensitic phases as well as structural compatibility between the parent and the martensitic phases. 187 alloys are found to survive after the screening. Some of the surviving alloys are constituted by the chemical elements already widely used in SMAs, but other various metallic elements are also found in the surviving alloys. The energetic stability of the surviving alloys is further analyzed by comparison with the data in Materials Project Database (MPD) to examine the alloys which may occur phase separation or transition. △ Less

Submitted 9 February, 2017; originally announced February 2017.

Comments: 8 figures. This will be submitted to PRB

Journal ref: npj Computational Materials 3, 52 (2017)

Temperature-dependent phonon spectra of magnetic random solid solutions

Authors: Yuji Ikeda, Fritz Körmann, Biswanath Dutta, Abel Carreras, Atsuto Seko, Jörg Neugebauer, Isao Tanaka

Abstract: A first-principles-based method for computing phonons of magnetic random solid solutions including thermal magnetic fluctuations is developed. The method takes fluctuations of force constants (FCs) due to magnetic excitations as well as due to chemical disorder into account. The developed approach correctly predicts the experimentally observed unusual phonon hardening of a transverse acoustic mode… ▽ More A first-principles-based method for computing phonons of magnetic random solid solutions including thermal magnetic fluctuations is developed. The method takes fluctuations of force constants (FCs) due to magnetic excitations as well as due to chemical disorder into account. The developed approach correctly predicts the experimentally observed unusual phonon hardening of a transverse acoustic mode in Fe-Pd an Fe-Pt Invar alloys with increasing temperature. This peculiar behavior, which cannot be explained within a conventional, harmonic picture, turns out to be a consequence of thermal magnetic fluctuations. △ Less

Submitted 8 February, 2017; originally announced February 2017.

First-principles calculation of the lattice thermal conductivities of αあるふぁ-, βべーた-, and γがんま-Si$_3$N$_4$

Authors: Kazuyoshi Tatsumi, Atsushi Togo, Isao Tanaka

Abstract: Lattice thermal conductivities (LTCs) of αあるふぁ-, βべーた-, and γがんま-Si$_3$N$_4$ single crystals are investigated from ab initio anharmonic lattice dynamics, within the single-mode relaxation-time approximation of the linearized phonon Boltzmann transport equation. At a temperature of 300 K, a κかっぱ$_{xx}$ of 70 and a κかっぱ$_{zz}$ of 98 (in units of Wm$^{-1}$K$^{-1}$) are obtained for αあるふぁ-Si$_3$N$_4$. For βべーた-Si$_3$N$_4$,… ▽ More Lattice thermal conductivities (LTCs) of αあるふぁ-, βべーた-, and γがんま-Si$_3$N$_4$ single crystals are investigated from ab initio anharmonic lattice dynamics, within the single-mode relaxation-time approximation of the linearized phonon Boltzmann transport equation. At a temperature of 300 K, a κかっぱ$_{xx}$ of 70 and a κかっぱ$_{zz}$ of 98 (in units of Wm$^{-1}$K$^{-1}$) are obtained for αあるふぁ-Si$_3$N$_4$. For βべーた-Si$_3$N$_4$, κかっぱ$_{xx}$ and κかっぱ$_{zz}$ are found to be 71 and 194, respectively, which are consistent with the reported experimental values of 69 and 180 for individual βべーた-Si$_3$N$_4$ grains in a ceramic. The theoretical κかっぱ$_{xx}$ values of αあるふぁ- and βべーた-Si$_3$N$_4$ are comparable, while the κかっぱ$_{zz}$ value of βべーた-Si$_3$N$_4$ is almost twice that of αあるふぁ-Si$_3$N$_4$, which demonstrates the very large anisotropy in the LTC of the βべーたphase. It is found that the large anisotropy in the LTC of βべーた-Si$_3$N$_4$ was caused by the elongated Brillouin zone along the c* axis, where the acoustic phonons mostly contribute to LTC and have large group velocities even near the Brillouin zone boundary. The LTC of γがんま-Si$_3$N$_4$ is 81, which is as small as that of αあるふぁ-Si$_3$N$_4$, although γがんま-Si$_3$N$_4$ has much larger elastic constants. This means that elastic constants are not always a good indicator of LTC. We show that knowing the detailed distributions of both the group velocities and phonon lifetimes in the Brillouin zones is important for characterizing the LTC of the three phases. △ Less

Submitted 26 December, 2016; originally announced December 2016.

Comments: 20 pages including 6 figures and 2 tables

MSC Class: 65Zxx

Representation of compounds for machine-learning prediction of physical properties

Authors: Atsuto Seko, Hiroyuki Hayashi, Keita Nakayama, Akira Takahashi, Isao Tanaka

Abstract: The representations of a compound, called "descriptors" or "features", play an essential role in constructing a machine-learning model of its physical properties. In this study, we adopt a procedure for generating a systematic set of descriptors from simple elemental and structural representations. First it is applied to a large dataset composed of the cohesive energy for about 18000 compounds com… ▽ More The representations of a compound, called "descriptors" or "features", play an essential role in constructing a machine-learning model of its physical properties. In this study, we adopt a procedure for generating a systematic set of descriptors from simple elemental and structural representations. First it is applied to a large dataset composed of the cohesive energy for about 18000 compounds computed by density functional theory (DFT) calculation. As a result, we obtain a kernel ridge prediction model with a prediction error of 0.041 eV/atom, which is close to the "chemical accuracy" of 1 kcal/mol (0.043 eV/atom). The procedure is also applied to two smaller datasets, i.e., a dataset of the lattice thermal conductivity (LTC) for 110 compounds computed by DFT calculation and a dataset of the experimental melting temperature for 248 compounds. We examine the performance of the descriptor sets on the efficiency of Bayesian optimization in addition to the accuracy of the kernel ridge regression models. They exhibit good predictive performances. △ Less

Submitted 21 December, 2016; v1 submitted 25 November, 2016; originally announced November 2016.

Comments: 11 pages, 5 figures

Journal ref: Phys. Rev. B 95, 144110 (2017)