-
Projected Augmented Wave (PAW) Method: extended resolution of unity method
Authors:
Garry Goldstein
Abstract:
The Projected Augmented Wave (PAW) Method is based on generating a transform between the pseudo wavefunction and all electron wavefunction. For the accuracy of the method, it is important that the local part of the transform (in each atomic sphere $\mathbf{R}$) be over a complete basis set (with deviations from completeness leading to corrections to the total energy not computed within current imp…
▽ More
The Projected Augmented Wave (PAW) Method is based on generating a transform between the pseudo wavefunction and all electron wavefunction. For the accuracy of the method, it is important that the local part of the transform (in each atomic sphere $\mathbf{R}$) be over a complete basis set (with deviations from completeness leading to corrections to the total energy not computed within current implementations of PAW). Here we show how to make this basis much closer to complete without significant additional computational work without modifying the transformation in any way by extending the resolution of unity used for the transform to include more wavefunctions and having them transform via the identity.
△ Less
Submitted 7 July, 2024; v1 submitted 30 May, 2024;
originally announced May 2024.
-
Energy Window Augmented Plane Waves (EWAPW)
Authors:
Garry Goldstein
Abstract:
In this work we present a new basis set for electronic structure calculations of crystalline solids using Density Functional theory (DFT) methods. In this construction we take advantage of the fact that most DFT calculations use a convergence loop in order to obtain the eigenstates of a final Khon Sham (KS) Hamiltonian matrix whose eigenstates also give the appropriate electron density needed to o…
▽ More
In this work we present a new basis set for electronic structure calculations of crystalline solids using Density Functional theory (DFT) methods. In this construction we take advantage of the fact that most DFT calculations use a convergence loop in order to obtain the eigenstates of a final Khon Sham (KS) Hamiltonian matrix whose eigenstates also give the appropriate electron density needed to obtain the KS potential needed for that KS Hamiltonian matrix. Here we propose that for the basis of each step of the iteration we use the previous eigenstate basis in the interstitial region but augmented inside the MT sphere with the solution to the spherically averaged KS Hamiltonian for the energy window of that eigenstate. To reduce the number of times the KS potential needs to be solved inside the MT spheres it is advantageous to use energy windows and solve the KS Hamiltonian inside the MT region only once per window (at some energy inside the window) so that the KS Hamiltonian needs only be solved a small number of times per iteration, for practical applications on the order of 10 to 100 windows. This method combines the energy dependence of methods such as Projected Augmented Wave functions (PAW) with the ability of the basis set to adjust to the solid state (rather then atomic) environment of basis sets such as Linearized Augmented Plane Waves (LAPW).
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
Four center integrals for Coulomb interactions in small molecules
Authors:
Garry Goldstein
Abstract:
In this work we make some progress on studying four center integrals for the Coulomb energy for both Hartree Fock (HF) and Density Functional Theory (DFT) calculations for small molecules. We consider basis wave functions of the form of an arbitrary radial wave function multiplied by a spherical harmonic and study four center Coulomb integrals for them. We reformulated these Coulomb four center in…
▽ More
In this work we make some progress on studying four center integrals for the Coulomb energy for both Hartree Fock (HF) and Density Functional Theory (DFT) calculations for small molecules. We consider basis wave functions of the form of an arbitrary radial wave function multiplied by a spherical harmonic and study four center Coulomb integrals for them. We reformulated these Coulomb four center integrals in terms of some derivatives of integrals of nearly factorable functions which then depend on the Bessel transform of the radial wave functions considered.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Projected Augmented Waves (PAW) inspired mixed basis sets for small molecules
Authors:
Garry Goldstein
Abstract:
The success behind many pseudopotential methods such as Projected Augmented Waves (PAW) and Phillips-Kleinman pseudopotential methods is that these methods are all electron methods in disguise. For the Phillips Kleinman and PAW methods we show that there is an all electron formulation. In the all electron formulation there are regular low wavevector wavefunctions and several, specially chosen, hig…
▽ More
The success behind many pseudopotential methods such as Projected Augmented Waves (PAW) and Phillips-Kleinman pseudopotential methods is that these methods are all electron methods in disguise. For the Phillips Kleinman and PAW methods we show that there is an all electron formulation. In the all electron formulation there are regular low wavevector wavefunctions and several, specially chosen, high wavevector wavefunctions that are specialized to the atomic environment. Using this as inspiration here we propose a new, PAW pseudopotential inspired, basis set method for small molecules, where we use the recently introduced LO or lo (Localized orbitals) basis set is paired with a Gaussian basis for a hybrid basis for small molecules.
△ Less
Submitted 16 April, 2024; v1 submitted 12 April, 2024;
originally announced April 2024.
-
Multi-radius Soler-Williams Augmented Plane Waves (SAPWMR), Multi-Radius Soler-Williams Linearized Augmented Plane Waves (SLAPWMR) and extensions
Authors:
Garry Goldstein
Abstract:
In this work we present a new basis set for electronic structures (Density Functional Theory (DFT)) calculations. This basis set extends Soler Williams Linearized Augmented Plane Wave (SLAPW) basis sets by allowing variable Muffin Tin (MT) sphere radii for the different angular momentum channels and for different magnitude wave vectors of the augmented plane waves. With the correct choice of MT ra…
▽ More
In this work we present a new basis set for electronic structures (Density Functional Theory (DFT)) calculations. This basis set extends Soler Williams Linearized Augmented Plane Wave (SLAPW) basis sets by allowing variable Muffin Tin (MT) sphere radii for the different angular momentum channels and for different magnitude wave vectors of the augmented plane waves. With the correct choice of MT radius, this allows us to match additional derivatives of the wave function at the MT radius without having to resort to additional higher derivative terms as part of the wave function expansion inside the MT sphere. This should lead to low wave vector basis set and low linearization energy errors, arguably as good as APW basis set size and LAPW basis linearization errors. We call these basis sets SAPWMR and SLAPWMR depending on the number of derivative like wave functions kept inside the MT radii. Similarly local orbital (LO and lo basis wave functions depending on the number of derivative terms in the MT radius) are suggested with a variable radius. that reduces the number of derivative like terms needed to make them continuous or continuously differentiable. As such semi-core states can also be well handled by our methods. Furthermore in the appendix Full Potential Hamiltonian calculations (FLAPW) are extended to FSLAPWMR full potential calculations for Hamiltonian matrix elements. In the Appendix we introduce some further ideas to improve the speed of DFT calculations which are relevant to the basis sets presented here and to other basis sets such as the LAPW basis set.
△ Less
Submitted 11 April, 2024; v1 submitted 23 March, 2024;
originally announced March 2024.
-
Higher Derivative Muffin Tin Orbitals (HDMTO) and Higher Derivative Koringa Khon and Rostoker (HDKKR) methods
Authors:
Garry Goldstein
Abstract:
In this work we introduce a Linearized version of the Koringa Khon and Rostoker method (LKKR) and show it to be equivalent to the Linearized Muffin Tin Orbitals method (LMTO). We then present higher derivative versions of both methods, e.g. HDKKR and HDMTO and show them to be partially distinct (not equivalent). In particular HDKKR basis set does not have an equivalent ground state for the Khon Sh…
▽ More
In this work we introduce a Linearized version of the Koringa Khon and Rostoker method (LKKR) and show it to be equivalent to the Linearized Muffin Tin Orbitals method (LMTO). We then present higher derivative versions of both methods, e.g. HDKKR and HDMTO and show them to be partially distinct (not equivalent). In particular HDKKR basis set does not have an equivalent ground state for the Khon Sham (KS) Hamiltonian as the HDKKR basis set and has greater variational power than the HDMTO one. Because the KS method, for Density Functional Theory (DFT), is variational HDKKR will give better ground state energies than HDMTO. However HDKKR is much harder to work with then HDMTO requiring much greater computer resources so HDMTO can often be preferred.
△ Less
Submitted 19 March, 2024;
originally announced March 2024.
-
The Wilson-Fisher Fixed point revisited: importance of the form of the cutoff
Authors:
Garry Goldstein
Abstract:
In this work we re-examine the Wilson Fisher fixed point. We study Wilsonian momentum space renormalization group (RG) flow for various forms of the cutoff. We show that already at order $\left(4-d\right)^{1}$, where $d$ is the dimension of the $φ^{4}$ theory, there are changes to the position of the fixed point and the direction of irrelevant coupling parameters. We also show in a multi-flavor…
▽ More
In this work we re-examine the Wilson Fisher fixed point. We study Wilsonian momentum space renormalization group (RG) flow for various forms of the cutoff. We show that already at order $\left(4-d\right)^{1}$, where $d$ is the dimension of the $φ^{4}$ theory, there are changes to the position of the fixed point and the direction of irrelevant coupling parameters. We also show in a multi-flavor $φ^{4}$ model that symmetries of the Lagrange function can be destroyed if the different flavors have different cutoffs (that is the Lagrangian can flow to a non-symmetric fixed point). Some related comments are made about a similar situation in parquet RG (pRG). In future works we will study Wilsonian RG to order $\left(4-d\right)^{2}$ and find non-universal critical exponents that depend on the cutoff.
△ Less
Submitted 28 March, 2024; v1 submitted 9 March, 2024;
originally announced March 2024.
-
The effect of Coulomb assisted hopping on STM signal: extended two site Hubbard model analysis
Authors:
Garry Goldstein
Abstract:
In this work we study STM signal in the presence of Coulomb assisted hopping. We perform an extended two site Hubbard model analysis between the atom on the tip and the atom in the sample nearest to each other. We show that in the presence of Coulomb assisted hopping the STM signal depends on several spectral functions thereby complicating its interpretation. Furthermore in the broadband tip limit…
▽ More
In this work we study STM signal in the presence of Coulomb assisted hopping. We perform an extended two site Hubbard model analysis between the atom on the tip and the atom in the sample nearest to each other. We show that in the presence of Coulomb assisted hopping the STM signal depends on several spectral functions thereby complicating its interpretation. Furthermore in the broadband tip limit there are now three different competing rates for the total current (instead of one for the usual two site Hubbard model analysis used in the literature so far). We find an exact (within the Fermi golden rule - that is in the limit of weak coupling between tip and sample) expression for the current as a function of the bias voltage. As an example we apply our calculations to the case of free fermions with a uniform density of states. Even in this simple case there are non-trivial corrections - where the $dI/d\mathcal{V}$ (the rate of change of the current with respect to bias voltage) is not uniform as a two site (non-extended) Hubbard model analysis would predict. We also show that for realistic conditions the corrections predicted here are order one.
△ Less
Submitted 26 April, 2024; v1 submitted 18 March, 2024;
originally announced March 2024.
-
Order parameters in quasi-1D spin systems
Authors:
Garry Goldstein
Abstract:
In this work we extend the notion of what is meant by a meanfield. Meanfields are approximately maps - through some self consistency relation - of a complex, usually manybody, problem to a simpler more readily solvable problem. This mapping can then be solved to represent properties of the complex many body problem using some self consistency relations. Prototypical examples of simpler meanfield p…
▽ More
In this work we extend the notion of what is meant by a meanfield. Meanfields are approximately maps - through some self consistency relation - of a complex, usually manybody, problem to a simpler more readily solvable problem. This mapping can then be solved to represent properties of the complex many body problem using some self consistency relations. Prototypical examples of simpler meanfield problems (meanfield systems) are the single site and free particle problems. Here we propose a new class of simple meanfield systems where the simple problem to be solved is a 1D spin chain. These meanfields are particularly useful for studying quasi-1D models, where there is a 3D system composed of weakly coupled 1D spin chains with the coupling in the transverse direction weaker than in the 1D direction. We illustrate this idea by considering meanfields for the Ising (of any coupling sign) and ferromagnetic Heisenberg models with one direction coupled much more strongly then the other directions (quasi-1D systems) which map at meanfield level onto the 1D Ising and 1D ferromagnetic Heisenberg models. We also consider more exotic models to illustrate other methods of solving 1D systems, namely the $N$-state Potts model. Magnetic phase transition temperatures and are obtained for all three models, we see that they significantly differ from the usual meanfield estimates. Indeed if the 1D direction has coupling $Γ$ and the transverse directions have coupling $J$ with $λ\sim\fracΓ{J}\gg1$ then regular meanfield would predict the transition temperature to be $k_{B}T_{c}\simΓ$ for all three models while 1D meanfield predicts temperatures of $k_{B}T_{c}\sim\fracΓ{\log\left(λ\right)}$ for the Ising and Potts models and $k_{B}T_{c}\sim\fracΓ{\sqrtλ}$ for the ferromagnetic Heisenberg model. Cluster 1D meanfield extensions are also proposed.
△ Less
Submitted 11 April, 2024; v1 submitted 13 March, 2024;
originally announced March 2024.
-
Classical action for the height within the Kardar-Parisi-Zhang (KPZ) equation for surface growth
Authors:
Garry Goldstein
Abstract:
In this work we write down a classical (not quantum) action for the surface height for the Kardar-Parisi-Zhang (KPZ) equation for surface growth. We do so starting with the regular Martin-Siggia-Rose (MSR) action (which is quantum - contains the constraint field) and integrate out the quantum constraint field exactly. We analyze the classical action, we thereby obtain, within the gaussian and one…
▽ More
In this work we write down a classical (not quantum) action for the surface height for the Kardar-Parisi-Zhang (KPZ) equation for surface growth. We do so starting with the regular Martin-Siggia-Rose (MSR) action (which is quantum - contains the constraint field) and integrate out the quantum constraint field exactly. We analyze the classical action, we thereby obtain, within the gaussian and one loop approximations various instabilities to rough surfaces. The gaussian analysis predicts instabilities to rough surfaces below two dimensions while one loop analysis predicts stronger surface stability to roughness and instabilities to roughness below one dimension rather than two dimensions. The one loop analysis shows that the KPZ action is incomplete and that we generate additional terms (not found in the initial classical action) in the action perturbatively. In the supplement we also modify the KPZ equation for growing surfaces to include the effect of surface tilt on the noise (that is have two sources of noise one of which is multiplicatively coupled noise).
△ Less
Submitted 3 May, 2024; v1 submitted 11 March, 2024;
originally announced March 2024.
-
Time reversal invariant topological 1D and 2D superconductors: doubling the Sau-Luchtin-Tewari-Sarma and Oreg-Refael-von Oppen proposals
Authors:
Garry Goldstein
Abstract:
In this work we present a doubled version of the Sau-Luchtin-Tewari-Sarma and Oreg-Refael-von Open proposals thereby obtaining time reversal invariant p-wave superconductivity in both 1D and 2D. This construction is much like the Kane-Mele spin Hall model which is a time reversal invariant doubling of the Haldane model. We show that the low energy effective action for these doubled versions of the…
▽ More
In this work we present a doubled version of the Sau-Luchtin-Tewari-Sarma and Oreg-Refael-von Open proposals thereby obtaining time reversal invariant p-wave superconductivity in both 1D and 2D. This construction is much like the Kane-Mele spin Hall model which is a time reversal invariant doubling of the Haldane model. We show that the low energy effective action for these doubled versions of the Sau-Luchtin-Tewari-Sarma and Oreg-Refael-von Open models corresponds to a single band p-wave time reversal invariant superconductors with pseudo spin degree of freedom instead of spin degree of freedom. There are Majorana fermions at the ends of wires or in vortex cores of these superconductors. Furthermore these are shown to be stable to small perturbations. In the supplement we present physical realization of the system with cold atoms and show a related "no-go" theorem given in Haim et. al. (2019) has too restrictive assumptions to apply to this proposal.
△ Less
Submitted 24 April, 2024; v1 submitted 3 March, 2024;
originally announced March 2024.
-
Logical Majorana zero modes in a nanowire network
Authors:
Sayandip Dhara,
Garry Goldstein,
Claudio Chamon,
Eduardo R. Mucciolo
Abstract:
We present a scheme to use physical Majorana quasi-zero modes at each junction of a two-dimensional nanowire network to build a logical Majorana zero mode, the location of which is controllable through gate voltages. The wire-network is a way to realize a proposal by Yang et al. arXiv:1808.04825 to imprint a Kekulé vortex pattern on a honeycomb lattice via gate voltages. We show that a specific ty…
▽ More
We present a scheme to use physical Majorana quasi-zero modes at each junction of a two-dimensional nanowire network to build a logical Majorana zero mode, the location of which is controllable through gate voltages. The wire-network is a way to realize a proposal by Yang et al. arXiv:1808.04825 to imprint a Kekulé vortex pattern on a honeycomb lattice via gate voltages. We show that a specific type of junction -- other than a naive Y- or T-junction -- is needed to realize, without breaking time-reversal symmetry, an artificial ``graphene'' system with Majorana fermions instead of complex ones at each site. The junction we propose (i) traps exactly one physical Majorana (quasi-)zero mode at each site of either a brick wall or honeycomb lattice and (ii) allows this mode to hybridize with all three neighboring sites. Using a lattice of these junctions and starting from an electronic, tight-binding model for the wires, we imprint the voltage patterns corresponding to the Kekulé vortex and observe the emergence of the logical Majorana zero mode at the vortex core. We also provide the range of parameters where this excitation could be realized experimentally.
△ Less
Submitted 14 October, 2022;
originally announced October 2022.
-
Lifetime of Excitations in Atomic and Molecular Bose-Einstein Condensates
Authors:
Matteo Bellitti,
Garry Goldstein,
Chris R. Laumann
Abstract:
Recent experimental progress has produced Molecular Superfluids (MSF) in thermal equilibrium; this opens the door to a new class of experiments investigating the associated thermodynamic and dynamical responses. We review the theoretical picture of the phase diagram and quasiparticle spectrum in the Atomic Superfluid (ASF) and MSF phases. We further compute the parametric dependence of the quasipa…
▽ More
Recent experimental progress has produced Molecular Superfluids (MSF) in thermal equilibrium; this opens the door to a new class of experiments investigating the associated thermodynamic and dynamical responses. We review the theoretical picture of the phase diagram and quasiparticle spectrum in the Atomic Superfluid (ASF) and MSF phases. We further compute the parametric dependence of the quasiparticle lifetimes at one-loop order. In the MSF phase, the $U(1)$ particle number symmetry breaks to $\mathbb{Z}_2$ and the spectrum exhibits a gapless Goldstone mode in addition to a gapped $\mathbb{Z}_2$-protected atom-like mode. In the ASF phase, the $U(1)$ symmetry breaks completely, leaving behind a Goldstone mode and an unprotected gapped mode. In both phases, the Goldstone mode decays with a rate given by the celebrated Belyaev result, as in a single component condensate. In the MSF phase, the gapped mode is sharp up to a critical Cherenkov momentum beyond which it emits phonons. In the ASF phase, the gapped mode decays with a constant rate even at small momenta. These decay rates govern the spectral response in microtrap tunneling experiments and lead to sharp features in the transmission spectrum of atoms fired through molecular clouds.
△ Less
Submitted 12 October, 2022;
originally announced October 2022.
-
Semiclassical transport with Berry curvature: Chambers formula and applications to systems with Fermi surface topological transitions
Authors:
Emmanouil K. Kokkinis,
Garry Goldstein,
Dmitry Efremov,
Joseph J. Betouras
Abstract:
Starting with general semiclassical equations of motion for electrons in the presence of electric and magnetic fields, we extend the Chambers formula to include in addition to a magnetic field, time-dependent electric fields and bands with Berry curvature. We thereby compute the conductivity tensor $σ_{αβ}(B,ω)$ in the presence of magnetic field for bands in two (2D) and three (3D) dimensions with…
▽ More
Starting with general semiclassical equations of motion for electrons in the presence of electric and magnetic fields, we extend the Chambers formula to include in addition to a magnetic field, time-dependent electric fields and bands with Berry curvature. We thereby compute the conductivity tensor $σ_{αβ}(B,ω)$ in the presence of magnetic field for bands in two (2D) and three (3D) dimensions with Berry curvature. We focus then on several applications to magnetotransport for metals with Fermi surface topological transitions in 2D. In particular, we consider a rectangular lattice and a model related to overdoped graphene, to investigate the signatures of different types of Fermi surface topological transitions in metals in the Hall coefficient, Hall conductivity $σ_{xy}$ and longitudinal conductivity $σ_{xx}$. The behavior of those quantities as a function of frequency, when the electric field is time dependent, is also investigated. As an example of non-zero Berry curvature, we study the magnetotransport of the Haldane model within this context. In addition, we provide the linear and nonlinear electric current formula to order $E^2$.
△ Less
Submitted 6 April, 2022; v1 submitted 17 November, 2021;
originally announced November 2021.
-
Extending the Gutzwiller approximation to intersite interactions
Authors:
Garry Goldstein,
Nicola Lanata,
Gabriel Kotliar
Abstract:
We develop an extension of the Gutzwiller Approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping and Coulomb assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multi-band Hubbard models in the presence of only local interactions. This is accomplished by combining the 1/z…
▽ More
We develop an extension of the Gutzwiller Approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping and Coulomb assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multi-band Hubbard models in the presence of only local interactions. This is accomplished by combining the 1/z expansion -- where z is the coordination number and only the leading order terms contribute in the limit of infinite dimensions -- with a P^2-I expansion, where P is the Gutzwiller projector. The method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single band Hubbard model. Similarly to the usual Brinkman-Rice mechanism we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wavefunction is enhanced with respect to the uncorrelated Slater determinant wavefunction.
△ Less
Submitted 1 June, 2020; v1 submitted 5 January, 2020;
originally announced January 2020.
-
Generalised route to effective field theories for quantum systems with local constraints
Authors:
Attila Szabó,
Garry Goldstein,
Claudio Castelnovo,
Alexei M. Tsvelik
Abstract:
Some of the exciting phenomena uncovered in strongly correlated systems in recent years - for instance quantum topological order, deconfined quantum criticality, and emergent gauge symmetries -- appear in systems in which the Hilbert space is effectively projected at low energies in a way that imposes local constraints on the original degrees of freedom. Cases in point include spin liquids, valenc…
▽ More
Some of the exciting phenomena uncovered in strongly correlated systems in recent years - for instance quantum topological order, deconfined quantum criticality, and emergent gauge symmetries -- appear in systems in which the Hilbert space is effectively projected at low energies in a way that imposes local constraints on the original degrees of freedom. Cases in point include spin liquids, valence bond systems, dimer models, and vertex models. In this work, we use a slave boson description coupled to a large-$S$ path integral formulation to devise a generalised route to obtain effective field theories for such systems. We demonstrate the validity and capability of our approach by studying quantum dimer models and by comparing our results with the existing literature. Field-theoretic approaches to date are limited to bipartite lattices, they depend on a gauge-symmetric understanding of the constraint, and they lack generic quantitative predictive power for the coefficients of the terms that appear in the Lagrangians of these systems. Our method overcomes all these shortcomings and we show how the results up to quadratic order compare with the known height description of the square lattice quantum dimer model, as well as with the numerical estimate of the speed of light of the photon excitations on the diamond lattice. Finally, instanton considerations allow us to infer properties of the finite-temperature behaviour in two dimensions.
△ Less
Submitted 5 August, 2019; v1 submitted 29 April, 2019;
originally announced April 2019.
-
Steady-state superconductivity in electronic materials with repulsive interactions
Authors:
Oliver Hart,
Garry Goldstein,
Claudio Chamon,
Claudio Castelnovo
Abstract:
We study the effect of laser driving on a minimal model for a hexagonal two-dimensional material with broken inversion symmetry. Through the application of circularly polarised light and coupling to a thermal free electron bath, the system is driven into a nonequilibrium steady state with asymmetric, nonthermal carrier populations in the two valleys. We show that, in this steady state, interband s…
▽ More
We study the effect of laser driving on a minimal model for a hexagonal two-dimensional material with broken inversion symmetry. Through the application of circularly polarised light and coupling to a thermal free electron bath, the system is driven into a nonequilibrium steady state with asymmetric, nonthermal carrier populations in the two valleys. We show that, in this steady state, interband superconducting correlations between electrons can develop independent of the sign of the electron-electron interactions. We discuss how our results apply, for example, to transition metal dichalcogenides. This work opens the door to technological applications of superconductivity in a range of materials that were hitherto precluded from it.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
Hall coefficient and magnetoresistance in boson+fermion dimer models for the pseudogap phase of high Tc superconductors
Authors:
Garry Goldstein
Abstract:
We show that the Hall coefficient of the boson+fermion dimer model for the pseudogap phase of high temperature superconductivity introduced in Punk et al. changes sign from negative at low temperatures to positive at high temperatures at a characteristic temperature scaling proportional to the the cyclotron frequency of the fermionic dimers (divided by k), (here k~0.7 fits the experimental data we…
▽ More
We show that the Hall coefficient of the boson+fermion dimer model for the pseudogap phase of high temperature superconductivity introduced in Punk et al. changes sign from negative at low temperatures to positive at high temperatures at a characteristic temperature scaling proportional to the the cyclotron frequency of the fermionic dimers (divided by k), (here k~0.7 fits the experimental data well). We show that this is an effect of the changing of the sign of the coupling between the fermionic dimer and the magnetic field from negative coupling -e at low temperatures to positive coupling +e at high temperature, with the Hall coefficient being proportional to R_H~e_B*e_E*e_J (the product of the magnetic charge, electric charge and current charge all of which we carefully define). We relate the Hall conductivity to the coefficient in Kohler's like rule for magnetoconductivity and calculate some corrections which are relevant near the intermediate temperature range ~50K (typical values for the cyclotron frequency of the dimers). Furthermore we make a sharp prediction that the magnetoresistance effect vanishes to order B^2 at the temperature and magnetic field where the Hall coefficient vanishes.
△ Less
Submitted 1 April, 2018;
originally announced April 2018.
-
Quantum oscillations and criticality in a fermionic and bosonic dimer model for the cuprates
Authors:
Garry Goldstein,
Claudio Chamon,
Claudio Castelnovo
Abstract:
We study quantum oscillations for a system of fermionic and bosonic dimers and compare the results to those experimentally observed in the cuprate superconductors in their underdoped regime. Based on gauge invariance, we argue that the charge carriers obey the Onsager quantization condition and quantum oscillations take on a Lifshitz-Kosevich form. We obtain the effective mass and find good qualit…
▽ More
We study quantum oscillations for a system of fermionic and bosonic dimers and compare the results to those experimentally observed in the cuprate superconductors in their underdoped regime. Based on gauge invariance, we argue that the charge carriers obey the Onsager quantization condition and quantum oscillations take on a Lifshitz-Kosevich form. We obtain the effective mass and find good qualitative agreement with experiments if we tune the model to the point where the observed mass divergence at optimum doping is associated to a van Hove singularity at which four free-dimer Fermi pockets touch pairwise in the interior of the Brillouin zone. The same van Hove singularity leads to a maximum in the d-wave superconducting pairing amplitude when anti-ferromagnetic interactions are included. Our combined results therefore suggest that a quantum critical point separating the underdoped and overdoped regimes is marked by the location of the van Hove saddle point in the fermionic dimer dispersion.
△ Less
Submitted 22 February, 2019; v1 submitted 6 March, 2018;
originally announced March 2018.
-
Superfluid density of a photo-induced superconducting state
Authors:
A. Shtyk,
G. Goldstein,
C. Aron,
C. Chamon
Abstract:
Nonequilibrium conditions offer novel routes to superconductivity that are not available at equilibrium. For example, by engineering nonequilibrium electronic populations, pairing may develop between electrons in different energy bands. A concrete proposal has been made to photo-induce superconductivity in a semiconductor, where pairing occurs between electrons in the conduction and valence bands,…
▽ More
Nonequilibrium conditions offer novel routes to superconductivity that are not available at equilibrium. For example, by engineering nonequilibrium electronic populations, pairing may develop between electrons in different energy bands. A concrete proposal has been made to photo-induce superconductivity in a semiconductor, where pairing occurs between electrons in the conduction and valence bands, even for repulsive interactions. Here, we calculate the superfluid density for such a nonequilibrium paired state, and find it to be positive for repulsive interactions and interband pairing. The positivity of the superfluid density implies the stability of the photo-induced superconducting state as well as the existence of the Meissner effect.
△ Less
Submitted 7 February, 2018;
originally announced February 2018.
-
Electrons at the monkey saddle: a multicritical Lifshitz point
Authors:
Alex Shtyk,
Garry Goldstein,
Claudio Chamon
Abstract:
We consider 2D interacting electrons at a monkey saddle with dispersion $\propto p_x^3-3p_xp_y^2$. Such a dispersion naturally arises at the multicritical Lifshitz point when three van Hove saddles merge in an elliptical umbilic elementary catastrophe, which we show can be realized in biased bilayer graphene. A multicritical Lifshitz point of this kind can be identified by its signature Landau lev…
▽ More
We consider 2D interacting electrons at a monkey saddle with dispersion $\propto p_x^3-3p_xp_y^2$. Such a dispersion naturally arises at the multicritical Lifshitz point when three van Hove saddles merge in an elliptical umbilic elementary catastrophe, which we show can be realized in biased bilayer graphene. A multicritical Lifshitz point of this kind can be identified by its signature Landau level behavior $E_m\propto (Bm)^{3/2}$ and related oscillations in thermodynamic and transport properties, such as de Haas-van Alphen and Shubnikov-de Haas oscillations, whose period triples as the system crosses the singularity. We show, in the case of a single monkey saddle, that the non-interacting electron fixed point is unstable to interactions under the renormalization group flow, developing either a superconducting instability or non-Fermi liquid features. Biased bilayer graphene, where there are two non-nested monkey saddles at the $K$ and $K^\prime$ points, exhibits an interplay of competing many-body instabilities, namely s-wave superconductivity, ferromagnetism, and spin- and charge-density wave.
△ Less
Submitted 15 June, 2016;
originally announced June 2016.
-
D-wave superconductivity in boson+fermion dimer models
Authors:
Garry Goldstein,
Claudio Chamon,
Claudio Castelnovo
Abstract:
We present a slave-particle mean-field study of the mixed boson+fermion quantum dimer model introduced by Punk, Allais, and Sachdev [PNAS 112, 9552 (2015)] to describe the physics of the pseudogap phase in cuprate superconductors. Our analysis naturally leads to four charge e fermion pockets whose total area is equal to the hole doping p, for a range of parameters consistent with the t-J model for…
▽ More
We present a slave-particle mean-field study of the mixed boson+fermion quantum dimer model introduced by Punk, Allais, and Sachdev [PNAS 112, 9552 (2015)] to describe the physics of the pseudogap phase in cuprate superconductors. Our analysis naturally leads to four charge e fermion pockets whose total area is equal to the hole doping p, for a range of parameters consistent with the t-J model for high temperature superconductivity. Here we find that the dimers are unstable to d-wave superconductivity at low temperatures. The region of the phase diagram with d-wave rather than s-wave superconductivity matches well with the appearance of the four fermion pockets. In the superconducting regime, the dispersion contains eight Dirac cones along the diagonals of the Brillouin zone.
△ Less
Submitted 13 June, 2016;
originally announced June 2016.
-
Correlation functions of the Lieb-Liniger gas and the LeClair-Mussardo formula
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
In this letter we derive formulas for multi point correlation functions, in the thermodynamic limit, for the Lieb Liniger gas taken with respect to arbitrary eigenstates. These results apply for the ground state, thermal states and GGE states. We obtain these correlation functions as a series of multiple integrals of progressively higher dimensions. These integrals converge rapidly for short dista…
▽ More
In this letter we derive formulas for multi point correlation functions, in the thermodynamic limit, for the Lieb Liniger gas taken with respect to arbitrary eigenstates. These results apply for the ground state, thermal states and GGE states. We obtain these correlation functions as a series of multiple integrals of progressively higher dimensions. These integrals converge rapidly for short distance correlation functions and low densities of particles. The series derived matches exactly the LeClair Mussardo formula for correlation functions of relativistic integrable models.
△ Less
Submitted 30 June, 2015;
originally announced July 2015.
-
Equilibration and GGE for hard wall boundary conditions
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
In this work we present an analysis of a quench for the repulsive Lieb-Liniger gas confined to a large box with hard wall boundary conditions. We study the time average of local correlation functions and show that both the quench action approach and the GGE formalism are applicable for the long time average of local correlation functions. We find that the time average of the system corresponds to…
▽ More
In this work we present an analysis of a quench for the repulsive Lieb-Liniger gas confined to a large box with hard wall boundary conditions. We study the time average of local correlation functions and show that both the quench action approach and the GGE formalism are applicable for the long time average of local correlation functions. We find that the time average of the system corresponds to an eigenstate of the Lieb-Liniger Hamiltonian and that this eigenstate is related to an eigenstate of a Lieb-Liniger Hamiltonian with periodic boundary conditions on an interval of twice the length and with twice as many particles (a doubled system). We further show that local operators with support far away from the boundaries of the hard wall have the same expectation values with respect to this eigenstate as corresponding operators for the doubled system. We present an example of a quench where the gas is initially confined in several moving traps and then released into a bigger container, an approximate description of the Newton cradle experiment. We calculate the time average of various correlation functions for long times after the quench.
△ Less
Submitted 11 May, 2015;
originally announced May 2015.
-
Driven-dissipative ising model: mean field solution
Authors:
Garry Goldstein,
Camille Aron,
Claudio Chamon
Abstract:
We study the fate of the Ising model and its universal properties when driven by a rapid periodic drive and weakly coupled to a bath at equilibrium. The far from equilibrium steady-state regime of the system is accessed by means of a Floquet mean field approach. We show that, depending on the details of the bath, the drive can strongly renormalize the critical temperature to higher temperatures, m…
▽ More
We study the fate of the Ising model and its universal properties when driven by a rapid periodic drive and weakly coupled to a bath at equilibrium. The far from equilibrium steady-state regime of the system is accessed by means of a Floquet mean field approach. We show that, depending on the details of the bath, the drive can strongly renormalize the critical temperature to higher temperatures, modify the critical exponents, or even change the nature of the phase transition from second to first order after the emergence of a tricritical point. Moreover, by judiciously selecting the frequency of the field and by engineering the spectrum of the bath, one can drive a ferromagnetic Hamiltonian to an antiferromagnetically ordered phase and vice-versa.
△ Less
Submitted 10 February, 2015;
originally announced February 2015.
-
Band-edge superconductivity
Authors:
Garry Goldstein,
Camille Aron,
Claudio Chamon
Abstract:
We show that superconductivity can arise in semiconductors with a band in the shape of a Mexican hat when the chemical potential is tuned close to the band edge, but not intersecting the band, as long as interactions are sufficiently strong. Hence, this is an example where superconductivity can emerge from a band insulator when interactions exceed a threshold. Semiconductors with simple cubic symm…
▽ More
We show that superconductivity can arise in semiconductors with a band in the shape of a Mexican hat when the chemical potential is tuned close to the band edge, but not intersecting the band, as long as interactions are sufficiently strong. Hence, this is an example where superconductivity can emerge from a band insulator when interactions exceed a threshold. Semiconductors with simple cubic symmetry point groups and with strong spin-orbit coupling provide an example of a system with such band dispersion.
△ Less
Submitted 9 May, 2016; v1 submitted 22 October, 2014;
originally announced October 2014.
-
Stron eigenstate thermalization hypothesis
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
We present a generalization of the ETH conjecture. Using this generalization we are able to derive the fact that an arbitrary eigenstate of a general many body system may be used to represent microcanonical ensemble in any many body experiment that involves only local operators and projectors onto eigenstates of the system Hamiltonian. In particular we extend the ETH to include some non-local oper…
▽ More
We present a generalization of the ETH conjecture. Using this generalization we are able to derive the fact that an arbitrary eigenstate of a general many body system may be used to represent microcanonical ensemble in any many body experiment that involves only local operators and projectors onto eigenstates of the system Hamiltonian. In particular we extend the ETH to include some non-local operators. We present a derivation of this conjecture in the case of a many body model whose Hamiltonian is composed of two parts: an integrable Hamiltonian and a small but finite Gaussian perturbation.
△ Less
Submitted 15 August, 2014;
originally announced August 2014.
-
Photo-induced superconductivity in semiconductors
Authors:
Garry Goldstein,
Camille Aron,
Claudio Chamon
Abstract:
We show that optically pumped semiconductors can exhibit superconductivity. We illustrate this phenomenon in the case of a two-band semiconductor tunnel-coupled to broad-band reservoirs and driven by a continuous wave laser. More realistically, we also show that superconductivity can be induced in a two-band semiconductor interacting with a broad-spectrum light source. We furthermore discuss the c…
▽ More
We show that optically pumped semiconductors can exhibit superconductivity. We illustrate this phenomenon in the case of a two-band semiconductor tunnel-coupled to broad-band reservoirs and driven by a continuous wave laser. More realistically, we also show that superconductivity can be induced in a two-band semiconductor interacting with a broad-spectrum light source. We furthermore discuss the case of a three-band model in which the middle band replaces the broad-band reservoirs as the source of dissipation. In all three cases, we derive the simple conditions on the band structure, electron-electron interaction, and hybridization to the reservoirs that enable superconductivity. We compute the finite superconducting gap and argue that the mechanism can be induced through both attractive and repulsive interactions and is robust to high temperatures.
△ Less
Submitted 19 May, 2015; v1 submitted 27 June, 2014;
originally announced June 2014.
-
Quench between a Mott insulator and a Lieb-Liniger liquid
Authors:
G. Goldstein,
N. Andrei
Abstract:
In this work we study a quench between a Mott insulator and a repulsive Lieb-Liniger liquid. We find explicitly the stationary state when a long time has passed after the quench. It is given by a GGE density matrix which we completely characterize, calculating the quasiparticle density describing the system after the quench. In the long time limit we find an explicit form for the local three body…
▽ More
In this work we study a quench between a Mott insulator and a repulsive Lieb-Liniger liquid. We find explicitly the stationary state when a long time has passed after the quench. It is given by a GGE density matrix which we completely characterize, calculating the quasiparticle density describing the system after the quench. In the long time limit we find an explicit form for the local three body density density density correlation function and the asymptotic long distance limit of the density density correlation function. The later is shown to have a Gaussian decay at large distances.
△ Less
Submitted 18 June, 2014;
originally announced June 2014.
-
How to experimentally detect a GGE? - Universal Spectroscopic Signatures of the GGE in the Tonks gas
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
In this work we study the properties of the density density correlation function of the 1-D Lieb-Liniger model with infinite repulsion in the GGE regime. The GGE describes the equilibrated system in the long time limit after a quench from a generic initial state. In the case that the initial and hence the final state has low entropy per particle we find that the density density correlation functio…
▽ More
In this work we study the properties of the density density correlation function of the 1-D Lieb-Liniger model with infinite repulsion in the GGE regime. The GGE describes the equilibrated system in the long time limit after a quench from a generic initial state. In the case that the initial and hence the final state has low entropy per particle we find that the density density correlation function has a universal form, in particular it depends on a few parameters corresponding to "key" momenta and has power law dependence on the distance. This provides an experimental signature of the GGE which may readily be identified through spectroscopy. These signatures are universal and robust to initial sate preparation.
△ Less
Submitted 17 June, 2014; v1 submitted 25 May, 2014;
originally announced May 2014.
-
Failure of the GGE hypothesis for integrable models with bound states
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
In this work we study the applicability of the GGE to integrable one dimensional systems with bound states. We find that the GGE fails to describe the long time dynamics for most initial states including eigenstates. We present our calculations studying the attractive Lieb-Liniger gas and the XXZ magnet, though similar results may be obtained for other models.
In this work we study the applicability of the GGE to integrable one dimensional systems with bound states. We find that the GGE fails to describe the long time dynamics for most initial states including eigenstates. We present our calculations studying the attractive Lieb-Liniger gas and the XXZ magnet, though similar results may be obtained for other models.
△ Less
Submitted 27 March, 2016; v1 submitted 16 May, 2014;
originally announced May 2014.
-
Equilibration and Generalized GGE in the Lieb Liniger gas
Authors:
Garry Goldstein,
Natan Andrei
Abstract:
We study the nonequilibrium properties of the one dimensional Lieb Liniger model in the finite repulsion regime. Introducing a new version of the Yudson representation applicable to finite size systems and appropriately taking the infinite volume limit we are able to study equilibration of the Lieb Liniger gas in the thermodynamic limit. We provide a formalism to compute various correlation functi…
▽ More
We study the nonequilibrium properties of the one dimensional Lieb Liniger model in the finite repulsion regime. Introducing a new version of the Yudson representation applicable to finite size systems and appropriately taking the infinite volume limit we are able to study equilibration of the Lieb Liniger gas in the thermodynamic limit. We provide a formalism to compute various correlation functions for highly non equilibrium initial states. We are able to find explicit analytic expressions for the long time limit of the expectation of the density, density density and related correlation functions. We show that the gas equilibriates to a diagonal ensemble which we show is equivalent to a generalized version of the GGE for sufficiently simple correlation functions, which in particular include density correlations.
△ Less
Submitted 15 May, 2014; v1 submitted 13 September, 2013;
originally announced September 2013.
-
Exact Zero Modes in Closed Systems of Interacting Fermions
Authors:
Garry Goldstein,
Claudio Chamon
Abstract:
We show that for closed finite sized systems with an odd number of real fermionic modes, even in the presence of interactions, there are always at least two fermionic operators that commute with the Hamiltonian.There is a zero mode corresponding to the fermion parity operator, as shown by Akhmerov, as well as additional linearly independent zero modes, one of which is 1) the one that is continuous…
▽ More
We show that for closed finite sized systems with an odd number of real fermionic modes, even in the presence of interactions, there are always at least two fermionic operators that commute with the Hamiltonian.There is a zero mode corresponding to the fermion parity operator, as shown by Akhmerov, as well as additional linearly independent zero modes, one of which is 1) the one that is continuously connected to the Majorana mode solution in the non-interacting limit, and 2) less prone to decoherence when the system is opened to contact with an infinite bath. We also show that in the idealized situation where there are two or more well separated zero modes each associated with a finite number of fermions at a localized vortex, these modes have non-Abelian Ising statistics under braiding. Furthermore the algebra of the zero mode operators makes them useful for fermionic quantum computation.
△ Less
Submitted 26 August, 2012; v1 submitted 8 August, 2011;
originally announced August 2011.
-
Decay rates for topological memories encoded with Majorana fermions
Authors:
Garry Goldstein,
Claudio Chamon
Abstract:
Recently there have been numerous proposals to create Majorana zero modes in solid state heterojunctions, superconducting wires and optical lattices. Putatively the information stored in qubits constructed from these modes is protected from various forms of decoherence. Here we show that this is not the case. We present a generic method to study the effect of external perturbations on these modes.…
▽ More
Recently there have been numerous proposals to create Majorana zero modes in solid state heterojunctions, superconducting wires and optical lattices. Putatively the information stored in qubits constructed from these modes is protected from various forms of decoherence. Here we show that this is not the case. We present a generic method to study the effect of external perturbations on these modes. We focus on the case where there are no interactions between different Majorana modes either directly or through intermediary fermions. To quantify the rate of loss of the information stored in the Majorana modes we study the two-time correlators for qubits built from them. We analyze a generic gapped fermionic environment (bath) interacting via tunneling with different components of the qubit (different Majorana modes). We find that static (time-independent) perturbations are not harmful, and just lead to a depletion of the information stored that saturates in the long-time limit. On the other hand we find markedly different results for dynamic (time dependent) systems, with no coherence surviving in the long time limit. In particular we find that dephasing and energy fluctuations of the modes in the fermionic bath have adverse effects on the Majorana memories.
△ Less
Submitted 17 November, 2011; v1 submitted 1 July, 2011;
originally announced July 2011.
-
Imaging mesoscopic nuclear spin noise with a diamond magnetometer
Authors:
C. A. Meriles,
L. Jiang,
G. Goldstein,
J. S. Hodges,
J. R. Maze,
M. D. Lukin,
P. Cappellaro
Abstract:
Magnetic Resonance Imaging (MRI) can characterize and discriminate among tissues using their diverse physical and biochemical properties. Unfortunately, submicrometer screening of biological specimens is presently not possible, mainly due to lack of detection sensitivity. Here we analyze the use of a nitrogen-vacancy center in diamond as a magnetic sensor for nanoscale nuclear spin imaging and spe…
▽ More
Magnetic Resonance Imaging (MRI) can characterize and discriminate among tissues using their diverse physical and biochemical properties. Unfortunately, submicrometer screening of biological specimens is presently not possible, mainly due to lack of detection sensitivity. Here we analyze the use of a nitrogen-vacancy center in diamond as a magnetic sensor for nanoscale nuclear spin imaging and spectroscopy. We examine the ability of such a sensor to probe the fluctuations of the "classical" dipolar field due to a large number of neighboring nuclear spins in a densely protonated sample. We identify detection protocols that appropriately take into account the quantum character of the sensor and find a signal-to-noise ratio compatible with realistic experimental parameters. Through various example calculations we illustrate different kinds of image contrast. In particular, we show how to exploit the comparatively long nuclear spin correlation times to reconstruct a local, high-resolution sample spectrum.
△ Less
Submitted 29 April, 2010;
originally announced April 2010.