Source author record

Koun Shirai

Koun Shirai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

First-principles study on the specific heat jump in the glass transition of silica glass and the Prigogine-Defay ratio

The most important characteristic of glass transition is a jump in specific heat $ΔC_{p}$. Despite its significance, no standard theory exists to describe it. In this study, first-principles molecular-dynamics (MD) simulations are used to describe the glass transition of silica glass, which presents many challenges. The novel view that state variables are extended to include the equilibrium positions of atoms $\{ \bar{\bf R}_{j} \}$ is fully used in analyzing the simulation results. Decomposing the internal energy into three components (structural, phonon, and thermal expansion energies) reveals that the jump $ΔC_{p}$ of silica glass is entirely determined by the component of structural energy. The reason for the small $ΔC_{p}$ is its high glass-transition temperature, which makes the fluctuation in the structural energy insensitive to temperature changes. This significantly affects how the Prigogine-Defay ratio $\itΠ$ is interpreted, which was previously unknown. The ratio $\itΠ$ represents the ratio of the total energy change to the contribution of thermal expansion energy at the glass transition. The general property, $\itΠ>1$, of glasses indicates that glass transitions occur mainly by changes in the structural energy. Silica glass is an extreme case in that the transition occurs entirely by the change in internal structure, such as the distribution of the bending angle of Si--O--Si bond.

preprint2022arXiv

First-principles study on the specific heat of glass at glass transition with a case study on glycerol

The standard method to determine the transition temperature (Tg) of glass transition is the jump in the specific heat. Despite this importance, standard theory for this jump is lacking. The difficulties encompass from lack of proper treatments of specific heat of liquids, hysteresis, to the timescale issue. The first part of this paper provides a non-empirical method to calculate specific heat of glass. The method consists of molecular dynamics (MD) simulations based on density-functional theory (DFT) and thermodynamics methods. The total-energy approach based on DFT enables us to calculate specific heat, irrespective of solids or liquids. A serious problem for glass-transition states is involvement of complicated energy dissipation processes, which is resolved by adiabatic MD simulations. The problems of hysteresis and the timescale issue are alleviated by restricting the scope of calculations to equilibrium states only. The second part of this paper describes an application of the theory to the specific-heat jump (Delta Cp) of glycerol in order to show the validity of the methods. Despite severe limitations due to the small size of supercells, a reasonable value for the specific-heat jump is obtained. By decomposing Delta Cp into contributions of the structural energy, phonon, and thermal expansion parts, we have a sound interpretation for the Delta Cp: the major contribution to Delta Cp comes from the change in the structural energy. From this, a neat energy diagram about the glass transition is obtained. An outcome of this study is verification of the empirical relationship between the fragility and specific-heat jump. These two energies are scaled by the ratio k=Tg/ Delta T, where Delta T is the width of the transition, through which the two quantities are interrelated. This is useful for organizing various relationships that are found empirically.

preprint2022arXiv

State Variables and Constraints in Thermodynamics of Solids and Their Implications

There is a common view in thermodynamics that the behavior of a macroscopic system can be described by only a few state variables. Although this is true for many cases, it is unclear whether it is meaningful to ask how many state variables are acceptable. This is indeed a problem when solids are investigated within the framework of thermodynamics, which is scarcely discussed in textbooks. The present study gives an answer to this question: the mean values of all the atom positions of a given solid together with the internal energy constitute a commensurate set of state variables (thermodynamic coordinates, TCs). The argument begins with constructing consistent definitions of equilibrium and TCs. TCs are created by the constraints which characterize the system under consideration. The values of TCs are uniquely determined in equilibrium and the mutual relationships between them constitute the fundamental relation of equilibrium (FRE). Specific heat can be deduced from the FRE. Therefore, the TCs of a solid must be to give a full expression of the specific heat in the entire range of temperature, from which the above conclusion is deduced. Contrary to the conventional view, an infinite number of the atom positions and their microscopic characters do not conflict with the principles of thermodynamics. The most important requirement for TCs to meet is the uniqueness of their values in equilibrium against random motions of the constituent particles. This conclusion is compatible with the principle of information theory that the information needed to determine the probability distribution of states is the expectation values of statistical variables. A few working examples of TCs in solids are given.

preprint2020arXiv

A thermodynamic description of the glass state and its application to glass transition

Many properties of solids such as the glass state are commonly treated as nonequilibrium phenomena, which involve many conceptual difficulties. However, few studies have addressed the problem of understanding equilibrium itself. Equilibrium is commonly assessed based on the assumption that its thermodynamic state should be determined solely by temperature and pressure. However, this assumption must be reappraised from the beginning through a rigorous definition of equilibrium because no rigorous proof for this assumption exists. Previous work showed that for solids, the equilibrium positions of all constituent atoms of the solid are state variables (TC). In this study, this conclusion is further elaborated starting from the principles of solid-state physics. This theory has been applied to glass materials. Results show that first, the glass state is an equilibrium state. Accordingly, the properties of a glass can be solely described by the present positions of the atom, irrespective of their previous history. Second, the glass transition, although a nonequilibrium phenomenon, is well described using TCs if the thermal part of the energy is be assumed to be well separated from the structural part. The only dynamic parameter involved in this approach is the relaxation time, which is uniquely determined using the present values of the TCs. Therefore, complicated functions describing the past history are unnecessary. This implies that the activation energy for the structural relaxation strongly depends on TCs. This finding provides a reasonable understanding of the large deviations from the Arrhenius law. The unrealistic values of the activation energy and of the pre-exponential factor of the relaxation time can be resolved on this basis. The theory is particularly suitable for experiments that do not involve hypothetical quantities and models; therefore, all quantities are measurable.

