Source author record

Koichi Kusakabe

Koichi Kusakabe 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

20works
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

20 published item(s)

preprint2022arXiv

Discovery of nanographene for hydrogen storage solving low reversibility issues

We found good reversibility of hydrogen uptake-release in vacancy-centered hexagonal armchair nanographene (VANG) based on density functional theory calculation. VANG has a triply hydrogenated vacancy (V$_{111}$) at the center and acts as a self-catalytic property to reduce an activation barrier of hydrogen uptake-release. We found remarkable features in an almost equal value of the activation energy barrier of 1.19 eV for hydrogen uptake and 1.25 eV for hydrogen release on V$_{111}$ of VANG. The dehydrogenation showed slightly exothermic and the hydrogenation became slightly endothermic, suggesting the efficiency of hydrogen uptake-release. In high hydrogen coverage, the quintuply hydrogenated vacancy (V$_{221}$) is formed with some hydrogenated located in the in-plane and armchair edges. This structure produces an exothermic hydrogen release from the in-plane with an energy barrier of not more than 2 eV. This finding potentially addresses the low reversibility issues in the organic chemical hydrides as hydrogen storage materials.

preprint2022arXiv

Theoretical Analysis on the Stability of 1-Pyrenebutanoic Acid Succinimidyl Ester Adsorbed on Graphene

The adsorbed structure of 1-pyrenebutanoic acid succinimidyl ester (PASE) on graphene was investigated based on density functional theory. We found two locally stable structures: a straight structure with the chainlike part of butanoic acid succinimidyl ester (BSE) lying down and a bent structure with the BSE part directed away from graphene, keeping the pyrene (Py) part adsorbed on graphene. Then, to elucidate the adsorption mechanism, we separately estimated the contributions of the Py and BSE parts to the entire PASE adsorption, and the adsorption effect of the BSE part was found to be secondary in comparison to the contribution of the Py. Next, the mobility of the BSE part at room temperature was confirmed by the activation energy barrier between straight and bent structures. To take account of the external environment, we considered the presence of amino acids and the hydration effect by a three-dimensional reference interaction site model. The contributions of glycine molecules and the solvent environment to stabilizing the bent PASE structure relative to the straight PASE structure were found. Therefore, the effect of the external environment around PASE is of importance when the standing-up process of the BSE part from graphene is considered.

preprint2019arXiv

Interplanar stiffness in defect-free monocrystalline graphite

The interplanar bond strength in graphite has been identified to be very low owing to the contribution of the van der Waals interaction. However, in this study, we use microscopic picosecond ultrasound to demonstrate that the elastic constant, $C_{33}$, along the $c$ axis of defect-free monocrystalline graphite exceeds 45 GPa, which is higher than reported values by 20\%. Existing theories fail to reproduce this strongly correlated interplanar system, and our results, thus, indicate the necessity for improvement. Since the LDA+U+RPA method, including both random phase approximation correlation and short-range correlation in $p$ Wannier orbitals, shows better agreement with the observation than LDA or even than ACFDT-RPA, the experimental results indicate non-negligible electron correlation effects with respect to both the short-range and long-range interactions.

preprint2019arXiv

Spin-Current Control by Induced Electric-Polarization Reversal in Ni/hBN/Ni: A Cross-Correlation Material

We undertook an ab-initio study of hexagonal boron nitride (hBN) sandwiched between Ni(111) layers to examine the interface of this material structure. We considered Ni(111)/hBN/Ni(111) with a slab with three Ni atomic layers to determine the exact atom arrangement at the interface. The density functional theory calculations for 36 stacking arrangements, which are doubled with respect to the magnetic alignment of slabs in an anti-parallel configuration (APC) and parallel configuration (PC), revealed that the number of formed weak chemical bonds, in the pd-hybridization between the N and Ni atoms, is decisive. A maximum of two pd-hybridization bonds stabilized the structure, with APC proving to be the most favorable magnetic alignment, in line with the results of previous experimental studies. In the lowest energy state, an induced magnetic moment at an N site appears when N is moved closer to one of the Ni atoms. Interestingly, the moment direction is switched by the position of the N layer in the resulting bi-stable state with electrical polarization when APC is chosen. The transmission probability calculation of Ni/hBN/Ni having the determined interface structure at the center of the junction exhibits a spin-filtering effect where the spin-polarized current is controlled by the electric field when a field-induced reversal of the polarization is realized.

preprint2016arXiv

Theoretical Analysis on Pseudo-Degenerate Zero-Energy Modes in Vacancy-Centered Hexagonal Armchair Nanographene

