Source author record

Y. Alhassid

Y. Alhassid appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2021arXiv

Nuclear level densities: from empirical models to microscopic methods

The level density is among the most important statistical nuclear properties. It appears in Fermi's golden rule for transition rates and is an important input to the Hauser-Feshbach theory of compound nucleus reactions. We discuss empirical models of level densities and summarize the main experimental methods used to determine them. The microscopic calculation of level densities in the presence of correlations is a challenging many-body problem. We review recent microscopic approaches to calculate level densities. Mean-field and combinatorial methods have been applied across the nuclear chart, but often need to be augmented with empirical collective enhancement factors. The moment method and the auxiliary-field quantum Monte Carlo (AFMC) method are formulated in the context of the configuration-interaction shell model approach, and include correlations beyond the mean-field approximation.

preprint2020arXiv

Derivation of K-matrix reaction theory in a discrete basis formalism

The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann-Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present an alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.

preprint2020arXiv

Pairing correlations across the superfluid phase transition in the unitary Fermi gas

In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature $T_c$ and a temperature $T^* > T_c$. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor $ν\simeq 0.06$ for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate $T_c$ at the corresponding filling factor. We find that the energy-staggering pairing gap vanishes above $T_c$, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.

preprint2020arXiv

Statistical-model description of $γ$ decay from compound-nucleus resonances

The statistical model of compound-nucleus reactions predicts that the fluctuations of the partial $γ$-decay widths for a compound-nucleus resonance are governed by the Porter-Thomas distribution (PTD), and that consequently the distribution of total $γ$-decay widths is very narrow. However, a recent experiment [Koehler, Larsen, Guttormsen, Siem, and Guber, Phys. Rev. C 88, 041305(R) (2013)] reported large fluctuations of the total $γ$-decay widths in the $^{95}$Mo$(n,γ)^{96}$Mo* reaction, contrary to this expectation. Furthermore, in recent theoretical works it was argued that sufficiently strong channel couplings can cause deviations of the partial width distributions from PTD. Here, we investigate whether the combined influence of a large number of nonequivalent $γ$-decay channels, each of which couples weakly to the compound-nucleus resonances, can modify the statistics of the partial widths. We study this effect in neutron scattering off $^{95}$Mo within a random-matrix model that includes coupling to the entrance neutron channel and to the large number of $γ$ channels. Using realistic coupling parameters obtained from empirical models for the level density and the $γ$ strength function, we find that the PTD describes well the distribution of partial widths for all decay channels, in agreement with the statistical-model expectation. Furthermore, we find that the width of the distribution of the total $γ$-decay widths is insensitive to wide variations in the parameters of the $γ$ strength function, as well as to deviations of the partial-width distributions from the PTD. Our results rule out an explanation of the recent experimental data within a statistical-model description of the compound nucleus.

preprint2019arXiv

The contact in the unitary Fermi gas across the superfluid phase transition

A quantity known as the contact plays a fundamental role in quantum many-body systems with short-range interactions. The determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy of the unitary gas as a function of temperature in the continuum limit.

preprint2016arXiv

Microsopic nuclear level densities by the shell model Monte Carlo method

The configuration-interaction shell model approach provides an attractive framework for the calculation of nuclear level densities in the presence of correlations, but the large dimensionality of the model space has hindered its application in mid-mass and heavy nuclei. The shell model Monte Carlo (SMMC) method permits calculations in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We discuss recent progress in the SMMC approach to level densities, and in particular the calculation of level densities in heavy nuclei. We calculate the distribution of the axial quadrupole operator in the laboratory frame at finite temperature and demonstrate that it is a model-independent signature of deformation in the rotational invariant framework of the shell model. We propose a method to use these distributions for calculating level densities as a function of intrinsic deformation.

preprint2016arXiv

The shell model Monte Carlo approach to level densities: recent developments and perspectives

We review recent advances in the shell model Monte Carlo approach for the microscopic calculation of statistical and collective properties of nuclei. We discuss applications to the calculation of (i) level densities in nickel isotopes, implementing a recent method to circumvent the odd-particle sign problem; (ii) state densities in heavy nuclei; (iii) spin distributions of nuclear levels; and (iv) finite-temperature quadrupole distributions.

preprint2015arXiv

Benchmarking mean-field approximations to level densities

