Source author record

F. Marsiglio

F. Marsiglio 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

42works
13topics
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

42 published item(s)

preprint2022arXiv

Clear evidence against superconductivity in hydrides under high pressure

The Meissner effect, magnetic field expulsion, is a hallmark of superconductivity. Associated with it, superconductors exclude applied magnetic fields. Recently Minkov et al. presented experimental results reportedly showing "definitive evidence of the Meissner effect" in sulfur hydride and lanthanum hydride under high pressure [1], and more recently Eremets et al. argued that "the arguments against superconductivity (in hydrides) can be either refuted or explained" [2]. Instead, we show here that the evidence presented in those papers does not support the case for superconductivity in these materials. Together with experimental evidence discussed in earlier papers, we argue that this strongly suggests that hydrides under pressure are not high-temperature superconductors.

preprint2022arXiv

Impact of Retardation in the Holstein-Hubbard Model: a two-site Calculation

Eliashberg theory provides a theoretical framework for understanding the phenomenon of superconductivity when pairing between two electrons is mediated by phonons, and retardation effects are fully accounted for. However, when a direct Coulomb interaction between two electrons is also present, this interaction is only partially accounted for. In this work we use a well-defined Hamiltonian, the Hubbard-Holstein model, to examine this competition more rigorously, using exact diagonalization on a two-site system. We find that the direct electron-electron repulsion between two electrons has a significantly more harmful effect on pairing than suggested through the standard treatment of this interaction.

preprint2021arXiv

Absence of evidence of superconductivity in sulfur hydride in optical reflectance experiments

Capitani and coworkers [1] reported that infrared optical reflectance measurements provided evidence for a superconducting transition in sulfur hydride [2] under 150 GPa pressure, and that the transition is driven by the electron-phonon interaction. Here we argue that the measured data did not provide evidence that the system undergoes a transition to a superconducting state, nor do the data support any role of phonons in driving a transition. Rather, the data are consistent with the system remaining in the normal state down to temperature 50K, the lowest temperature measured in the experiment. This calls into further question [3,4] the generally accepted view [5] that sulfur hydride under pressure is a high temperature superconductor.

preprint2021arXiv

Intrinsic hysteresis in the presumed superconducting transition of hydrides under high pressure

Superconducting transitions in the absence of magnetic field should be non-hysteretic. Here we address the fact that the drops in electrical resistance that have been interpreted as evidence of superconductivity in several hydrides under high pressure (so-called "superhydrides") show hysteresis. We argue that the experimental evidence shows that the observed hysteresis cannot be attributed to experimental artifacts but is intrinsic to the samples. Assuming that the drops in resistance signal a thermodynamic phase transition, we argue that the presence of intrinsic thermal hysteresis indicates that these are first order transitions, whereas for standard superconductors the transition in the absence of applied magnetic field is always second order. We conclude that this is another feature that qualitatively distinguishes superhydrides from standard superconductors, in addition to the ones that have been pointed out earlier [1,2], $assuming$ these materials are superconductors. Alternatively and more likely, whether or not the drops in resistance signal a thermodynamic phase transition, our analysis indicates that superhydrides are not superconductors.

preprint2021arXiv

Phonon Self-Energy Effects in Migdal-Eliashberg Theory

Recent work on an electron-phonon model in two-dimensions is reviewed and some new results are presented. We utilize both Quantum Monte Carlo simulations and Migdal-Eliashberg theory. The type of electron-ion coupling considered is on-site. Competing instabilities are charge density wave (CDW) and singlet superconductivity. When the electron band is half-filled, the charge density wave dominates. Away from half-filling CDW correlations actually suppress superconductivity.

preprint2021arXiv

Scattering Problems via Real-time Wave Packet Scattering

