Routes to probe Bismuth induced strong-coupling superconductivity in bimetallic BiIn alloys

We report the observation of strong electron-phonon coupling in intergranular linked BiIn superconductors over an infinite range mediated by low-lying phonons. An enhanced superconducting transition temperature was observed from the magnetization, revealing a main diamagnetic Meissner state below TC(0) = 5.86(1) K and a critical field HC(0) = 1355(15) Oe with an In2Bi phase of the composite sample. The electron-phonon coupling to low lying phonons is found to be the leading mechanism for observed strong-coupling superconductivity in the BiIn system. Our findings suggest that In2Bi is in the strong-coupling region with TC(0) = 5.62(1) K, λep = 1.45, ωln = 45.92 K and α = 2.23. The estimated upper critical field can be well-described by a power law with α value higher than 2, consistent with the strong electron-phonon coupling.

that of virgin In. Since the atomic radii of Bi (156 pm) is smaller than In (167 pm), Bi can therefore occupy interstitial sites in the In lattice and can form an interstitial solid solution 8 . The In phase becomes unstable beyond y = 0.1 and shows the presence of two immiscible phases in each alloy up to y = 0.4. A horizontal red-colored dashed line drawn at y = 0.5 marks the two nuclear peaks (110) and (102) indexed based on P6 3 /mmc (No. 194), indicating the formation of a pure In 2 Bi crystalline phase. For y = 0.6, two additional nuclear peaks (202) and (111) become visible indexed based on I4/mcm (No. 140) and P4/nmm (No. 129), respectively indicating the formation of immiscible In 5 Bi 3 and InBi crystalline phases. For y from 0.7 to 0.99, along with InBi one additional nuclear peak (10-2) becomes visible indexed based on R-3m (No. 166), and trigonal-hexagonal scalenohedral Bi indicates the presence of two immiscible phases in each alloy. In bimetallic alloys, because of the different atomic and structural properties of the constituent's elements and the state of the alloy (pure or mixed), lattice expansion or contraction is unavoidable. The detailed structural analysis of the x-ray diffraction pattern of Bi y In 1-y bimetallic alloys was further refined by Rietveld analysis 9 using the GSAS software package 10 . The refined patterns and the fitting parameters are as shown in Fig. 2. From the refined XRD spectra of Bi y In 1-y bimetallic alloys, it was observed that depending on their initial composition, formed alloys are either in pure or mixed form. A single phase was obtained at y = 0.01, as shown in Fig. 2(a), exhibiting the expected In (S) solid solution structure. Since the atomic radii of Bi is smaller than In, Bi can therefore occupy interstitial sites in the In lattice and form an interstitial solid solution. The major detectable phase was In (S), along with In 2 Bi (S) as a secondary phase, for the region of y = 0.1 to 0.4, presenting two immiscible phases in this concentration region, as shown in Fig. 2(b-e). After this stage, the reactant phases have converted fully into In 2 Bi (S) for x = 0.5, as shown in Fig. 2(f). It is to be noted that there is a retardation in the growth of In 2 Bi (S) as the percentage of Bi increases. Structural transformation was observed as y is larger than 0.5, and then a new immiscible phase including In 5 Bi 3 (S) and InBi (N) was detected at y = 0.6, forming a co-existence of a superconductor (S) and a normal (N) bimetallic compound, and the pairing correlation between electrons is delivered from the superconducting phase to the normal phase through proximity effect (PE), as shown in Fig. 2(g). Figure S1(a) represents the observed state of the alloy (i.e. pure or mixed phases) with respect to initial composition y and the corresponding crystal structure of In, In 2 Bi, In 5 Bi 3 , InBi and Bi, respectively (see the supporting information: Figure S1). The formed pure In 2 Bi phase at y = 0.5, is marked by a red circle. The observed newly formed compounds with respect to the initial composition contradict with previous finding. Giessen et al. 11 reported the equilibrium phase diagram of In-Bi according to which pure In 2 Bi, In 5 Bi 3 and InBi are formed at 33.5(0.5), 37.5 (5) and 50 atomic % of Bi. The reason for the observed discrepancy is not clear yet, but for our alloys the estimated atomic % and structural properties are in good agreement with the initial In-Bi composition. The effect of Bi doping in In (solid solution) and the effect of strain on the lattice constants of the In-solid solution, In 2 Bi, In 5 Bi 3 and InBi, is plotted with respect to initial composition y, as shown in Figure S1(b-e) (the dotted line represents the standard values 5 ). Lattice expansion was observed for the In-solid solution, InBi and Bi phases, and contraction was observed for the In 5 Bi 3 phase of the mixed phase alloys. For In 2 Bi, phase lattice contraction was observed along the axial plane, and expansion along the basal plane. The observed lattice expansion for the solid solution is the combined effects of electron doping, as Bi has one extra electron compared to In and the strain. On the other hand, the observed lattice expansion/ contraction in immiscible mixed alloys could be due to strain effect, since thermal coefficients vary from material to material. The atomic position and structural properties of the parent elements, newly formed In 2 Bi (y = 0.5), In 5 Bi 3 (y = 0.6), and InBi (y = 0.6) crystalline phases obtained after refinement are as tabulated in the supplementary information of Table S1, along with the atomic radii of In and Bi. The observed lattice constants of the different crystalline phases are in excellent agreement with previously published results. As tabulated in Table S1, it can be observed that the c/a of In 2 Bi, In 5 Bi 3 and InBi, 1.1992, 1.4830 and 0.9548, respectively, is smaller than that of both In and Bi, 1.5208 and 2.6025, respectively. Such a contracted lattice, compared to that of the constituent elements, may result in softening of the vibrational mode, that is, phonon softening, thus giving rise to the enhanced electron-phonon coupling strength. However, it has been reported that In 2 Bi and In 5 Bi 3 are superconductors below 5.1 and 4.1 K, respectively, whereas InBi is non-superconducting down to 0.5 K 6 . Allen et al. 2 observed strong-coupling behavior for the In 2 Bi phase mediated by low energy phonons. Therefore, studying the superconducting properties of In-Bi bimetallic alloys will provide more insight into these newly formed compounds.