We assess the accuracy of finite-temperature mean-field theory using as a standard the Hamiltonian and model space of the shell model Monte Carlo calculations. Two examples are considered: the nucleus $^{162}$Dy, representing a heavy deformed nucleus, and $^{148}$Sm, representing a nearby heavy spherical nucleus with strong pairing correlations. The errors inherent in the finite-temperature Hartree-Fock and Hartree-Fock-Bogoliubov approximations are analyzed by comparing the entropies of the grand canonical and canonical ensembles, as well as the level density at the neutron resonance threshold, with shell model Monte Carlo (SMMC) calculations, which are accurate up to well-controlled statistical errors. The main weak points in the mean-field treatments are seen to be: (i) the extraction of number-projected densities from the grand canonical ensembles, and (ii) the symmetry breaking by deformation or by the pairing condensate. In the absence of a pairing condensate, we confirm that the usual saddle-point approximation to extract the number-projected densities is not a significant source of error compared to other errors inherent to the mean-field theory. We also present an alternative formulation of the saddle-point approximation that makes direct use of an approximate particle-number projection and avoids computing the usual three-dimensional Jacobian of the saddle-point integration. We find that the pairing condensate is less amenable to approximate particle-number projection methods due to the explicit violation of particle-number conservation in the pairing condensate. Nevertheless, the Hartree-Fock-Bogoliubov theory is accurate to less than one unit of entropy for $^{148}$Sm at the neutron threshold energy, which is above the pairing phase transition.

preprint2015arXiv

Collective enhancement of nuclear state densities by the shell model Monte Carlo approach

The shell model Monte Carlo (SMMC) approach allows for the microscopic calculation of statistical and collective properties of heavy nuclei using the framework of the configuration-interaction shell model in very large model spaces. We present recent applications of the SMMC method to the calculation of state densities and their collective enhancement factors in rare-earth nuclei.

preprint2015arXiv

Magnetic response of energy levels of superconducting nanoparticles with spin-orbit scattering

Discrete energy levels of ultrasmall metallic grains are extracted in single-electron-tunneling-spectroscopy experiments. We study the response of these energy levels to an external magnetic field in the presence of both spin-orbit scattering and pairing correlations. In particular, we investigate $g$-factors and level curvatures that parametrize, respectively, the linear and quadratic terms in the magnetic-field dependence of the many-particle energy levels of the grain. Both of these quantities exhibit level-to-level fluctuations in the presence of spin-orbit scattering. We show that the distribution of $g$-factors is not affected by the pairing interaction and that the distribution of level curvatures is sensitive to pairing correlations even in the smallest grains in which the pairing gap is smaller than the mean single-particle level spacing. We propose the level curvature in a magnetic field as a tool to probe pairing correlations in tunneling spectroscopy experiments.

preprint2015arXiv

Recent Advances in the Application of the Shell Model Monte Carlo Approach to Nuclei

The shell model Monte Carlo (SMMC) method is a powerful technique for calculating the statistical and collective properties of nuclei in the presence of correlations in model spaces that are many orders of magnitude larger than those that can be treated by conventional diagonalization methods. We review recent advances in the development and application of SMMC to mid-mass and heavy nuclei.

preprint2014arXiv

Mesoscopic superconductivity in ultrasmall metallic grains

A nano-scale metallic grain (nanoparticle) with irregular boundaries in which the single-particle dynamics are chaotic is a zero-dimensional system described by the so-called universal Hamiltonian in the limit of a large number of electrons. The interaction part of this Hamiltonian includes a superconducting pairing term and a ferromagnetic exchange term. Spin-orbit scattering breaks spin symmetry and suppresses the exchange interaction term. Of particular interest is the fluctuation-dominated regime, typical of the smallest grains in the experiments, in which the bulk pairing gap is comparable to or smaller than the single-particle mean-level spacing, and the Bardeen-Cooper-Schrieffer (BCS) mean-field theory of superconductivity is no longer valid. Here we study the crossover between the BCS and fluctuation-dominated regimes in two limits. In the absence of spin-orbit scattering, the pairing and exchange interaction terms compete with each other. We describe the signatures of this competition in thermodynamic observables, the heat capacity and spin susceptibility. In the presence of strong spin-orbit scattering, the exchange interaction term can be ignored. We discuss how the magnetic-field response of discrete energy levels in such a nanoparticle is affected by pairing correlations. We identify signatures of pairing correlations in this response, which are detectable even in the fluctuation-dominated regime.

