Source author record

Yuki Nagai

Yuki Nagai 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

53works
14topics
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

53 published item(s)

preprint2022arXiv

Extensively parallelizable chiral fermion

Chiral symmetry is a key to investigating quantum physics, from condensed matter to particle physics. We propose a novel way of realizing a chiral fermion, known as the overlap-Dirac operator, without explicitly calculating the low modes of the Wilson-Dirac operator. We introduce a projection operator inspired by the Sakurai-Sugiura method and formulate the exact sign function and overlap-Dirac operator with a contour-integral form. Like the Sakurai-Sugiura method, the proposing method is multi-scale parallelizable, which fits the multi-core/multi-GPGPU paradigm. We confirm that the quality of chiral symmetry realized with the proposed method is sufficient for double precision. We evaluate the strong scaling of the proposing method.

preprint2022arXiv

Intrinsic vortex pinning in superconducting quasicrystals

We numerically show that a vortex pinning occurs in superconducting quasicrystal without impurities and defects. This vortex pinning is intrinsic since the superconducting order parameter in quasicrystal is always inhomogeneous due to the lack of the translational symmetry. We propose that experiments influenced by vortex pinning effects can detect the atomic-scale inhomogeneous superconducting order parameter in quasicrystals. We develop a numerical method to solve the Bogoliubov-de Gennes equations and gap equations in large systems, which is based on the localized-Krylov subspace and a sparse modeling technique. Two two-dimensional quasicrystals, the Penrose and Amman-Beenker tiling, are considered.

preprint2021arXiv

Effective Ruderman-Kittel-Kasuya-Yosida-like interaction in diluted double-exchange model: self-learning Monte Carlo approach

We study the site-diluted double exchange (DE) model and its effective Ruderman-Kittel-Kasuya-Yosida-like interactions, where localized spins are randomly distributed, with the use of the Self-learning Monte Carlo (SLMC) method. The SLMC method is an accelerating technique for Markov chain Monte Carlo simulation using trainable effective models. We apply the SLMC method to the site-diluted DE model to explore the utility of the SLMC method for random systems. We check the acceptance ratios and investigate the properties of the effective models in the strong coupling regime. The effective two-body spin-spin interaction in the site-diluted DE model can describe the original DE model with a high acceptance ratio, which depends on temperatures and spin concentration. These results support a possibility that the SLMC method could obtain independent configurations in systems with a critical slowing down near a critical temperature or in random systems where a freezing problem occurs in lower temperatures.

preprint2020arXiv

First- and Second-Order Topological Superconductivity and Temperature-Driven Topological Phase Transitions in the Extended Hubbard Model with Spin-Orbit Coupling

The combination of spin-orbit coupling with interactions results in many exotic phases of matter. In this Letter, we investigate the superconducting pairing instability of the two-dimensional extended Hubbard model with both Rashba and Dresselhaus spin-orbit coupling within the mean-field level at both zero and finite temperature. We find that both first- and second-order time-reversal symmetry breaking topological gapped phases can be achieved under appropriate parameters and temperature regimes due to the presence of a favored even-parity $s+id$-wave pairing even in the absence of an external magnetic field or intrinsic magnetism. This results in two branches of chiral Majorana edge states on each edge or a single zero-energy Majorana corner state at each corner of the sample. Interestingly, we also find that not only does tuning the doping level lead to a direct topological phase transition between these two distinct topological gapped phases, but also using the temperature as a highly controllable and reversible tuning knob leads to different direct temperature-driven topological phase transitions between gapped and gapless topological superconducting phases. Our findings suggest new possibilities in interacting spin-orbit coupled systems by unifying both first- and higher-order topological superconductors in a simple but realistic microscopic model.

preprint2020arXiv

Majorana corner flat bands in two-dimensional second-order topological superconductors

In this paper we find that confining a second-order topological superconductor with a harmonic potential leads to a proliferation of Majorana corner modes. As a consequence, this results in the formation of Majorana corner flat bands which have a fundamentally different origin from that of the conventional mechanism. This is due to the fact that they arise solely from the one-dimensional gapped boundary states of the hybrid system that become gapless without the bulk gap closing under the increase of the trapping potential magnitude. The Majorana corner states are found to be robust against the strength of the harmonic trap and the transition from Majorana corner states to Majorana flat bands is merely a smooth crossover. As a harmonic trap can potentially be realized in heterostructures, this proposal paves a way to observe these Majorana corner flat bands in an experimental context.

preprint2020arXiv

N-independent Localized Krylov Bogoliubov-de Gennes Method: Ultra-fast Numerical Approach to Large-scale Inhomogeneous Superconductors

We propose the ultra-fast numerical approach to large-scale inhomogeneous superconductors, which we call the Localized Krylov Bogoliubov-de Gennes method (LK-BdG). In the LK-BdG method, the computational complexity of the local Green's function, which is used to calculate the local density of states and the mean-fields, does not depend on the system size $N$. The calculation cost of self-consistent calculations is ${\cal O}(N)$, which enables us to open a new avenue for treating extremely large systems with millions of lattice sites. To show the power of the LK-BdG method, we demonstrate a self-consistent calculation on the 143806-site Penrose quasicrystal lattice with a vortex and a calculation on 1016064-site two-dimensional nearest-neighbor square-lattice tight-binding model with many vortices. We also demonstrate that it takes less than 30 seconds with one CPU core to calculate the local density of states with whole energy range in 100-millions-site tight-binding model.

preprint2020arXiv

Self-learning Monte Carlo method with Behler-Parrinello neural networks

