Source author record

Masahiko Machida

Masahiko Machida 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

46works
11topics
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

46 published item(s)

preprint2016arXiv

A screened automated structural search with semiempirical methods

We developed an interface program between a program suite for an automated search of chemical reaction pathways, GRRM, and a program package of semiempirical methods, MOPAC. A two-step structural search is proposed as an application of this interface program. A screening test is first performed by semiempirical calculations. Subsequently, a reoptimization procedure is done by ab initio or density functional calculations. We apply this approach to ion adsorption on cellulose. The computational efficiency is also shown for a GRRM search. The interface program is suitable for the structural search of large molecular systems for which semiempirical methods are applicable.

preprint2016arXiv

Emergence of $η$-pairing ground-state in population-imbalanced attractive Fermi-gases filling $p$ orbitals on 1-D optical lattice

We explore the ground states in population-imbalanced attractive 1-D fermionic optical lattice filling $p$ orbitals over the lowest $s$ one by using the density-matrix-renormalization-group (DMRG) method. The DMRG calculations find the occurrence of spatially non-uniform off-diagonal long-range order. In contrast to Fulde-Ferrel Larkin-Ovchinikov pair as observed in the single-band Hubbard model. The spatial oscillation period of the pair correlation function is widely fixed to be $π$ irrespective of the mismatch between spin-split Fermi surfaces. The ground-state $π$ order corresponds to $η$-pair condensate predicted by Yang [Phys. Rev. Lett. \textbf{63}, 2144 (1989)]. Taking account of the effects of harmonic traps, we confirm that the $η$-pair state distinctly emerges at the center of the trap potential surrounded by perfectly-polarized states even in the trapped cases.

preprint2016arXiv

Superconductivity in repulsively interacting fermions on a diamond chain: flat-band induced pairing

To explore whether a flat-band system can accommodate superconductivity, we consider repulsively interacting fermions on the diamond chain, a simplest quasi-one-dimensional system that contains a flat band. Exact diagonalization and the density-matrix renormalization group (DMRG) are used to show that we have a significant binding energy of a Cooper pair with a long-tailed pair-pair correlation in real space when the total band filling is slightly below $1/3$, where the dispersive band interacts with the flat band that is empty but close to $E_F$. Pairs selectively formed across the outer sites of the diamond chain are responsible for the pairing correlation. At exactly $1/3$-filling an insulating phase emerges, where the entanglement spectrum indicates the particles on the outer sites are highly entangled and topological. These come from a peculiarity of the flat band in which "Wannier orbits" are not orthogonalizable.

preprint2015arXiv

Evaluation of ambient dose equivalent rates influenced by vertical and horizontal distribution of radioactive cesium in soil in Fukushima Prefecture

The air dose rate in an environment contaminated with 134Cs and 137Cs depends on the amount, depth profile and horizontal distribution of these contaminants within the ground. This paper introduces and verifies a tool that models these variables and calculates ambient dose equivalent rates at 1 m above the ground. Good correlation is found between predicted dose rates and dose rates measured with survey meters in Fukushima Prefecture in areas contaminated with radiocesium from the Fukushima Dai-ichi Nuclear Power Plant accident. This finding is insensitive to the choice for modelling the activity depth distribution in the ground using activity measurements of collected soil layers, or by using exponential and hyperbolic secant fits to the measurement data. Better predictions are obtained by modelling the horizontal distribution of radioactive cesium across an area if multiple soil samples are available, as opposed to assuming a spatially homogeneous contamination distribution. Reductions seen in air dose rates above flat, undisturbed fields in Fukushima Prefecture are consistent with decrement by radioactive decay and downward migration of cesium into soil. Analysis of remediation strategies for farmland soils confirmed that topsoil removal and interchanging a topsoil layer with a subsoil layer result in similar reductions in the air dose rate. These two strategies are more effective than reverse tillage to invert and mix the topsoil.

preprint2015arXiv

Fields of View for Environmental Radioactivity