preprint2014arXiv

Nuclear deformation at finite temperature

Deformation, a key concept in our understanding of heavy nuclei, is based on a mean-field description that breaks the rotational invariance of the nuclear many-body Hamiltonian. We present a method to analyze nuclear deformations at finite temperature in a framework that preserves rotational invariance. The auxiliary-field Monte-Carlo method is used to generate the statistical ensemble and calculate the probability distribution associated with the quadrupole operator. Applying the technique to nuclei in the rare-earth region, we identify model-independent signatures of deformation and find that deformation effects persist to higher temperatures than the spherical-to-deformed shape phase-transition temperature of mean-field theory.

preprint2014arXiv

Recent Advances in the Microscopic Calculations of Level Densities by the Shell Model Monte Carlo Method

The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (i) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (ii) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes $^{59-64}$Ni and of a heavy deformed rare-earth nucleus $^{162}$Dy and found them to be in close agreement with various experimental data sets.

preprint2014arXiv

Stabilizing Canonical-Ensemble Calculations in the Auxiliary-Field Monte Carlo Method

Quantum Monte Carlo methods are powerful techniques for studying strongly interacting Fermi systems. However, implementing these methods on computers with finite-precision arithmetic requires careful attention to numerical stability. In the auxiliary-field Monte Carlo (AFMC) method, low-temperature or large-model-space calculations require numerically stabilized matrix multiplication. When adapting methods used in the grand-canonical ensemble to the canonical ensemble of fixed particle number, the numerical stabilization increases the number of required floating-point operations for computing observables by a factor of the size of the single-particle model space, and thus can greatly limit the systems that can be studied. We describe an improved method for stabilizing canonical-ensemble calculations in AFMC that exhibits better scaling, and present numerical tests that demonstrate the accuracy and improved performance of the method.

preprint2013arXiv

Collectivity in Heavy Nuclei in the Shell Model Monte Carlo Approach

The microscopic description of collectivity in heavy nuclei in the framework of the configuration-interaction shell model has been a major challenge. The size of the model space required for the description of heavy nuclei prohibits the use of conventional diagonalization methods. We have overcome this difficulty by using the shell model Monte Carlo (SMMC) method, which can treat model spaces that are many orders of magnitude larger than those that can be treated by conventional methods. We identify a thermal observable that can distinguish between vibrational and rotational collectivity and use it to describe the crossover from vibrational to rotational collectivity in families of even-even rare-earth isotopes. We calculate the state densities in these nuclei and find them to be in close agreement with experimental data. We also calculate the collective enhancement factors of the corresponding level densities and find that their decay with excitation energy is correlated with the pairing and shape phase transitions.

preprint2013arXiv

Configuration-interaction Monte Carlo method and its application to the trapped unitary Fermi gas

We develop a quantum Monte Carlo method to estimate the ground-state energy of a fermionic many-particle system in the configuration-interaction shell model approach. The fermionic sign problem is circumvented by using a guiding wave function in Fock space. The method provides an upper bound on the ground-state energy whose tightness depends on the choice of the guiding wave function. We argue that the antisymmetric geminal product class of wave functions is a good choice for guiding wave functions. We demonstrate our method for the trapped two-species fermionic cold atom system in the unitary regime of infinite scattering length using the particle-number projected Hartree-Fock-Bogoliubov wave function as the guiding wave function. We estimate the ground-state energy and energy-staggering pairing gap as a function of the number of particles. Our results compare favorably with exact numerical diagonalization results and with previous coordinate-space Monte Carlo calculations.

preprint2013arXiv

Crossover from vibrational to rotational collectivity in heavy nuclei in the shell-model Monte Carlo approach

Heavy nuclei exhibit a crossover from vibrational to rotational collectivity as the number of neutrons or protons increases from shell closure towards midshell, but the microscopic description of this crossover has been a major challenge. We apply the shell model Monte Carlo approach to families of even-even samarium and neodymium isotopes and identify a microscopic signature of the crossover from vibrational to rotational collectivity in the low-temperature behavior of $<\mathbf{J}^2>_T$, where $\bf J$ is the total spin and $T$ is the temperature. This signature agrees well with its values extracted from experimental data. We also calculate the state densities of these nuclei and find them to be in very good agreement with experimental data. Finally, we define a collective enhancement factor from the ratio of the total state density to the intrinsic state density as calculated in the finite-temperature Hartree-Fock-Bogoliubov approximation. The decay of this enhancement factor with excitation energy is found to correlate with the pairing and shape phase transitions in these nuclei.