We propose a general way to construct an effective Hamiltonian in the Self-learning Monte Carlo method (SLMC), which speeds up Monte Carlo simulations by training an effective model to propose uncorrelated configurations in the Markov chain. Its applications are, however, limited. This is because it is not obvious to find the explicit form of the effective Hamiltonians. Particularly, it is difficult to make trainable effective Hamiltonians including many-body interactions. In order to overcome this critical difficulty, we introduce the Behler-Parrinello neural networks (BPNNs) as "effective Hamiltonian'' without any prior knowledge, which is used to construct the potential-energy surfaces in interacting many particle systems for molecular dynamics. We combine SLMC with BPNN by focusing on a divisibility of Hamiltonian and propose how to construct the element-wise configurations. We apply it to quantum impurity models. We observed significant improvement of the acceptance ratio from 0.01 (the effective Hamiltonian with the explicit form) to 0.76 (BPNN). This drastic improvement implies that the BPNN effective Hamiltonian includes many body interaction, which is omitted in the effective Hamiltonian with the explicit forms. The BPNNs make SLMC more promising.

preprint2020arXiv

Sparse modeling approach to obtaining the shear viscosity from smeared correlation functions

We propose the sparse modeling method to estimate the spectral function from the smeared correlation functions. We give a description of how to obtain the shear viscosity from the correlation function of the renormalized energy-momentum tensor (EMT) measured by the gradient flow method ($C(t,τ)$) for the quenched QCD at finite temperature. The measurement of the renormalized EMT in the gradient flow method reduces a statistical uncertainty thanks to its property of the smearing. However, the smearing breaks the sum rule of the spectral function and the over-smeared data in the correlation function may have to be eliminated from the analyzing process of physical observables. In this work, we demonstrate that the sparse modeling analysis in the intermediate-representation basis (IR basis), which connects between the Matsubara frequency data and real frequency data. It works well even using very limited data of $C(t,τ)$ only in the fiducial window of the gradient flow. We utilize the ADMM algorithm which is useful to solve the LASSO problem under some constraints. We show that the obtained spectral function reproduces the input smeared correlation function at finite flow-time. Several systematic and statistical errors and the flow-time dependence are also discussed.

preprint2020arXiv

Sparse modeling of large-scale quantum impurity models with low symmetries

Quantum embedding theories provide a feasible route for obtaining quantitative descriptions of correlated materials. However, a critical challenge is solving an effective impurity model of correlated orbitals embedded in an electron bath. Many advanced impurity solvers require the approximation of a bath continuum using a finite number of bath levels, producing a highly nonconvex, ill-conditioned inverse problem. To address this drawback, this study proposes an efficient fitting algorithm for matrix-valued hybridization functions based on a data-science approach, sparse modeling, and a compact representation of Matsubara Green's functions. The efficiency of the proposed method is demonstrated by fitting random hybridization functions with large off-diagonal elements as well as those of a 20-orbital impurity model for a high-Tc compound, LaAsFeO, at low temperatures (T). The results set quantitative goals for the future development of impurity solvers toward quantum embedding simulations of complex correlated materials.

preprint2019arXiv

Self-learning Hybrid Monte Carlo: A First-principles Approach

We propose a novel approach called Self-Learning Hybrid Monte Carlo (SLHMC) which is a general method to make use of machine learning potentials to accelerate the statistical sampling of first-principles density-functional-theory (DFT) simulations. The trajectories are generated on an approximate machine learning (ML) potential energy surface. The trajectories are then accepted or rejected by the Metropolis algorithm based on DFT energies. In this way the statistical ensemble is sampled exactly at the DFT level for a given thermodynamic condition. Meanwhile the ML potential is improved on the fly by training to enhance the sampling, whereby the training data set, which is sampled from the exact ensemble, is created automatically. Using the examples of $α$-quartz crystal SiO$_2^{}$ and phonon-mediated unconventional superconductor YNi$_2^{}$B$_2^{}$C systems, we show that SLHMC with artificial neural networks (ANN) is capable of very efficient sampling, while at the same time enabling the optimization of the ANN potential to within meV/atom accuracy. The ANN potential thus obtained is transferable to ANN molecular dynamics simulations to explore dynamics as well as thermodynamics. This makes the SLHMC approach widely applicable for studies on materials in physics and chemistry.

preprint2016arXiv

Critical Temperature Enhancement of Topological Superconductors: A Dynamical Mean Field Study

We show that a critical temperature Tc for spin-singlet two-dimensional superconductivity is enhanced by a cooperation between the Zeeman magnetic field and the Rashba spin-orbit coupling, where a superconductivity becomes topologically non-trivial below Tc. The dynamical mean field theory (DMFT) with the segment-based hybridization-expansion continuous-time quantum Monte Carlo impurity solver (ct-HYB) is used for accurately evaluating a critical temperature, without any Fermion sign problem. A strong-coupling approach shows that spin-flip driven local pair hopping leads to part of this enhancement, especially effects of the magnetic field. We propose physical settings suitable for verifying the present calculations, one-atom-layer system on Si(111) and ionic-liquid based electric double-layer transistors (EDLTs).

preprint2016arXiv

Excitation spectra and wave functions of quasiparticle bound states in bilayer Rashba superconductors

We study the excitation spectra and the wave functions of quasiparticle bound states at a vortex and an edge in bilayer Rashba superconductors under a magnetic field. In particular, we focus on the quasiparticle states at the zero energy in the pair-density wave state in a topologically non-trivial phase. We numerically demonstrate that the quasiparticle wave functions with zero energy are localized at both the edge and the vortex core if the magnetic field exceed the critical value.

preprint2016arXiv

Manifestation of chirality in the vortex lattice in a two-dimensional topological superconductor

