Source author record

J. Dukelsky

J. Dukelsky 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

33works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

33 published item(s)

preprint2016arXiv

Structure of the number projected BCS wave function

We study the structure of the number projected BCS (PBCS) wave function in the particle-hole basis, displaying its similarities with coupled clusters theory (CCT). The analysis of PBCS together with several modifications suggested by the CCT wave function is carried out for the exactly solvable Richardson model involving a pure pairing hamiltonian acting in a space of equally-spaced doubly-degenerate levels. We point out the limitations of PBCS to describe the non-superconducting regime and suggest possible avenues for improvement.

preprint2015arXiv

Advanced density matrix renormalization group method for nuclear structure calculations

We present an efficient implementation of the Density Matrix Renormalization Group (DMRG) algorithm that includes an optimal ordering of the proton and neutron orbitals and an efficient expansion of the active space utilizing various concepts of quantum information theory. We first show how this new DMRG methodology could solve a previous $400$ KeV discrepancy in the ground state energy of $^{56}$Ni. We then report the first DMRG results in the $pf+g9/2$ shell model space for the ground $0^+$ and first $2^+$ states of $^{64}$Ge which are benchmarked with reference data obtained from Monte Carlo shell model. The corresponding correlation structure among the proton and neutron orbitals is determined in terms of the two-orbital mutual information. Based on such correlation graphs we propose several further algorithmic improvement possibilities that can be utilized in a new generation of tensor network based algorithms.

preprint2014arXiv

Disentangling phase transitions and critical points in the proton-neutron interacting boson model by catastrophe theory

We introduce the basic concepts of catastrophe theory needed to derive analytically the phase diagram of the proton-neutron interacting boson model (IBM-2). Previous studies [1,2,3] were based on numerical solutions. We here explain the whole IBM-2 phase diagram including the precise order of the phase transitions in terms of the cusp catastrophe.

preprint2012arXiv

Alpha-like quartet condensation and isovector pairing correlations in N=Z nuclei

We propose a simple quartet condensation model (QCM) which describes with very high accuracy the isovector pairing correlations in self-conjugate nuclei. The quartets have an alpha-like structure and are formed by collective isovector pairs. The accuracy of the QCM is tested for N=Z nuclei for which exact shell model diagonalizations can be performed. The calculations are done with two isovector pairing forces, one extracted from standard shell model interactions and the other of seniority type, acting, respectively, upon spherical and axially-deformed single-particle states. It is shown that for all calculated nuclei the QCM gives very accurate values for the pairing correlations energies, with errors which do not exceed 1%. These results show clearly that the correlations induced by the isovector pairing in self-conjugate nuclei are of quartet type and also indicate that QCM is the proper tool to calculate the isovector proton-neutron correlations in mean field pairing models.

preprint2012arXiv

The Lipkin-Meshkov-Glick model as a particular limit of the SU(1,1) Richardson-Gaudin integrable models

The Lipkin-Meshkov-Glick (LMG) model has a Schwinger boson realization in terms of a two-level boson pairing Hamiltonian. Through this realization, it has been shown that the LMG model is a particular case of the SU (1, 1) Richardson-Gaudin (RG) integrable models. We exploit the exact solvability of the model tostudy the behavior of the spectral parameters (pairons) that completely determine the wave function in the different phases, and across the phase transitions. Based on the relation between the Richardson equations and the Lamé differential equations we develop a method to obtain numerically the pairons. The dynamics of pairons in the ground and excited states provides new insights into the first, second and third order phase transitions, as well as into the crossings taking place in the LMG spectrum.

preprint2011arXiv

Comment on "Fermi-Bose Mixtures near Broad Interspecies Feshbach Resonances"

In a recent Letter (Phys. Rev. lett. 105, 195301 (2010)), Song et al. introduced a new variational approach to treat strong attractive boson-fermion (BF) correlations in BF atomic mixtures. The proposed theory predicts a first order phase transition to a condensate of composite BF pairs with center of mass momentum Q=0 as opposed to a composite fermionic molecular Fermi gas. We show in this comment that their approach is incorrect and moreover, by resorting to an exactly solvable model we demonstrate that there cannot be more than one correlated Q=0 BF pair in complete contradiction with their conclusions.

preprint2011arXiv

Exactly solvable pairing Hamiltonian for heavy nuclei

We present a new exactly solvable Hamiltonian with a separable pairing interaction and non-degenerate single-particle energies. It is derived from the hyperbolic family of Richardson-Gaudin models and possesses two free parameters, one related to an interaction cutoff and the other to the pairing strength. These two parameters can be adjusted to give an excellent reproduction of Gogny self-consistent mean-field calculations in the canonical basis.