Temperature Dependency of ZFC and FC Magnetization.
A superconducting material when cooled below unique superconducting transition temperature T C shows perfect Meissner effect. However in Bi y In 1-y bimetallic alloys, the effect of electron doping, strain effect, and superconducting proximity (SPE) 12 structural and vibrational properties of newly formed In 2 Bi, In 5 Bi 3 and InBi compounds on the superconducting properties are unavoidable. To study the superconducting properties of these bimetallic alloys we carried out detailed field and temperature dependent magnetization measurements. Magnetic properties were measured using a Quantum Design MPMS VSM SQUID magnetometer. The transition temperature below which the diamagnetic Meissner effect becomes dominant was confirmed by measuring the zero-field cooled (ZFC) and field-cooled (FC) magnetization. The ZFC is a screening measurement, performed by first cooling the sample in a zero magnetic field to 2 K, subsequently applying a small magnetic field, H a = 10 Oe, and then measuring the resultant magnetization while warming the sample. The FC is the Meissner flux expulsion measurement, performed by first cooling the sample to 2 K in a small magnetic field, H a = 10 Oe, and then measuring the resultant magnetization while warming the sample. Four selected samples of field cooled and zero-field cooled M(T) measurements recorded in the 1.8-6 K range with H a = 10 Oe, as represented in Fig. 3(a-d). As expected, two step-like transitions were observed in Fig. 3(b) due to two immiscible phases that can be described using a superposition of the London equation 13 for In 2 Bi(M1) and α-In(M2), showing two step-like transitions, T C1 and T C2 , which are consistent with the multiphase bismuth-indium system obtained by Hutcherson et al. 6 and Currie et al. 14 . The step-like magnetization can be well fitted using a modified London equation: Where <d> is the mean size, ρ is mass density, and λ L is the London penetration depth which varies according to . Conventionally, P = 4 has been utilized to define bulk-like superconductivity for spherical particles. Li et al. 15 used the above London equation to define the susceptibility in well-separated spherical particles of lead with fix P = 4. However, the value of P defines only the distribution of the transition temperature, i.e., the higher the value of P ( ≥ 1) the steeper the distribution of transition temperature will be. In such a scenario, the bulk-like magnetization M(T), can be well fitted using P > 4. The solid red and blue lines shown in Fig. 3(a-h) represent the fitted curves for ZFC and FC, respectively, with M(T) = M 1 + M 2 , and the steepness of M(T) represented by the value of P and the penetration depth at H a = 10 Oe. The fitted values of T C1 (10 Oe) and T C2 (10 Oe) with respect to the initial composition y are as shown in Fig. 3(i). The highest value of T C1 = 5.85 K was observed for the y = 0.6 sample (and the T C1 of y = 0.7 is close to that of the In 2 Bi phase), indicating the presence of an In 2 Bi phase. However, from the XRD data (see Fig. 2(g and h)), we did not observe any sign for In 2 Bi phase, which could be due to the relatively low superconducting volume fraction (SVF) of the In 2 Bi phase, as can be seen from the M(T) data. Therefore, the above finding shows the potentiality of superconducting measurement even over synchrotron radiation XRD, particularly for superconducting materials. Apparently, enhanced T C was observed for all superconducting phases. The observed enhanced T C for In 2 Bi and In 5 Bi 3 , as compared to that of the reported value of the respective single crystal, is because of the stress induced strain effect, as can be seen from the refined SR-XRD data. Hutcherson et al. 6   for In 2 Bi and In 5 Bi 3 phases, respectively. These reported values are somewhat lower, but still in good agreement with our finding. However, for higher concentrations, i.e. 40 to 50 atomic % of Bi, only transition temperature, 4.1 K, corresponding to In 5 Bi 3 has been reported which contradicts our observations. Furthermore, the observed enhanced T C for an In-Bi solid solution, α-In could be the combined effect of enhanced density of state on the Fermi surface and the stress induced strain effect. Since, T C is proportional to density of state and Bi has two extra valence electrons compared to In, this results in an increase of density of state N(E F ) on the Fermi surface with an increase of Bi concentration. The effect of electron and hole doping has been studied in a Tl-Pb-Bi system (which retains the same fcc-Pb structure for 0.8 Tl to 0.2 Bi concentrations) in which Tl has one fewer, and Bi has one more, valence electron than Pb. It has been reported that an increase of N(E F ) from Tl to Bi resulted in strong electron-phonon coupling mediated by low lying phonon and showed an enhancement of T C from 2.4 to 8 K 16 . However, the effect of stress induced strain on T C cannot be neglected, as samples from y = 0.1 to 0.4 are composite in nature. A normal to superconducting transition, namely T C1 , was observed for the y = 0.6 sample due to Josephson coupled with an anisotropic superconducting granular system via inter-granular weak links 7 , as shown in Fig. 3(g). The weak link could be achieved through a normal phase (e.g. InBi and Bi) or a superconductor with low T C or the grain boundaries. The calculated superconducting volume fraction (SVF) from the ZFC of the pure In 2 Bi phase ( Fig. 3(f)) is 94%, close to the perfect diamagnetism (−1/4π, according to the CGS system), indicating bulk like superconductivity without considering the demagnetization factor. The calculated low value of flux expulsion (from FC) of 2.7%, indicated strong flux pinning, possibly on the grain boundaries of In 2 Bi. The values of SVF and flux expulsion at 2 K with respect to initial composition are as shown in the inset of Fig. 4. The high value of flux expulsion and relatively low value of SVF in superconducting-normal phase composites, y = 0.6 and 0.7, suggest relatively strong-linked grains in these samples. The SVF exceeding 100% in composite y = 0.4 to 0.1 samples is attributed to the magnetization contribution from the two superconducting phases and the demagnetization factor of the sample. Due to the high density of defects and strong flux pinning, a very low value of flux expulsion was observed for these composites. As can be seen from Fig. 3(a), the FC and ZFC curve of y = 0.01, α-In, overlaps below T C (10 Oe) = 3.46 K down to lowest available temperature, which reflects reversible behavior due to weaker pinning and/or a lower density of defects. The calculated high value of SVF, 179% is attributed to the demagnetization effect of the sample. As expected, 101% of flux expulsion indicated an ideal, pinning-free bulk sample, with the sample always remaining in thermodynamic equilibrium. The similar reversible region was also observed for the In 2 Bi phase of y = 0.1 and y = 0.6 samples, suggesting that both grain coupling and intra-grain pinning are relatively poor. As we can see, the magnetization SVF for the In 2 Bi phase in the former and later samples is only ~9% and ~0.09%, respectively. Furthermore, the reversibility in the y = 0.6 composite could be due to strain caused by normal vortex cores embedded in the superconducting phase. According to Kogan et al. 17 the strains in this type of composite system arise due to the difference in density of the normal (i.e. InBi) and superconducting (i.e. In 2 Bi and In 5 Bi 3 ) phases, which are related to the stress dependency of the critical temperature.
Isothermal Magnetization and theoretical analysis. The observed overlapping and strong flux pinning from M(T) was further confirmed from the isothermal magnetization M(H a ) measured at 2 K, as shown in Fig. 5(a-h). The observed pronounced reversibility for y = 0.01, α-In is in excellent agreement with the overlapping of FC and ZFC. The rounded shape of the M(H a ) near penetration field Hp is plainly due to large local demagnetization at the corners 18 . Similar reversible behavior was also spotted for y = 0.1-0.2 and 0.6 composite samples above irreversible field H irr , as shown in Fig. 5  The M(H C1 ) is the magnetization at lower critical field H C1 , below which the superconductor retains the perfect Meissner effect, and H C2 is the upper critical field above which the superconductor behaves as normal conductor. The deviation of magnetization in the low field region of linearity gives the value of H C1 20 . The solid red lines in Fig. 5(a-h) represent the linear fit to M(H a ). The goodness of linearity fit in the low field region can be represented by differential 21 , as shown in Fig. 5(i). In general, the obtained H C1 in M(H a ) measurement suffers not only from demagnetization factor but also surface barrier effect. The type-II like hysteresis loop M(H a ) then reflects the characteristic of the penetration of magnetic flux into the alloy. As can be seen in Fig. 5(e), the M(H a ) loop displays an asymmetric profile in the field-increasing and -decreasing branches and a large penetration field Hp ∼ 230 Oe was obtained at T = 2 K for y = 0.4 sample. The irreversibility of the entry and exit of magnetic flux lines through the alloys surfaces shows that it is mainly the geometric edge barriers that control the movements of flux lines at the surface; a Bean-Livingston 22 surface barrier effect has been reported in high T C superconductors 23 . It should be noted that the determination of H C1 is of primary importance since it allows one to extract the magnetic penetration depth fundamental parameter characterizing the superconducting condensate and carrying information about the underlying pairing mechanism. A popular approach to measure H C1 consists of measuring the magnetization M(H a ) as a function of H a and then identify the deviation of the linear Meissner response which would correspond to the vortex penetration. One may argue that the H C1 values obtained from the deviation of magnetization in the low field region may not reflect the true H C1 but the flux entry field because of the Bean-Livingston surface barrier 22 . However, it is clear that the influence of surface barrier is important in our investigated alloys system since: (i) the magnetic hysteresis loops are asymmetric close to T C . (ii) Therefore, if the surface barrier should be taken into account, the true H C1 would be much smaller than the obtained value from the deviation of magnetization in the low field region. The upper critical field H C2 is defined as the field H a at which the magnetization M(H a ) shows zero value. Except for y = 0.01, α-In, the value of N is close to ~0.95(1), indicating superconducting lasting over an infinite range with H a applied perpendicular to the surface. This accordance reveals that the coupling between the observed spherical grains from the SEM images extends over infinity. Thus, the demagnetization effect is governed by the outward appearance of the sample, not by the size and shape of each grain. In other words, as discussed before, the sample shows bulk-like superconductivity due to weakly-connected grains, even though the grain boundaries exist. Furthermore, Fig. 6 Table S2.
To estimate the thermal effect of the superconducting properties of a Bi y In 1-y system, M(H a ) loop measurement of all the samples (y = 0.01 to 0.7) was carried out at various temperatures below the transition temperature, as can be seen in Figure S2 (see the supporting information: Figure S2). Following a similar procedure as mentioned before, the values of H C1 and H C2 were estimated at various temperatures and the observed values plotted as shown in Fig. 7(a and b). The H C1 and H C2 can be fitted with a power law, in accordance with , where β and δ are the fitting parameters with a value close to ~1(0.1), and H C1 (0) and H C2 (0) are the lower and upper critical field at zero temperature, respectively (see in the supplementary information of Table S2). The value of α represents the coupling strength obtained from the α-model (which will be discussed later in the text) and it was kept fixed for respective samples. The resulting fitted values of H C1 (0) and H C2 (0) allow us to calculate the coherence length ξ GL (0) and the penetration depth λ(0) using a Ginzburg-Landau (GL) formulation 24 2 is the quantum flux 25 , which yields minimum ξ GL (0) ~ 40.1 nm (y = 0.4) and maximum λ(0) ~ 220 nm (y = 0.7) at 0 K. Using the relation ξ κ = λ (0) (0)/ (0) GL we calculated the value κ(0), which varied from 1.6 (y = 0.01) to 5.0 (y = 0.7), indicating type II superconducting behavior. In addition, the thermodynamic critical field, was calculated, showing a maximum value of ~534 Oe for the In 2 Bi superconductor (y = 0.5). The estimated superconducting parameters for all the BiIn samples are as tabulated in the supplementary information of Table S3.
Applied magnetic field dependency of ZFC and FC Magnetization. The cause behind an enhanced transition temperature of In 2 Bi and In 5 Bi 3 has not been clearly identified yet. There are a few possibilities, such as internal strain or short range disorder, as reported for an InSn system, or the enhanced electron-phonon coupling mediated by low lying phonons, as observed for a BiPb system. The effect due to a short range disorder can be avoided, as the calculated mean free path at 2 K is a multiple of lattice constants. From GL theory, the carrier mean free path, , at low temperature T = 2 K can be estimated using where ξ(0) = 3640 Å is the BCS coherence length of bulk In 26 . To estimate the superconducting coupling strength, we carried out ZFC measurement at various external applied field H a , as shown in Figure S3 (see the supporting information: Figure S3). The Meissner diamagnetic screening signal appears to shift consistently to a lower temperature as H a is increased. The relatively sharp distribution of the transition temperature from the y = 0.  to what was observed in the InBi system. However, they correlated the broad distribution of transition temperatures with weak-pinning behavior, which is contradictory to a BiIn system. The observed broad distribution could be the sign of a weakly-linked granular superconductor, as these BiIn superconductors include both weak and strong-pinned superconductors. The fitted value of P to the ZFC of all the superconducting samples within the Meissner state is higher than 4, whereas in the flux pinning state it lies between 1 to 2. Therefore, values of P can help in defining the state of the linked intergranular superconducting system. The fitted values of P(H a ) and the London penetration depth λ(H a ) are as plotted in Fig. 8(a and b), respectively. Within a Meissner state λ(H a ) remains constant or shows slightly decreasing behavior with an increase of H a (≤H C1 (0)). The reduction of the diamagnetic screening in the flux pinning state is known to be caused by the increase of λ with H a , which in turn reduces the SVF, hence the shielding capability of a superconductor 29 . Furthermore, the suppression of critical temperature with an increase of field H a can be well described using , where T C (0) is the zero field critical temperature, H C (0) is the zero temperature critical field, and γ is a fitting parameter, as shown in Fig. 9(a). All the fitting parameters are as tabulated in the supplementary information of Table S4. The fitted value of H C (0) for each sample is close to the H C2 (0) and the maximum value is ~6.5 times that of the pure Indium phase. Padamsee et al. 30 introduced the so-called α-model, by which the relative coupling strength, α = Δ(0)/k B T C (where, Δ(0) is the superconducting energy gap at 0 K and k B is the Boltzmann constant) of a superconducting system can be revealed in the deviation of H C (T)/H C (0) from the parabolic dependency of 1 − (T/T C ) 2 . It is known that weak-coupled systems yield negative deviations, while strong-coupled systems yield positive deviations 31 . The observed negative deviation for α-In (y = 0.01) and the positive deviation for the remaining BiIn superconductor (y = 0.1 to 0.7) can be fitted by , as shown in Fig. 9(b). The solid line represents the fitted curve, and the fitted values of α are as shown in the supplementary information of Table S4. A minimum of α = 1.835(2) (close to weak-coupled bulk In, 1.9) from y = 0.01 and a maximum of 2.231(9) (larger than strong-coupled bulk Pb) from y = 0.5, indicate the pure In 2 Bi phase is larger than BCS prediction of 1.764. The fitted values of the coupling strength for a In 2 Bi superconductor could be mediated by low lying phonons resulting in an enhanced T C . The average value of logarithmic phonon energy ω ln and the dimensionless electron-phonon coupling constant λ ep for a sample of BiIn can be calculated from the Eliashberg theory-based Allen-Dynes formulation 2 . The corrections of the BCS values by strong electron-phonon interactions have been deduced in the following approximate analytic formulas that link ω ln /T C to experimental thermodynamic quantities: The calculated values of ω ln are as shown in Table S4, and we see that a minimum of 45.92 K was calculated for the pure In 2 Bi sample (y = 0.5). Next, the dimensionless λ ep can be estimated from McMillan's formulation 32 . However, as can be seen from Fig. 10, the McMillan formulation is useful only for intermediate coupled superconductors, λ ep < 1. As shown in Fig. 10, most of the phonon mediated superconductors can be well described using Allen and Dynes' modified McMillan formula 2 , and ϖ ω = < > ln 2 2 is the square root average logarithm of the phonon frequency introduced by Allen and Dynes, and μ * is the Coulomb pseudo potential. The μ * describes the Coulomb pseudopotential to represent the repulsive part of the pairing interaction, and to calculate λ ep we set μ * = 0.11 and ϖ ω = .
1 3 ln 2 . The plot of T C /ω ln vs. λ ep of the BiIn system (red filled circle) along with different phonon mediated superconductors is as shown in Fig. 10. The blue solid line represents the calculated T C /ω ln using Allen and Dyne's formula. For comparison, we have also shown the results from the BCS and McMillan formulation (black and orange solid lines, respectively). From the above Fig. 10, for λ ep < 1, for the Allen and Dyne formulation reduced to McMillan formulation 32 and for λ ep < 0.2, both the Allen-Dyne and McMillan formulations are reduced to the BCS prediction 33,34 . The calculated value of λ ep = 1.453, ω ln = 45.92 K for BiIn is in good agreement with the earlier reported value of λ ep = 1.40, ω ln = 46 K 2 . However, using McMillan's formula with μ * = 0.11 and assuming ω ln to be same as a Debye frequency, we found that the ratio T C /ω ln for In 2 Bi is reproduced with relatively high λ ep ~1.68, which contradicts previous findings. In order to match the experimental value of T C the required λ would be 1.2-1.4 (with reasonable Coloumb pseudo potential μ* of 0.06-0.09), reflecting the uncertainty in the estimation of Coloumb pseudo potential. It is known that disorder/substitution can affect the Coulomb pseudo potential in superconductor 35 . Hence we could infer that the local inhomogeneity introduced by disorder/substitution in bimetallic BiIn alloys may play nontrivial  role in the influence of Coulomb pseudo potential. Indeed, the intimate relationship between local inhomogeneity and superconductivity has been extensively studied in the context of high T C superconductors, and the possible enhancement of superconductivity by local inhomogeneities has been discussed recently by Martin et al. 36 in the weak coupling BCS regime. From Fig. 10, we can conclude that enhanced T C is the combined effect of phonon energy and the electron-phonon coupling. For example, ε-Pb 0.7 Bi 0.3 (λ ep ~ 2.01, ω ln = 47 K) and Pb 0.8 Bi 0.2 (λ ep ~ 1.88, ω ln = 46 K) have the same phonon energy as In 2 Bi, but higher values of λ ep , resulting in enhanced Tc of 8.45 K and 7.95 K, respectively 2 . Similarly, Mo 6 Se 8 (λ ep = 1.27, ω ln = 70 K) and SrPt 3 P (λ ep = 1.33, ω ln = 77 K) which overlap InBi superconductors, as shown in the figure, also show enhanced T C of 8.5 and 6.34 K, respectively 3,4 . Since, λ ep for these two superconductors is lower than In 2 Bi, but with high value of phonon energy. Moreover, there are a few superconductors, as shown in the figure, which are beyond Allen and Dynes' prediction, and therefore more theoretical work needs to be carried out. It is known that disorder, substitution in weak coupling superconductivity could enhance the electron-phonon coupling constant. A simplified physical model of electron-phonon coupling has been developed by Gao et al. 37 to allow heat transfer from phonons to electrons, and applied to study defect or disorder as a function of the strength of electron-phonon coupling. The number of point defects produced in the primary damage state increases with increasing strength of electron-phonon coupling, implying the defect or disorder/substitution plays a role in the increasing strength of electron-phonon coupling constant. As reported in Cu x Bi 2 Se 3 F superconductor by Zhang et al. 35 , first principle calculation was introduced to investigate the role of electron-phonon coupling for the superconducting pairing in the prime candidate Cu x Bi 2 Se 3 . In their comprehensive analysis reveal that electron-phonon scattering process in this system is dominated by zone center and boundary optical modes, with coexistence of phonon stiffening and softening, suggesting that superconductivity may not only come from pure electron-phonon coupling but also disorder effect. In the present work, we also cannot ignore the possible enhancement of electron-phonon coupling by local inhomogeneity introduced by doping defect or disorder/substitution in bimetallic BiIn alloys. From the above finding, it's clear that the observed enhanced T C for a Bi y In 1-y system (i.e. In 2 Bi and In 5 Bi 3 ) is the effect of strong electron-phonon coupling to low-lying phonons. The effect of low lying phonons can be further proved by carrying out electron-tunneling spectroscopy measurements.