We study the vortex lattice in a two-dimensional s-wave topological superconductor with Rashba spin-orbit coupling and Zeeman field by solving the Bogoliubov-de Gennes equations self-consistently for the superconducting order parameter. We find that when spin-orbit coupling is relatively weak, one of the two underlying chiralities in the topological superconducting state can be strongly manifest in the vortex core structure and govern the response of the system to vorticity and a nonmagnetic impurity where the vortex is pinned. The Majorana zero mode in the vortex core is found to be robust against the nonmagnetic impurity in that it remains effectively a zero-energy bound state regardless of the impurity potential strength and the major chirality. The spin polarization of the Majorana bound state depends on the major chirality for weak spin-orbit coupling, while it is determined simply by the vorticity when spin-orbit coupling is relatively strong.

preprint2016arXiv

Multi-band Eilenberger theory of superconductivity: Systematic low-energy projection

We propose the general multi-band quasiclassical Eilenberger theory of superconductivity to describe quasiparticle excitations in inhomogeneous systems. With the use of low-energy projection matrix, the $M$-band quasiclassical Eilenberger equations are systematically obtained from $N$-band Gor'kov equations. Here $M$ is the internal degrees of freedom in the bands crossing the Fermi energy and $N$ is the degree of freedom in a model. Our framework naturally includes inter-band off-diagonal elements of Green's functions, which have usually been neglected in previous multi-band quasiclassical frameworks. The resultant multi-band Eilenberger and Andreev equations are similar to the single-band ones, except for multi-band effects. The multi-band effects can exhibit the non-locality and the anisotropy in the mapped systems. Our framework can be applied to an arbitrary Hamiltonian (e.g. a tight-binding Hamiltonian derived by the first-principle calculation). As examples, we use our framework in various kinds of systems, such as noncentrosymmetric superconductor CePt$_{3}$Si, three-orbital model for Sr$_{2}$RuO$_{4}$, heavy fermion CeCoIn$_{5}$/YbCoIn$_{5}$ superlattice, a topological superconductor with the strong spin-orbit coupling Cu$_{x}$Bi$_{2}$Se$_{3}$, and a surface system on a topological insulator.

preprint2016arXiv

Nuclear magnetic relaxation rates of unconventional superconductivity in doped topological insulators

We study the temperature dependence of nuclear magnetic relaxation (NMR) rates to detect a sign of topological superconductivity in doped topological insulators, such as $M$($=$Cu,Nb,Sr)$_{x}$Bi$_{2}$Se$_{3}$ and Sn$_{1-x}$In$_{x}$Te. The Hebel-Slichter coherence effect below a critical temperature Tc depends on the superconducting states predicted by a minimal model of doped topological insulators. In a nodal anisotropic topological state similar to the ABM-phase in $^{3}$He, the NMR rate has a conventional $s$-wave like coherence peak below Tc. In contrast, in a fully-gapped isotropic topological superconducting state, this rate below Tc exhibits an anti-peak profile. Moreover, in a two-fold in-plane anisotropic topological superconducting state, there is no coherence effect, which is similar to that in a chiral $p$-wave state. Thus, we reveal that the NMR rates shed light on unconventional superconductivity in doped topological insulators.

preprint2016arXiv

Orbital angular momentum in a topological superconductor with Chern number higher than $1$

We investigate the bulk orbital angular momentum (AM) in a two-dimensional hole-doped topological superconductor (SC) which is composed of a hole-doped semiconductor thin film, a magnetic insulator, and an $s$-wave SC and is characterized by the Chern number $C = -3$. In the topological phase, $L_z/N$ is strongly reduced from the intrinsic value by the non-particle-hole-symmetric edge states as in the corresponding chiral $f$-wave SCs when the spin-orbit interactions (SOIs) are small, while this reduction of $L_z/N$ does not work for the large SOIs. Here $L_z$ and $N$ are the bulk orbital AM and the total number of particles at zero temperature, respectively. As a result, $L_z/N$ is discontinuous or continuous at the topological phase transition depending on the strengths of the SOIs. We also discuss the effects of the edge states by calculating the radial distributions of the orbital AM.

preprint2016arXiv

Robust zero-energy bound states around a pair-density-wave vortex core in locally noncentrosymmetric superconductors

We numerically investigate the electronic structures around a vortex core in a bilayer superconducting system, with s-wave pairing, Rashba spin-orbit coupling and Zeeman magnetic field, with use of the quasiclassical Green's function method. The Bardeen-Cooper-Schrieffer (BCS) phase and the so-called pair-density wave (PDW) phase appear in the temperature-magnetic-field phase diagram in a bulk uniform system [Phys. Rev. B 86, 134514 (2012)]. In the low magnetic field perpendicular to the layers, the zero-energy vortex bound states in the BCS phase are split by the Zeeman magnetic field. On the other hand, the PDW state appears in the high magnetic field, and sign of the order parameter is opposite between the layers. We find that the vortex core suddenly shrinks and the zero-energy bound states appear by increasing the magnetic field through the BCS-PDW transition. We discuss the origin of the change in vortex core structure between the BCS and PDW states by clarifying the relation between the vortex bound states and the bulk energy spectra. In the high magnetic field region, the PDW state and vortex bound states are protected by the spin-orbit coupling. These characteristic behaviors in the PDW state can be observed by scanning tunneling microscopy/spectroscopy.

preprint2015arXiv

First-principle study of antimony doping effects on the iron-based superconductor CaFe(Sb$_{x}$As$_{1-x}$)$_{2}$