preprint2011arXiv

Integrable two-channel p_x+ip_y-wave superfluid model

We present a new two-channel integrable model describing a system of spinless fermions interacting through a p-wave Feshbach resonance. Unlike the BCS-BEC crossover of the s-wave case, the p-wave model has a third order quantum phase transition. The critical point coincides with the deconfinement of a single molecule within a BEC of bound dipolar molecules. The exact many-body wavefunction provides a unique perspective of the quantum critical region suggesting that the size of the condensate wavefunction, that diverges logarithmically with the chemical potential, could be used as an experimental indicator of the phase transition.

preprint2010arXiv

Comment on "Quantum phase transition in the four-spin exchange antiferromagnet" by V. N. Kotov, D.-X. Yao, A. H. Castro-Neto, and D. K. Campbell

In a recent paper [Phys. Rev. {\bf B80}, 174403 (2009)] Kotov {\it et al.} studied the paramagnetic-to-antiferromagnetic transition in the $J$-$Q$ model. Their findings were claimed to be in "fairly good agreement" with previous quantum Monte-Carlo (QMC) results. In this Comment we show that the above claim is misleading and in reality their phase transition point is not only far from the corresponding QMC value but also lies in a region of parameter space not yet explored in the literature. We also show that their reference dimer state is unstable against formation of a plaquette condensate, which could in part explain the large fluctuations they found.

preprint2010arXiv

Competition between Normal Superfluidity and Larkin-Ovchinnikov Phases of Polarized Fermi Gases in Elongated Traps

By applying the recently proposed antisymmetric superfluid local density approximation (ASLDA) to strongly interacting polarized atomic gases at unitarity in very elongated traps, we find families of Larkin-Ovchinnikov (LO) type of solutions with prominent transversal oscillation of pairing potential. These LO states coexist with a superfluid state having a smooth pairing potential. We suggest that the LO phase could be accessible experimentally by increasing adiabatically the trap aspect ratio. We show that the local asymmetry effects contained in ASLDA do not support a deformed superfluid core predicted by previous Bogoliubov-de Gennes treatments.

preprint2010arXiv

Decoherence as a signature of an excited state quantum phase transition

We analyze the decoherence induced on a single qubit by the interaction with a two-level boson system with critical internal dynamics. We explore how the decoherence process is affected by the presence of quantum phase transitions in the environment. We conclude that the dynamics of the qubit changes dramatically when the environment passes through a continuous excited state quantum phase transition. If the system-environment coupling energy equals the energy at which the environment has a critical behavior, the decoherence induced on the qubit is maximal and the fidelity tends to zero with finite size scaling obeying a power-law.

preprint2010arXiv

Decoherence due to an excited state quantum phase transition in a two-level boson model

The decoherence induced on a single qubit by its interaction with the environment is studied. The environment is modelled as a scalar two-level boson system that can go through either first order or continuous excited state quantum phase transitions, depending on the values of the control parameters. A mean field method based on the Tamm-Damkoff approximation is worked out in order to understand the observed behaviour of the decoherence. Only the continuous excited state phase transition produces a noticeable effect in the decoherence of the qubit. This is maximal when the system-environment coupling brings the environment to the critical point for the continuous phase transition. In this situation, the decoherence factor (or the fidelity) goes to zero with a finite size scaling power law.

preprint2010arXiv

Quantum quench influenced by an excited-state phase transition

We analyze excited-state quantum phase transitions (ESQPTs) in three schematic (integrable and nonintegrable) models describing a single-mode bosonic field coupled to a collection of atoms. It is shown that the presence of the ESQPT in these models affects the quantum relaxation processes following an abrupt quench in the control parameter. Clear cut evidence of the ESQPT effects is presented in integrable models, while in the nonintegrable model the evidence is blurred due to chaotic behavior of the system in the region around the critical energy.

preprint2010arXiv

Richardson-Gaudin Algebras and the Exact Solutions of the Proton-Neutron Pairing

Many exactly solvable models are based on Lie algebras. The pairing interaction is important in nuclear physics and its exact solution for identical particles in non-degenerate single-particle levels was first given by Richardson in 1963. His solution and its generalization to Richardson-Gaudin quasi-exactly solvable models have attracted the attention of many contemporary researchers and resulted in the exact solution of the isovector pn-pairing within the so(5) RG-model and the equal strength spin-isospin pn-pairing within the so(8) RG-model. Basic properties of the RG-models are summarized and possible applications to nuclear physics are emphasized. MSC2010 Classification: 81U15 Exactly and quasi-solvable systems, 17B81 Applications to physics, 81V35 Nuclear physics, 81R40 Symmetry breaking.