In this paper, we use a straightforward numerical method to solve scattering models in one-dimensional lattices based on a tight-binding band structure. We do this by using the wave packet approach to scattering, which presents a more intuitive physical picture than the traditional plane wave approach. Moreover, a general matrix diagonalization method that is easily accessible to undergraduate students taking a first course in quantum mechanics is used. Beginning with a brief review of wave packet transport in the continuum limit, comparisons are made with its counterpart in a lattice. The numerical results obtained through the diagonalization method are then benchmarked against analytic results. The case of a resonant dimer is investigated in the lattice, and several resonant values of the mean wave packet momentum are identified. The transmission coefficients obtained for a plane wave incident on a step potential and rectangular barrier are compared by investigating an equivalent scenario in a lattice. Lastly, we present several short simulations of the scattering process which emphasize how a simple methodology can be used to visualize some remarkable phenomena.

preprint2020arXiv

Thermodynamics of Eliashberg theory in the weak-coupling limit

The weak-coupling limits of the gap and critical temperature computed within Eliashberg theory surprisingly deviate from the BCS theory predictions by a factor of $1/\sqrt{e}$. Interestingly, however, the ratio of these two quantities agrees for both theories. Motivated by this result, here we investigate the weak-coupling thermodynamics of Eliashberg theory, with a central focus on the free energy, specific heat, and the critical magnetic field. In particular, we numerically calculate the difference between the superconducting and normal-state specific heats, and we find that this quantity differs from its BCS counterpart by a factor of $1/\sqrt{e}$, for all temperatures below $T_{c}$. We find that the dimensionless ratio of the specific-heat discontinuity to the normal-state specific heat reduces to the BCS prediction given by $ΔC_{V}(T_{c})/C_{V,n}(T_c)\approx1.43$. This gives further evidence to the expectation that all dimensionless ratios tend to their "universal values" in the weak-coupling limit.

preprint2019arXiv

Eliashberg theory in the weak-coupling limit: results on the real frequency axis

We formulate and solve the Eliashberg equations on the imaginary frequency axis at temperatures below $T_c$ in the weak-coupling limit. We find an excellent scaling at all temperatures, for a given coupling strength, and the normalized order parameter exhibits a BCS-like temperature dependence. The hybrid real-imaginary axis equations are also solved to obtain numerically exact analytic continuations from the imaginary frequency axis to the real frequency axis. This provides a determination of the gap edge, which, in the weak-coupling limit, is identical to the order parameter from the imaginary axis. The analytical result for the zero-temperature gap edge deviates from the BCS result by a factor of $1/\sqrt{e}$, which was also obtained for the transition temperature $T_c$. We show that the normalized gap function on both the real and imaginary frequency axes, for an electron-phonon Einstein spectrum ($δ$-function) of a given strength, is a universal function of frequency, independent of temperature. The $1/\sqrt{e}$ correction is a result of this non-trivial frequency dependence in the gap function. This modification, in the gap edge and in $T_{c}$, serves to preserve various dimensionless ratios to their BCS values.

preprint2019arXiv

Eliashberg Theory: a short review

Eliashberg theory is a theory of superconductivity that describes the role of phonons in providing the attractive interaction between two electrons. Phonon dynamics are taken into account, thus giving rise to retardation effects that impact the electrons, in the form of a frequency-dependent electron self-energy. In the superconducting state, this means that the order parameter, generally considered to be a static quantity in the Bardeen-Cooper-Schrieffer (BCS) theory, also becomes frequency-dependent. Here we review the finite temperature formulation of Eliashberg theory, both on the imaginary and real frequency axis, and briefly display some examples of the consequences of a dynamical, as opposed to static, interaction. Along the way we point out where further work is required, concerning the validity of some of the assumptions used.

preprint2016arXiv

A systematic study of the superconducting critical temperature in two and three dimensional tight-binding models: a possible scenario for superconducting H$_3$S ?