preprint2013arXiv

Direct microscopic calculations of nuclear level densities in the shell model Monte Carlo approach

Nuclear level densities are crucial for estimating statistical nuclear reaction rates. The shell model Monte Carlo method is a powerful approach for microscopic calculation of state densities in very large model spaces. However, these state densities include the spin degeneracy of each energy level, whereas experiments often measure level densities in which each level is counted just once. To enable the direct comparison of theory with experiments, we introduce a method to calculate directly the level density in the shell model Monte Carlo approach. The method employs a projection on the minimal absolute value of the magnetic quantum number. We apply the method to nuclei in the iron region as well as the strongly deformed rare-earth nucleus $^{162}$Dy. We find very good agreement with experimental data including level counting at low energies, charged particle spectra and Oslo method at intermediate energies, neutron and proton resonance data, and Ericson's fluctuation analysis at higher excitation energies.

preprint2013arXiv

Level Densities of Heavy Nuclei by the Shell Model Monte Carlo Method

The microscopic calculation of nuclear level densities in the presence of correlations is a difficult many-body problem. The shell model Monte Carlo method provides a powerful technique to carry out such calculations using the framework of the configuration-interaction shell model in spaces that are many orders of magnitude larger than spaces that can be treated by conventional methods. We present recent applications of the method to the calculation of level densities and their collective enhancement factors in heavy nuclei. The calculated level densities are in close agreement with experimental data.

preprint2013arXiv

Level densities of nickel isotopes: microscopic theory versus experiment

We apply a spin-projection method to calculate microscopically the level densities of a family of nickel isotopes $^{59-64}$Ni using the shell model Monte Carlo approach in the complete $pfg_{9/2}$ shell. Accurate ground-state energies of the odd-mass nickel isotopes, required for the determination of excitation energies, are determined using the Green's function method recently introduced to circumvent the odd particle-number sign problem. Our results are in excellent agreement with recent measurements based on proton evaporation spectra and with level counting data at low excitation energies. We also compare our results with neutron resonance data, assuming equilibration of parity and a spin-cutoff model for the spin distribution at the neutron binding energy, and find good agreement with the exception of $^{63}$Ni.

preprint2013arXiv

Nuclear state densities of odd-mass heavy nuclei in the shell model Monte Carlo approach

While the shell model Monte Carlo approach has been successful in the microscopic calculation of nuclear state densities, it has been difficult to calculate accurately state densities of odd-even heavy nuclei. This is because the projection on an odd number of particles in the shell model Monte Carlo method leads to a sign problem at low temperatures, making it impractical to extract the ground-state energy in direct Monte Carlo calculations. We show that the ground-state energy can be extracted to a good precision by using level counting data at low excitation energies and the neutron resonance data at the neutron threshold energy. This allows us to extend recent applications of the shell-model Monte Carlo method in even-even rare-earth nuclei to the odd-even isotopic chains of $^{149-155}$Sm and $^{143-149}$Nd. We calculate the state densities of the odd-even samarium and neodymium isotopes and find close agreement with the state densities extracted from experimental data.

preprint2013arXiv

Pair condensation in a Finite Trapped Fermi Gas

Superfluidity in the cold atomic two-species Fermi gas system in the unitary limit of infinite scattering length remains incompletely understood. In particular, a pseudogap phase has been proposed to exist above the superfluid critical temperature. Here we apply the auxiliary-field quantum Monte Carlo method to perform the first ab initio calculations of the temperature dependence of three quantities -- the energy-staggering pairing gap, the condensate fraction and the heat capacity -- in a trapped finite-size cold atom system. As the calculations of the energy-staggering pairing gap require the use of the canonical ensemble, we employ a novel algorithm for the stabilization of particle-number projection that is essential for reaching convergence in the size of the model space. We observe clear signatures of the superfluid phase transition in all three quantities, including a signature of the recently measured lambda peak in the heat capacity, but find no evidence of a pseudogap effect in the energy-staggering pairing gap.

preprint2013arXiv

The coexistence of superconductivity and ferromagnetism in nano-scale metallic grains