preprint2009arXiv

Local physics of magnetization plateaux in the Shastry-Sutherland model

We address the physical mechanism responsible for the emergence of magnetization plateaux in the Shastry-Sutherland model. By using a hierarchical mean-field approach we demonstrate that a plateau is stabilized in a certain {\it spin pattern}, satisfying {\it local} commensurability conditions derived from our formalism. Our results provide evidence in favor of a robust local physics nature of the plateaux states, and are in agreement with recent NMR experiments on $\scbo$.

preprint2009arXiv

Pairing in 4-component fermion systems: the bulk limit of SU(4)-symmetric Hamiltonians

Fermion systems with more than two components can exhibit pairing condensates of much more complex structure than the well-known single BCS condensate of spin-up and spin-down fermions. In the framework of the exactly solvable SO(8) Richardson-Gaudin model with SU(4)-symmetric Hamiltonians, we show that the BCS approximation remains valid in the thermodynamic limit of large systems for describing the ground state energy and the canonical and quasiparticle excitation gaps. Correlations beyond BCS pairing give rise to a spectrum of collective excitations, but these do not affect the bulk energy and quasiparticle gaps.

preprint2009arXiv

Phase diagram of the Heisenberg antiferromagnet with four-spin interactions

We study the quantum phase diagram of the Heisenberg planar antiferromagnet with a subset of four-spin ring exchange interactions, using the recently proposed heirarchical mean-field approach. By identifying relevant degrees of freedom, we are able to use a single variational anzatz to map the entire phase diagram of the model and uncover the nature of its various phases. It is shown that there exists a transition between a Néel state and a quantum paramagnetic phase, characterized by broken translational invariance. The non-magnetic phase preserves the lattice rotational symmetry, and has a correlated plaquette nature. Our results also suggest that this phase transition can be properly described within the Landau paradigm.

preprint2008arXiv

Unexpected features of quantum degeneracies in a pairing model with two integrable limits

The evolution pattern of level crossings and exceptional points is studied in a non-integrable pairing model with two different integrable limits. One of the integrable limits has two independent parameter-dependent integrals of motion. We demonstrate, and illustrate in our model, that quantum integrability of a system with more than one parameter-dependent integral of motion is always signaled by level crossings of a complex-extended Hamiltonian. We also find that integrability implies a reduced number of exceptional points. Both properties could uniquely characterize quantum integrability in small Hilbert spaces.

preprint2007arXiv

Exact solution of the spin-isospin proton-neutron pairing Hamiltonian

The exact solution of proton-neutron isoscalar-isovector (T=0,1) pairing Hamiltonian with non-degenerate single-particle orbits and equal pairing strengths (g_{T=1}= g_{T=0}) is presented for the first time. The Hamiltonian is a particular case of a family of integrable SO(8) Richardson-Gaudin (RG) models. The exact solution of the T=0,1 pairing Hamiltonian is reduced to a problem of 4 sets of coupled non linear equations that determine the spectral parameters of the complete set of eigenstates. The microscopic structure of individual eigenstates is analyzed in terms of evolution of the spectral parameters in the complex plane for system of A=80 nucleons. The spectroscopic trends of the exact solutions are discussed in terms of generalized rotations in isospace.

preprint2006arXiv

Exact Solution of the Isovector Proton Neutron Pairing Hamiltonian

The complete exact solution of the T=1 neutron-proton pairing Hamiltonian is presented in the context of the SO(5) Richardson-Gaudin model with non-degenerate single-particle levels and including isospin-symmetry breaking terms. The power of the method is illustrated with a numerical calculation for $^{64}$Ge for a $pf+g_{9/2}$ model space which is out of reach of modern shell-model codes.

preprint2006arXiv

Exactly-solvable models of proton and neutron interacting bosons

We describe a class of exactly-solvable models of interacting bosons based on the algebra SO(3,2). Each copy of the algebra represents a system of neutron and proton bosons in a given bosonic level interacting via a pairing interaction. The model that includes s and d bosons is a specific realization of the IBM2, restricted to the transition regime between vibrational and gamma-soft nuclei. By including additional copies of the algebra, we can generate proton-neutron boson models involving other boson degrees of freedom, while still maintaining exact solvability. In each of these models, we can study not only the states of maximal symmetry, but also those of mixed symmetry, albeit still in the vibrational to gamma-soft transition regime. Furthermore, in each of these models we can study some features of F-spin symmetry breaking. We report systematic calculations as a function of the pairing strength for models based on s, d, and g bosons and on s, d, and f bosons. The formalism of exactly-solvable models based on the SO(3,2) algebra is not limited to systems of proton and neutron bosons, however, but can also be applied to other scenarios that involve two species of interacting bosons.