Ever since BCS theory was first formulated it was recognized that a large electronic density of states at the Fermi level was beneficial to enhancing $T_c$. The A15 compounds and the high temperature cuprate materials both have had an enormous amount of effort devoted to studying the possibility that such peaks play an important role in the high critical temperatures existing in these compounds. Here we provide a systematic study of the effect of these peaks on the superconducting transition temperature for a variety of tight-binding models of simple structures, both in two and three dimensions. In three dimensions large enhancements in $T_c$ can occur, due to van Hove singularities that result in divergences in the density of states. Furthermore, even in more realistic structures, where the van Hove singularity disappears, large enhancements in $T_c$ continue due to the presence of `robust' peaks in the densities of states. Such a peak, recently identified in the bcc structure of H$_3$S, is likely the result of such a van Hove singularity. In certain regimes, anomalies in the isotope coefficient are also expected.

preprint2016arXiv

The Effect of Next-Nearest Neighbour Hopping in the One, Two, and Three Dimensional Holstein Model

Allowing a single electron to hop to next-nearest neighbours (NNN) in addition to the closest atomic sites in the Holstein model, a modified Trugman method is applied to exactly calculate the effect on the polaronic effective mass in one, two, and three dimensions, building on the previous study of the one-dimensional NNN Holstein model. We also present perturbative calculations and a heuristic scaling factor for the coupling strength and ion frequency to nearly map the NNN Holstein model back onto the original Holstein model. When account is taken of the modified electronic bandwidth near the electron energy, we find that including NNN hopping effectively increases the polaron effective mass.

preprint2015arXiv

Asymmetric wave functions from tiny perturbations

The quantum mechanical behavior of a particle in a double well defies our intuition based on classical reasoning. Not surprisingly, an asymmetry in the double well will restore results more consistent with the classical picture. What is surprising, however, is how a very small asymmetry can lead to essentially classical behavior. In this paper we use the simplest version of a double well potential to demonstrate these statements. We also show how this system accurately maps onto a two-state system, which we refer to as a `toy model'.

preprint2015arXiv

Pairing effects in the normal phase of a two-dimensional Fermi gas

In a recent experiment [M. Feld et al., Nature 480, 75 (2011); B. Froehlich et al., Phys. Rev. Lett. 109,130403 (2012)], a pairing gap was detected in a two-dimensional (2D) Fermi gas with attractive interaction at temperatures where superfluidity does not occur. The question remains open as to whether this gap is a pseudogap phenomenon or is due to a molecular state. In this paper, by using a t-matrix approach, we reproduce quite well the experimental data for a 2D Fermi gas, and set the boundary between the pseudogap and molecular regimes. We also show that pseudogap phenomena occurring in 2D and 3D can be related through a variable spanning the BCS-BEC crossover in a universal way.

preprint2015arXiv

Superconductivity in the elements, alloys and simple compounds