Deriving mathematical expressions of two zero modes for a $π$-band tight-binding model, we identify a class of bipartite graphs having the same number of subgraph sites, where each graph represents one of the quasi-hexagonal nanographene molecule with a center vacancy (VANG). Indeed, in a VANG molecule, C$_{60}$H$_{24}$, showing stability in a density-functional simulation at the highest occupied level, there appear two pseudo-degenerate zero modes, a vacancy-centered quasi-localized zero mode, and extending zero mode with a $\sqrt{3} \times \sqrt{3}$ structure. Since there is a finite energy gap between these two zero-energy modes and the other modes, low-lying states composed of quasi-degenerate zero modes appear as magnetic multiplets. Thus, the unique magnetic characteristics derived in our theory are expected to hold for synthesized VANG molecules in reality.

preprint2014arXiv

Direct imaging of monovacancy-hydrogen complexes in single graphitic layer

Understanding how foreign chemical species bond to atomic vacancies in graphene layers can advance our ability to tailor the electronic and magnetic properties of defective graphenic materials. Here we use ultra-high vacuum scanning tunneling microscopy (UHV-STM) and density functional theory to identify the precise structure of hydrogenated single atomic vacancies in a topmost graphene layer of graphite and establish a connection between the details of hydrogen passivation and the electronic properties of a single atomic vacancy. Monovacancy-hydrogen complexes are prepared by sputtering of the graphite surface layer with low energy ions and then exposing it briefly to an atomic hydrogen environment. High-resolution experimental UHV-STM imaging allows us to determine unambiguously the positions of single missing atoms in the defective graphene lattice and, in combination with the ab initio calculations, provides detailed information about the distribution of low-energy electronic states on the periphery of the monovacancy-hydrogen complexes. We found that a single atomic vacancy where each sigma-dangling bond is passivated with one hydrogen atom shows a well-defined signal from the non-bonding pi-state which penetrates into the bulk with a (\sqrt 3 \times \sqrt 3)R30^ \circ periodicity. However, a single atomic vacancy with full hydrogen termination of sigma-dangling bonds and additional hydrogen passivation of the extended pi-state at one of the vacancy's monohydrogenated carbon atoms is characterized by complete quenching of low-energy localized states. In addition, we discuss the migration of hydrogen atoms at the periphery of the monovacancy-hydrogen complexes which dramatically change the vacancy's low-energy electronic properties, as observed in our low-bias high-resolution STM imaging.

preprint2014arXiv

Magnetic correlation effects by the topological zero mode in a hydrogenated graphene vacancy $V_{111}$

Electron correlation effects caused by the topological zero mode of a hydrogenated graphene vacancy, $V_{111}$, with three adsorbed hydrogen atoms is discussed theoretically. A Kondo model is derived from the multi-reference representation of the density functional theory, where exchange scattering processes between the zero mode and low-energy modes in the Dirac cones are estimated. Even when the Dirac cone is slightly off from the charge neutral point, a finite on-site correlation energy, $U_0$, for the zero mode of an isolated $V_{111}$ allows the half-filling of the localized level giving a spin $s=1/2$. The anti-ferromagnetic Kondo screening mediated by higher order scattering processes becomes dominant in the dilute limit of the vacancies. Our estimation of relevant two body interactions certifies appearance of the Kondo effect at low temperatures.

preprint2014arXiv

Orbital mixture effect on the Fermi surface-$T_c$ correlation in the cuprate superconductors --- bilayer vs single layer

By constructing $d_{x^2-y^2}-d_{z^2}$ two-orbital models from first principles, we have obtained a systematic correlation between the Fermi surface warping and the evaluated $T_c$ for various bilayer as well as single-layer cuprates. This reveals that smaller mixture of the $d_{z^2}$ orbital component on the Fermi surface leads to both of larger Fermi surface warping and higher $T_c$. The theoretical correlation strikingly resembles a systematic plot for the experimentally observed $T_c$ against the Fermi surface warping due to Pavarini {\it et al.} [Phys. Rev. Lett. {\bf 87}, 047003 (2001)], and the present result unambiguously indicates that the $d_{z^2}$ mixture is a key factor that determines $T_c$ in the cuprates.

preprint2013arXiv

Material Optimization for Fermi Surface Shape Control of Tl-based Cuprate Superconductors