preprint2006arXiv

Two-level interacting boson models beyond the mean field

The phase diagram of two-level boson Hamiltonians, including the Interacting Boson Model (IBM), is studied beyond the standard mean field approximation using the Holstein-Primakoff mapping. The limitations of the usual intrinsic state (mean field) formalism concerning finite-size effects are pointed out. The analytic results are compared to numerics obtained from exact diagonalizations. Excitation energies and occupation numbers are studied in different model space regions (Casten triangle for IBM) and especially at the critical points.

preprint2004arXiv

Exactly-solvable models for atom-molecule hamiltonians

We present a family of exactly-solvable generalizations of the Jaynes-Cummings model involving the interaction of an ensemble of SU(2) or SU(1,1) quasi-spins with a single boson field. They are obtained from the trigonometric Richardson-Gaudin models by replacing one of the SU(2) or SU(1,1) degrees of freedom by an ideal boson. Application to a system of bosonic atoms and molecules is reported.

preprint2003arXiv

Pair Fluctuations in Ultra-small Fermi Systems within Self-Consistent RPA at Finite Temperature

A self-consistent version of the Thermal Random Phase Approximation (TSCRPA) is developed within the Matsubara Green's Function (GF) formalism. The TSCRPA is applied to the many level pairing model. The normal phase of the system is considered. The TSCRPA results are compared with the exact ones calculated for the Grand Canonical Ensemble. Advantages of the TSCRPA over the Thermal Mean Field Approximation (TMFA) and the standard Thermal Random Phase Approximation (TRPA) are demonstrated. Results for correlation functions, excitation energies, single particle level densities, etc., as a function of temperature are presented.

preprint2002arXiv

Large N limit of the exactly solvable BCS model: analytics versus numerics

We have studied the numerical solutions of Richardson equations of the BCS model in the limit of large number of energy levels at half-filling, and compare them with the analytic results derived by Gaudin and Richardson, which in turn leads to the standard BCS solution. We focus on the location and density of the roots, the eigenvalues of the conserved quantities, and the scaling properties of the total energy for the equally spaced and the two-level models.

preprint2002arXiv

The Density Matrix Renormalization Group Method and Large-Scale Nuclear Shell-Model Calculations

The particle-hole Density Matrix Renormalization Group (p-h DMRG) method is discussed as a possible new approach to large-scale nuclear shell-model calculations. Following a general description of the method, we apply it to a class of problems involving many identical nucleons constrained to move in a single large j-shell and to interact via a pairing plus quadrupole interaction. A single-particle term that splits the shell into degenerate doublets is included so as to accommodate the physics of a Fermi surface in the problem. We apply the p-h DMRG method to this test problem for two $j$ values, one for which the shell model can be solved exactly and one for which the size of the hamiltonian is much too large for exact treatment. In the former case, the method is able to reproduce the exact results for the ground state energy, the energies of low-lying excited states, and other observables with extreme precision. In the latter case, the results exhibit rapid exponential convergence, suggesting the great promise of this new methodology even for more realistic nuclear systems. We also compare the results of the test calculation with those from Hartree-Fock-Bogolyubov approximation and address several other questions about the p-h DMRG method of relevance to its usefulness when treating more realistic nuclear systems.

preprint2001arXiv

A New Approach to Large-Scale Nuclear Structure Calculations

A new approach to large-scale nuclear structure calculations, based on the Density Matrix Renormalization Group (DMRG), is described. The method is tested in the context of a problem involving many identical nucleons constrained to move in a single large-j shell and interacting via a pairing plus quadrupole interaction. In cases in which exact diagonalization of the hamiltonian is possible, the method is able to reproduce the exact results for the ground state energy and the energies of low-lying excited states with extreme precision. Results are also presented for a model problem in which exact solution is not feasible.

preprint2001arXiv

New mechanism for the enhancement of $sd$ dominance in interacting boson models

We introduce an exactly solvable model for interacting bosons that extend up to high spin and interact through a repulsive pairing force. The model exhibits a phase transition to a state with almost complete $sd$ dominance. The repulsive pairing interaction that underlies the model has a natural microscopic origin in the Pauli exclusion principle between contituent nucleons. As such, repulsive pairing between bosons seems to provide a new mechanism for the enhancement of $sd$ dominance, giving further support for the validity of the $sd$ Interacting Boson Model.