preprint2020arXiv

Interpretation of the apparent activation energy of glass transition

The temperature dependence of the viscosity of glass is a major concern in the field of glass research. Strong deviations from the Arrhenius law make the interpretation of the activation energy difficult. In the present study, a reasonable interpretation of the apparent activation energy is demonstrated along similar lines as those adopted in solid-state physics and chemistry. In contrast to the widely held view that phase transition in glass occurs at the reference temperature $T_{0}$ according to the Vogel--Fulcher--Tammann formula, in the present work the transition observed at the glass-transition temperature $T_{g}$ is regarded as a phase transition from the liquid to solid phases. A distinct feature of glass is that the energy barrier significantly changes in the transition range with width $ΔT_{g}$. This change in the energy barrier alters the manner in which the apparent activation energy constitutes the Arrhenius form. Analysis of available experimental data showed that the actual energy barrier is significantly smaller than the apparent activation energy, and importantly, the values obtained were in the reasonable range of energy expected for chemical bonds. The overestimation of the apparent activation energy depends on the ratio $T_{g}/ΔT_{g}$, which explains the existence of two types of glasses strong and fragile glasses. The fragility can be re-interpreted as an indication of the degree of increase in the energy barrier when approaching $T_{g}$ from high temperatures. Since no divergence in viscosity was observed below $T_{g}$, it is unlikely that a transition occurs at $T_{0}$.

preprint2020arXiv

Revisiting the stable structure of the Cu$_{4}$ complex in silicon

The photoluminescence (PL) spectrum of Cu-containing silicon has a sharp zero-phonon (ZP) band at 1.014 eV. The luminescence center corresponding to this band is called Cu$_{PL}$ and is known to have the local $C_{3v}$ symmetry. A recent measurement by ultrahigh-resolution PL spectroscopy revealed that the Cu$_{PL}$ center is a Cu$_{4}$ complex. Later, it was shown, by first-principles calculations, that the structure was Cu$_{(s)}$Cu$_{3(i)}$, that is, a complex composed of three interstitial Cu$_{(i)}$ atoms around a substitutional Cu$_{(s)}$ atom. This complex (called $C$-type) has the desired symmetry. However, in this study, we show that the lowest-energy structure is different. The tetrahedral structure Cu$_{4}$, called $T$-type, has the lowest energy, with the value being 0.26 eV lower than that of $C$-type. Between these two types, there is an energy barrier of 0.14 eV, which allows $C$-type to exist in a metastable state. Details of the electronic properties of the $T$-type complex are given, by comparing with $C$-type and other isovalent complexes such as Li$_{4}$. Whereas the Cu$_{4}$ tetrahedron is incorporated in silicon in a manner compatible with the tetrahedral network, it also has its own molecular orbitals that exhibit metallic characteristics, in contrast to other complexes. The ZP of the PL spectrum is very likely ascribed to the backflow mode of the Cu$_{4}$ tetrahedron.

preprint2016arXiv

Structure, non-stoichiometry, and geometrical frustration of alpha-tetragonal boron

Recent discoveries of supposedly pure alpha-tetragonal boron require to revisit its structure. The system is also interesting with respect to a new type of geometrical frustration in elemental crystals, which was found in beta-rhombohedral boron. Based on density functional theory calculations, the present study has resolved the structural and thermodynamic characteristics of pure alpha-tetragonal boron. Different from beta-rhombohedral boron, the conditions for stable covalent bonding (a band gap and completely filled valence bands) are almost fulfilled at a composition B_52 with two 4c interstitial sites occupied. This indicates that the ground state of pure alpha-tetragonal boron is stoichiometric. However, the covalent condition is not perfectly fulfilled because non-bonding in-gap states exist that cannot be eliminated. The half occupation of the 4c sites yields a macroscopic amount of residual entropy, which is as large as that of beta-rhombohedral boron. Therefore, alpha-tetragonal boron can be classified as an elemental crystal with geometrical frustration. Deviations from stoichiometry can occur only at finite temperatures. Thermodynamic considerations show that deviations delta from the stoichiometric composition B_(52+delta) are small and positive. For reported high-pressure syntheses conditions delta is predicted to be about 0.1 to 0.2. An important difference between pure and C- or N-containing alpha-tetragonal boron is found in the occupation of interstitial sites: the pure form prefers to occupy the 4c sites, whereas in C- or N-containing forms a mixture of 2a, 8h, and 8i sites are occupied. The present article provides relations of site occupation, delta values, and lattice parameters, which enable us to identify pure alpha-tetragonal and distinguish the pure form from other ones.

preprint2015arXiv

Impact of Ferroelectric Distortion on Thermopower in BaTiO3

We present a strategy for increasing thermoelectric performance by invoking ferroelectric distortion in a transition-metal oxide. By using an ab initio calculation approach with the Boltzmann transport equation, a large Seebeck coefficient is calculated in n-type BaTiO3. The polar structural distortion causes the peculiar dispersion of the lowest conduction band, which enhances the Seebeck coefficient along the polar direction. A microscopic mechanism of carrier concentration dependence of the Seebeck coefficient is discussed in terms of the band velocity distribution in the Brillouin zone.