We give a brief review of superconductivity at ambient pressure in elements, alloys, and simple three-dimensional compounds. Historically these were the first superconducting materials studied, and based on the experimental knowledge gained from them the BCS theory of superconductivity was developed in 1957. Extended to include the effect of phonon retardation, the theory is believed to describe the subset of superconducting materials known as `conventional superconductors', where superconductivity is caused by the electron-phonon interaction. These include the elements, alloys and simple compounds discussed in this article and several other classes of materials discussed in other articles in this Special Issue.

preprint2014arXiv

BLF-SSH Polarons coupled to Acoustic Phonons in the Adiabatic limit

We survey polaron formation in the BLF-SSH model using acoustic phonons in the adiabatic limit. Multiple different numerical optimization routines and strong coupling analytical calculations are used to find a robust ground state energy for a wide range of coupling strengths. The electronic configuration and accompanying ionic distortions of the polaron were determined, as well as a non-zero critical coupling strength for polaron formation in two and three dimensions.

preprint2014arXiv

Hole superconductivity in $H_2S$ and other sulfides under high pressure

Superconductivity at temperatures up to 190 K at high pressures has recently been observed in $H_2S$ and interpreted as conventional BCS-electron-phonon-driven superconductivity.\cite{h2s} Instead we propose that it is another example of the mechanism of hole superconductivity at work. Within this mechanism high temperature superconductivity arises when holes conduct through negatively charged anions in close proximity. We propose that electron transfer from $H$ to $S$ leads to conduction by holes in a nearly full band arising from direct overlap of $S^=$ $p$ orbitals in a planar structure. The superconductivity is non-phononic and is driven by pairing of heavily dressed hole carriers to lower their kinetic energy. Possible explanations for the observed lower critical temperature of $D_2S$ are discussed. We predict that high temperature superconductivity will also be found in other sulfides under high pressure such as $Li_2S$, $Na_2S$ and $K_2S$.

preprint2014arXiv

The extended vs standard Holstein model; results in two and three dimensions

We present numerically exact solutions to the problem of a single electron interacting through a long range interaction with optical phonons in two and three dimensions. Comparisons are made with results for the standard Holstein model, and with perturbative approaches from both the weak coupling and strong coupling sides. We find, in agreement with earlier work, that the polaron effective mass increases (decreases) in the weak (strong) coupling regime, respectively. However, in two dimensions, the decrease in effective mass still results in too large an effective mass to be relevant in realistic models of normal metals. In three dimensions the decrease can be more relevant, but exists only over a very limited range of coupling strengths.

preprint2014arXiv

The Kronig-Penney model extended to arbitrary potentials via numerical matrix mechanics

The Kronig-Penney model describes what happens to electron states when a confining potential is repeated indefinitely. This model uses a square well potential; the energies and eigenstates can be obtained analytically for a the single well, and then Bloch's Theorem allows one to extend these solutions to the periodically repeating square well potential. In this work we describe how to obtain simple numerical solutions for the eigenvalues and eigenstates for any confining potential within a unit cell, and then extend this procedure, with virtually no extra effort, to the case of periodically repeating potentials. In this way one can study the band structure effects which arise from differently-shaped potentials. One of these effects is the electron-hole mass asymmetry. More realistic unit cell potentials generally give rise to higher electron-hole mass asymmetries.

preprint2013arXiv

Solving for three-dimensional central potentials using matrix mechanics

Matrix mechanics is an important component of an undergraduate education in quantum mechanics. In this paper we present several examples of the use of matrix mechanics to solve for a number of three dimensional problems involving central forces. These include examples with which the student is familiar, such as the Coulomb interaction. In this case we obtain excellent agreement with exact analytical methods. More importantly, other interesting `non-solvable' examples, such as the Yukawa potential, can be solved as well. Much less mathematical expertise is required for these methods, while some minimal familiarity with the usage of numerical diagonalization software is necessary.

preprint2013arXiv

Vanishing of interband light absorption in a persistent spin helix state

Spin-orbit coupling plays an important role in various properties of very different materials. Moreover efforts are underway to control the degree and quality of spin-orbit coupling in materials with a concomitant control of transport properties. We calculate the frequency dependent optical conductivity in systems with both Rashba and Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling is tuned to be equal to the Rashba spin-orbit coupling, the interband optical conductivity disappears. This is taken to be the signature of the recovery of SU(2) symmetry. The presence of the cubic Dresselhaus spin-orbit coupling modifies the dispersion relation of the charge carriers and the velocity operator. Thus the conductivity is modified, but the interband contribution remains suppressed at most but not all photon energies for a cubic coupling of reasonable magnitude. Hence, such a measurement can serve as a diagnostic probe of engineered spin-orbit coupling.

preprint2012arXiv

Impact of Dresselhaus vs. Rashba spin-orbit coupling on the Holstein polaron

We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van-Hove singularity in the density of states will be shifted away from the bottom of the band and finally disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be suppressed; the trend will become more significant for low phonon frequency. The presence of two dominant spin-orbit couplings will make it possible to tune the effective mass with more varied observables.

preprint2012arXiv

Model of the electron-phonon interaction and optical conductivity of Ba$_{1-x}$K$_x$BiO$_3$ superconductors

We investigate the physical properties of the Ba$_{1-x}$K$_x$BiO$_3$ compounds with a focus on the optical properties. Results from the simple Holstein model, describing a single band coupled to an oxygen breathing mode with parameters derived from first principles calculations, are in excellent agreement with a broad range of experimental information. It accounts for an insulating parent compound at $x=0$ %, with a direct- (optical) and an indirect-gap, and a metal insulator transition around $x = 0.38$. Strong electron-phonon coupling leads to spectral weight redistribution over a frequency scale much larger than the characteristic phonon frequency and to strongly anharmonic phonons. We find that the metallic phase in the vicinity of phase boundary is close to the polaronic regime.