We study antimony doping effects in the iron-based superconductor CaFe(Sb$_{x}$As$_{1-x}$)$_{2}$ by using the first-principle calculation. The calculations reveal that the substitution of the doped antimony atom into As of the chainlike As layers is more stable than that in FeAs layers. This prediction can be checked by experiments. Our results suggest that doping homologous elements into the chainlike As layers existing only in novel 112 system is responsible for rising up the critical temperature. We discuss antimony doping effects on the electronic structure. It is found that the calculated band structures with and without the antimony doping are similar to each other within our framework.

preprint2015arXiv

Inverse coherence effects in nuclear magnetic relaxation rates as a sign of topological superconductivity

We reveal that three-dimensional multi-orbital topological superconductivity can be identified by a bulk measurement, i.e., the temperature dependence of nuclear magnetic relaxation (NMR) rates. Below a critical temperature $T_{\rm c}$, the NMR rate in the topological state exhibits an anti-peak profile, which is opposite to the conventional $s$-wave state. This inversion coherence effect comes from a twist of order parameters with respect to orbital and spin degrees of freedom. Our self-consistent calculations in the model for Cu$_{x}$Bi$_{2}$Se$_{3}$ prove that the inverse coherence effect appears as a concave temperature dependence of the NMR rates. We propose that a time-reversal-invariant orbital-singlet spin-triplet topological superconductivity is characterized by the temperature dependence of the NMR rate.

preprint2015arXiv

Non-fragile superconductivity with nodes in the superconducting topological insulator CuxBi2Se3: Zeeman orbital field and non-magnetic impurities

We study the robustness against non-magnetic impurities in the topological superconductor with point nodes, focusing on an effective model of Cu$_{x}$Bi$_{2}$Se$_{3}$. We find that the topological superconductivity with point-nodes is not fragile against non-magnetic impurities, although the superconductivity with nodes in past studies is usually fragile. Exchanging the role of spin with the one of orbital, and vice versa, we find that in the "dual" space the topological superconductor with point-nodes is regarded as the intra-orbital spin-singlet $s$-wave one. From the viewpoint of the dual space, we deduce that the point-node state is not fragile against non-magnetic impurity, when the orbital imbalance in the normal states is small. Since the spin imbalance is induced by the Zeeman magnetic field, we shall name this key quantity for the impurity effects Zeeman "orbital" field. The numerical calculations support that the deduction is correct. If the Zeeman orbital field is small, the topological superconductivity is not fragile in dirty materials, even with nodes. Thus, the topological superconductors can not be simply regarded as one of the conventional unconventional superconductors.

preprint2015arXiv

Orbital angular momentum in a nonchiral topological superconductor

We investigate the bulk orbital angular momentum in a two-dimensional time-reversal broken topological superconductor with the Rashba spin-orbit interaction, the Zeeman interaction, and the $s$-wave pairing potential. Prior to the topological phase transition, we find the crossover from $s$ wave to $p$ wave. For the large spin-orbit interaction, even in the topological phase, $L_z/N$ does not reach $-1/2$, which is the intrinsic value in chiral $p$-wave superconductors. Here $L_z$ and $N$ are the bulk orbital angular momentum and the total number of electrons at zero temperature, respectively. Finally, we discuss the effects of nonmagnetic impurities.

preprint2015arXiv

Surface states around a vortex in topological superconductors: Intersection of a surface and a vortex

We numerically show that the zero-energy Majorana surface states are suppressed around a vortex in the three-dimensional topological superconductors such as Cu$_{x}$Bi$_{2}$Se$_{3}$ and Sn$_{1-x}$In$_{x}$Te. On the other hand, the zero-energy Majorana bound states along the vortex line are robust against cut by the surface. The suppression of the surface bound states is similar to that with a magnetic impurity on the surface of the topological insulator. The suppression of the surface-bound states around a vortex can be observed as the unconventional energy dependence of imaging of the scanning tunneling microscopy/spectroscopy.

preprint2015arXiv

Theory of low-energy behaviors in topological $s$-wave pairing superconductors

We construct a low-energy effective theory of topological $s$-wave pairing superconductors, focusing on the mean-field model of superconductor $\mbox{Cu}_{x}\mbox{Bi}_{2}\mbox{Se}_{3}$. Our approach is second-order perturbation with respect to the inverse of the mass (i.e., large-mass expansion) in the Dirac-type electron dispersion from topological insulator $\mbox{Bi}_{2}\mbox{Se}_{3}$. Since the Dirac-type dispersion with a large mass describes non-relativistic electrons, the large-mass expansion corresponds to a low-energy theory with respect to the original setup. We show that the effective gap function has not only a $p$-wave-like component as the primary contribution, but also an $s$-wave-like one as higher-order corrections. The mixture of $p$- and $s$-wave explains the numerical results [Phys. Rev. B 89 (2014) 214506] of the non-magnetic impurity effects.

preprint2015arXiv

Topological s-wave pairing superconductivity with spatial inhomogeneity: Mid-gap-state appearance and robustness of superconductivity

We study the quasiparticle spectrum of 2D topological $s$-wave superconductors with the Zeeman magnetic field and the Rashba spin-orbit coupling in the presence of spatial inhomogeneity. Solving the real-space Bogoliubov-de Gennes equations, we focus on the excitations within the superconducting gap amplitude, i.e., the appearance of mid-gap states. Two kinds of potential functions, line-type (a chain of impurities) and point-type (a single impurity) ones are examined to take spatial inhomogeneity into account. The line setting shows a link of the mid-gap states with the gapless surface modes indicated by the bulk-boundary correspondence in topological superfluid. The point one shows that the quasiparticles with mid-gap energy are much easily excited by an impurity when the Zeeman magnetic field increases within the topological number to be unchanged. Thus, we obtain insights into the robustness of a topological superconductor against non-magnetic impurities. Moreover, we derive an effective theory applicable to high magnetic fields. The effective gap is the mixture of the chiral $p$-wave and $s$-wave characters. The former is predominant when the magnetic field increases. Therefore, we claim that a chiral $p$-wave character of the effective gap function creates the mid-gap states.