A nano-scale metallic grain in which the single-particle dynamics are chaotic is described by the so-called universal Hamiltonian. This Hamiltonian includes a superconducting pairing term and a ferromagnetic exchange term that compete with each other: pairing correlations favor minimal ground-state spin, while the exchange interaction favors maximal spin polarization. Of particular interest is the fluctuation-dominated regime where the bulk pairing gap is comparable to or smaller than the single-particle mean level spacing and the Bardeen-Cooper-Schrieffer theory of superconductivity breaks down. Superconductivity and ferromagnetism can coexist in this regime. We identify signatures of the competition between superconductivity and ferromagnetism in a number of quantities: ground-state spin, conductance fluctuations when the grain is weakly coupled to external leads and the thermodynamic properties of the grain, such as heat capacity and spin susceptibility.

preprint2013arXiv

Thermodynamics of ultrasmall metallic grains in the presence of pairing and exchange correlations: Mesoscopic fluctuations

We study the mesoscopic fluctuations of thermodynamic observables in a nanosized metallic grain in which the single-particle dynamics are chaotic and the dimensionless Thouless conductance is large. Such a grain is modeled by the universal Hamiltonian describing the competition between exchange and pairing correlations. The exchange term is taken into account exactly by a spin-projection method, and the pairing term is treated in the static-path approximation together with small-amplitude quantal fluctuations around each static fluctuation of the pairing field. Odd-even particle-number effects induced by pairing correlations are included using a number-parity projection. We find that the exchange interaction shifts the number-parity effects in the heat capacity and spin susceptibility to lower temperatures. In the regime where the pairing gap is similar to or smaller than the single-particle mean level spacing, these number-parity effects are suppressed by exchange correlations, and the fluctuations of the spin susceptibility may be particularly large. However, for larger values of the pairing gap, the number-parity effects may be enhanced by exchange correlations.

preprint2012arXiv

Odd-particle systems in the shell model Monte Carlo: circumventing a sign problem

We introduce a novel method within the shell model Monte Carlo approach to calculate the ground-state energy of a finite-size system with an odd number of particles by using the asymptotic behavior of the imaginary-time single-particle Green's functions. The method circumvents the sign problem that originates from the projection on an odd number of particles and has hampered direct application of the shell model Monte Carlo method to odd-particle systems. We apply this method to calculate pairing gaps of nuclei in the iron region. Our results are in good agreement with experimental pairing gaps.

preprint2012arXiv

Recent developments in the shell model Monte Carlo approach to nuclei

The shell model Monte Carlo (SMMC) approach provides a powerful method for the microscopic calculation of statistical and collective nuclear properties in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods. We discuss recent applications of the method to describe the emergence of collectivity in the framework of the configuration-interaction shell model and the crossover from vibrational to rotational collectivity in families of rare-earth nuclei. We have calculated state densities of these rare-earth nuclei and find their collective enhancement factors to be correlated with the pairing and shape phase transitions. We also discuss an accurate method to calculate the ground-state energy of odd-even and odd-odd nuclei, circumventing the sign problem that originates in the projection on an odd number of particles. We have applied this method to calculate pairing gaps in families of isotopes in the iron region.

preprint2012arXiv

Signatures of phase transitions in nuclei at finite excitation energies

The mean-field approximation predicts pairing and shape phase transitions in nuclei as a function of temperature or excitation energy. However, in the finite nucleus the singularities of these phase transitions are smoothed out by quantal and thermal fluctuations. An interesting question is whether signatures of these transitions survive despite the large fluctuations. The shell model Monte Carlo (SMMC) approach enables us to calculate the statistical properties of nuclei beyond the mean-field approximation in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We have extended the SMMC method to heavy nuclei and used it to study the transition from vibrational (spherical) to rotational (deformed) nuclei in families of rare-earth isotopes. We have calculated collective enhancement factors of level densities as a function of excitation energy and found that the decay of the vibrational and rotational enhancements is well correlated with the pairing and shape phase transitions, respectively.

preprint2012arXiv

Thermal signatures of pairing correlations in nuclei and nano-scale metallic grains

Atomic nuclei and nano-scale metallic grains are in the crossover regime of pairing correlations between the bulk limit, where the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity is valid, and the fluctuation-dominated regime, where BCS theory breaks down. In this fluctuation-dominated regime, the pairing gap is comparable to or smaller than the single-particle mean level spacing. We discuss thermal signatures of pairing correlations in nuclei and ultra-small metallic grains that survive despite the large fluctuations of the pairing field.