preprint2012arXiv

Optical conductivity for a dimer in the Dynamic Hubbard model

The Dynamic Hubbard Model represents the physics of a multi-band Hubbard model by using a pseudo-spin degree of freedom to dynamically modify the on-site Coulomb interaction. Here we use a dimer system to obtain analytical results for this model. The spectral function and the optical conductivity are calculated analytically for any number of electrons, and the distribution of optical spectral weight is analyzed in great detail. The impact of polaron-like effects due to overlaps between pseudo-spin states on the optical spectral weight distribution is derived analytically. Our conclusions support results obtained previously with different models and techniques: holes are less mobile than electrons.

preprint2012arXiv

The double well potential in quantum mechanics: a simple, numerically exact formulation

The double well potential is arguably one of the most important potentials in quantum mechanics, because the solution contains the notion of a state as a linear superposition of `classical' states, a concept which has become very important in quantum information theory. It is therefore desirable to have solutions to simple double well potentials that are accessible to the undergraduate student. We describe a method for obtaining the numerically exact eigenenergies and eigenstates for such a model, along with the energies obtained through the Wentzel-Kramers-Brillouin (WKB) approximation. The exact solution is accessible with elementary mathematics, though numerical solutions are required. We also find that the WKB approximation is remarkably accurate, not just for the ground state, but for the excited states as well.

preprint2012arXiv

The polaron-like nature of an electron coupled to phonons

When an electron interacts with phonons, the electron can exhibit either free electron-like or polaron-like properties. The latter tends to occur for very strong coupling, and results in a phonon cloud accompanying the electron as it moves, thus raising its mass considerably. We summarize this behaviour for the Holstein model in one, two and three dimensions, and note that the crossover occurs for fairly low coupling strengths compared to those attributed to real materials exhibiting conventional superconductivity.

preprint2012arXiv

The Spectral Decomposition of the Helium atom two-electron configuration in terms of Hydrogenic orbitals

The two electron configuration in the Helium atom is known to very high precision. Yet, we tend to refer to this configuration as a $1s\uparrow 1s\downarrow$ singlet, where the designations refer to Hydrogen orbitals. The high precision calculations utilize basis sets that are suited for high accuracy and ease of calculation, but do not really aid in our understanding of the electron configuration in terms of product states of Hydrogen orbitals. Since undergraduate students are generally taught to think of Helium, and indeed, the rest of the periodic table, in terms of Hydrogenic orbitals, we present in this paper a detailed spectral decomposition of the two electron ground state for Helium in terms of these basis states. The $1s\uparrow 1s\downarrow$ singlet contributes less than 93% to the ground state configuration, with other contributions coming from both bound and continuum Hydrogenic states.

preprint2012arXiv

Why is the ground state electron configuration for Lithium $1s^22s$ ?

The electronic ground state for Lithium is $1s^22s$, and not $1s^22p$. The traditional argument for why this is so is based on a screening argument that claims that the $2p$ electron is better shielded by the $1s$ electrons, and therefore higher in energy then the configuration that includes the $2s$ electron. We show that this argument is flawed, and in fact the actual reason for the ordering is because the electron-electron interaction energy is higher for the $2p-1s$ repulsion than it is for the $2s-1s$ repulsion.

preprint2011arXiv

Impact of spin-orbit coupling on the Holstein polaron

We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of a Rashba-like spin orbit interaction. Our results corroborate with previous work performed with the Momentum Average approximation and with weak coupling perturbation theory. We find that spin orbit coupling increases the effective mass in the regime with weak electron phonon coupling, and decreases the effective mass in the intermediate and strong electron phonon coupling regime. Analytical strong coupling perturbation theory results confirm our numerical results in the small polaron regime. A large amount of spin orbit coupling can lead to a significant lowering of the polaron effective mass.

preprint2010arXiv

Ground state properties of the Holstein model near the adiabatic limit

We adapt a variational procedure to calculate ground state properties of the Holstein model in the adiabatic limit. At strong coupling, this adaption leads to rapid convergence of results. The intermediate coupling regime is further handled with an adaptive algorithm. We also use semi-classically derived results for the adiabatic end-point, along with weak coupling perturbation theory. These establish weak and strong coupling (or large and small polaron, respectively) regimes in two dimensions or higher. As is well known, these are connected smoothly, but the cross-over becomes increasingly abrupt as the phonon frequency decreases.

preprint2010arXiv

Perturbation theory of the mass enhancement for a polaron coupled to acoustic phonons

We use both a perturbative Green's function analysis and standard perturbative quantum mechanics to calculate the decrease in energy and the effective mass for an electron interacting with acoustic phonons. The interaction is between the difference in lattice displacements for neighbouring ions, and the hopping amplitude for an electron between those two sites. The calculations are performed in one, two, and three dimensions, and comparisons are made with results from other electron-phonon models. We also compute the spectral function and quasiparticle residue, as a function of characteristic phonon frequency. There are strong indications that this model is always polaronic on one dimension, where an unusual relation between the effective mass and the quasiparticle residue is also found.

preprint2010arXiv

Two-site dynamical mean field theory for the dynamic Hubbard model

At zero temperature, two-site dynamical mean field theory is applied to the Dynamic Hubbard model. The Dynamic Hubbard model describes the orbital relaxation that occurs when two electrons occupy the same site, by using a two-level boson field at each site. At finite boson frequency, the appearance of a Mott gap is found to be enhanced even though it shows a metallic phase with the same bare on-site interaction $U$ in the conventional Hubbard model. The lack of electron-hole symmetry is highlighted through the quasi-particle weight and the single particle density of states at different fillings, which qualitatively differentiates the dynamic Hubbard model from other conventional Hubbard-like models.

preprint2008arXiv

Hall Conductivity in a Spin-Triplet Superconductor

We calculate the Hall conductivity for a spin-triplet superconductor, using a generalized pairing symmetry dependent on an arbitrary phase, $\vphi$. A promising candidate for such an order parameter is Sr$_{2}$RuO$_{4}$, whose superconducting order parameter symmetry is still subject to investigation. The value of this phase can be determined through Kerr rotation and DC Hall conductivity measurements. Our calculations impose significant constraints on $\vphi$.

preprint2008arXiv

The Harmonic Oscillator in Quantum Mechanics: A Third Way

Courses on undergraduate quantum mechanics usually focus on solutions of the Schrödinger equation for several simple one-dimensional examples. When the notion of a Hilbert space is introduced only academic examples are used, such as the matrix representation of Dirac's raising and lowering operators or the angular momentum operators. We introduce some of the same one-dimensional examples as matrix diagonalization problems, with a basis that consists of the infinite set of square well eigenfunctions. Undergraduate students are well equipped to handle such problems in familiar contexts. We pay special attention to the one-dimensional harmonic oscillator. This paper should equip students to obtain the low lying bound states of any one-dimensional short range potential.

preprint2007arXiv

Optical properties and electronic structure of ZrB$_{12}$

We report optical (6 meV - 4 eV) properties of a boride superconductor ZrB$_{12}$ ($T_c$ = 6 K) in the normal state from 20 to 300 K measured on high-quality single crystals by a combination of reflectivity and ellipsometry. The Drude plasma frequency and interband optical conductivity calculated by self-consistent full-potential LMTO method agree well with experimental data. The Eliashberg function $α_{tr}^2F(ω)$ extracted from optical spectra features two peaks at about 25 and 80 meV, in agreement with specific heat data. The total coupling constant is $λ_{tr}=1.0\pm0.35$. The low energy peak presumably corresponds to the displacement mode of Zr inside $B_{24}$ cages, while the second one involves largely boron atoms. In addition to the usual narrowing of the Drude peak with cooling down, we observe an unexpected removal of about 10 % of the Drude spectral weight which is partially transferred to the region of the lowest-energy interband transition ($\sim$ 1 eV). This effect may be caused by the delocalization of the metal ion from the center of the $B_{24}$ cluster.

preprint2004arXiv

Pairing Fluctuations in The Attractive Hubbard Model in the Atomic Limit

BCS theory accounts for the pairing instability in the weak coupling limit, but fails to describe pairing fluctuations above $T_c$. One possibility for describing these fluctuations in the dilute limit is the T-matrix approximation. We critically examine various degrees of self-consistency in the T-matrix formalism, along with a non-diagrammatic two-particle self-consistent (TPSC) formulation, in the strong coupling regime, where an exact solution is readily available. We find that one particular degree of self-consistency is quite accurate, particularly at low temperature as evidenced by examining both static and dynamic properties.

preprint2003arXiv

Dynamic Hubbard Model: Effect of Finite Boson Frequency

Dynamic Hubbard models describe coupling of a boson degree of freedom to the on-site electronic double occupancy. In the limit of infinite boson frequency this coupling gives rise to a correlated hopping term in the effective Hamiltonian and to superconductivity when the Fermi level is near the top of the band. Here we study the effect of finite boson frequency through a generalized Lang-Firsov transformation and a high frequency expansion. It is found that finite frequency enhances the tendency to superconductivity in this model.

preprint2000arXiv

Anderson's "Theorem" and Bogoliubov-de Gennes Equations for Surfaces and Impurities

In order to incorporate spatial inhomogeneity due to nonmagnetic impurities, Anderson [1] proposed a BCS-type theory in which single-particle states in such an inhomogeneous system are used. We examine Anderson's proposal, in comparison with the Bogoliubov-de Gennes equations, for the attractive Hubbard model on a system with surfaces and impurities. [1] P. W. Anderson, J. Phys. Chem. Solids {\bf 11}, 26 (1959).

preprint1999arXiv

Fixed Number and Quantum Size Effects in Nanoscale Superconductors

In recent experiments on nanoscale Al particles, whose electron number was fixed by charging effects, a ``negative gap'' was observed in particles with an odd number of electrons. This observation has called into question the use of a grand canonical ensemble in describing superconductivity in such ultrasmall particles. We have studied the effects of fixed electron number and finite size in nanoscale superconductors, by applying the canonical BCS theory for the attractive Hubbard model. The ground state energy and the energy gap are compared with the conventional and parity-projected grand canonical BCS results, and in one dimension also with the exact solutions by the Bethe ansatz. The crossover from the bulk to quantum limit is studied for various regimes of electron density and coupling strength. The effects of boundary conditions and different lattice structures are also examined. A ``negative gap'' for odd electron number emerges most naturally in the canonical scheme. For even electron number, the gap is particularly large for ``magic numbers'' of electrons for a given system size or of atoms for a fixed electron density. These features are in accordance with the exact solutions, but are essentially missed in the grand canonical results.

preprint1999arXiv

New Solutions of the T-Matrix Theory of the Attractive Hubbard Model

This short paper summarizes a calculational method for obtaining the dynamical properties of many-body theories formulated in terms of (unrenormalized) bare propagators (and more generally, in terms of meromorphic functions, or convolutions over meromorphic functions) to a very high accuracy. We demonstrate the method by applying it to a T-matrix theory of the attractive Hubbard model in two dimensions. We expand the pair propagator using a partial fraction decomposition, and then solve for the residues and pole locations of such a decomposition using a computer algebra system to an arbitrarily high accuracy (we used MapleV and obtained these quantities to a relative error of 10^(-80)). Thus, this method allows us to bypass all inaccuracies associated with the traditional analytical continuation problem. Our results for the density of states make clear the pronounced development of a pseudogap as the temperature is lowered in this formulation of the attractive Hubbard model.

preprint1997arXiv

Quasiparticle Lifetimes and the Conductivity Scattering Rate

We compute the single-particle inverse lifetime, along with the conductivity-derived scattering rate, for a metallic system in an s-wave superconducting state. When both electron-phonon and electron-impurity scattering are included, we find that while these scattering rates are in qualitative agreement, in general quantitative agreement is lacking. We also derive results for the quasiparticle lifetime within the BCS framework with impurity scattering, which makes it clear that impurity scattering is suppressed for electrons near the Fermi surface in the superconducting state.