preprint2014arXiv

Field-angle Resolved Flux-flow Resistivity as a Phase-sensitive Probe of Unconventional Cooper Pairing

We theoretically investigate the applied magnetic field-angle dependence of the flux-flow resistivity $ρ_{\rm f}(α_{\rm M})$ for an uniaxially anisotropic Fermi surface. $ρ_{\rm f}$ is related to the quasiparticle scattering rate $\varGamma$ inside a vortex core, which reflects the sign change in the superconducting pair potential. We find that $ρ_{\rm f}(α_{\rm M})$ is sensitive to the sign-change in the pair potential and has its maximum when the magnetic field is parallel to the gap-node direction. We propose the measurement of the field-angle dependent oscillation of $ρ_{\rm f}(α_{\rm M})$ as a phase-sensitive field-angle resolved experiment.

preprint2014arXiv

Field-angle-dependent low-energy excitations around a vortex in the superconducting topological insulator CuxBi2Se3

We study the quasiparticle excitations around a single vortex in the superconducting topological insulator Cu$_{x}$Bi$_{2}$Se$_{3}$, focusing on a superconducting state with point nodes. Inspired by the recent Knight shift measurements, we propose two ways to detect the positions of point nodes, using an explicit formula of the density of states with Kramer-Pesch approximation in the quasiclassical treatment. The zero-energy local density of states around a vortex parallel to the $c$-axis has a twofold shape and splits along the nodal direction with increasing energy; these behaviors can be detected by the scanning tunneling microscopy. An angular dependence of the density of states with a rotating magnetic field on the $a$-$b$ plane has deep minima when the magnetic field is parallel to the directions of point nodes, which can be detected by angular-resolved heat capacity and thermal conductivity measurements. All the theoretical predictions are detectable via standard experimental techniques in magnetic fields.

preprint2014arXiv

Imaging Josephson Vortices on the Surface Superconductor Si(111)-(root7xroot3)-In using a Scanning Tunneling Microscope

We have studied the superconducting Si(111)-(root7xroot3)-In surface using a 3He-based low-temperature scanning tunneling microscope (STM). Zero-bias conductance (ZBC) images taken over a large surface area reveal that vortices are trapped at atomic steps after magnetic fields are applied. The crossover behavior from Pearl to Josephson vortices is clearly identified from their elongated shapes along the steps and significant recovery of superconductivity within the cores. Our numerical calculations combined with experiments clarify that these characteristic features are determined by the relative strength of the interterrace Josephson coupling at the atomic step.

preprint2014arXiv

Impurity effects in a two-dimensional topological superconductor: A link of Tc-robustness with a topological number

Impurity effects are probes for revealing an unconventional property in superconductivity. We study effects of non-magnetic impurities, in a 2D topological superconductor with s-wave pairing, the Rashba spin-orbit coupling, and the Zeeman term. Using a self-consistent T-matrix approach, we calculate a phenomenological formula for the Thouless-Kohmoto-Nightingale-Nijs (TKNN) invariant in interacting systems, as well as density of states, with different magnetic fields. This quantity weakly depends on the magnetic field, when a spectral gap opens, whereas this changes drastically, when in-gap states occurs. Furthermore, in the latter case, we find that the Anderson's theorem (robustness of s-wave superconductivity against non-magnetic impurities) breaks down. We discuss the origin, from the viewpoints of both unconventional superconductivity and the TKNN invariant.

preprint2014arXiv

Non-magnetic impurity effects in a three-dimensional topological superconductor: From p- to s-wave behaviors

Unconventional features in superconductivity are revealed by responses to impurity scattering. We study non-magnetic impurity effects in a three-dimensional topological superconductor, focusing on an effective model (massive Dirac Hamiltonian with s-wave on-site pairing) of Copper-doped bismuth-selenium compounds. Using a self-consistent T-matrix approach for impurity scattering, we examine in-gap states in density of states. We find that the results are well characterized by a single material variable, which measures relativistic effects in the Dirac Hamiltonian. In non-relativistic regime, an odd-parity superconducting state is fragile against non-magnetic impurities. We show that this behavior is caused by a p-wave character involved in the topological superconducting state. In contrast, we show that in relativistic regime the superconductivity is robust against non-magnetic impurities, owing to an s-wave character. To summarize, the system has two aspects, p- and s-wave features, depending on the weight of relativistic effects.

preprint2014arXiv

Quasiclassical Treatment and Odd-parity/Triplet Correspondence in Topological Superconductors

We construct a quasiclassical framework for topological superconductors with the strong spin-orbit coupling such as CuxBi2Se3. In the manner of the quasiclassical treatment, decomposing the slowly varying component from the total quasi-particle wave function, the original massive Dirac Bogoliubov-de Gennes (BdG) Hamiltonian derived from the tight-binding model represented by 8 x 8 matrix is reduced to 4 x 4 one. The resultant equations are equivalent to Andreev-type equations of singlet or triplet superconductors, in which the apparent spin-orbit coupling vanishes. Using this formalism, we find a fact that the odd-parity superconductivity in topological superconductors turns to the spin-triplet one. % without the spin-orbit coupling through the quasiclassical treatment. Moreover, in terms of the quasiclassical treatment, we show that the topologically-protected zero-energy states in topological superconductors has the correspondence to the Andreev bound states established in a long history of studies for the unconventional superconductors. This clearly indicates that low-energy non-trivial superconducting properties in the topological superconductors can be analyzed by the established theoretical descriptions on the spin-triplet superconductors.