preprint2011arXiv

A new effective interaction for the trapped fermi gas: the BEC-BCS crossover

We extend a recently introduced separable interaction for the unitary trapped Fermi gas to all values of the scattering length. We derive closed expressions for the interaction matrix elements and the two-particle eigenvectors and analytically demonstrate the convergence of this interaction to the zero-range two-body pseudopotential for s-wave scattering. We apply this effective interaction to the three- and four-particle systems along the BEC-BCS crossover, and find that their low-lying energies exhibit convergence in the regularization parameter that is much faster than for the conventional renormalized contact interaction. We find similar convergence properties of the three-particle free energy at unitarity.

preprint2010arXiv

Number-conserving theory of nuclear pairing gaps: a global assessment

We study odd-even mass staggering of nuclei, also called pairing gaps, using a Skyrme self-consistent mean-field theory and a numerically exact treatment of the pairing Hamiltonian. We find that the configuration-space Monte Carlo method proposed by Cerf and Martin offers a practical computational procedure to carry out the numerical solutions in large-dimensional model spaces. Refitting the global strength of the pairing interaction for 443 neutron pairing gaps in our number-conserving treatment, we find the correction to the pairing correlation energies and pairing gaps to have rms values of 0.6 MeV and 0.12 MeV, respectively. The exact treatment provides a significant improvement in the fit to experimental gaps, although it is partially masked by a larger rms error due to deficiencies in other aspects of the theory such as the mean-field energy functional.

preprint2010arXiv

Signatures of exchange correlations in the thermopower of quantum dots

We use a many-body rate-equation approach to calculate the thermopower of a quantum dot in the presence of an exchange interaction. At temperatures much smaller than the single-particle level spacing, the known quantum jumps (discontinuities) in the thermopower are split by the exchange interaction. The origin and nature of the splitting are elucidated with a simple physical argument based on the nature of the intermediate excited state in the sequential tunneling approach. We show that this splitting is sensitive to the number parity of electrons in the dot and the dot's ground-state spin. These effects are suppressed when cotunneling dominates the electrical and thermal conductances. We calculate the thermopower in the presence of elastic cotunneling, and show that some signatures of exchange correlations should still be observed with current experimental methods. In particular, we propose a method to determine the strength of the exchange interaction from measurements of the thermopower.

preprint2010arXiv

The competition between superconductivity and ferromagnetism in small metallic grains: thermodynamic properties

We study the thermodynamic properties of a small superconducting metallic grain using a quantum Monte Carlo method. The grain is described by the universal Hamiltonian, containing pairing and ferromagnetic exchange correlations. In particular, we study how the thermodynamic signatures of pairing correlations are affected by the spin exchange interaction. We find the exchange interaction effects to be qualitatively different in the BCS and fluctuation-dominated regimes of pairing correlations.

preprint2009arXiv

Interaction Matrix Element Fluctuations in Ballistic Quantum Dots: Dynamical Effects

We study matrix element fluctuations of the two-body screened Coulomb interaction and of the one-body surface charge potential in ballistic quantum dots, comparing behavior in actual chaotic billiards with analytic results previously obtained in a normalized random wave model. We find that the matrix element variances in actual chaotic billiards typically exceed by a factor of 3 or 4 the predictions of the random wave model, for dot sizes commonly used in experiments. We discuss dynamical effects that are responsible for this enhancement. These dynamical effects have an even more striking effect on the covariance, which changes sign when compared with random wave predictions. In billiards that do not display hard chaos, an even larger enhancement of matrix element fluctuations is possible. These enhanced fluctuations have implications for peak spacing statistics and spectral scrambling for quantum dots in the Coulomb blockade regime.

preprint2009arXiv

Mesoscopic interplay of superconductivity and ferromagnetism in ultra-small metallic grains

We review the effects of electron-electron interactions on the ground-state spin and the transport properties of ultra-small chaotic metallic grains. Our studies are based on an effective Hamiltonian that combines a superconducting BCS-like term and a ferromagnetic Stoner-like term. Such terms originate in pairing and spin exchange correlations, respectively. This description is valid in the limit of a large dimensionless Thouless conductance. We present the ground-state phase diagram in the fluctuation-dominated regime where the single-particle mean level spacing is comparable to the bulk BCS pairing gap. This phase diagram contains a regime in which pairing and spin exchange correlations coexist in the ground-state wave function. We discuss the calculation of the tunneling conductance for an almost-isolated grain in the Coulomb-blockade regime, and present measurable signatures of the competition between superconductivity and ferromagnetism in the mesoscopic fluctuations of the conductance.

