Showing 1–50 of 51 results for author: Santos, R R d

Four interacting spins: addition of angular momenta, spin-spin correlation functions, and entanglement

Authors: Raimundo R. dos Santos, Lucas Alves Oliveira, Natanael C. Costa

Abstract: We study four spins on a ring coupled through competing Heisenberg interactions between nearest neighbors, $J$, and next-nearest neighbors, $J_2\equivαあるふぁJ>0$. The spectrum is obtained in a simple way by using the rules for addition of 4 angular momenta. This allows us to follow the evolution of the ground state with $αあるふぁ$, characterized by level crossings and by analyses of spin-spin correlation funct… ▽ More We study four spins on a ring coupled through competing Heisenberg interactions between nearest neighbors, $J$, and next-nearest neighbors, $J_2\equivαあるふぁJ>0$. The spectrum is obtained in a simple way by using the rules for addition of 4 angular momenta. This allows us to follow the evolution of the ground state with $αあるふぁ$, characterized by level crossings and by analyses of spin-spin correlation functions. Further insight is obtained by examining the entanglement between different parts of the system: we observe that the entanglement entropy is strongly dependent on how the system is partitioned. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 10 pages, 5 figures

Magnetism and metal-insulator transitions in the anisotropic kagome lattice

Authors: Lucas O. Lima, Andressa R. Medeiros-Silva, Raimundo R. dos Santos, Thereza Paiva, Natanael C. Costa

Abstract: The interest in the physical properties of kagome lattices has risen considerably. In addition to the synthesis of new materials, the possibility of realizing ultracold atoms on an optical kagome lattice (KL) raises interesting issues. For instance, by considering the Hubbard model on an anisotropic KL, with a hopping $t^\prime$ along one of the directions, one is able to interpolate between the L… ▽ More The interest in the physical properties of kagome lattices has risen considerably. In addition to the synthesis of new materials, the possibility of realizing ultracold atoms on an optical kagome lattice (KL) raises interesting issues. For instance, by considering the Hubbard model on an anisotropic KL, with a hopping $t^\prime$ along one of the directions, one is able to interpolate between the Lieb lattice ($t^\prime=0$) and the isotropic KL ($t^\prime=t$). The ground state of the former is a ferrimagnetic insulator for any on-site repulsion, $U$, while the latter displays a transition between a paramagnetic metal and a Mott insulator. One may thus consider $t^\prime$ as a parameter controlling the degree of magnetic frustration in the system. By means of extensive quantum Monte Carlo simulations, we have examined magnetic and transport properties as $t^\prime$ varies between these limits in order to set up a phase diagram in the $(U/t, t^\prime/t)$ parameter space. As an auxiliary response, analysis of the average sign of the fermionic determinant provides consistent predictions for critical points in the phase diagram. We observe a metal-insulator transition occurring at some critical point $U_c^\text{M}(t^\prime)$, which increases monotonically with $ t^\prime $, from the unfrustrated lattice limit. In addition, we have found that the boundary between the ferrimagnetic insulator and the Mott insulator rises sharply with $t^\prime$. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 8 pages, 7 figures

The half-filled extended Hubbard model on a square lattice: Phase boundaries from determinant quantum Monte Carlo simulations

Authors: Sebastião dos Anjos Sousa-Júnior, Natanael C. Costa, Raimundo R. dos Santos

Abstract: The extended Hubbard model (EHM) describes fermions on a lattice coupled through on-site, $U$, and first-neighbor, $V$, interactions. In the context of high-$T_c$ cuprates, antiferromagnetic fluctuations may lead to an attractive channel, hence to superconductivity. Despite interest in the two-dimensional version of the model, the current knowledge about the phase diagram is still far from complet… ▽ More The extended Hubbard model (EHM) describes fermions on a lattice coupled through on-site, $U$, and first-neighbor, $V$, interactions. In the context of high-$T_c$ cuprates, antiferromagnetic fluctuations may lead to an attractive channel, hence to superconductivity. Despite interest in the two-dimensional version of the model, the current knowledge about the phase diagram is still far from complete. Here, we report on the results of extensive determinant quantum Monte Carlo simulations for this model at half filling, in which we have used the average sign of the product of fermionic determinants as an additional observable to locate critical points. We arrive at a ground state phase diagram in the $U$-$V$ plane in which the boundaries involving antiferromagnetic, charge-ordered, $s$- and $d$-wave superconductivity, and phase-separated phases are quantitatively set with good accuracy. We have also proposed a partial phase diagram, $T_c(U,V)$, featuring critical temperatures for the CDW and $s$-wave superconducting phases. △ Less

Submitted 11 February, 2024; v1 submitted 17 April, 2023; originally announced April 2023.

Magnetic impurities in a charge-ordered background

Authors: Sebastião dos Anjos Sousa-Júnior, Raimundo R. dos Santos, Natanael C. Costa

Abstract: We investigate how magnetic impurities may affect a system exhibiting charge-density wave (CDW) in its ground state. We consider a disordered Hubbard-Holstein model with a homogeneous electron-phonon interaction, but with a (randomly chosen) fraction of sites displaying a non-zero Coulomb repulsion, $U$, and perform state-of-the-art finite-temperature quantum Monte Carlo simulations. For a single… ▽ More We investigate how magnetic impurities may affect a system exhibiting charge-density wave (CDW) in its ground state. We consider a disordered Hubbard-Holstein model with a homogeneous electron-phonon interaction, but with a (randomly chosen) fraction of sites displaying a non-zero Coulomb repulsion, $U$, and perform state-of-the-art finite-temperature quantum Monte Carlo simulations. For a single magnetic impurity, charge-charge correlations hamper the spin-spin ones around the repulsive site, thus requiring a strong enough value of $U$ to create non-negligible antiferromagnetic (AFM) correlations. As the number of magnetic impurities increases, these AFM correlations become deleterious to CDW order and its features. First, the critical temperature is drastically reduced, and seems to vanish around 40$\%$ of impurities (for fixed $U/λらむだ=2$), which we correlate with the classical percolation threshold. We also notice that just a small amount of disorder suffices to create a \textit{bad insulating} state, with the suppression of both Peierls and spin gaps, even within the charge-ordered phase. Finally, we have also found that pairing correlations are enhanced at large doping, driven by the competition between CDW and AFM tendencies. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Magnetism and metal-insulator transitions in the Rashba-Hubbard model

Authors: Welberth Kennedy, Sebastião dos Anjos Sousa Júnior, Natanael C. Costa, Raimundo R. dos Santos

Abstract: The nature of metal-insulator and magnetic transitions is still a subject under intense debate in condensed matter physics. Amongst the many possible mechanisms, the interplay between electronic correlations and spin-orbit couplings is an issue of a great deal of interest, in particular when dealing with quasi-2D compounds. In view of this, here we use a Hartree-Fock approach to investigate how th… ▽ More The nature of metal-insulator and magnetic transitions is still a subject under intense debate in condensed matter physics. Amongst the many possible mechanisms, the interplay between electronic correlations and spin-orbit couplings is an issue of a great deal of interest, in particular when dealing with quasi-2D compounds. In view of this, here we use a Hartree-Fock approach to investigate how the Rashba spin-orbit coupling, $V_\text{SO}$, affects the magnetic ordering provided by a Hubbard interaction, $U$, on a square lattice. At half-filling, we have found a sequence of transitions for increasing $V_\text{SO}$: from a Mott insulator to a metallic antiferromagnet, and then to a paramagnetic Rashba metal. Also, our results indicate that the Rashba coupling favors magnetic striped phases in the doped regime. By analyzing spectral properties, we associate the rearrangement of the magnetic ordering with the emerging chirality created by the spin-orbit coupling. Our findings provide insights towards clarifying the competition between these tendencies. △ Less

Submitted 17 May, 2022; originally announced May 2022.

The 2D attractive Hubbard model and the BCS-BEC crossover

Authors: Rodrigo A. Fontenele, Natanael C. Costa, Raimundo R. dos Santos, Thereza Paiva