The gamma component of air radiation dose rates is a function of the amount and spread of radioactive nuclides in the environment. These radionuclides can be natural or anthropogenic in origin. The field of view describes the area of radionuclides on, or below, the ground that is responsible for determining the air dose rate, and hence correspondingly the external radiation exposure. This work describes Monte Carlo radiation transport calculations for the field of view under a variety of situations. Presented first are results for natural 40K and thorium and uranium series radionuclides distributed homogeneously within the ground. Results are then described for atmospheric radioactive caesium fallout, such as from the Fukushima Daiichi Nuclear Power Plant accident. Various stages of fallout evolution are considered through the depth distribution of 134Cs and 137Cs in soil. The fields of view for the natural radionuclides and radiocaesium are different. This can affect the responses of radiation monitors to these nuclides if the detector is partially shielded from the ground within its field of view. The field of view also sets the maximum reduction in air dose rates that can be achieved through local decontamination or remediation measures. This maximum efficiency can be determined quickly from the data presented here for the air dose rate versus the spatial extent of radioactive source on the ground.

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

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

Topographic Effects on Ambient Dose Equivalent Rates from Radiocesium Fallout

Land topography can affect air radiation dose rates by locating radiation sources closer to, or further from, detector locations when compared to perfectly flat terrain. Hills and slopes can also shield against the propagation of gamma rays. To understand the possible magnitude of topographic effects on air dose rates, this study presents calculations for ambient dose equivalent rates at a range of heights above the ground for varying land topographies. The geometries considered were angled ground at the intersection of two planar surfaces, which is a model for slopes neighboring flat land, and a simple conical geometry, representing settings from hilltops to valley bottoms. In each case the radiation source was radioactive cesium fallout, and the slope angle was varied systematically to determine the effect of topography on the air dose rate. Under the assumption of homogeneous fallout across the land surface, and for these geometries and detector locations, the dose rates at high altitudes are more strongly affected by the underlying land topography than those close to ground level. At a height of 300m, uneven topographies can lead to a 50% change in air dose rates compared to if the ground were uniformly flat. However, in practice the effect will more often than not be smaller than this, and heterogeneity in the source distribution is likely to be a more significant factor in determining local air dose rates.

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

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

Quantum phases in $p$-orbital degenerated attractive 1D fermionic optical lattices

We examine quantum phases emerged by double degeneracy of $p$-orbital bands in attractive atomic Fermi gases loaded on a 1D optical lattice. Our numerical simulations by the density-matrix renormalization group predict the emergence of a state with a charge excitation gap, the Haldane insulator phase. A mapping onto an effective spin-$1$ model reveals its physical origin. Moreover, we show that population imbalance leads to richer diversity of the quantum phases, including a phase-separated polarized state. Finally, we study the effects of harmonic trap potential in this 1D chain.

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

Theory of macroscopic quantum tunneling with Josephson-Leggett collective excitations in multi-band superconducting Josephson junctions

Collective excitations reveal fundamental properties and potential applications of superconducting states. We theoretically study macroscopic quantum tunneling (MQT) in a Josephson junction composed of multi-band superconductors, focusing on a phase mode induced by inter-band fluctuations: the Josephson-Leggett (JL) collective excitation mode. Using the imaginary-time path-integral method, we derive a formula for the MQT escape rate for low-temperature switching events. We clarify that the JL mode has two major effects on the MQT: (i) the zero-point fluctuations enhance the escape rate, and (ii) the quantum dissipation induced by the couplings to the gauge-invariant phase difference suppresses the quantum tunneling. We show that the enhancement exceeds the suppression for a wide range of junction parameters. This enhancement originates from the single-mode interaction between the tunneling variable and the inter-band fluctuations.

preprint2014arXiv

Two-dimensional macroscopic quantum tunneling in multi-gap superconductor Josephson junctions

Low-temperature characters of superconducting devices yield definite probes for different superconducting phenomena. We study the macroscopic quantum tunneling (MQT) in a Josephson junction, composed of a single-gap superconductor and a two-gap superconductor. Since this junction has two kinds to the superconducting phase differences, calculating the MQT escape rate requires the analysis of quantum tunneling in a multi-dimensional configuration space. Our approach is the semi-classical approximation along a 1D curve in a 2D potential- energy landscape, connecting two adjacent potential (local) minimums through a saddle point. We find that this system has two plausible tunneling paths; an in-phase path and an out-of-phase path. The former is characterized by the Josephson-plasma frequency, whereas the latter is by the frequency of the characteristic collective mode in a two-band superconductor, Josephson- Leggett mode. Depending on external bias current and inter-band Josephson-coupling energy, one of them mainly contributes to the MQT. Our numerical calculations show that the difference between the in-phase path and the out-of-phase path is manifest, with respect to the bias- current-dependence of the MQT escape rate. This result suggests that our MQT setting be an indicator of the Josephson-Leggett mode.

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

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

