Source author record

Toshiaki Iitaka

Toshiaki Iitaka 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

12works
8topics
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

12 published item(s)

preprint2026arXiv

The Impact of Ionic Anharmonicity on Superconductivity in Metal-Stuffed B-C Clathrates

Metal-stuffed B$-$C compounds with sodalite clathrate structure have captured increasing attention due to their predicted exceptional superconductivity above liquid nitrogen temperature at ambient pressure. However, by neglecting the quantum lattice anharmonicity, the existing studies may result in an incomplete understanding of such a lightweight system. Here, using state-of-the-art ab initio methods incorporating quantum effects and machine learning potentials, we revisit the properties of a series of $XY$$\text{B}_{6}\text{C}_{6}$ clathrates where $X$ and $Y$ are metals. Our findings show that ionic quantum and anharmonic effects can harden the $E_g$ and $E_u$ vibrational modes, enabling the dynamical stability of 15 materials previously considered unstable in the harmonic approximation, including materials with previously unreported ($XY$)$^{1+}$ state, which is demonstrated here to be crucial to reach high critical temperatures. Further calculations based on the anisotropic Migdal-Eliashberg equation demonstrate that the $T_\text{c}$ values for KRb$\text{B}_{6}\text{C}_{6}$ and Rb$\text{B}_{3}\text{C}_{3}$ among these stabilized compounds are 102 and 115 K at 0 and 15 GPa, respectively, both being higher than $T_\text{c}$ of 92 K of KPb$\text{B}_{6}\text{C}_{6}$ at the anharmonic level. These record-high $T_\text{c}$ values, surpassing liquid nitrogen temperatures, emphasize the importance of anharmonic effects in stabilizing B-C clathrates with large electron-phonon coupling strength and advancing the search for high-$T_\text{c}$ superconductivity at (near) ambient pressure.

preprint2020arXiv

Nonlinear optical response from quantum kinetic equation

Motivated by the nonlinear Hall effect observed in topological semimetals, we studied the photocurrent by the quantum kinetic equation. We recovered the shift current and injection current discovered by Sipe et al., and the nonlinear Hall current induced by Berry curvature dipole (BCD) proposed by Inti Sodemann and Liang Fu. Especially, we further proposed that 3-form tensor can also induce photocurrent, in addition to the Berry curvature and BCD. This work will supplement the existing mechanisms for photocurrent. In contrast to the shift current induced by shift vector, all photocurrents induced by gradient/curl of Berry curvature, and high rank tensor require circularly polarized light and topologically non-trivial band structure, viz. non-vanishing Berry curvature.

preprint2020arXiv

Random Phase Product Sate for Canonical Ensemble

Method of random phase product state (RPPS) is proposed to calculate canonical ensemble average of quantum systems described with matrix product states and also with tensor network states in general. The RPPS method is an extension of the method of random phase state for full Hilbert space representation. The validity of the method is confirmed by comparing the average energy of N-site antiferromagnetic spin-1/2 Heisenberg chain model with open boundary conditions with the result of direct method (for up to N=14) and minimally entangled typical thermal state (METTS) method (for N=100). Numerical advantages of the RPPS such as parallelization, combined calculation of thermal averages at different temperatures, parameters for controlling error are discussed. View point of self-averaging for the super-convergence of random state method is emphasized in addition to that of typicality.

preprint2015arXiv

GdN thin-film: Chern Insulating State on Square Lattice

Using first-principles calculations, we predict a Chern insulating phase in thin films of the ferromagnetic semi-metal GdN. In contrast to previously proposed Chern insulator candidates, which mostly rely on honeycomb lattices, this system affords a great chance to realize the quantum anomalous Hall Effect on a square lattice without either a magnetic substrate or transition metal doping, making synthesis easier. The band inversion between 5d-orbitals of Gd and 2p-orbitals of N is verified by first-principles calculation based on density functional theory, and the band gap can be as large as 100 meV within GdN trilayer. With further increase of film thickness, the band gap tends to close and the metallic bulk property becomes obvious.

preprint2012arXiv

First Principles Calculation of Helical Spin Order in Iron Perovskite SrFeO3 and BaFeO3

Motivated by recent discovery of ferromagnetism in cubic perovskite BaFeO3 under small magnetic field, we investigate spin order in BaFeO3 and isostructual SrFeO3 by the first principles calculation. The on-site Coulomb and exchange interactions are necessary for the helical spin order consistent with experiments. SrFeO3 exhibits stable G-type helical order, while A- and G-type helical orders in BaFeO3 are almost degenerate at short propagating vector with tiny energetic barrier with respect to ferromagnetic spin order, explaining ferromagnetism under small field. The results are consistent with model calculation where negative charge-transfer energy is explicitly taken into account.

preprint2012arXiv

First principles molecular dynamics study of filled ice hydrogen hydrate