To show an optimization method of element substitution for cuprate superconductors, we investigate Fermi surface shape of TlR2A2Cu3O9 with R=La, Y, and A=Li, Na, K, Rb, Cs. We adopt the generalized gradient approximation in the density-functional theory (DFT-GGA) for the study of over-doped phases of these unknown cuprates. The electronic structures of crystals optimized by DFT-GGA show systematic element dependence in a Fermi surface shape controlling parameter, r, of Cu dx2-y2 bands, where nearly absent dz2 component at the Fermi level and smaller r keeping t1 suggest enhancement of the superconducting transition temperature within the spin-fluctuation mechanism. For TlYRb2Cu3O9, smaller r by a reduction factor larger than 10% compared to a reference system of TlBa2Ca2Cu3O9 (TBCCO) appears in the outer CuO2 plane, but with 12 % reduction in t1. For TlR2Li2Cu3O9, (R=Y, La), smaller r by a factor larger than 1% appears keeping t1 as large as TBCCO in the inner plane. Our method may be used to predict an optimized material structure referencing known cuprate superconductors.

preprint2011arXiv

Determination of Boundary Scattering, Intermagnon Scattering, and the Haldane Gap in Heisenberg Chains

Low-lying magnon dispersion in a S=1 Heisenberg antiferromagnetic (AF) chain is analyzed using the non-Abelian DMRG method. The scattering length $a_{\rm b}$ of the boundary coupling and the inter-magnon scattering length $a$ are determined. The scattering length $a_{\rm b}$ is found to exhibit a characteristic diverging behavior at the crossover point. In contrast, the Haldane gap $Δ$, the magnon velocity $v$, and $a$ remain constant at the crossover. Our method allowed estimation of the gap of the S=2 AF chain to be $Δ= 0.0891623(9)$ using a chain length longer than the correlation length $ξ$.

preprint2011arXiv

Hyperbolic Deformation Applied to S = 1 Spin Chains - Scaling Relation in Excitation Energy -

We investigate excitation energies of hyperbolically deformed S = 1 spin chains, which are specified by the local energy scale f_j^{~} = \cosh j λ, where j is the lattice index and λis the deformation parameter. The elementary excitation is well described by a quasiparticle hopping model, which is also expressed in the form of hyperbolic deformation. It is possible to estimate the excitation gap Δin the uniform limit λ\rightarrow 0, by means of a finite size scaling with respect to the system size N and the deformation parameter λ.

preprint2010arXiv

A Theoretical Simulation of Deformed Carbon Nanotubes with Adsorbed Metal Atoms: Enhanced Reactivity by Deformation

First-principles simulations were performed to investigate reaction of carbon nanotubes with adsorbed metal atoms. Mechanical modification of their structures enhances chemical reactivity of carbon nanotubes. Adsorption of a tungsten, tantalum, or niobium atom on a (5, 0) nanotube with a Stone-Wales defect was shown to have characteristically strong chemisorption. Bond-breaking in the carbon-carbon network and formation of a local metal-carbon complex were observed during the simulation. Adsorption of W, Ta, Ni or Mo on a twisted (5, 0) nanotube showed a preferred breaking of several bonds, even creating an opening in the wall. The enhanced chemical reactivity of deformed nanotubes is characterized by formation of a metal-carbon complex. Applications of the reaction are suggested.

preprint2010arXiv

Scaling Relation for Excitation Energy Under Hyperbolic Deformation

We introduce a one-parameter deformation for one-dimensional (1D) quantum lattice models, the hyperbolic deformation, where the scale of the local energy is proportional to cosh lambda j at the j-th site. Corresponding to a 2D classical system, the deformation does not strongly modify the ground state. In this situation, the effective Hamiltonian of the quantum system shows that the quasi particle is weakly bounded around the center of the system. By analyzing this binding effect, we derive scaling relations for the mean-square width <w^2> of confinement, the energy correction with respect to the excitation gap Δ, and the deformation parameter $λ$. This finite-size scaling allows us to investigate excitation gap of 1D non-deformed bulk quantum systems.

preprint2010arXiv

Theorems on ground-state phase transitions in Kohn-Sham models given by the Coulomb density functional

Some theorems on derivatives of the Coulomb density functional with respect to the coupling constant $λ$ are given. Consider an electron density $n_{GS}({\bf r})$ given by a ground state. A model Fermion system with the reduced coupling constant, $λ<1$, is defined to reproduce $n_{GS}({\bf r})$ and the ground state energy. Fixing the charge density, possible phase transitions as level crossings detected in a value of the reduced density functional happen only at discrete points along the $λ$ axis. If the density is $v$-representable also for $λ<1$, accumulation of phase transition points is forbidden when $λ\rightarrow 1$. Relevance of the theorems for the multi-reference density functional theory is discussed.

preprint2009arXiv

Pair-Hopping Mechanism for Layered Superconductors