Massless collective excitations in frustrated multi-band superconductors

We study collective excitations in three- and four-band superconductors with inter-band frustration, which causes neither 0 nor $π$ inter-band phases in the superconducting state. Using a low-energy spin-Hamiltonian originating from a multi-band tight-binding model, we find that a massless Leggett mode occurs in a wide parameter region of this four-band system. This massless mode is related to the fact that the mean-field energy does not depend on a continuous superconducting phase. In other words, we find a link of the massless mode with a degeneracy between a time-reversal-symmetry-breaking state (neither 0 nor $π$ phases) and a time-reversal-symmetric state (either 0 or $π$ phases). Therefore, the mass of this collective mode characterizes well the time-reversal symmetry in frustrated multi-band superconductors.

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}$.

preprint2013arXiv

Phonons and Spin Excitations in Fe-Based Superconductor Ca10Pt4As8(Fe1-xPtxAs)10 (x ~ 0.2)

By means of neutron inelastic scattering, magnetic excitations and phonons were measured for a single crystal of slightly overdoped superconductor Ca10Pt4As8(Fe1-xPtxAs)10 (x ~ 0.2) with the transition temperature Tc of ~33 K. Below Tc, magnetic excitation spectra \{chi}"(Q, ω) measured at Q = QM [magnetic Γ points] are gapped, and in the relatively higher ω region, the \{chi}"(QM, ω)-increase was observed with decreasing T, where the maximum of the increase was found at ~18 meV at 3 K (<< Tc). These characteristics are favorable to orbital-fluctuation-mediated superconductivity with the so-called S++symmetry of the order parameter. The energy dependence of δQ-temperature(T) curve of the magnetic excitations seems to have anomalous behavior in rather wide T region above Tc, δQ being the width of the Q scan profile. In the phonon measurements, we observed softening of the in-plane TA mode, which corresponds to the elastic constant C66. This softening seems to start at rather high temperature T, as T is lowered. Additionally, anomalous increase in spectral weights of the TO-phonons at around QM in the region 35 < ω < 40 meV was found even above Tc, as T is lowered from ambient T . Because the spectral weights in this ω region mainly correspond to the in-plane motions of Fe atoms and because orbital fluctuations are expected to be strong at around QM, the result may present clues to investigate a possible coupling between the fluctuations of the orbitals and lattice system.

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

Nontrivial Haldane phase of an atomic two-component Fermi gas trapped in a 1d optical lattice

We propose how to create a non-trivial Haldane phase in atomic two-component Fermi-gas loaded on one-dimensional (1-D) optical lattice with trap potential. The Haldane phase is naturally formed on $p$-band Mott core in a wide range of the strong on-site repulsive interaction. The present proposal is composed of two steps, one of which is theoretical derivation of an effective 1-D S=1 interacting-chain model from the original tight-binding Hamiltonian handling the two $p$-orbitals, and the other of which is numerical demonstration employing the density-matrix renormalization-group for the formation of the Haldane phase on $p$-band Mott core and its associated features in the original tight-binding model with the harmonic trap potential.

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

Collective modes in multiband superfluids and superconductors: Multiple dynamical classes

One important way to characterize the states having a gauge symmetry spontaneously broken over multibands should be to look at their collective excitation modes. We find that a three-band system has multiple Leggett modes with significantly different masses, which can be classified into different dynamical classes according to whether multiple inter-band Josephson currents add or cancel. This provides a way to dynamically characterize multiband superconductivity while the pairing symmetry is a static property.

preprint2011arXiv

Energy Spectra of Quantum Turbulence: Large-scale Simulation and Modeling