We investigated structural changes, phase diagram, and vibrational properties of hydrogen hydrate in filled-ice phase C2 by using first principles molecular dynamics simulation. It was found that the experimentally reported 'cubic' structure is unstable at low temperature and/or high pressure. The 'cubic' structure reflects the symmetry at high (room) temperature where the hydrogen bond network is disordered and the hydrogen molecules are orientationally disordered due to thermal rotation. In this sense, the 'cubic' symmetry would definitely be lowered at low temperature where the hydrogen bond network and the hydrogen molecules are expected to be ordered. At room temperature and below 30 GPa, it is the thermal effects that play an essential role in stabilizing the structure in 'cubic' symmetry. Above 60 GPa, the hydrogen bonds in the framework would be symmetrized and the hydrogen bond order-disorder transition would disappear. These results also suggest the phase behavior of other filled-ice hydrates. In the case of rare gas hydrate, there would be no guest molecues rotation-nonrotation transition since the guest molecules keep their spherical symmetry at any temperature. On the contrary methane hydrate MH-III would show complex transitions due to the lower symmetry of the guest molecule. These results would encourage further experimental studies, especially NMR spectroscopy and neutron scattering, on the phases of filled-ice hydrates at high pressures and/or low temperatures.

preprint2012arXiv

Pressure Induced Ferromagnetism in Cubic Perovskite SrFeO3 and BaFeO3

The spin order in cubic perovskite SrFeO3 and BaFeO3 under high pressure is studied by density functional theory (DFT) calculation with local spin density approximation plus Hubbard U (LSDA+U). At ambient pressure, A-type and G-type helical spin orders are almost degenerate in BaFeO3 whose lattice constant is 3.97Å . When the lattice constant is reduced to 3.85 Å which is same as the lattice constant of SrFeO3 at ambient pressure, G-type helical spin order becomes stable, being consistent with SrFeO3. This is because superexchange interaction is enhanced as compared with double exchange interaction. Phase transition from helical spin state to ferromagnetic state in both SrFeO3 and BaFeO3 takes place if the lattice constant is further reduced to 3.70 Å. This is because reduced local spin moment weakens the contribution from superexchange interaction. Our result agrees with recent experimental result of BaFeO3 under high pressure. Additionally, our calculation predicts that half-metal BaFeO3 at ambient pressure will become a good metal under high pressure

preprint2012arXiv

Superconductive "sodalite"-like clathrate calcium hydride at high pressures

Hydrogen-rich compounds hold promise as high-temperature superconductors under high pressures. Recent theoretical hydride structures on achieving high-pressure superconductivity are composed mainly of H2 fragments. Through a systematic investigation of Ca hydrides with different hydrogen contents using particle-swam optimization structural search, we show that in the stoichiometry CaH6 a body-centred cubic structure with hydrogen that forms unusual "sodalite" cages containing enclathrated Ca stabilizes above pressure 150 GPa. The stability of this structure is derived from the acceptance by two H2 of electrons donated by Ca forming a "H4" unit as the building block in the construction of the 3-dimensional sodalite cage. This unique structure has a partial occupation of the degenerated orbitals at the zone centre. The resultant dynamic Jahn-Teller effect helps to enhance electron-phonon coupling and leads to superconductivity of CaH6. A superconducting critical temperature (Tc) of 220-235 K at 150 GPa obtained from the solution of the Eliashberg equations is the highest among all hydrides studied thus far.

preprint2010arXiv

Stability of ferroelectric ice

We theoretically study the stability conditions of the ferroelectric ice of the Cmc21 structure, which has been considered, for decades, one of the most promising candidates of the low temperature proton-ordered phase of pure ice Ih. It turned out that the Cmc21 structure is stable only with a certain amount of dopant and the true proton-ordered phase of pure ice Ih remains to be found at lower temperature. Implication for spin ice is mentioned.

preprint2009arXiv

Large scale simulation of quantum-mechanical molecular dynamics for nano-polycrystalline diamond

Quantum-mechanical molecular-dynamics simulations are carried out to explore possible precursor states of nano-polycrystalline diamond, a novel ultra-hard material produced directly from graphite. Large-scale simulation with 10^5 atoms is realized by using the ' order-N' simulation code 'ELSES' (http://www.elses.jp). The simulation starts with a diamond structure that contains initial structural defects and results in a mixture of graphite(sp^2)-like and diamond(sp^3)-like regions as nano-meter-scale domains. We speculate that the domains are metastable and are possible candidates of the precursor structures.

preprint2005arXiv

Temperature dependence of ESR intensity for the nanoscale molecular magnet V15

The electron spin resonance (ESR) of nanoscale molecular magnet ${\rm V}_{15}$ is studied. Since the Hamiltonian of ${\rm V}_{15}$ has a large Hilbert space and numerical calculations of the ESR signal evaluating the Kubo formula with exact diagonalization method is difficult, we implement the formula with the help of the random vector technique and the Chebyshev polynominal expansion, which we name the double Chebyshev expansion method. We calculate the temperature dependence of the ESR intensity of ${\rm V}_{15}$ and compare it with the data obtained in experiment. As another complementary approach, we also implement the Kubo formula with the subspace iteration method taking only important low-lying states into account. We study the ESR absorption curve below $100{\rm K}$ by means of both methods. We find that side peaks appear due to the Dzyaloshinsky-Moriya interaction and these peaks grows as temperature decreases.