preprint2014arXiv

Spin-polarized Majorana Bound States inside a Vortex Core in Topological Superconductors

We reveal that Majorana bound states inside the vortex core in an odd-parity topological superconductivity classified as "pseudo-scalar" type in the gap function are distinctly spin-polarized by solving the massive Dirac Bogoliubov-de Gennes (BdG) equation considering the spin-orbit coupling. This result is universal for "Dirac superconductivity" whose rotational degree of freedom is characterized by the total angular momentum J = S + L and in marked contrast to the spin-degeneracy of the core bound states as the consequence of the conventional BdG equation. The spin-polarized vortex core can be easily detected by spin-sensitive probes such as the neutron scattering and other measurements well above the first critical magnetic field H_{c1}.

preprint2014arXiv

Vortex Core Structure in Multilayered Rashba Superconductors

We numerically study the electronic structure of a single vortex in two dimensional superconducting bilayer systems within the range of the mean-field theory. The lack of local inversion symmetry in the system is taken into account through the layer dependent Rashba spin-orbit coupling. The spatial profiles of the pair potential and the local quasiparticle density of states are calculated in the clean spin-singlet superconductor on the basis of the quasiclassical theory. In particular, we discuss the characteristic core structure in the pair-density wave state, which is spatially modulated exotic superconducting phase in a high magnetic field.

preprint2013arXiv

Effect of anisotropic Fermi surface on the flux-flow resistivity under rotating magnetic field

We numerically investigate the effect of in-plane anisotropic Fermi surface (FS) on the flux-flow resistivity $ρ_{\rm f}$ under rotating magnetic field on the basis of the quasiclassical Green's function method. We demonstrate that one can detect the phase in pairing potential of Cooper pair through the field-angular dependence of $ρ_{\rm f}$ even if the FS has in-plane anisotropy. In addition, we point out one can detect the gap-node directions irrespective of the FS anisotropy by measuring $ρ_{\rm f}$ under rotating field.

preprint2013arXiv

Effect of uniaxially anisotropic Fermi surface on the quasiparticle scattering inside a vortex core in unconventional superconductors

We theoretically study the dependence of the quasiparticle (QP) scattering rate $\varGamma$ on the uniaxial anisotropy of a Fermi surface with changing the magnetic field angle $α_{\rm M}$. We consider the QP scattering due to the non-magnetic impurities inside a single vortex core. The field-angle dependence of the quasiparticle scattering rate $\varGamma(α_{\rm M})$ is sensitive to the sign-change of the pair potential. We show that with increasing the two dimensionality of the system, $\varGamma(α_{\rm M})$ reflects more clearly whether there is the sign-change in the pair potential.

preprint2013arXiv

Field-angle dependence of the quasiparticle scattering inside a vortex core in unconventional superconductors

We theoretically investigate the quasiparticle scattering rate $\varGamma$ inside a vortex core in the existence of non-magnetic impurities distributed randomly in a superconductor. We show that the dependence of $\varGamma$ on the magnetic field direction is sensitive to the sign of the pair potential. The behavior of $\varGamma$ is quite different between an s-wave and a d-wave pair potential, where these are assumed to have the same amplitude anisotropy, but a sign change only for the d-wave one. It is suggested that measurements of the microwave surface impedance with changing applied-field directions would be used for the phase-sensitive identification of pairing symmetry.

preprint2013arXiv

Impurity Effect on the Local Density of States around a Vortex in Noncentrosymmetric Superconductors

We numerically study the effect of non-magnetic impurities on the vortex bound states in noncentrosymmetric systems. The local density of states (LDOS) around a vortex is calculated by means of the quasiclassical Green's function method. We find that the zero energy peak of the LDOS splits off with increasing the impurity scattering rate.

preprint2013arXiv

Impurity scattering effect on the zero-energy peak of the local density of states in a multi-quantum vortex core

We theoretically study a non-magnetic impurity effect on the vortex bound states of a multi-quantum vortex. The zero-energy peak of the local density of states is investigated for vortex cores with the winding numbers 2 and 4 within the framework of the quasiclassical theory of superconductivity. We find that the zero-energy peaks, which appear away from the vortex center in the clean limit, move towards the vortex center with increasing the impurity scattering rate, resolving a contradiction between an experimental result and previous theoretical predictions.

preprint2013arXiv

Inhomogeneity Effects in Topological Superconductors

We have constructed a quasiclassical framework on superconductors with strong spin-orbit couplings, applicable to CuxBi2Se3[Y. Nagai, H. Nakamura, and M. Machida: arXiv:1305.3025]. The notable point is that in this framework the Bogoliubov-de Gennes Hamiltonians with suggested odd-parity pairing states turn to quasiclassical ones with usual spin-triplet Cooper pairs. Using this quasiclassical theory, we can investigate inhomogeneity effects such as the phenomena with vortices and surfaces in this superconductors and shed light on the pairing state of topological superconductors. In this paper, we apply the quasiclassical framework to the surface bound states with the Dirac-cone energy dispersion originated from the topological invariant in the parent compound Bi2Se3 in order to investigate the robustness of these bound states under the superconducting order parameter. The odd-parity gap functions can not open on the Dirac-cone-dispersion band in the Cu-doped Bi$_{2}$Se$_{3}$ superconductor. We show that the massless Dirac quasiparticles originated from the normal-state topological invariant and the Majorana quasiparticles coexist with each other on the surface in the odd-parity topological superconductivity. Inhomogeneity effects can be easily investigated with the use of our quasiclassical framework in topological superconductors.

