Researcher profile

Koun Shirai

Koun Shirai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
4topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 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.