Discussion
The strong-coupling superconducting behavior mediated by low lying phonons has been reported to be the cause for enhanced T C in Bi y In 1-y superconductors. The maximum strong-coupling superconducting strength, α = 2.231(9), phonon energy, ω ln = 45.9 K, electron-phonon coupling constant, λ ep = 1.453, and critical field, H C (0) = 1909(3) Oe has been estimated from pure In 2 Bi (y = 0.5) samples using Allen and Dynes' formulation. The different superconducting parameters estimated using GL formulation show Bi y In 1-y is a type-II superconductor with GL constant, 1.6 ≤ k(0) ≤ 5 for 0.01 ≤ y ≤ 0.7. The room temperature SR-XRD refinement carried out by using weighted analysis and superconducting measurement reveals the formation of both pure α-In (y = 0.01), In 2 Bi (y = 0.5) and composite α-In + In 2 Bi (y = 0.1-0.4), In 2 Bi + In 5 Bi 3 + InBi (y = 0.6-0.7) samples, respectively. The ZFC-FC temperature dependent measurement revealed two-gap superconductivity in y = 0.1-0.4 (α-In + In 2 Bi) and 0.6-0.7 (In 2 Bi + In 5 Bi 3 ) and a single gap in y = 0.01(α-In) and y = 0.5 (In 2 Bi) superconductors, respectively. The highest transition temperature was determined to be T C1 (10 Oe) = 5.85 K and T C2 (10 Oe) = 4.27 K (y = 0.6) when fitted to the London equation with free fitting parameter P. The value of P gives an idea about the distribution of transition temperatures and can help in naming the Meissner/flux pinning state. The SVF estimated from ZFC revealed bulk-like superconductivity for all the superconducting samples. The high and low values of flux expulsion calculated from FC revealed that y = 0.01, 0.6-0.7 are strong-linked, whereas y = 0.1 to 0.5 are relatively weak-linked intergranular superconductors. The estimated demagnetization factor N ~ 0.95(1) (y = 0.1 to 0.7) from field dependent magnetization measurement signals that linkage within the grains, as observed from SEM images, extends over infinity. We argue that the observed broad distribution of transition temperatures in the flux pinning state is the special feature of strong/weak-linked granular superconductors. Finally, in this report, we have demonstrated the potentiality of superconducting measurement over the SR-XRD for identifying very small volume fractions of a superconducting phase.