preprint2013arXiv

Numerical construction of a low-energy effective Hamiltonian in a self-consistent Bogoliubov-de Gennes approach of superconductivity

We propose a fast and efficient approach for solving the Bogoliubov-de Gennes (BdG) equations in superconductivity, with a numerical matrix-size reduction procedure proposed by Sakurai and Sugiura [J. Comput. Appl. Math. 159, 119 (2003)]. The resultant small-size Hamiltonian contains the information of the original BdG Hamiltonian in a given energy domain. In other words, the present approach leads to a numerical construction of a low-energy effective theory in superconductivity. The combination with the polynomial expansion method allows a self-consistent calculation of the BdG equations. Through numerical calculations of quasi-particle excitations in a vortex lattice, thermal conductivity, and nuclear magnetic relaxation rate, we show that our approach is suitable for evaluating physical quantities in a large-size superconductor and a nano-scale superconducting device, with the mean-field superconducting theory.

preprint2013arXiv

Phase-Sensitive Flux-Flow resistivity in Unconventional Superconductors

We theoretically investigate the magnetic-field-angle dependence of the flux-flow resistivity $ρ_{\rm f}$ in unconventional superconductors. Two contributions to $ρ_{\rm f}$ are considered: one is the quasiparticle (QP) relaxation time $τ(\bm{k}_{\rm F})$ and the other is $ω_0(\bm{k}_{\rm F})$, which is a counterpart to the interlevel spacing of the QP bound states in the quasiclassical approach. Here, $\bm{k}_{\rm F}$ denotes the position on a Fermi surface. Numerical calculations are conducted for a line-node s-wave and a d-wave pair potential with the same anisotropy of their amplitudes, but with a sign change only for a d-wave one. We show that the field-angle dependence of $ρ_{\rm f}$ differs prominently between s-wave and d-wave pairs, reflecting the phase of the pair potentials. We also discuss the case where $τ$ is constant and compare it with the more general case where $τ$ depends on $\bm{k}_{\rm F}$.

preprint2012arXiv

Direct Numerical Demonstration of Sign-preserving Quasiparticle Interference via Impurity inside Vortex Core in Unconventional Superconductors

We perform large-scale numerical calculations self-consistently solving the Bogoliubov-de Gennes (BdG) equations in the magnetic field together with random impurities to directly demonstrate the typical quasi-particle interference (QPI) in the presence of vortices as observed by scanning tunneling microscopy/spectroscopy experiments in unconventional superconductors. The calculations reveal that vortex itself never works as a scatter causing the QPI pattern but vortex core containing impurity brings about the enhancement of the sign-preserving QPI peaks. Its origin is Andreev bound-states distorted by impurity, and all the measurement findings are consistently explained by the scenario based on the numerical results.

preprint2012arXiv

Efficient Numerical Self-consistent Mean-field Approach for Fermionic Many-body Systems by Polynomial Expansion on Spectral Density

We propose an efficient numerical algorithm to solve Bogoliubov de Gennes equations self-consistently for inhomogeneous superconducting systems with a reformulated polynomial expansion scheme. This proposed method is applied to typical issues such as a vortex under randomly distributed impurities and a normal conducting junction sandwiched between superconductors. With various technical remarks, we show that its efficiency becomes remarkable in large-scale parallel performance.

preprint2012arXiv

Q-scan-analysis of the neutron scattering in iron-based superconductors

We propose a way to determine the pairing state of the iron pnictide superconductors exploiting the momentum (Q) scan of the neutron scattering data. We investigate the spin susceptibility in the s+- and s++ superconducting states for various doping levels using the effective five-orbital model and considering the quasiparticle damping. The peak position of the intensity shifts from the position on the line Qx = pi to that on the line Qy = 0 as the doping level is decreased from electron doping to hole doping. We find that the Q-dependence of the ratio of the intensity in the superconducting state to that in the normal state is qualitatively different between the s+- -wave and s++ -wave pairings. We propose to investigate experimentally this ratio in Q-space to distinguish the two pairing states.

preprint2012arXiv

Quasiclassical numerical method for mesoscopic superconductors: bound states in a circular d-wave island with a single vortex

We demonstrate an efficient numerical method for obtaining unique solutions to the Eilenberger equation for a mesoscopic or nanoscale superconductor. In particular, we calculate the local density of states of a circular d-wave island containing a single vortex. The "vortex shadow" effect is found to strongly depend on the quasiparticle energy in such small systems. We show how to construct by geometry quasiparticle trajectories confined in a finite-size system with specular reflections at the boundary, and discuss the stability of the numerical solutions even in the case of vanishing order parameter as for nodal quasiparticles in a d-wave superconductor, or for quasiparticles passing through the vortex center with zero energy.

preprint2012arXiv

Rotational Isotropy Breaking as Proof for Spin-polarized Cooper Pairs in the Topological Superconductor CuxBi2Se3

In a promising candidate of topological superconductors, CuxBi2Se3, we propose a way to exclusively determine the pairing symmetry. The proposal suggests that the angle dependence of the thermal conductivity in the basal ab-plane shows a distinct strong anisotropy only when the pairing symmetry is an odd-parity spin-polarized triplet below the superconducting transition temperature (Tc). Such striking isotropy breaking below Tc is explicitly involved in Dirac formalism for superconductors, in which the spin-orbit coupling is essential. We classify possible gap functions based on the Dirac formalism and clarify an origin of the isotropy breaking.

preprint2011arXiv