In $2048^3$ simulation of quantum turbulence within the Gross-Pitaevskii equation we demonstrate that the large scale motions have a classical Kolmogorov-1941 energy spectrum E(k) ~ k^{-5/3}, followed by an energy accumulation with E(k) ~ const at k about the reciprocal mean intervortex distance. This behavior was predicted by the L'vov-Nazarenko-Rudenko bottleneck model of gradual eddy-wave crossover [J. Low Temp. Phys. 153, 140-161 (2008)], further developed in the paper.

preprint2011arXiv

First-principles Studies for the Hydrogen Doping Effects on Iron-based Superconductors

We study hydrogen doping effects in an iron-based superconductor LaFeAsO_(1-y) by using the first-principles calculation and explore the reason why the superconducting transition temperature is remarkably enhanced by the hydrogen doping. The present calculations reveal that a hydrogen cation stably locating close to an iron atom attracts a negatively-charged FeAs layer and results in structural distortion favorable for further high temperature transition. In fact, the lattice constant a averaged over the employed supercell shrinks and then the averaged As-Fe-As angle approaches 109.74 degrees with increasing the hydrogen doping amount. Moreover, the calculations clarify electron doping effects of the solute hydrogen and resultant Fermi-level shift. These insights are useful for design of high transition-temperature iron-based superconductors.

preprint2011arXiv

Ground-state properties of the one-dimensional attractive Hubbard model with confinement: a comparative study

We revisit the one-dimensional attractive Hubbard model by using the Bethe-ansatz based density-functional theory and density-matrix renormalization method. The ground-state properties of this model are discussed in details for different fillings and different confining conditions in weak-to-intermediate coupling regime. We investigate the ground-state energy, energy gap, and pair-binding energy and compare them with those calculated from the canonical Bardeen-Cooper-Schrieffer approximation. We find that the Bethe-ansatz based density-functional theory is computationally easy and yields an accurate description of the ground-state properties for weak-to-intermediate interaction strength, different fillings, and confinements. In order to characterize the quantum phase transition in the presence of a harmonic confinement, we calculate the thermodynamic stiffness, the density-functional fidelity, and fidelity susceptibility, respectively. It is shown that with the increase of the number of particles or attractive interaction strength, the system can be driven from the Luther-Emery-type phase to the composite phase of Luther-Emery-like in the wings and insulating-like in the center.

preprint2011arXiv

Macroscopic quantum tunneling in multigap superconducting Josephson junctions: Escape rate enhancement via quantum fluctuations of Josephson-Leggett mode

We theoretically study the macroscopic quantum tunneling (MQT) in a hetero Josephson junction formed by a conventional single-gap superconductor and a multi-gap one such as ${MgB}_{2}$ and iron-based superconductors. In such a Josephson junction more than one phase differences are defined. We clarify the quantum dynamics of the phase differences and construct a theory for the MQT in the multi-gap Josephson junctions. The dynamics of the phase differences are strongly affected by the Josephson-Leggett mode, which is the out-of-phase oscillation mode of the phase differences. The escape rate is calculated in terms of the effective action renormalized by the Josephson-Leggett mode. It is shown that the escape rate is drastically enhanced when the frequency of the Josephson-Leggett mode is less than the Josephson-plasma frequency.

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.

preprint2011arXiv

Variety of c-axis collective excitations in layered multigap superconductors