preprint2007arXiv

A new effective interaction for the trapped Fermi gas

We apply the configuration-interaction method to calculate the spectra of two-component Fermi systems in a harmonic trap, studying the convergence of the method at the unitary interaction limit. We find that for a fixed regularization of the two-body interaction the convergence is exponential or better in the truncation parameter of the many-body space. However, the conventional regularization is found to have poor convergence in the regularization parameter, with an error that scales as a low negative power of this parameter. We propose a new regularization of the two-body interaction that produces exponential convergence for systems of three and four particles. From the systematics, we estimate the ground-state energy of the four-particle system to be (5.05 +- 0.024)hbar omega.

preprint2007arXiv

Effect of a Zeeman field on the superconductor-ferromagnet transition in metallic grains

We investigate the competition between pairing correlations and ferromagnetism in small metallic grains in the presence of a Zeeman field. Our analysis is based on the universal Hamiltonian, valid in the limit of large Thouless conductance. We show that the coexistence regime of superconducting and ferromagnetic correlations can be made experimentally accessible by tuning an external Zeeman field. We compare the exact solution of the model with a mean-field theory and find that the latter cannot describe pairing correlations in the intermediate regime. We also study the occurrence of spin jumps across the phase boundary separating the superconducting and coexistence regimes.

preprint2003arXiv

Fano interference and cross-section fluctuations in molecular photodissociation

We derive an expression for the total photodissociation cross section of a molecule incorporating both indirect processes that proceed through excited resonances, and direct processes. We show that this cross section exhibits generalized Beutler-Fano line shapes in the limit of isolated resonances. Assuming that the closed system can be modeled by random matrix theory, we derive the statistical properties of the photodissociation cross section and find that they are significantly affected by the direct processes. We identify a unique signature of the direct processes in the cross-section distribution in the limit of isolated resonances.

preprint2003arXiv

The Importance of Parity-Dependence of the Nuclear Level Density in the Prediction of Astrophysical Reaction Rates

A simple description for obtaining the parity distribution of nuclear levels in the pf + g9/2 shell as a function of excitation energy was recently derived. We implement this in a global nuclear level density model. In the framework of the statistical model, cross sections and astrophysical reaction rates are calculated in the Fe region and compared to rates obtained with the common assumption of an equal distribution of parities. We find considerable differences, especially for reactions involving particles in the exit channel.

preprint2002arXiv

Statistical Fluctuations of Electromagnetic Transition Intensities and Electromagnetic Moments in pf-Shell Nuclei

We study the fluctuation properties of $ΔT=0$ electromagnetic transition intensities and electromagnetic moments in $A \sim 60$ nuclei within the framework of the interacting shell model, using a realistic effective interaction for $pf$-shell nuclei with a $^{56}$Ni core. The distributions of the transition intensities and of the electromagnetic moments are well described by the Gaussian orthogonal ensemble of random matrices. In particular, the transition intensity distributions follow a Porter-Thomas distribution. When diagonal matrix elements (i.e., moments) are included in the analysis of transition intensities, we find that the distributions remain Porter-Thomas except for the isoscalar $M1$. The latter deviation is explained in terms of the structure of the isoscalar $M1$ operator.

preprint1998arXiv

The Spectral Autocorrelation Function in Weakly Open Chaotic Systems: Indirect Photodissociation of Molecules

We derive the statistical limit of the spectral autocorrelation function and of the survival probability for the indirect photodissociation of molecules in the regime of non-overlapping resonances. The results are derived in the framework of random matrix theory, and hold more generally for any chaotic quantum system that is weakly coupled to the continuum. The "correlation hole" that characterizes the spectral autocorrelation in the bound molecule diminishes as the typical average total width of a resonance increases.

preprint1995arXiv

Universal Parametric Correlations of Eigenfunctions in Chaotic and Disordered Systems

This paper establishes the universality of parametric correlations of eigenfunctions in chaotic and weakly disordered systems. We demonstrate this universality in the framework of the gaussian random matrix process and obtain predictions for a number of parametric correlators, one of them analytically. We present numerical evidence from different models that verifies our predictions.