Comment on arXiv:1105.6233 entitled "Neutron-Inelastic-Scattering Peak by Dissipationless Mechanism in the s++-wave State in Iron-based Superconductors" by S. Onari and H. Kontani

Recently, Onari and Kontani submitted a paper [arXiv:1105.6233] which criticizes our recent theoretical study [arXiv:1103.0586] on the neutron scattering experiment as a probe for determining the superconducting gap in the iron pnictides. In their paper, Onari and Kontani have developed a formalism in which the imaginary part of the dynamical spin susceptibility (Im$χ_s(\Vec{q},ω)$) in the superconducting state can be more accurately calculated especially in the $ω<2Δ$ regime, where $Δ$ is the superconducting gap. In section IIIC of their paper, they mention that the conclusions of our paper are "incorrect based on inaccurate numerical calculation". In the present Comment, we show that this in fact is not correct.

preprint2011arXiv

Determination of the pairing state in iron-based superconductors through neutron scattering

We calculate the spin susceptibility in the s_{+-} and s_{++} superconducting states of the iron pnictides using the effective five orbital model and considering the quasiparticle damping. For the experimentally evaluated magnitude of the quasiparticle damping and the superconducting gap, the results at the wave vector ~ (pi,0) show that the s_{+-} state is more consistent with the neutron scattering experiments, while for larger quasiparticle damping and the superconducting gap, the s_{++} state can be more consistent. To distinguish between two cases that reproduce the experiments at the wave vector ~ (pi,0), we propose to investigate experimentally the wave vector ~ (pi,pi).

preprint2011arXiv

Superconducting gap function in the organic superconductor (TMTSF)2ClO4 with anion ordering; First-principles calculations and quasiclassical analysis for angle-resolved heat capacity

We calculate angle-dependent heat capacity in a low magnetic field range on the basis of Kramer-Pesch approximation together with an electronic structure obtained by first-principles calculations to determine a superconducting gap function of (TMTSF)2ClO4 through its comparisons with experiments. The present comparative studies reveal that a nodal d-wave gap function consistently explains the experimental results for (TMTSF)2ClO4. Especially, it is emphasized that the observed unusual axis-asymmetry of the angle-dependence eliminates the possibility of s-wave and node-less d-wave functions. It is also found that the directional ordering of ClO4 anions does not have any significant effects on the Fermi surface structure contrary to the previous modelings since the two Fermi surfaces obtained by the band calculations almost cross within the present full accuracy in first-principles calculations.

preprint2010arXiv

Impurity effect in multi-orbital, sign-reversing s-wave superconductors

We study the impurity effects on the transition temperature Tc with use of the T-matrix approximation. We propose a way to visualize the multi-orbital effect by introducing the hybridization function characterizing the multi-orbital effect for the impurity scattering. Characterizing function does not depend on the superconducting pairing symmetry, since this function is defined by the eigenvectors in normal states. The result indicates that an impurity-robust superconductivity does not necessarily imply a sign- preserving pairing. Visualizing the hybridization effect in the effective five-band model for LaFeAsO, we show that the impurity effect on Tc is relatively weaker than that in single-band models.

preprint2010arXiv

Impurity Scattering Rate and Coherence factor in Vortex Core of Sign-reversing s-wave Superconductors

We investigate the impurity scattering rates for quasi-particles in vortex cores of sign-reversing s-wave superconductors as a probe to detect the internal phase difference of the order parameters among different Fermi surfaces. The impurity scattering rates and coherence factors are related to quasiparticle interference effect by the scanning tunneling microscopy and spectroscopy technique. With use of the Born and Kramer-Pesch approximations for the Andreev bound states, we show that the sign-reversed forward scatterings are dominant in vortex cores. Owing to the coherence factor in vortex cores of \pm s-wave superconductors, the impurity scattering rate of the Andreev bound states has a characteristic distribution on the Fermi surfaces. For comparison, the impurity scattering rates in vortex cores of s-wave and d-wave superconductors are also discussed.

preprint2010arXiv

Microwave Surface Impedance Measurements of LiFeAs Single Crystals

We report results of microwave surface impedance measurements of LiFeAs single crystals. The in-plane penetration depth depends on temperature exponentially at low temperatures, which strongly suggests that this material has the nodeless superconducting gap. The temperature dependence of the superfluid density indicates that LiFeAs is a multi-gap superconductor with at least two isotropic gaps. In addtion, the real part of the microwave conductivity exhibits a large enhancement below $T_\mathrm{c}$, indicating that the quasi-particle relaxation time, $τ$, increases rapidly below $T_\mathrm{c}$. We believe that this enhancement is rather common to all superconductors where an inelastic scattering is dominant above $T_\mathrm{c}$, irrespective of the strength of the electron correlation.

preprint2009arXiv

Surface-angle dependence of the tunneling spectroscopy in iron-based superconductors: sign-reversing s-wave scenarios

We discuss the surface Andreev bound states in Fe-based superconductors with the use of an effective five-band model and investigate the surface-angle dependence of the tunneling spectroscopy by a quasiclassical approach for an isotropic and an anisotropic /pm s-wave gap superconductivity. We show that information on the normal state is important for the Andreev bound state and its peak positions do not depend on the gap amplitude anisotropy.

preprint2009arXiv

Uniform current in graphene strip with zigzag edges

Graphene exhibits zero-gap massless-Dirac fermion and zero density of states at E = 0. These particles form localized states called edge states on finite width strip with zigzag edges at E = 0. Naively thinking, one may expect that current is also concentrated at the edge, but Zarbo and Nikolic numerically obtained a result that the current density shows maximum at the center of the strip. We derive a rigorous relation for the current density, and clarify the reason why the current density of edge state has a maximum at the center.