We present a dynamical theory for the phase differences along a stacked direction of intrinsic Josephson junctions (IJJ's) in layered multigap superconductors, motivated by the discovery of highly-anisotropic iron-based superconductors with thick perovskite-type blocking layers. The dynamical equations describing AC and DC intrinsic Josephson effects peculiar to multigap IJJ's are derived, and collective Leggett mode excitations in addition to the Josephson plasma established in single-gap IJJ's are predicted. The dispersion relations of their collective modes are explicitly displayed, and the remarkable peculiarity of the Leggett mode is demonstrated.

preprint2010arXiv

Ambegaokar-Baratoff relations of Josephson critical current in heterojunctions with multi-gap superconductors

An extension of the Ambegaokar-Baratoff relation to a superconductor-insulator-superconductor (SIS) Josephson junction with multiple tunneling channels is derived. Appling the resultant relation to a SIS Josephson junction formed by an iron-based (five-band) and a single-band Bardeen-Cooper-Schrieffer (BCS) type superconductors, a theoretical bound of the Josephson critical current ($I_{\rm c}$) multiplied by the resistance of the junction ($R_{\rm n}$) is given. We reveal that such a bound is useful for identifying the pairing symmetry of iron-pnictide superconductors. One finds that if a measured value of $I_{\rm c}R_{\rm n}$ is smaller than the bound then the symmetry is $\pm s$-wave, and otherwise $s$-wave without any sign changes. In addition, we stress that temperature dependence of $I_{\rm c}R_{\rm n}$ is sensitive to the difference of the gap functions from the BCS type gap formula in the above heterojunction.

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

Magnetic Ordering in Blocking Layer and Highly Anisotropic Electronic Structure of High-Tc Iron-based Superconductor Sr2VFeAsO3: LDA+U Studies

We calculate electronic structures of a high-Tc iron-based superconductor Sr2VFeAsO3 by LDA+U method. We assume a checker-board antiferromagnetic order on blocking layers including vanadium and strong correlation in d-orbits of vanadium through the Hubbard U. While the standard LDA brings about metallic blocking layers and complicated Fermi surface as in the previous literatures, our calculation changes the blocking layer into insulating one and the Fermi surface becomes quite similar to those of other iron-based superconductors. Moreover, the appearance of the insulating blocking layers predicts high anisotropy on quasi-particle transports and new types of intrinsic Josephson effects.

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.

preprint2010arXiv

Shapiro Steps as a Direct Probe of $\pm s$-wave Symmetry in Multi-gap Superconducting Josephson Junctions

We theoretically study the Shapiro steps in a hetero Josephson junction made of a single-gap superconductor and a two-gap one. We find that an anomalous DC Josephson current is induced by tuning the frequency of an applied microwave to the Josephson-Leggett mode frequency, which creates an extra step structure in the $I$-$V$ characteristics besides the conventional Shapiro steps. The step heights at the resonance voltages exhibit an alternate structure of a large and small value reflecting the gap symmetry of the two-gap superconductor. In the $\pm s-wave case in which the two gaps have opposite signs in the two-gap superconductor the steps with odd index are enhanced, whereas in the s-wave case the ones with even index have larger values. The existence of the fractional Shapiro steps is also predicted.

preprint2010arXiv

Theory of vortex structure in Josephson junctions with multiple tunneling channels: Vortex enlargement as a probe of $\pm s$-wave superconductivity

We theoretically study Josephson vortex structures in Josephson junctions which have multiple tunneling channels caused by multiple superconducting gaps. Deriving "coupled sine-Gordon equations" from the free-energy taking account of the multiple tunneling channels, we examine two typical situations, a heterotic junction composed of multigap-superconductor, insulator, and single-gap superconductor, and a grain-boundary junction formed by two identical multigap superconductors. Then, we reveal in both situations that the magnetic field distribution of the Josephson vortex for $\pm s$-wave superconductivity is more enlarged than that for s-wave without sign change between the order parameters. Its mechanism is ascribed to a cancellation of the multiple Josephson currents. We display such an anomalous Josephson vortex and suggest how to evaluate the enlargement.

preprint2009arXiv

First-Principles Study for the Anisotropy of Iron-based Superconductors toward Power and Device Applications

Performing the first-principles calculations, we investigate the anisotropy in the superconducting state of iron-based superconductors to gain an insight into their potential applications. The anisotropy ratio $γ_λ$ of the c-axis penetration depth to the ab-plane one is relatively small in BaFe2As2 and LiFeAs, i.e., $γ_λ\sim 3$, indicating that the transport applications are promising in these superconductors. On the other hand, in those having perovskite type blocking layers such as Sr2ScFePO3 we find a very large value, $γ_λ\sim 200$, comparable to that in strongly anisotropic high-Tc cuprate Bi2Sr2CaCu2O{8-δ}. Thus, the intrinsic Josephson junction stacks are expected to be formed along the c-axis, and novel Josephson effects due to the multi-gap nature are also suggested in these superconductors.

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.