We propose a possible charge fluctuation effect expected in layered superconducting materials. In the multireference density functional theory, relevant fluctuation channels for the Josephson coupling between superconducting layers include the interlayer pair hopping derived from the Coulomb repulsion. When interlayer single-electron tunneling processes are irrelevant in the Kohn-Sham electronic band structure calculation, the two-body effective interactions stabilize a superconducting phase. This state is also regarded as a valence-bond solid in a bulk electronic state. The hidden order parameters coexist with the superconducting order parameter when the charging effect of a layer is comparable to the pair hopping. Relevant materials structures favorable for the pair-hopping mechanism are discussed.

preprint2008arXiv

Finite-size scaling of string order parameters characterizing the Haldane phase

We have developed a numerical procedure to clarify the critical behavior near a quantum phase transition by analyzing a multi-point correlation function characterizing the ground state. This work presents a successful application of this procedure to the string order parameter of the S = 1 XXZ chain with uniaxial single-ion anisotropy. The finite-size string correlation function is estimated by the density matrix renormalization group method. We focus on the gradient of the inversedsystem- size dependence of the correlation function on a logarithmic plot. This quantity shows that the finite-size scaling sensitively changes at the critical point. The behavior of the gradient with increasing system size is divergent, stable at a finite value, or rapidly decreases to zero when the system is in the disordered phase, at the critical point, or in the ordered phase, respectively. The analysis of the finite-size string correlation functions allows precise determination of the boundary of the Haldane phase and estimation of the critical exponent of the correlation length. Our estimates of the transition point and the critical exponents, which are determined only by the ground-state quantities, are consistent with results obtained from the analysis of the energy-level structure. Our analysis requires only the correlation functions of several finite sizes under the same condition as a candidate for the long-range order. The quantity is treated in the same manner irrespective of the kind of elements which destroy the order concerned. This work will assist in the development of a method to directly observe quantum phase transitions.

preprint2008arXiv

General Rule and Materials Design of Negative Effective U System for High-T_c Superconductivity

Based on the microscopic mechanisms of (1) charge-excitation-induced negative effective U in s^1 or d^9 electronic configurations, and (2) exchange-correlation-induced negative effective U in d^4 or d^6 electronic configurations, we propose a general rule and materials design of negative effective U system in itinerant (ionic and metallic) system for the realization of high-T_c superconductors. We design a T_c-enhancing layer (or clusters) of charge-excitation-induced negative effective $U$ connecting the superconducting layers for the realistic systems.

preprint2008arXiv

Generation of Nano-Catalyst Particles by Spinodal Nano-Decomposition in Perovskite

A new mechanism of nano-catalyst generation based on the spinodal nano-decomposition in self-regenerating perovskite catalysts for automotive-emissions control is proposed. To demonstrate existence of the spinodal nano-decomposition in real perovskite catalysts, we performed first-principles calculations to evaluate the free energy of La(Fe$_{1-x}$Pd$_x$)O$_3$ and La(Fe$_{1-x}$Rh$_x$)O$_3$. The result indicates appearance of a spinodal region in the phase diagram of each material. Formation of nano-catalyst particles in the perovskite host matrix is crucial for the self-regeneration of perovskite catalyst. Based on the spinodal nano-decomposition model, possible materials are designed for new three-way catalyst with no contents of precious metal.

preprint2008arXiv

Two-Site shift Product Wave Function Renormalization Group Method Applied to Quantum Systems

We report a way of wave function estimation for the density matrix renormalization group (DMRG) method applied to quantum systems, which has 2-site modulation, when the system size extension is necessary in both the finite and the infinite algorithms. The estimation is performed by renormalization group (RG) transformation applied to the ground state wave function, which is represented as the matrix product. This RG scheme is known as the product wave function renormalization group (PWFRG) method. In order to treat the 2-site modulation, the operation of the RG transformation is shifted by amount of 2 lattice sites. It turns out that this 2-site shift algorithm provides better wave function estimation in the thermodynamic limit, compared with the previously known PWFRG algorithm.

preprint2004arXiv

Zinc-blende CaP, CaAs and CaSb as half-metals: A new route to magnetism in calcium compounds

Existence of ferromagnetism in bulk calcium compounds is discovered theoretically. First-principles calculations of calcium phosphide, calcium arsenide and calcium antimonide in the zinc-blende structure have been performed to show the half-metallic ground state in each optimized stable structure. Magnetism comes from spin-polarization of electrons in $p$-orbitals of P, As or Sb and $d$-orbitals of calcium atoms. The half-metallicity is analogous to the half-metallic zinc-blende compounds, {\it e.g.} CrAs or CrSb, but the predicted compounds become ferromagnetic without transition metals. In (In$_{1-x}$Ca$_x$)Sb, the magnetism remains to be stable in a range of the doping rate ($x>0.8$).