Abstract: Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Mo… ▽ More Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature, $T_c$, in terms of the band filling, $\ave{n}$, and interaction strength, $U$, from which we pinpoint a somewhat wide region $|U|/t \approx 5 \pm 1$ ($t$ is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum $T_c \approx 0.16 t$. Two additional temperature scales, namely pairing, $T_p$, and degeneracy, $T_d$, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Comments: 9 pages, 12 figures

Unsupervised machine learning approaches to the $q$-state Potts model

Authors: Andrea Tirelli, Danyella O. Carvalho, Lucas A. Oliveira, J. P. Lima, Natanael C. Costa, Raimundo R. dos Santos

Abstract: In this paper with study phase transitions of the $q$-state Potts model, through a number of unsupervised machine learning techniques, namely Principal Component Analysis (PCA), $k$-means clustering, Uniform Manifold Approximation and Projection (UMAP), and Topological Data Analysis (TDA). Even though in all cases we are able to retrieve the correct critical temperatures $T_c(q)$, for $q = 3, 4$ a… ▽ More In this paper with study phase transitions of the $q$-state Potts model, through a number of unsupervised machine learning techniques, namely Principal Component Analysis (PCA), $k$-means clustering, Uniform Manifold Approximation and Projection (UMAP), and Topological Data Analysis (TDA). Even though in all cases we are able to retrieve the correct critical temperatures $T_c(q)$, for $q = 3, 4$ and $5$, results show that non-linear methods as UMAP and TDA are less dependent on finite size effects, while still being able to distinguish between first and second order phase transitions. This study may be considered as a benchmark for the use of different unsupervised machine learning algorithms in the investigation of phase transitions. △ Less

Submitted 18 March, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Added computation of critical exponents; exposition improved

Percolation on Lieb lattices

Authors: W. S. Oliveira, J. Pimentel de Lima, Natanael C. Costa, R. R. dos Santos

Abstract: We study site- and bond-percolation on a class of lattices referred to as Lieb lattices. In two dimensions the Lieb lattice (LL) is also known as the decorated square lattice, or as the CuO$_2$ lattice; in three dimensions it can be generalized to a layered Lieb lattice (LLL) or to a perovskite lattice (PL). Emergent electronic phenomena, such as topological states and ferrimagnetism, have been pr… ▽ More We study site- and bond-percolation on a class of lattices referred to as Lieb lattices. In two dimensions the Lieb lattice (LL) is also known as the decorated square lattice, or as the CuO$_2$ lattice; in three dimensions it can be generalized to a layered Lieb lattice (LLL) or to a perovskite lattice (PL). Emergent electronic phenomena, such as topological states and ferrimagnetism, have been predicted to occur in these systems, which may be realized in optical lattices as well as in solid state. Since the study of the interplay between quantum fluctuations and disorder in these systems requires the availability of accurate estimates of geometrical critical parameters, such as percolation thresholds and correlation length exponents, here we use Monte Carlo simulations to obtain these data for Lieb lattices when a site (or bond) is present with probability $p$. We have found that the thresholds satisfy a mean-field (Bethe lattice) trend, namely that the critical concentration, $p_c$, increases as the average coordination number decreases; our estimates for the correlation length exponent are in line with the expectation that there is no change in the universality class. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 5 pages, 5 figures

Superconducting Kondo phase in an orbitally-separated bilayer

Authors: Sebastião dos Anjos Sousa-Júnior, José P. de Lima, Natanael C. Costa, Raimundo R. dos Santos

Abstract: The nature of superconductivity in heavy-fermion materials is a subject under intense debate, and controlling this many-body state is central for its eventual understanding. Here, we examine how proximity effects may change this phenomenon, by investigating the effects of an additional metallic layer on the top of a Kondo-lattice, and allowing for pairing in the former. We analyze a bilayer Kondo… ▽ More The nature of superconductivity in heavy-fermion materials is a subject under intense debate, and controlling this many-body state is central for its eventual understanding. Here, we examine how proximity effects may change this phenomenon, by investigating the effects of an additional metallic layer on the top of a Kondo-lattice, and allowing for pairing in the former. We analyze a bilayer Kondo Lattice Model with an on-site Hubbard interaction, $-U$, on the additional layer, using a mean-field approach. For $U=0$, we notice a drastic change in the density-of-states due to multiple-orbital singlet resonating combinations. It destroys the well-known Kondo insulator at half filling, leading to a metallic ground state, which, in turn, enhances antiferromagnetism through the polarization of the conduction electrons. For $U\neq 0$, a superconducting Kondo state sets in at zero temperature, with the occurrence of unconventional pairing amplitudes involving $f$-electrons. We establish that this remarkable feature is only possible due to the proximity effects of the additional layer. At finite temperatures we find that the critical superconducting temperature, $T_c$, decreases with the interlayer hybridization. We have also established that a zero temperature superconducting amplitude tracks $T_c$, which reminisces the BCS proportionality between the superconducting gap and $T_c$. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: 11 pages, 10 figures

Journal ref: Phys. Rev. Research 2, 033168 (2020)

Dynamical resilience to disorder: the dilute Hubbard model on the Lieb lattice

Authors: L. Oliveira-Lima, N. C. Costa, J. Pimentel de Lima, R. T. Scalettar, R. R. dos Santos

Abstract: In itinerant systems, electron-electron interactions may lead to the formation of local magnetic moments and their effective exchange coupling, which in turn gives rise to long-range magnetic order. Therefore, when moment formation is weakened, such as in the single-band Hubbard model on a square lattice with the on-site repulsion being randomly switched off on a fraction $x$ of sites, magnetic or… ▽ More In itinerant systems, electron-electron interactions may lead to the formation of local magnetic moments and their effective exchange coupling, which in turn gives rise to long-range magnetic order. Therefore, when moment formation is weakened, such as in the single-band Hubbard model on a square lattice with the on-site repulsion being randomly switched off on a fraction $x$ of sites, magnetic order is suppressed beyond some critical $x_c$, which was found to lie below the classical percolation threshold, $x_c^\text{(perc,sq)}$. Here we study dilute magnetism in flat band systems, namely in the Hubbard model on a `Lieb' lattice. Interestingly, we show that magnetic order persists to $x$ almost twice as large as the classical percolation threshold for the lattice, thus emphasizing the central role of electron itinerancy to the magnetic response. The analysis of the orbital-resolved order parameters reveals that the contribution of the four-fold coordinated `d' sites to magnetism is dramatically affected by dilution, while the localized `p' states of the flat band provide the dominant contribution to long-range correlations. We also examine the transport properties, which suggest the existence of an insulator-to-metal transition in the same range of the critical magnetic dilution. △ Less

Submitted 9 April, 2020; v1 submitted 10 December, 2019; originally announced December 2019.

Comments: 9 pages, 10 figures

Journal ref: Phys. Rev. B 101, 165109 (2020)

Determinant Quantum Monte Carlo Study of Exhaustion in the Periodic Anderson Model

Authors: Lufeng Zhang, Tianxing Ma, Natanael C. Costa, Raimundo R. dos Santos, Richard T. Scalettar

Abstract: The Kondo and Periodic Anderson models describe many of the qualitative features of local moments coupled to a conduction band, and thereby the physics of materials such as the heavy fermions. In particular, when the exchange coupling $J$ or hybridization $V$ between the moments and the electrons of the metallic band is large, singlets form, quenching the magnetism. In the opposite, small $J$ or… ▽ More The Kondo and Periodic Anderson models describe many of the qualitative features of local moments coupled to a conduction band, and thereby the physics of materials such as the heavy fermions. In particular, when the exchange coupling $J$ or hybridization $V$ between the moments and the electrons of the metallic band is large, singlets form, quenching the magnetism. In the opposite, small $J$ or $V$, limit, the moments survive, and the conduction electrons mediate an effective interaction which can trigger long range, often antiferromagnetic, order. In the case of the Kondo model, where the moments are described by local spins, Nozières considered the possibility that the available conduction electrons within the Kondo temperature of the Fermi surface would be insufficient in number to accomplish the screening. Much effort in the literature has been devoted to the study of the temperature scales in the resulting `exhaustion' problem, and how the `coherence temperature' where a heavy Fermi liquid forms is related to the Kondo temperature. In this paper, we study a version of the Periodic Anderson model in which some of the conduction electrons are removed in a way which avoids the fermion sign problem and hence allows low temperature Quantum Monte Carlo simulations which can access both singlet formation and magnetic ordering temperature scales. We are then able to focus on a somewhat different aspect of exhaustion physics than previously considered: the effect of dilution on the critical $V$ for the singlet-antiferromagnetic transition. △ Less

Submitted 23 May, 2019; v1 submitted 24 March, 2019; originally announced March 2019.

Comments: 10 figures and 9 pages. Accepted for publication in Physical Review B

Journal ref: Phys. Rev. B 99, 195147 (2019)

A kaleidoscope of phases in the dipolar Hubbard model

Authors: Tiago Mendes-Santos, Rubem Mondaini, Thereza Paiva, Raimundo R. dos Santos

Abstract: We investigate the emergence of a myriad of phases in the strong coupling regime of the dipolar Hubbard model in two dimensions. By using a combination of numerically unbiased methods in finite systems with analytical perturbative arguments, we show the versatility that trapped dipolar atoms possess in displaying a wide variety of many-body phases, which can be tuned simply by changing the collect… ▽ More We investigate the emergence of a myriad of phases in the strong coupling regime of the dipolar Hubbard model in two dimensions. By using a combination of numerically unbiased methods in finite systems with analytical perturbative arguments, we show the versatility that trapped dipolar atoms possess in displaying a wide variety of many-body phases, which can be tuned simply by changing the collective orientation of the atomic dipoles. We further investigate the stability of these phases to thermal fluctuations in the strong coupling regime, highlighting that they can be accessed with current techniques employed in cold atoms experiments on optical lattices. Interestingly, both quantum and thermal phase transitions are signalled by peaks or discontinuities in local moment-local moment correlations, which have been recently measured in some of these experiments, so that they can be used as probes for the onset of different phases. △ Less

Submitted 2 March, 2019; originally announced March 2019.

Coherence Temperature in the Diluted Periodic Anderson Model

Authors: N. C. Costa, T. Mendes-Santos, T. Paiva, N. J. Curro, R. R. dos Santos, R. T. Scalettar

Abstract: The Kondo and Periodic Anderson Model (PAM) are known to provide a microscopic picture of many of the fundamental properties of heavy fermion materials and, more generally, a variety of strong correlation phenomena in $4f$ and $5f$ systems. In this paper, we apply the Determinant Quantum Monte Carlo (DQMC) method to include disorder in the PAM, specifically the removal of a fraction $x$ of the loc… ▽ More The Kondo and Periodic Anderson Model (PAM) are known to provide a microscopic picture of many of the fundamental properties of heavy fermion materials and, more generally, a variety of strong correlation phenomena in $4f$ and $5f$ systems. In this paper, we apply the Determinant Quantum Monte Carlo (DQMC) method to include disorder in the PAM, specifically the removal of a fraction $x$ of the localized orbitals. We determine the evolution of the coherence temperature $T^*$, where the local moments and conduction electrons become entwined in a heavy fermion fluid, with $x$ and with the hybridization $V$ between localized and conduction orbitals. We recover several of the principal observed trends in $T^*$ of doped heavy fermions, and also show that, within this theoretical framework, the calculated Nuclear Magnetic Resonance (NMR) relaxation rate tracks the experimentally measured behavior in pure and doped CeCoIn$_5$. Our results contribute to important issues in the interpretation of local probes of disordered, strongly correlated systems. △ Less

Submitted 9 May, 2019; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: 8 pages, 4 figures

Journal ref: Phys. Rev. B 99, 195116 (2019)

Compressible Ferrimagnetism in the depleted Periodic Anderson Model

Authors: N. C. Costa, M. V. Araújo, J. P. Lima, T. Paiva, R. R. dos Santos, R. T. Scalettar

Abstract: Tight-binding Hamiltonians with single and multiple orbitals exhibit an intriguing array of magnetic phase transitions. In most cases the spin ordered phases are insulating, while the disordered phases may be either metallic or insulating. In this paper we report a Determinant Quantum Monte Carlo study of interacting electrons in a geometry which can be regarded as a two-dimensional Periodic Ander… ▽ More Tight-binding Hamiltonians with single and multiple orbitals exhibit an intriguing array of magnetic phase transitions. In most cases the spin ordered phases are insulating, while the disordered phases may be either metallic or insulating. In this paper we report a Determinant Quantum Monte Carlo study of interacting electrons in a geometry which can be regarded as a two-dimensional Periodic Anderson Model with depleted interacting ($f$) orbitals. For a single depletion, we observe an enhancement of antiferromagnetic correlations and formation of localized states. For half of the $f$-orbitals regularly depleted, the system exhibits a ferrimagnetic ground state. We obtain a quantitative determination of the nature of magnetic order, which we discuss in the context of Tsunetsugu's theorem, and show that, although the dc conductivity indicates insulating behavior at half-filling, the compressibility remains finite. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: 12 pages, 13 figures

Journal ref: Phys. Rev. B 97, 085123 (2018)

A mean-field approach to Kondo-attractive-Hubbard model

Authors: N. C. Costa, J. P. Lima, T. C. L. Paiva, M. ElMassalami, R. R. dos Santos

Abstract: With the purpose of investigating coexistence between magnetic order and superconductivity, we consider a model in which conduction electrons interact with each other, via an attractive Hubbard on-site coupling $U$, and with local moments on every site, via a Kondo-like coupling, $J$. The model is solved on a simple cubic lattice through a Hartree-Fock approximation, within a `semi-classical' fram… ▽ More With the purpose of investigating coexistence between magnetic order and superconductivity, we consider a model in which conduction electrons interact with each other, via an attractive Hubbard on-site coupling $U$, and with local moments on every site, via a Kondo-like coupling, $J$. The model is solved on a simple cubic lattice through a Hartree-Fock approximation, within a `semi-classical' framework which allows spiral magnetic modes to be stabilized. For a fixed electronic density, $n_c$, the small $J$ region of the ground state ($T=0$) phase diagram displays spiral antiferromagnetic (SAFM) states for small $U$. Upon increasing $U$, a state with coexistence between superconductivity (SC) and SAFM sets in; further increase in $U$ turns the spiral mode into a Néel antiferromagnet. The large $J$ region is a (singlet) Kondo phase. At finite temperatures, and in the region of coexistence, thermal fluctuations suppress the different ordered phases in succession: the SAFM phase at lower temperatures and SC at higher temperatures; also, reentrant behaviour is found to be induced by temperature. Our results provide a qualitative description of the competition between local moment magnetism and superconductivity in the borocarbides family. △ Less

Submitted 17 December, 2018; v1 submitted 10 February, 2017; originally announced February 2017.

Comments: 12 pages, 12 figures

Journal ref: J. Phys.: Cond. Matter 30, 045602 (2018)

Ferromagnetism beyond Lieb's theorem

Authors: Natanael C. Costa, Tiago Mendes-Santos, Thereza Paiva, Raimundo R. dos Santos, Richard T. Scalettar

Abstract: The noninteracting electronic structures of tight binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a \emph{uniform} on-site Hubbard interaction $U$ is turned on, Lieb proved rigorously that at half filling ($ρろー=1$) the ground state has a non-zero… ▽ More The noninteracting electronic structures of tight binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a \emph{uniform} on-site Hubbard interaction $U$ is turned on, Lieb proved rigorously that at half filling ($ρろー=1$) the ground state has a non-zero spin. In this paper we consider a `CuO$_2$ lattice (also known as `Lieb lattice', or as a decorated square lattice), in which `$d$-orbitals' occupy the vertices of the squares, while `$p$-orbitals' lie halfway between two $d$-orbitals. We use exact Determinant Quantum Monte Carlo (DQMC) simulations to quantify the nature of magnetic order through the behavior of correlation functions and sublattice magnetizations in the different orbitals as a function of $U$ and temperature. We study both the homogeneous (H) case, $U_d= U_p$, originally considered by Lieb, and the inhomogeneous (IH) case, $U_d\neq U_p$. For the H case at half filling, we found that the global magnetization rises sharply at weak coupling, and then stabilizes towards the strong-coupling (Heisenberg) value, as a result of the interplay between the ferromagnetism of like sites and the antiferromagnetism between unlike sites; we verified that the system is an insulator for all $U$. For the IH system at half filling, we argue that the case $U_p\neq U_d$ falls under Lieb's theorem, provided they are positive definite, so we used DQMC to probe the cases $U_p=0,U_d=U$ and $U_p=U, U_d=0$. We found that the different environments of $d$ and $p$ sites lead to a ferromagnetic insulator when $U_d=0$; by contrast, $U_p=0$ leads to to a metal without any magnetic ordering. In addition, we have also established that at density $ρろー=1/3$, strong antiferromagnetic correlations set in, caused by the presence of one fermion on each $d$ site. △ Less

Submitted 11 October, 2016; originally announced October 2016.

Comments: 10 pages, 14 figures

Journal ref: Phys.Rev.B 94, 155107 (2016)

Impurities near an Antiferromagnetic-Singlet Quantum Critical Point

Authors: T. Mendes, N. Costa, G. Batrouni, N. Curro, R. R. dos Santos, T. Paiva, R. T. Scalettar

Abstract: Heavy fermion systems, and other strongly correlated electron materials, often exhibit a competition between antiferromagnetic (AF) and singlet ground states. Using exact Quantum Monte Carlo (QMC) simulations, we examine the effect of impurities in the vicinity of such AF- singlet quantum critical points, through an appropriately defined impurity susceptibility, $χかい_{imp}$. Our key finding is a con… ▽ More Heavy fermion systems, and other strongly correlated electron materials, often exhibit a competition between antiferromagnetic (AF) and singlet ground states. Using exact Quantum Monte Carlo (QMC) simulations, we examine the effect of impurities in the vicinity of such AF- singlet quantum critical points, through an appropriately defined impurity susceptibility, $χかい_{imp}$. Our key finding is a connection, within a single calculational framework, between AF domains induced on the singlet side of the transition, and the behavior of the nuclear magnetic resonance (NMR) relaxation rate $1/T_1$. We show that local NMR measurements provide a diagnostic for the location of the QCP which agrees remarkably well with the vanishing of the AF order parameter and large values of $χかい_{imp}$. We connect our results with experiments on Cd-doped CeCoIn$_5$. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Journal ref: Phys. Rev. B 95, 054419 (2017)

Spiral magnetic phases on the Kondo Lattice Model: A Hartree-Fock approach

Authors: Natanael de Carvalho Costa, José Pimentel de Lima, Raimundo R. dos Santos

Abstract: We study the Kondo Lattice Model (KLM) on a square lattice through a Hartree-Fock approximation in which the local spins are treated semi-classically, in the sense that their average values are modulated by a magnetic wavevector $\mathbf{Q}$ while they couple with the conduction electrons through fermion operators. In this way, we obtain a ground state phase diagram in which spiral magnetic phases… ▽ More We study the Kondo Lattice Model (KLM) on a square lattice through a Hartree-Fock approximation in which the local spins are treated semi-classically, in the sense that their average values are modulated by a magnetic wavevector $\mathbf{Q}$ while they couple with the conduction electrons through fermion operators. In this way, we obtain a ground state phase diagram in which spiral magnetic phases (in which the wavevector depends on the coupling constants and on the density) interpolate between the low-density ferromagnetic phase and the antiferromagnetic phase at half filling; within small regions of the phase diagram commensurate magnetic phases can coexist with Kondo screening. We have also obtained `Doniach-like' diagrams, showing the effect of temperature on the ground state phases, and established that for some ranges of the model parameters (the exchange coupling and conduction electron density) the magnetic wavevector changes with temperature, either continuously or abruptly (e.g., from spiral to ferromagnetic). △ Less

Submitted 11 October, 2016; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 13 pages, 15 figures

Journal ref: J. Magn. Magn. Mater., 423, 74 (2017)

Size and shape of Mott regions for fermionic atoms in a two-dimensional optical lattice

Authors: Tiago Mendes-Santos, Thereza Paiva, Raimundo R. dos Santos

Abstract: We investigate the harmonic-trap control of size and shape of Mott regions in the Fermi Hubbard model on a square optical lattice. The use of Lanczos diagonalization on clusters with twisted boundary conditions, followed by an average over 50-80 samples, drastically reduce finite-size effects in some ground state properties; calculations in the grand canonical ensemble together with a local-densit… ▽ More We investigate the harmonic-trap control of size and shape of Mott regions in the Fermi Hubbard model on a square optical lattice. The use of Lanczos diagonalization on clusters with twisted boundary conditions, followed by an average over 50-80 samples, drastically reduce finite-size effects in some ground state properties; calculations in the grand canonical ensemble together with a local-density approximation (LDA) allow us to simulate the radial density distribution. We have found that as the trap closes, the atomic cloud goes from a metallic state, to a Mott core, and to a Mott ring; the coverage of Mott atoms reaches a maximum at the core-ring transition. A `phase diagram' in terms of an effective density and the on-site repulsion is proposed, as a guide to maximize the Mott coverage. We also predict that the usual experimentally accessible quantities, the global compressibility and the average double occupancy (rather, its density derivative) display detectable signatures of the core-ring transition. Some spin correlation functions are also calculated, and predict the existence Néel ordering within Mott cores and rings. △ Less

Submitted 6 November, 2014; originally announced November 2014.

Comments: 5 pages, 6 figures

Journal ref: Physical Review A 91, 023632 (2015)

Evolution of the magnetism of Tb(Co_{x}Ni_{1-x})_{2}B_{2}C

Authors: M. ElMassalami, H. Takeya, A. M. Gomes, T. Paiva, R. R. dos Santos

Abstract: The magnetic properties of polycrystalline Tb(Co_{x}Ni_{1-x})_{2}B_{2}C (x=0.2,0.4,0.6,0.8) samples were probed by magnetization, specific heat, ac susceptibility, and resistivity techniques. For x{\neq}0.4, the obtained curves are consistent with the features expected for the corresponding magnetic modes, namely k_{1}=(0.55,0,0) at x=0; k_{2}=([nicefrac]\nicefrac{1}{2}</LaTeX>,0,[nicefrac]<LaTeX>… ▽ More The magnetic properties of polycrystalline Tb(Co_{x}Ni_{1-x})_{2}B_{2}C (x=0.2,0.4,0.6,0.8) samples were probed by magnetization, specific heat, ac susceptibility, and resistivity techniques. For x{\neq}0.4, the obtained curves are consistent with the features expected for the corresponding magnetic modes, namely k_{1}=(0.55,0,0) at x=0; k_{2}=([nicefrac]\nicefrac{1}{2}</LaTeX>,0,[nicefrac]<LaTeX>\nicefrac{1}{2}) at x= 0.2; k_{3}=(0,0,[nicefrac]\nicefrac{1}{3}) at x= 0.6, and k_{4}=(0,0,0) at x= 0.8 and 1. For x=0.4, even though the neutron diffraction indicates a k_{2} mode, but with a reduced magnetic moment, the magnetization, the ac susceptibility, and resistivity indicate two magnetic events; furthermore, deviation from Curie-Weiss behavior is observed below 150 K for this sample. These features, together with the evolution of both magnetic moment and critical temperature, are attributed to an interplay between competing magnetic couplings; for the particular x=0.4 case, additional factors such as crystalline electric field effects may be in operation. △ Less

Submitted 2 July, 2012; originally announced July 2012.

Tuning in magnetic modes in Tb(Co_{x}Ni_{1-x})_{2}B_{2}C: from longitudinal spin-density waves to simple ferromagnetism

Authors: M. ElMassalami, H. Takeya, B. Ouladdiaf, R. Maia Filho, A. M. Gomes, T. Paiva, R. R. dos Santos

Abstract: Neutron diffraction and thermodynamics techniques were used to probe the evolution of the magnetic properties of Tb(Co_{x}Ni_{1-x})_{2}B_{2}C. A succession of magnetic modes was observed as x is varied: the longitudinal modulated k=(0.55,0,0) state at x=0 is transformed into a collinear k=([nicefrac]<LaTeX>\nicefrac{1}{2}</LaTeX>,0,[nicefrac]<LaTeX>\nicefrac{1}{2}</LaTeX>) antiferromagnetic state… ▽ More Neutron diffraction and thermodynamics techniques were used to probe the evolution of the magnetic properties of Tb(Co_{x}Ni_{1-x})_{2}B_{2}C. A succession of magnetic modes was observed as x is varied: the longitudinal modulated k=(0.55,0,0) state at x=0 is transformed into a collinear k=([nicefrac]<LaTeX>\nicefrac{1}{2}</LaTeX>,0,[nicefrac]<LaTeX>\nicefrac{1}{2}</LaTeX>) antiferromagnetic state at x= 0.2, 0.4; then into a transverse c-axis modulated k=(0,0,[nicefrac]<LaTeX>\nicefrac{1}{3}</LaTeX>) mode at x= 0.6, and finally into a simple ferromagnetic structure at x= 0.8 and 1. Concomitantly, the low-temperature orthorhombic distortion of the tetragonal unit cell at x=0 is reduced smoothly such that for x >= 0.4 only a tetragonal unit cell is manifested. Though predicted theoretically earlier, this is the first observation of the k=(0,0,[nicefrac]<LaTeX>\nicefrac{1}{3}</LaTeX>) mode in borocarbides; our findings of a succession of magnetic modes upon increasing x also find support from a recently proposed theoretical model. The implication of these findings and their interpretation on the magnetic structure of the RM_{2}B_{2}C series are also discussed. △ Less

Submitted 7 May, 2012; originally announced May 2012.

Finite-size effects in transport data from Quantum Monte Carlo simulations

Authors: Rubem Mondaini, K. Bouadim, Thereza Paiva, Raimundo R. dos Santos

Abstract: We have examined the behavior of the compressibility, the dc-conductivity, the single-particle gap, and the Drude weight as probes of the density-driven metal-insulator transition in the Hubbard model on a square lattice. These quantities have been obtained through determinantal quantum Monte Carlo simulations at finite temperatures on lattices up to 16 X 16 sites. While the compressibility, the d… ▽ More We have examined the behavior of the compressibility, the dc-conductivity, the single-particle gap, and the Drude weight as probes of the density-driven metal-insulator transition in the Hubbard model on a square lattice. These quantities have been obtained through determinantal quantum Monte Carlo simulations at finite temperatures on lattices up to 16 X 16 sites. While the compressibility, the dc-conductivity, and the gap are known to suffer from `closed-shell' effects due to the presence of artificial gaps in the spectrum (caused by the finiteness of the lattices), we have established that the former tracks the average sign of the fermionic determinant (<sign>), and that a shortcut often used to calculate the conductivity may neglect important corrections. Our systematic analyses also show that, by contrast, the Drude weight is not too sensitive to finite-size effects, being much more reliable as a probe to the insulating state. We have also investigated the influence of the discrete imaginary-time interval (Δでるたτたう) on <sign>, on the average density (ρろー), and on the double occupancy (d): we have found that <sign> and ρろーare more strongly dependent on Δでるたτたうaway from closed-shell configurations, but d follows the Δでるたτたう^2 dependence in both closed- and open-shell cases. △ Less

Submitted 16 March, 2012; v1 submitted 1 July, 2011; originally announced July 2011.

Comments: 11 pages and 16 figures. To appear in Phys Rev B

Journal ref: Phys. Rev. B 85, 125127 (2012)

Magnetic exchange mechanism for electronic gap opening in graphene

Authors: T. G. Rappoport, M. Godoy, B. Uchoa, R. R dos Santos, A. H. Castro Neto

Abstract: We show within a local self-consistent mean-field treatment that a random distribution of magnetic adatoms can open a robust gap in the electronic spectrum of graphene. The electronic gap results from the interplay between the nature of the graphene sublattice structure and the exchange interaction between adatoms.The size of the gap depends on the strength of the exchange interaction between carr… ▽ More We show within a local self-consistent mean-field treatment that a random distribution of magnetic adatoms can open a robust gap in the electronic spectrum of graphene. The electronic gap results from the interplay between the nature of the graphene sublattice structure and the exchange interaction between adatoms.The size of the gap depends on the strength of the exchange interaction between carriers and localized spins and can be controlled by both temperature and external magnetic field. Furthermore, we show that an external magnetic field creates an imbalance of spin-up and spin-down carriers at the Fermi level, making doped graphene suitable for spin injection and other spintronic applications. △ Less

Submitted 8 October, 2011; v1 submitted 18 August, 2010; originally announced August 2010.

Comments: 5 pages, 5 figures

Journal ref: EPL 96, 27010 (2011)

Suppression of Anderson localization of light and Brewster anomalies in disordered superlattices containing a dispersive metamaterial

Authors: D. Mogilevtsev, F. A. Pinheiro, R. R. dos Santos, S. B. Cavalcanti, L. E. Oliveira

Abstract: Light propagation through 1D disordered structures composed of alternating layers, with random thicknesses, of air and a dispersive metamaterial is theoretically investigated. Both normal and oblique incidences are considered. By means of numerical simulations and an analytical theory, we have established that Anderson localization of light may be suppressed: (i) in the long wavelength limit, for… ▽ More Light propagation through 1D disordered structures composed of alternating layers, with random thicknesses, of air and a dispersive metamaterial is theoretically investigated. Both normal and oblique incidences are considered. By means of numerical simulations and an analytical theory, we have established that Anderson localization of light may be suppressed: (i) in the long wavelength limit, for a finite angle of incidence which depends on the parameters of the dispersive metamaterial; (ii) for isolated frequencies and for specific angles of incidence, corresponding to Brewster anomalies in both positive- and negative-refraction regimes of the dispersive metamaterial. These results suggest that Anderson localization of light could be explored to control and tune light propagation in disordered metamaterials. △ Less

Submitted 31 May, 2010; originally announced June 2010.

Comments: 4 two-column pages, 3 figures

Finite-size scaling behavior in trapped systems

Authors: S. L. A. de Queiroz, R. R. dos Santos, R. B. Stinchcombe

Abstract: Numerical transfer-matrix methods are applied to two-dimensional Ising spin systems, in presence of a confining magnetic field which varies with distance $|{\vec x}|$ to a "trap center", proportionally to $(|{\vec x}|/\ell)^p$, $p>0$. On a strip geometry, the competition between the "trap size" $\ell$ and the strip width, $L$, is analysed in the context of a generalized finite-size scaling {\em an… ▽ More Numerical transfer-matrix methods are applied to two-dimensional Ising spin systems, in presence of a confining magnetic field which varies with distance $|{\vec x}|$ to a "trap center", proportionally to $(|{\vec x}|/\ell)^p$, $p>0$. On a strip geometry, the competition between the "trap size" $\ell$ and the strip width, $L$, is analysed in the context of a generalized finite-size scaling {\em ansatz}. In the low-field regime $\ell \gg L$, we use conformal-invariance concepts in conjunction with a linear-response approach to derive the appropriate ($p$-dependent) limit of the theory, which agrees very well with numerical results for magnetization profiles. For high fields $\ell \lesssim L$, correlation-length scaling data broadly confirms an existing picture of $p$-dependent characteristic exponents. Standard spin-1/2 and spin-1 Ising systems are considered, as well as the Blume-Capel model. △ Less

Submitted 27 May, 2010; v1 submitted 4 March, 2010; originally announced March 2010.

Comments: RevTeX, 8 pages, 7 .eps figures (published version

Journal ref: Physical Review E 81, 051122 (2010)

Kondo--attractive-Hubbard model for the ordering of local magnetic moments in superconductors

Authors: P. R. Bertussi, A. L. Malvezzi, T. Paiva, R. R. dos Santos

Abstract: We consider local magnetic moments coupled to conduction electrons with on-site attraction, in order to discuss the interplay between pairing and magnetic order. We probe the ground state properties of this model on a one-dimensional lattice through pair binding energies and several correlation functions, calculated by means of density-matrix renormalization group. A phase diagram is obtained (f… ▽ More We consider local magnetic moments coupled to conduction electrons with on-site attraction, in order to discuss the interplay between pairing and magnetic order. We probe the ground state properties of this model on a one-dimensional lattice through pair binding energies and several correlation functions, calculated by means of density-matrix renormalization group. A phase diagram is obtained (for fixed electron density 1/3), from which we infer that coexistence between magnetic order and superconductivity is robust, at the expense of a continuous distortion of the magnetic arrangement of the local moments, as evidenced by a strong dependence of the characteristic wave vector k* with the coupling constants. This allows us to understand some trends of the coexistence, such as the influence of the rare earth on k*, as observed experimentally in the borocarbides. △ Less

Submitted 15 June, 2009; originally announced June 2009.

Comments: 4 two-column pages, 7 figures, RevTeX, to appear in Phys Rev B (RC)

Disordered two-dimensional superconductors: roles of temperature and interaction strength

Authors: Felipe Mondaini, Thereza Paiva, Raimundo R. dos Santos, R. T. Scalettar

Abstract: We have considered the half-filled disordered attractive Hubbard model on a square lattice, in which the on-site attraction is switched off on a fraction $f$ of sites, while keeping a finite $U$ on the remaining ones. Through Quantum Monte Carlo (QMC) simulations for several values of $f$ and $U$, and for system sizes ranging from $8\times 8$ to $16\times 16$, we have calculated the configuratio… ▽ More We have considered the half-filled disordered attractive Hubbard model on a square lattice, in which the on-site attraction is switched off on a fraction $f$ of sites, while keeping a finite $U$ on the remaining ones. Through Quantum Monte Carlo (QMC) simulations for several values of $f$ and $U$, and for system sizes ranging from $8\times 8$ to $16\times 16$, we have calculated the configurational averages of the equal-time pair structure factor $P_s$, and, for a more restricted set of variables, the helicity modulus, $ρろー_s$, as functions of temperature. Two finite-size scaling {\it ansatze} for $P_s$ have been used, one for zero-temperature and the other for finite temperatures. We have found that the system sustains superconductivity in the ground state up to a critical impurity concentration, $f_c$, which increases with $U$, at least up to U=4 (in units of the hopping energy). Also, the normalized zero-temperature gap as a function of $f$ shows a maximum near $f\sim 0.07$, for $2\lesssim U\lesssim 6$. Analyses of the helicity modulus and of the pair structure factor led to the determination of the critical temperature as a function of $f$, for $U=3,$ 4 and 6: they also show maxima near $f\sim 0.07$, with the highest $T_c$ increasing with $U$ in this range. We argue that, overall, the observed behavior results from both the breakdown of CDW-superconductivity degeneracy and the fact that free sites tend to "push" electrons towards attractive sites, the latter effect being more drastic at weak couplings. △ Less

Submitted 24 July, 2008; originally announced July 2008.

Comments: 9 two-column pages, 14 figures, RevTeX

Modulation of charge-density waves by superlattice structures

Authors: Andre L Malvezzi, Thereza Paiva, Raimundo R dos Santos

Abstract: We discuss the interplay between electronic correlations and an underlying superlattice structure in determining the period of charge density waves (CDW's), by considering a one-dimensional Hubbard model with a repeated (non-random) pattern of repulsive (U>0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems (up to 120 sites) is used to calculate the ch… ▽ More We discuss the interplay between electronic correlations and an underlying superlattice structure in determining the period of charge density waves (CDW's), by considering a one-dimensional Hubbard model with a repeated (non-random) pattern of repulsive (U>0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems (up to 120 sites) is used to calculate the charge-density correlation function and structure factor in the ground state. The modulation period can still be predicted through effective Fermi wavevectors, k_F*, and densities, and we have found that it is much more sensitive to electron (or hole) doping, both because of the narrow range of densities needed to go from q*=0 to πぱい, but also due to sharp 2k_F*-4k_F* transitions; these features render CDW's more versatile for actual applications in heterostructures than in homogeneous systems. △ Less

Submitted 21 April, 2006; originally announced April 2006.

Comments: 4 pages, 5 figures, to appear in Phys Rev B

Journal ref: Phys. Rev. B 73, 193407 (2006)

Destruction of Superconductivity by Impurities in the Attractive Hubbard Model

Authors: Daniel Hurt, Evan Odabashian, Warren Pickett, Richard Scalettar, Felipe Mondaini, Thereza Paiva, Raimundo R. dos Santos

Abstract: We study the effect of U=0 impurities on the superconducting and thermodynamic properties of the attractive Hubbard model on a square lattice. Removal of the interaction on a critical fraction of $f_{\rm crit} \approx 0.30$ of the sites results in the destruction of off-diagonal long range order in the ground state. This critical fraction is roughly independent of filling in the range… ▽ More We study the effect of U=0 impurities on the superconducting and thermodynamic properties of the attractive Hubbard model on a square lattice. Removal of the interaction on a critical fraction of $f_{\rm crit} \approx 0.30$ of the sites results in the destruction of off-diagonal long range order in the ground state. This critical fraction is roughly independent of filling in the range $0.75 < ρろー< 1.00$, although our data suggest that $f_{\rm crit}$ might be somewhat larger below half-filling than at $ρろー=1$. We also find that the two peak structure in the specific heat is present at $f$ both below and above the value which destroys long range pairing order. It is expected that the high $T$ peak associated with local pair formation should be robust, but apparently local pairing fluctuations are sufficient to generate a low temperature peak. △ Less

Submitted 8 July, 2005; originally announced July 2005.

Journal ref: Phys. Rev B 72, 144513 (2005)

Disorder and the effective Mn-Mn exchange interaction in Ga$_{1-x}$Mn$_x$As diluted magnetic semiconductors

Authors: Antonio J. R. da Silva, A. Fazzio, Raimundo R. dos Santos, Luiz E. Oliveira

Abstract: We perform a theoretical study, using {\it ab initio} total energy density-functional calculations, of the effects of disorder on the $Mn-Mn$ exchange interactions for $Ga_{1-x}Mn_xAs$ diluted semiconductors. For a 128 atoms supercell, we consider a variety of configurations with 2, 3 and 4 Mn atoms, which correspond to concentrations of 3.1%, 4.7%, and 6.3%, respectively. In this way, the disor… ▽ More We perform a theoretical study, using {\it ab initio} total energy density-functional calculations, of the effects of disorder on the $Mn-Mn$ exchange interactions for $Ga_{1-x}Mn_xAs$ diluted semiconductors. For a 128 atoms supercell, we consider a variety of configurations with 2, 3 and 4 Mn atoms, which correspond to concentrations of 3.1%, 4.7%, and 6.3%, respectively. In this way, the disorder is intrinsically considered in the calculations. Using a Heisenberg Hamiltonian to map the magnetic excitations, and {\it ab initio} total energy calculations, we obtain the effective $\JMn$, from first ($n=1$) all the way up to sixth ($n=6$) neighbors. Calculated results show a clear dependence in the magnitudes of the $\JMn$ with the Mn concentration $x$. Also, configurational disorder and/or clustering effects lead to large dispersions in the Mn-Mn exchange interactions, in the case of fixed Mn concentration. Moreover, theoretical results for the ground-state total energies for several configurations indicate the importance of a proper consideration of disorder in treating temperature and annealing effects. △ Less

Submitted 18 July, 2005; v1 submitted 20 May, 2005; originally announced May 2005.

Critical temperature for the two-dimensional attractive Hubbard Model

Authors: Thereza Paiva, Raimundo R. dos Santos, R. T. Scalettar, P. J. H. Denteneer

Abstract: The critical temperature for the attractive Hubbard model on a square lattice is determined from the analysis of two independent quantities, the helicity modulus, $ρろー_s$, and the pairing correlation function, $P_s$. These quantities have been calculated through Quantum Monte Carlo simulations for lattices up to $18\times 18$, and for several densities, in the intermediate-coupling regime. Imposin… ▽ More The critical temperature for the attractive Hubbard model on a square lattice is determined from the analysis of two independent quantities, the helicity modulus, $ρろー_s$, and the pairing correlation function, $P_s$. These quantities have been calculated through Quantum Monte Carlo simulations for lattices up to $18\times 18$, and for several densities, in the intermediate-coupling regime. Imposing the universal-jump condition for an accurately calculated $ρろー_s$, together with thorough finite-size scaling analyses (in the spirit of the phenomenological renormalization group) of $P_s$, suggests that $T_c$ is considerably higher than hitherto assumed. △ Less

Submitted 16 March, 2004; originally announced March 2004.

Comments: 5 pages, 6 figures. Accepted for publication in Phys. Rev. B

Journal ref: Physical Review B 69, 184501 (2004)

A layering model for superconductivity in the borocarbides

Authors: Thereza Paiva, Mohammed El-Massalami, Raimundo R. dos Santos

Abstract: We propose a superlattice model to describe superconductivity in layered materials, such as the borocarbide families with the chemical formulæ $RT_2$B$_2$C and $RT$BC, with $R$ being (essentially) a rare earth, and $T$ a transition metal. We assume a single band in which electrons feel a local attractive interaction (negative Hubbard-$U$) on sites representing the $T$B layers, while U=0 on sites… ▽ More We propose a superlattice model to describe superconductivity in layered materials, such as the borocarbide families with the chemical formulæ $RT_2$B$_2$C and $RT$BC, with $R$ being (essentially) a rare earth, and $T$ a transition metal. We assume a single band in which electrons feel a local attractive interaction (negative Hubbard-$U$) on sites representing the $T$B layers, while U=0 on sites representing the $R$C layers; the multi-band structure is taken into account minimally through a band offset $εいぷしろん$. The one-dimensional model is studied numerically through the calculation of the charge gap, the Drude weight, and of the pairing correlation function. A comparison with the available information on the nature of the electronic ground state (metallic or superconducting) indicates that the model provides a systematic parametrization of the whole borocarbide family. △ Less

Submitted 7 November, 2003; originally announced November 2003.

Comments: 4 figures

Journal ref: J. Phys.: Condens. Matter 15, 7917 (2003)

On the nature of ferromagnetism in Ga_{1-x}Mn_xAs semiconductors

Authors: A. J. R. da Silva, A. Fazzio, R. R. dos Santos, L. E. Oliveira

Abstract: We have performed ab initio calculations within the density-functional theory for Ga_{1-x}Mn_xAs diluted semiconductors. Total energy results unambiguously show that a quasi-localized down-spin hole, with strong p-like character, surrounds the fully polarized Mn up-spin d^5-electrons. Calculations indicate that the holes form an essentially dispersionless impurity band, thus rendering effective-… ▽ More We have performed ab initio calculations within the density-functional theory for Ga_{1-x}Mn_xAs diluted semiconductors. Total energy results unambiguously show that a quasi-localized down-spin hole, with strong p-like character, surrounds the fully polarized Mn up-spin d^5-electrons. Calculations indicate that the holes form an essentially dispersionless impurity band, thus rendering effective-mass descriptions of hole states open to challenge. We obtain estimates both for the s = 1/2 hole and S = 5/2 Mn exchange coupling, and for the distance dependence of the effective Mn-Mn exchange interaction. Results demonstrate that the effective Mn-Mn coupling is always ferromagnetic, thus non-RKKY, and intermediated by the antiferromagnetic coupling of each Mn spin to the holes. △ Less

Submitted 29 July, 2003; originally announced July 2003.

Comments: RevTeX 4, 4 two-column pages, 3 .gif figures; submitted to Phys Rev Lett

Introduction to Quantum Monte Carlo simulations for fermionic systems

Authors: Raimundo R. dos Santos

Abstract: We tutorially review the determinantal Quantum Monte Carlo method for fermionic systems, using the Hubbard model as a case study. Starting with the basic ingredients of Monte Carlo simulations for classical systems, we introduce aspects such as importance sampling, sources of errors, and finite-size scaling analyses. We then set up the preliminary steps to prepare for the simulations, showing th… ▽ More We tutorially review the determinantal Quantum Monte Carlo method for fermionic systems, using the Hubbard model as a case study. Starting with the basic ingredients of Monte Carlo simulations for classical systems, we introduce aspects such as importance sampling, sources of errors, and finite-size scaling analyses. We then set up the preliminary steps to prepare for the simulations, showing that they are actually carried out by sampling discrete Hubbard-Stratonovich auxiliary fields. In this process the Green's function emerges as a fundamental tool, since it is used in the updating process, and, at the same time, it is directly related to the quantities probing magnetic, charge, metallic, and superconducting behaviours. We also discuss the as yet unresolved "minus-sign problem", and two ways to stabilize the algorithm at low temperatures. △ Less

Submitted 26 March, 2003; originally announced March 2003.

Comments: 18 two-column pages, 3 figures, RevTeX 4. Short course delivered at the Brazilian School on Statistical Mechanics, held at the Universidade de Sao Paulo at Sao Carlos, in February 2002; other courses can be found at http://pcsbf1.if.usp.br/bjp/Vol33/Num1

Journal ref: Braz. J. Phys. 33, 36 (2003)

Magnetic-field effects in defect-controlled ferromagnetic Ga_{1-x}Mn_xAs semiconductors

Authors: Raimundo R. dos Santos, J. d'Albuquerque e Castro, L. E. Oliveira

Abstract: We have studied the magnetic-field and concentration dependences of the magnetizations of the hole and Mn subsystems in diluted ferromagnetic semiconductor Ga_{1-x}Mn_xAs. A mean-field approximation to the hole-mediated interaction is used, in which the hole concentration p(x) is parametrized in terms of a fitting (of the hole effective mass and hole/local moment coupling) to experimental data o… ▽ More We have studied the magnetic-field and concentration dependences of the magnetizations of the hole and Mn subsystems in diluted ferromagnetic semiconductor Ga_{1-x}Mn_xAs. A mean-field approximation to the hole-mediated interaction is used, in which the hole concentration p(x) is parametrized in terms of a fitting (of the hole effective mass and hole/local moment coupling) to experimental data on the Tc critical temperature. The dependence of the magnetizations with x, for a given temperature, presents a sharply peaked structure, with maxima increasing with applied magnetic field, which indicates that application to diluted-magnetic-semiconductor devices would require quality-control of the Mn-doping composition. We also compare various experimental data for Tc(x) and p(x) on different Ga_{1-x}Mn_xAs samples and stress the need of further detailed experimental work to assure that the experimental measurements are reproducible. △ Less

Submitted 13 January, 2003; originally announced January 2003.

Comments: RevTeX 4, 3 two-column pages, 4 colour figures; to appear in J Appl Phys

Superconductivity Induced by Negative Centers

Authors: G. Litak, T. Paiva, R. T. Scalettar, C. Huscroft, R. R. dos Santos

Abstract: We study the effect of random dilution of attraction centers on electron pairing in a generalization of the negative $U$ Hubbard model where the interaction $U_i$, defined on the two dimensional lattice, can be $-|U|$ and 0 with probability $c$ and $1-c$ depending on the lattice site $i$. Using the determinant quantum Monte Carlo approach, we find the critical concentration $c_0$ of negative cen… ▽ More We study the effect of random dilution of attraction centers on electron pairing in a generalization of the negative $U$ Hubbard model where the interaction $U_i$, defined on the two dimensional lattice, can be $-|U|$ and 0 with probability $c$ and $1-c$ depending on the lattice site $i$. Using the determinant quantum Monte Carlo approach, we find the critical concentration $c_0$ of negative centers which leads to superconductivity, and show how the evolution of the local density on the two types of sites lends insight into the formation of global phase coherence. △ Less

Submitted 8 June, 2002; originally announced June 2002.

Comments: 10 pages, 3 figures, (submitted toPhysica C)

Multiperiodic magnetic structures in Hubbard superlattices

Authors: Andre L Malvezzi, Thereza Paiva, Raimundo R dos Santos

Abstract: We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard-U couplings arranged in a repeated pattern of repulsive (i.e., U>0) and free (U=0) sites. Density Matrix Renormalization Group (DMRG) diagonalization of finite systems is used to calculate the local moment and the magnetic structure factor in the ground state. We have found four regimes for magnetic b… ▽ More We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard-U couplings arranged in a repeated pattern of repulsive (i.e., U>0) and free (U=0) sites. Density Matrix Renormalization Group (DMRG) diagonalization of finite systems is used to calculate the local moment and the magnetic structure factor in the ground state. We have found four regimes for magnetic behavior: uniform local moments forming a spin-density wave (SDW), `floppy' local moments with short-ranged correlations, local moments on repulsive sites forming long-period SDW's superimposed with short-ranged correlations, and local moments on repulsive sites solely with long-period SDW's; the boundaries between these regimes depend on the range of electronic densities, rho, and on the SL aspect ratio. Above a critical electronic density, rho_{uparrow downarrow}, the SDW period oscillates both with rho and with the spacer thickness. The former oscillation allows one to reproduce all SDW wave vectors within a small range of electronic densities, unlike the homogeneous system. The latter oscillation is related to the exchange oscillation observed in magnetic multilayers. A crossover between regimes of `thin' to `thick' layers has also been observed. △ Less

Submitted 13 May, 2002; originally announced May 2002.

Comments: 9 two-column pages, 10 figures

Journal ref: Phys. Rev. B 66, 064430 (2002)

Luttinger liquid superlattices: realization of gapless insulating phases

Authors: J. Silva-Valencia, E. Miranda, Raimundo R. dos Santos

Abstract: We investigate Luttinger Liquid superlattices, a periodic structure composed of two kinds of one-dimensional systems of interacting electrons. We calculate several properties of the low-energy sector: the effective charge and spin velocities, the compressibility, various correlation functions, the Landauer conductance and the Drude weight. The low-energy properties are subsumed into effective pa… ▽ More We investigate Luttinger Liquid superlattices, a periodic structure composed of two kinds of one-dimensional systems of interacting electrons. We calculate several properties of the low-energy sector: the effective charge and spin velocities, the compressibility, various correlation functions, the Landauer conductance and the Drude weight. The low-energy properties are subsumed into effective parameters, much like homogeneous one-dimensional systems. A generic result is the weighted average nature of these parameters, in proportion to the spatial extent of the underlying subunits, pointing to the possibility of ``engineered'' structures. As a specific realization, we consider a one-dimensional Hubbard superlattice, which consists of a periodic arrangement of two long Hubbard chains with different coupling constants and different hopping amplitudes. This system exhibits a rich phase diagram with several phases, both metallic and insulating. We have found that gapless insulating phases are present over a wide range of parameters. △ Less

Submitted 24 April, 2002; v1 submitted 2 January, 2002; originally announced January 2002.

Comments: 16 pages, 15 figures, RevTeX4

Journal ref: Phys. Rev. B 65, 115115 (2002)

Charge-density waves in one-dimensional Hubbard superlattices

Authors: Thereza Paiva, Raimundo R. dos Santos

Abstract: We study the formation of charge density waves (CDW's) in one-dimensional Hubbard superlattices, modeled by a repeated pattern of repulsive (U>0) and free (U=0) sites. By means of Lanczos diagonalizations for the ground state, we calculate the charge structure factor. Our results show that while the superlattice structure affects the modulation of the charge density waves, the periodicity can st… ▽ More We study the formation of charge density waves (CDW's) in one-dimensional Hubbard superlattices, modeled by a repeated pattern of repulsive (U>0) and free (U=0) sites. By means of Lanczos diagonalizations for the ground state, we calculate the charge structure factor. Our results show that while the superlattice structure affects the modulation of the charge density waves, the periodicity can still be predicted through an effective density. We also show that, for a fixed repulsive layer thickness, the periodicity of the CDW is an oscillatory function of the free layer thickness. △ Less

Submitted 25 September, 2001; originally announced September 2001.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. B 65, 153101 (2002).

Hole concentration in a diluted ferromagnetic semiconductor

Authors: Raimundo R. dos Santos, Luiz E. Oliveira, J. d'Albuquerque e Castro

Abstract: We consider a mean-field approach to the hole-mediated ferromagnetism in III-V Mn-based semiconductor compounds to discuss the dependence of the hole density on that of Mn sites in Ga_{1-x}Mn_xAs. The hole concentration, p, as a function of the fraction of Mn sites, x, is parametrized in terms of the product m*J_{pd}^2 (where m* is the hole effective mass and J_{pd} is the Kondo-like hole/local-… ▽ More We consider a mean-field approach to the hole-mediated ferromagnetism in III-V Mn-based semiconductor compounds to discuss the dependence of the hole density on that of Mn sites in Ga_{1-x}Mn_xAs. The hole concentration, p, as a function of the fraction of Mn sites, x, is parametrized in terms of the product m*J_{pd}^2 (where m* is the hole effective mass and J_{pd} is the Kondo-like hole/local-moment coupling), and the critical temperature Tc. By using experimental data for these quantities, we have established the dependence of the hole concentration with x, which can be associated with the occurrence of a reentrant metal-insulator transition taking place in the hole gas. We also calculate the dependence of the Mn magnetization with x, for different temperatures (T), and found that as T increases, the width of the composition-dependent magnetization decreases drammatically, and that the magnetization maxima also decreases, indicating the need for quality-control of Mn-doping composition in diluted magnetic semiconductor devices. △ Less

Submitted 11 January, 2002; v1 submitted 15 August, 2001; originally announced August 2001.

Comments: 4 pages, 3 figures, RevTeX 3; Fig. 1 changed, new references added

Journal ref: J Phys: Condensed Matter 14, 3751 (2002)

Luttinger liquid superlattices

Authors: J. Silva-Valencia, E. Miranda, Raimundo R. dos Santos

Abstract: We calculate the correlation functions and the DC conductivity of Luttinger liquid superlattices, modeled by a repeated pattern of interacting and free Luttinger liquids. In a specific realization, where the interacting subsystem is a Hubbard chain, the system exhibits a rich phase diagram with four different phases: two metals and two compressible insulators. In general, we find that the effect… ▽ More We calculate the correlation functions and the DC conductivity of Luttinger liquid superlattices, modeled by a repeated pattern of interacting and free Luttinger liquids. In a specific realization, where the interacting subsystem is a Hubbard chain, the system exhibits a rich phase diagram with four different phases: two metals and two compressible insulators. In general, we find that the effective low energy description amalgamates features of both types of liquids in proportion to their spatial extent, suggesting the interesting possibility of `engineered' Luttinger liquids. △ Less

Submitted 5 July, 2001; originally announced July 2001.

Comments: RevTeX, 5 pages, 3 figures

Journal ref: J. Phys.:Condens. Matter 13, L619 (2001)

Charge-density waves in the Hubbard chain: evidence for 4k_F instability

Authors: Thereza Paiva, Raimundo R dos Santos

Abstract: Charge density waves in the Hubbard chain are studied by means of finite-temperature Quantum Monte Carlo simulations and Lanczos diagonalizations for the ground state. We present results both for the charge susceptibilities and for the charge structure factor at densities ρろー=1/6 and 1/3; for ρろー=1/2 (quarter filled) we only present results for the charge structure factor. The data are consistent wi… ▽ More Charge density waves in the Hubbard chain are studied by means of finite-temperature Quantum Monte Carlo simulations and Lanczos diagonalizations for the ground state. We present results both for the charge susceptibilities and for the charge structure factor at densities ρろー=1/6 and 1/3; for ρろー=1/2 (quarter filled) we only present results for the charge structure factor. The data are consistent with a 4k_F instability dominating over the 2k_F one, at least for sufficiently large values of the Coulomb repulsion, U. This can only be reconciled with the Luttinger liquid analyses if the amplitude of the 2k_F contribution vanishes above some U^*(ρろー). △ Less

Submitted 2 February, 2000; originally announced February 2000.

Comments: RevTeX, 4 two-column pages with 7 colour figures embedded in text

Journal ref: Phys Rev B 61, 13480 (2000)

Universality, frustration and conformal invariance in two-dimensional random Ising magnets

Authors: F. D. A. Aarao Reis, S. L. A. de Queiroz, Raimundo R. dos Santos

Abstract: We consider long, finite-width strips of Ising spins with randomly distributed couplings. Frustration is introduced by allowing both ferro- and antiferromagnetic interactions. Free energy and spin-spin correlation functions are calculated by transfer-matrix methods. Numerical derivatives and finite-size scaling concepts allow estimates of the usual critical exponents $γがんま/νにゅー$, $αあるふぁ/νにゅー$ and $νにゅー$ to be o… ▽ More We consider long, finite-width strips of Ising spins with randomly distributed couplings. Frustration is introduced by allowing both ferro- and antiferromagnetic interactions. Free energy and spin-spin correlation functions are calculated by transfer-matrix methods. Numerical derivatives and finite-size scaling concepts allow estimates of the usual critical exponents $γがんま/νにゅー$, $αあるふぁ/νにゅー$ and $νにゅー$ to be obtained, whenever a second-order transition is present. Low-temperature ordering persists for suitably small concentrations of frustrated bonds, with a transition governed by pure--Ising exponents. Contrary to the unfrustrated case, subdominant terms do not fit a simple, logarithmic-enhancement form. Our analysis also suggests a vertical critical line at and below the Nishimori point. Approaching this point along either the temperature axis or the Nishimori line, one finds non-diverging specific heats. A percolation-like ratio $γがんま/νにゅー$ is found upon analysis of the uniform susceptibility at the Nishimori point. Our data are also consistent with frustration inducing a breakdown of the relationship between correlation-length amplitude and critical exponents, predicted by conformal invariance for pure systems. △ Less

Submitted 8 June, 1999; originally announced June 1999.

Comments: RevTeX code for 10 pages, 9 eps figures, to appear in Physical Review B (September 1999)

Journal ref: Physical Review B 60, 6740 (1999)

Field-induced Ordering in Critical Antiferromagnets

Authors: S. L. A. de Queiroz, Thereza Paiva, Jorge S. de Sá Martins, Raimundo R. dos Santos

Abstract: Transfer-matrix scaling methods have been used to study critical properties of field-induced phase transitions of two distinct two-dimensional antiferromagnets with discrete-symmetry order parameters: triangular-lattice Ising systems (TIAF) and the square-lattice three-state Potts model (SPAF-3). Our main findings are summarised as follows. For TIAF, we have shown that the critical line leaves t… ▽ More Transfer-matrix scaling methods have been used to study critical properties of field-induced phase transitions of two distinct two-dimensional antiferromagnets with discrete-symmetry order parameters: triangular-lattice Ising systems (TIAF) and the square-lattice three-state Potts model (SPAF-3). Our main findings are summarised as follows. For TIAF, we have shown that the critical line leaves the zero-temperature, zero -field fixed point at a finite angle. Our best estimate of the slope at the origin is $(dT_c/dH)_{T=H=0} = 4.74 \pm 0.15$. For SPAF-3 we provided evidence that the zero-field correlation length diverges as $ξくしー(T \to 0, H=0) \simeq \exp (a/T^{x})$, with $x=1.08 \pm 0.13$, through analysis of the critical curve at $H \neq 0$ plus crossover arguments. For SPAF-3 we have also ascertained that the conformal anomaly and decay-of-correlations exponent behave as: (a) H=0: $c=1, ηいーた=1/3$; (b) $H \neq 0: c=1/2, ηいーた=1/4$. △ Less

Submitted 21 December, 1998; originally announced December 1998.

Comments: RevTex, 7 pages, 4 eps figures, to be published in Phys. Rev. E

Journal ref: Physical Review E, vol. 59, pp. 2772-2777 (1999)

Metal-Insulator transition in one-dimensional Hubbard superlattices

Authors: Thereza Paiva, Raimundo R. dos Santos

Abstract: We study the Metal-Insulator transition in one-dimensional Hubbard superlattices (SL's), modelled by a repeated pattern of repulsive (i.e., positive on-site coupling) and free sites. The evolution of the local moment and of the charge gap (calculated from Lanczos diagonalization of chains up to 18 sites), together with a strong coupling analysis, show that the electron density at which the syste… ▽ More We study the Metal-Insulator transition in one-dimensional Hubbard superlattices (SL's), modelled by a repeated pattern of repulsive (i.e., positive on-site coupling) and free sites. The evolution of the local moment and of the charge gap (calculated from Lanczos diagonalization of chains up to 18 sites), together with a strong coupling analysis, show that the electron density at which the system is insulating increases with the size of the free layer, relative to the repulsive one. In the insulating state, the mechanism of interaction between fermions separated by a free layer is the analog of superexchange, and the charge gap display universal features. △ Less

Submitted 22 July, 1998; originally announced July 1998.

Comments: 4 pages, 3 colour figures; Phys Rev B 15 Sept (tent)

Journal ref: Phys. Rev B 58, 9607 (1998)

Universality and logarithmic corrections in two-dimensional random Ising ferromagnets

Authors: F. D. A. Aarao Reis, S. L. A. de Queiroz, Raimundo R. dos Santos

Abstract: We address the question of weak versus strong universality scenarios for the random-bond Ising model in two dimensions. A finite-size scaling theory is proposed, which explicitly incorporates $\ln L$ corrections ($L$ is the linear finite size of the system) to the temperature derivative of the correlation length. The predictions are tested by considering long, finite-width strips of Ising spins… ▽ More We address the question of weak versus strong universality scenarios for the random-bond Ising model in two dimensions. A finite-size scaling theory is proposed, which explicitly incorporates $\ln L$ corrections ($L$ is the linear finite size of the system) to the temperature derivative of the correlation length. The predictions are tested by considering long, finite-width strips of Ising spins with randomly distributed ferromagnetic couplings, along which free energy, spin-spin correlation functions and specific heats are calculated by transfer-matrix methods. The ratio $γがんま/νにゅー$ is calculated and has the same value as in the pure case; consequently conformal invariance predictions remain valid for this type of disorder. Semilogarithmic plots of correlation functions against distance yield average correlation lengths $ξくしー^{av}$, whose size dependence agrees very well with the proposed theory. We also examine the size dependence of the specific heat, which clearly suggestsa divergency in the thermodynamic limit. Thus our data consistently favour the Dotsenko-Shalaev picture of logarithmic corrections (enhancements) to pure system singularities, as opposed to the weak universality scenario. △ Less

Submitted 19 June, 1997; originally announced June 1997.

Comments: RevTeX code for 8 pages, plus 3 Postscript figures; to appear in Physical Review B (1 September 1997)

Report number: IF/UFF preprint 017/97 (CM 012/97)

On weak vs. strong universality in the two-dimensional random-bond Ising ferromagnet

Authors: F. D. A. Aarão Reis, S. L. A. de Queiroz, Raimundo R. dos Santos

Abstract: We address the issue of universality in two-dimensional disordered Ising systems, by considering long, finite-width strips of ferromagnetic Ising spins with randomly distributed couplings. We calculate the free energy and spin-spin correlation functions (from which averaged correlation lengths, $ξくしー^{ave}$, are computed) by transfer-matrix methods. An {\it ansatz} for the size-dependence of logari… ▽ More We address the issue of universality in two-dimensional disordered Ising systems, by considering long, finite-width strips of ferromagnetic Ising spins with randomly distributed couplings. We calculate the free energy and spin-spin correlation functions (from which averaged correlation lengths, $ξくしー^{ave}$, are computed) by transfer-matrix methods. An {\it ansatz} for the size-dependence of logarithmic corrections to $ξくしー^{ave}$ is proposed. Data for both random-bond and site-diluted systems show that pure system behaviour (with $νにゅー=1$) is recovered if these corrections are incorporated, discarding the weak--universality scenario. △ Less

Submitted 20 August, 1996; originally announced August 1996.

Comments: RevTeX code, 4 pages plus 2 Postscript figures; to appear in Physical Review B Rapid Communications

Report number: IF/UFF preprint 021/96 (CM 014/96)

Bond-charge Interaction in the extended Hubbard chain

Authors: Daisy M. Luz, Raimundo R. dos Santos

Abstract: We study the effects of bond-charge interaction (or correlated hopping) on the properties of the extended ({\it i.e.,} with both on-site ($U$) and nearest-neighbor ($V$) repulsions) Hubbard model in one dimension at half-filling. Energy gaps and correlation functions are calculated by Lanczos diagonalization on finite systems. We find that, irrespective of the sign of the bond-charge interaction… ▽ More We study the effects of bond-charge interaction (or correlated hopping) on the properties of the extended ({\it i.e.,} with both on-site ($U$) and nearest-neighbor ($V$) repulsions) Hubbard model in one dimension at half-filling. Energy gaps and correlation functions are calculated by Lanczos diagonalization on finite systems. We find that, irrespective of the sign of the bond-charge interaction, $X$, the charge--density-wave (CDW) state is more robust than the spin--density-wave (SDW) state. A small bond-charge interaction term is enough to make the differences between the CDW and SDW correlation functions much less dramatic than when $X=0$. For $X=t$ and fixed $V<2t$ ($t$ is the uncorrelated hopping integral), there is an intermediate phase between a charge ordered phase and a phase corresponding to singly-occupied sites, the nature of which we clarify: it is characterized by a succession of critical points, each of which corresponding to a different density of doubly-occupied sites. We also find an unusual slowly decaying staggered spin-density correlation function, which is suggestive of some degree of ordering. No enhancement of pairing correlations was found for any $X$ in the range examined. △ Less

Submitted 26 March, 1996; originally announced March 1996.

Comments: 10 pages, 7 PostScript figures, RevTeX 3; to appear in Phys Rev B

Magnetism and Pairing in Hubbard Bilayers.

Authors: Raimundo R dos Santos

Abstract: We study the Hubbard model on a bilayer with repulsive on-site interactions, $U$, in which fermions undergo both intra-plane ($t$) and inter-plane ($t_z$) hopping. This situation is what one would expect in high-temperature superconductors such as YBCO, with two adjacent CuO$_2$ planes. Magnetic and pairing properties of the system are investigated through Quantum Monte Carlo simulations for both… ▽ More We study the Hubbard model on a bilayer with repulsive on-site interactions, $U$, in which fermions undergo both intra-plane ($t$) and inter-plane ($t_z$) hopping. This situation is what one would expect in high-temperature superconductors such as YBCO, with two adjacent CuO$_2$ planes. Magnetic and pairing properties of the system are investigated through Quantum Monte Carlo simulations for both half- and quarter-filled bands. We find that in all cases inter-planar pairing with $d_{x^2-z^2}$ symmetry is dominant over planar pairing with $d_{x^2-y^2}$ symmetry, and that for $t_z$ large enough pair formation is possible through antiferromagnetic correlations. However, another mechanism is needed to make these pairs condense into a superconducting state at lower temperatures. We identify the temperature for pair formation with the spin gap crossover temperature. [Submitted to Phys. Rev. B] △ Less

Submitted 10 February, 1995; originally announced February 1995.

Comments: 7 pages, uuencoded self-unpacking PS file with text and figures.

Report number: IF/UFF-01/95.

Spin gap and superconductivity in the three-dimensional attractive Hubbard model

Authors: Raimundo R dos Santos

Abstract: We study the phase diagram for the attractive ({\it i.e.,} negative-$U$) Hubbard model on a simple cubic lattice, through Monte Carlo simulations. We obtain the critical temperature, $T_c$, for superconductivity from a finite-size scaling analysis of the data for the pairing correlations. For fixed on-site attraction, $U$, $T_c$ displays a maximum near the filling factor 0.9, roughly independent… ▽ More We study the phase diagram for the attractive ({\it i.e.,} negative-$U$) Hubbard model on a simple cubic lattice, through Monte Carlo simulations. We obtain the critical temperature, $T_c$, for superconductivity from a finite-size scaling analysis of the data for the pairing correlations. For fixed on-site attraction, $U$, $T_c$ displays a maximum near the filling factor 0.9, roughly independent of $U$. For fixed filling we estimate the crossover temperature $T^\times(U)$, separating the normal states: metallic and spin-gap. There is also a critical value $U_p$ for pair formation, the magnitude of which seems to be independent of doping. The relevance of these results to the high-$T_c$ oxides is discussed. △ Less

Submitted 26 May, 1994; originally announced May 1994.

Comments: 5 pages, uuencoded PS text and figures, NC-FIS-PUC-09/93