Conference on Quantum-Many-Body Correlations in memory of Peter Schuck (QMBC 2023)

21 mars 2023, 08:30 → 23 mars 2023, 18:30 Europe/Paris

IJCLab, Orsay, France

This conference is dedicated to the memory of Peter Schuck, who was an outstanding personality not only in the field of nuclear theory, but in many-body physics in general. Thanks to his broad interests, he was able to initiate discussions between different communities with the aim to transfer methods from one field to another and to overcome barriers due to different terminologies. Following these ideas, the objective of this conference is to bring together experts from different fields of many-body physics and to discuss recent progress, methods and results covering the large spectrum of subjects in which Peter Schuck was interested, from nuclei and neutron stars to condensed-matter systems and theoretical chemistry. Examples of topics that will be discussed are:

• clustering phenomena in nuclei,
• quartet condensation in nuclear and condensed-matter systems,
• equation of state and density-functional theory,
• BCS-BEC crossover in ultracold Fermi gases,
• pairing and superfluidity in nuclear systems and neutron stars,
• RPA and its extensions in nuclear physics, solid-state physics, and theoretical chemistry,
• semiclassical methods for many-body systems.

Confirmed invited speakers: M. Baldo (Catania), B. Borderie (IJCLab Orsay), L. Canet (Grenoble), G. Chanfray (IP2I Lyon), D. Delion (Bucharest), Y. Funaki (Yokohama), F. Gulminelli (LPC Caen), E. Litvinova (Western Michigan Univ.), U. Lombardo (Catania), A. Minguzzi (Grenoble), V. Olevano (Grenoble), P. Pieri (Bologna), V. Piselli (Camerino), C. Providência (Coimbra), A. Recati (Trento), P. Ring (Munich), C. Salomon (Paris), J. Toulouse (Paris), X. Viñas (Barcelona), B. Zhou (Shanghai)

Local Organizing committee: Guillaume Hupin, Denis Lacroix, Michael Urban (chair)

Scientific committee: Markus Holzmann (Grenoble), Elias Khan (Orsay), Gerd Röpke (Rostock), Giancarlo Strinati (Camerino)

Conference secretary: Émilie Bonnardel, Valérie Brouillard

From_Orsay_and_Bures_IJCLab_com_2022.pdf

IJCLab-map-zoom_IJCLab-A4-RectoVerso-0822.pdf

mardi 21 mars

Tuesday 09:00-10:45

Président de session: Michael Urban (IPN Orsay)

1 Welcome

Orateur: Michael Urban (IJCLab)

2 Opening address

Orateur: Marcella Grasso (IN2P3)

Grasso-QMBC.pdf

3 Relativistic Brueckner-Hartree-Fock Theory: an ab initio Approach for Finite Nuclei

Most investigations of collective vibrations in medium-heavy and heavy nuclei are based on nuclear density functional theory. Very successful relativistic and non-relativistic functionals are available nowadays. However, most of them are phenomenological functionals. Therefore it is important to study the connection of such functionals to ab-initio nucleon-nucleon forces. Non-relativistic Brueckner-Hartree-Fock theory was a starting point of ab-initio investigations in nuclear structure in the fifties and sixties. However, it failed because three-body forces were not included at that time. Later it was found that relativistic Brueckner-Hartree-Fock (RBHF) theory can reproduce the saturation properties of nuclear matter. In this contribution, we discuss recent developments of RBHF theory for infinite nuclear matter finite nuclei, particularly applications for neutron stars.

Orateur: Peter Ring (Physics Department Technical University Munich)

230321-Ring-Orsay.ppt

4 Studies of odd-A nuclei and high-K isomeric states with the BCPM functional

The BCPM (Barcelona-Catania-Paris-Madrid) is a nuclear energy density functional with a central term obtained by using the density of finite nuclei in a polynomial fit to nuclear matter equation of state (both symmetric and neutron matter). As formulated, the BCPM functional does not contain time-odd densities and therefore it should not be used to describe odd-A nuclei and/or multiquasiparticle excitations like the ones making high-K isomeric states. However, the equal filling approximation to full blocking has proved to be very accurate in describing the above mentioned systems, with the advantage of not requiring time odd fields. We have implemented the EFA along with the BCPM functional to obtain properties of odd-A nuclei like ground state spin and parity and the excitation energy spectrum. The idea can easily be generalized to describe high-K isomeric states.

Orateur: Luis Robledo (Universidad Autonoma de Madrid)

Talk_1.pdf

5 Relativistic equation of motion framework for nuclear physics

Orateur: Dr Elena Litvinova (Western Michigan University)

Litvinova-QMBC.pdf

10:45 Coffee break

Tuesday 11:15-12:45

Président de session: Elias Khan (IPN Orsay)

6 Superfluidity in nuclear systems

Peter Schuck devoted much attention to the study of nuclear superfluidity and to many-body effects that renormalize the pairing interaction in nuclear systems, going from his works with the Catania group to his recent papers with E. Litvinova.
I will resume the research done on this topic by the Milano group, starting from a paper we wrote together with him [1]. Such research deals not only with atomic nuclei but also with superfluidity in the inner crust of neutron stars and with the structure of vortices.

[1] F. Barranco, P.F.Bortignon, R.A. Broglia, G. Colò, P. Schuck, E. Vigezzi and X. Vinas, Pairing matrix elements and pairing gaps with bare, effective and induced interactions, Phys. Rev. C 73 (2005) 054314

Orateur: Dr Enrico Vigezzi (INFN Milano)

vigezzi_schuck_2023_a.pdf

7 Anti-Screening Effects in the Nuclear Pairing

Orateur: Dr Umberto Lombardo (INFN Catania)

Lombardo.pdf

8 Suppression of superfluidity in neutron-star crusts

Formed in the aftermath of gravitational core-collapse supernova explosions, neutron stars contain matter crushed at densities exceeding that found inside the heaviest atomic nuclei under conditions so extreme that they cannot be reproduced on Earth. Neutron stars are therefore unique laboratories for exploring phases of matter not observed in any other celestial bodies. As early as 1959, the inner crust and the outer core of a neutron star were predicted by Arkhady Migdal to exhibit neutron superfluidity at temperatures below $\sim 10^{10}$ K. Since then, neutron superfluidity has found strong support from radio-timing observations of pulsar frequency glitches, and more recently from the rapid decline of luminosity of the youngest known neutron star in the supernova remnant of Cassiopeia A [1].

In a similar way to laboratory superfluidity of atomic gases in optical lattices, neutron superfluidity in the inner crust of neutron stars is partially suppressed due to the presence of neutron-proton clusters [2]. The reduction of the neutron superfluid density implies that the neutron superfluid cannot flow completely freely despite the absence of viscous drag. Due to this entrainment effect, clusters move with an effective mass and collective excitations of the superfluid are mixed with those of the crust. The extent to which neutron superfluidity is suppressed depends on whether the crust is crystalline or disordered [3]. Recent fully 3D neutron band-structure calculations taking BCS pairing into account will be presented. Astrophysical implications will be discussed.

[1] N. Chamel, J. Astrophys. Astr. 38, 43 (2017).
[2] N. Chamel, J. Low Temp. Phys. 189, 328 (2017).
[3] J. A. Sauls, N. Chamel, M. A. Alpar, eprint arXiv:2001.09959

Orateur: Nicolas Chamel

Chamel_QMBC.pdf

13:00 Lunch

Tuesday 14:00-15:30

Président de session: Pierbiagio Pieri (University of Bologna)

9 Probing the BCS-BEC crossover with persistent currents

Orateur: Dr Anna Minguzzi (LPMMC, Grenoble)

Minguzzi-QMBC.pdf

10 Vortices in ultracold Fermi gases: peculiarity of their structure and impact on dynamics

It will be shown that spin polarized vortices in Fermi superfluid acquire a peculiar structure with a reversed circulation inside the core. Their structure admits the vanishing minigap with a characteristic pattern of single-quasiparticle level crossings at the Fermi surface. It is also predicted that the dynamics along the vortex line of spatially localized polarization inside the core will be suppressed. The impact of the vortex core structure on dynamics will be analyzed in the context of recent experiment involving colliding vortex dipoles [Nature (London) 600, 64 (2021)] which revealed nonuniversal dissipative dynamics. Moreover consequences for neutron star crust and the decay of turbulent state will be discussed.
References
[1] Piotr Magierski, Gabriel Wlazłowski, Andrzej Makowski, Konrad Kobuszewski, Phys. Rev. A 106, 033322 (2022)
[2] Andrea Barresi, Antoine Boulet, Piotr Magierski, Gabriel Wlazłowski, Phys. Rev. Lett. 130, 043001 (2023)
[3] Khalid Hossain, Konrad Kobuszewski, Michael McNeil Forbes, Piotr Magierski, Kazuyuki Sekizawa, Gabriel Wlazłowski, Phys. Rev. A 105, 013304 (2022)
[4] Gabriel Wlazłowski, Klejdja Xhani, Marek Tylutki, Nikolaos P. Proukakis, Piotr Magierski, Phys. Rev. Lett. 130, 023003 (2023)
[5] Daniel Pęcak, Nicolas Chamel, Piotr Magierski, Gabriel Wlazłowski, Phys. Rev. C 104, 055801 (2021)

Orateur: Prof. Piotr Magierski (Warsaw University of Technology)

qmbc2023.pdf

11 Space-time correlations in turbulence

The problem of solving a strongly-coupled many-body system arises in classical fluid dynamics when turbulence fully develops. The statistical properties of the turbulent fluid are encoded in an infinite hierachy of coupled equations which leads to what is called the "closure problem" of turbulence. The functional renormalisation group offers an efficient method to tackle this problem and achieve a controlled closure in the limit of large wavenumbers. I will present the principles of this approach and illustrate the results through comparisons with data from both direct numerical simulations and experiments.

Ref: L. Canet, Journal of Fluid Mechanics, Perspectives, 950, 1 (2022)

Orateur: Léonie CANET (Université Grenoble Alpes)

Canet-QMBC.pdf

15:30 Coffee break

Tuesday 16:00-17:30

Président de session: guy chanfray (IPN Lyon)

12 Neutron star properties from semiclassical methods in nuclear physics

Semiclassical methods in nuclear physics are a very good tool to study the behavior and composition of neutron stars. Using the Vlasov equation, the neutron star crust-core transition will be discussed, also in the presence of a strong magnetic field as expected inside magnetars. The description of the inner crust within a Thomas Fermi approach will be presented.

Orateur: Constança Providência (University of Coimbra)

PSchuck_QMBC2023_CProvidencia.pdf

13 An attempt at semiclassical methods for electronic-structure theory

Orateur: Julien Toulouse (Laboratoire de Chimie Théorique, Sorbonne Université)

Toulouse-QMBC.pdf

14 In memory of Peter Schuck: "Semiclassical approximation to pairing in the weak coupling regime: nuclei, cold atoms, and neutron stars"

A novel Thomas-Fermi (TF) approach to inhomogeneous superfluid Fermi-systems is presented and shown that it works well also in cases where the Local Density Approximation (LDA) breaks down. The novelty lies in the fact that the semiclassical approximation is applied to the pairing matrix elements not implying a local version of the chemical potential as with LDA. Applications will be given to the generic fact that if a fermionic superfluid in the BCS regime overflows from a narrow container into a much wider one, pairing is substantially reduced at the overflow point. Two examples pertinent to the physics of the outer crust of neutron stars and superfluid fermionic atoms in traps will be presented. The TF results will be compared to quantal and LDA ones.
[1] P. Schuck and X. Viñas, Phys. Rev. Lett. 107, 205301 (2011).
[2] X. Viñas, P. Schuck and M. Farine, J. Phys. Conf. Ser. 321, 012024 (2011); J. Phys. Conf. Ser. 338, 012016 (2012).
[3] Fifty years of Nuclear BCS: Pairing and Finite Systems, 212, Ed. R. Broglia, World Scientific, 2013.

Orateur: Prof. Xavier Viñas (Universitat de Barcelona and Institut Menorquı́ d’Estudis, Maó)

Vinas-QMBC.pdf

Cocktail & Poster session

15 Core level chemical shift by ab initio mehods: from mean-field to many-body theory

The binding energy of core electrons may not only provide information on the chemical composition, but also some additional information such as the type of bonding, which could be inferred from the shift of the binding energy, (also known as chemical shift). We present the study of the chemical shift using different theories, from Hartree-Fock and density-functional theory to many-body perturbation theory (COHSEX, GW). We benchmarked the accuracy of the chemical shift of the carbon 1s electron in a set of molecules against experiments. Besides, our study reveals the physical origin of the chemical shift.

Orateur: Valerio Olevano (CNRS, Institut NEEL)

16 Disordered structures in ultracold spin-imbalanced Fermi gas

We investigate the properties of spin-imbalanced ultracold Fermi gas in a large range of spin polarizations at low temperatures. We present the results of microscopic calculations based on mean-field and density functional theory approaches, with no symmetry constraints. At low polarization values we predict the structure of the system as consisting of several spin-polarized droplets. As the polarization increases, the system self-organizes into a disordered structure similar to liquid crystals, and energetically they can compete with ordered structures such as grid-like domain walls. At higher polarizations the system starts to develop regularities that, in principle, can be called supersolid, where periodic density modulation and pairing correlations coexist. The robustness of the results has been checked with respect to temperature effects, dimensionality, and the presence of a trapping potential. Dynamical stability has also been investigated.

Orateur: Piotr Magierski (Warsaw University of Technology)

17 Electron removal energies in noble gas atoms up to 100 keV: ab initio GW vs XPS

X-ray photoelectron spectroscopy (XPS) measures electron removal (quasiparticle) energies, providing direct access to core and valence electron binding energies, hence probing the electronic structure. We present the benchmark of the ab initio many-body GW approximation on the complete electron binding energies of noble gas atoms (He-Rn), which spans 100 keV. Our results demonstrate that GW achieves an accuracy within 1.2% in XPS binding energies, by systematically restoring the underestimation from density-functional theory (DFT, error of 14%) or the overestimation from Hartree-Fock (HF, error of 4.7%). Such results also imply the correlations of d electrons are very well described by GW.

Orateur: Valerio Olevano (CNRS, Institut NEEL)

18 Equation of state of superfluid neutron matter with low-momentum interactions

In this work, we calculate the ground state energy of pure neutron matter using the renormalization group based low-momentum effective interaction $V_{\text{low-}k}$ in Bogoliubov many-body perturbation theory (BMBPT), which is a perturbative expansion around the Hartree-Fock-Bogoliubov (HFB) ground state. In order to capture the low-density behavior of neutron matter, it turns out to be better to use a density dependent cutoff in the $V_{\text{low-}k}$ interaction. Perturbative corrections to the HFB energy up to third order are included. We find that at low densities corresponding to the inner crust of neutron stars, the HFB state that includes pairing is a better starting point for perturbation expansion. It is observed that including the higher order perturbative corrections, the cutoff dependence of the ground state energy is reduced.

Reference: V. Palaniappan, S. Ramanan, M. Urban, Phys. Rev. C 107, 025804 (2023)

Orateur: Viswanathan Palaniappan (IJCLab)

19 Neutrino Scattering Rates of Neutron-star and Supernova Matter within Skyrme RPA

Supernova explosions, which will leave behind a neutron star, are the most powerful neutrino sources. The neutrino emission is also the dominating cooling mechanism for a (proto-)neutron star, whose interior is mainly composed of extremely dense and hot nuclear matter. Neutrinos can be scattered frequently inside stars before they escape. We study the neutrino scattering rates of neutron-star and supernova matter within Skyrme RPA response functions. The neutrino scattering rates in neutron-star matter depend sensitively on the adopted interaction. The minimum scattering angle is different for different interactions because it depends on the Fermi velocity. We also find that many Skyrme interactions present the unphysical feature that the Fermi velocity of neutrons in neutron-star matter exceeds the speed of light at a density below the maximum central density of the neutron star predicted by the Skyrme interactions.

Orateur: Mingya Duan (IJCLab)

20 Neutron stars within the Bogoliubov quark-meson coupling model

A quark-meson coupling model based on the quark model proposed by Bogoliubov for the description of the quark dynamics is developed and applied to the description of neutron stars. Starting from a SU(3) symmetry approach, it is shown that this symmetry has to be broken in order to satisfy the constraints set by the hypernuclei and by neutron stars. The model is able to describe observations such as two-solar-mass stars or the radius of canonical neutron stars within the uncertainties presently accepted. If the optical potentials for Λ and Ξ hyperons in symmetric nuclear matter at saturation obtained from laboratory measurements of hypernuclei properties are imposed, the model predicts no strangeness inside neutron stars.

Orateur: Aziz Rabhi (Université de Carthage)

21 Spin-polarized droplets in ultracold Fermi gas

We demonstrate the existence of a new type of spatially localized excitations in the unitary Fermi gas: spin polarized droplets with a peculiar internal structure involving the abrupt change of the pairing phase at the surface of the droplet. It resembles the structure of the Josephson-π junction occurring when a slice of a ferromagnet is sandwiched between two superconductors. The stability of the impurity is enhanced by the mutual interplay between the polarization effects and the pairing field, resulting in an exceptionally long-lived state. We show that the motion of spin-polarized impurity (ferron) in ultracold atomic gas is characterized by a certain critical velocity which can be traced back to the amount of spin imbalance inside the impurity. We have calculated the effective mass of ferron in two dimensions. We show that the effective mass scales with the surface of the ferron. We discuss the impact of these findings; in particular, we demonstrate that ferrons become unstable in the vicinity of a vortex.

Orateur: Piotr Magierski (Warsaw University of Technology)

22 Static self energy and effecive mass of the homogeneous electron gas from Quantum Monte Carlo calculations

Landau’s Fermi liquid theory has provided a paradigmatic frame for the phenomenological description of equilibrium and transport properties of degenerate fermions in terms of very few characteristic parameters. Silin has set out the path to generalize for long-range forces, such to extend it for normal metals in condensed matter. Although the formal structure of the underlying microscopic theory is known for a long time, most explicit calculations of the Fermi liquid parameters basically rely on approximative, perturbative schemes. As diagrammatic perturbation theory is not expected to converge for typical electronic densities, basic Fermi liquid parameters of the homogeneous electron gas (jellium), like the effective mass $m^∗$ and the renormalization factor Z, considerably vary according to the underlying approximation. Recently, diagrammatic Monte Carlo calculations (DiaQMC) [1] have been performed to include and control higher order terms of the perturbation series. Those calculations found an overall agreement for Z with previous quantum Monte Carlo calculations (QMC) [2]. However, DiaQMC results on $m^∗$ have been strongly questioned by QMC calculations yielding substantial different values [3]. Here, we rewise the methodology of zero temperature QMC calculations of the effective mass, in order to resolve the discrepancy between QMC and perturbative/DiaQMC results.

[1] K. Haule and K. Chen, Scientific Reports 12, 2294 (2022).
[2] M. Holzmann, B. Bernu, C. Pierleoni, J. McMinis, D. M. Ceperley, V. Olevano, and L. Delle Site, Phys. Rev. Lett.
107, 110402 (2011).
[3] S. Azadi, N.D. Drummond, and W.M.C. Foulkes, Phys. Rev. Lett. 127, 086401 (2021).

Orateur: Dr Markus Holzmann (LPMMC, Grenoble)

23 Thermal phase transitions and giant resonances in atomic nuclei with chiral effective interactions

Ab initio approaches to the nuclear many-body problem have seen their reach considerably extended over the last decade. However, collective excitations have seldom been addressed in that context, due to the prohibitive cost of solving the corresponding equations of motion.

We adapt a novel method originally proposed in the framework of the nuclear energy density functional method, and derive its extension to non-zero temperatures.
We thus obtain the first fully self-consistent ab initio study of multipolar responses of atomic nuclei using two- and three-body chiral interactions, at zero and finite temperature.

The method is applied to the mid-mass nucleus $^{56}$Fe, which is challenging for ab initio interactions due to its mid-mass, deformed and superfluid character. After reproducing the temperature-induced shape transition at the mean-field level, we analyse the evolution of the isovector dipole strength function within a temperature range $k_BT \in [0;4]$ MeV. A strong dependence is observed, characterised by a downwards shift of the resonance centroid when the system is heated up.

[1] T. Nakatsukasa, T. Inakura, and K. Yabana. Finite amplitude method for the solution of the random-phase approximation. Phys. Rev. C, 76:024318, Aug 2007. doi: 10.1103/PhysRevC.76.024318.
[2] Y. Beaujeault-Taudière. Study of the multipolar excitations in cold and hot, deformed and superfluid systems with the method of finite amplitudes. Theses, Université Paris-Saclay, Oct 2021.
[3] Y. Beaujeault-Taudière et al. Ab initio description of multipolar responses in superfluid and deformed nuclei at finite temperature: application to dipole modes in 56Fe. Mar 2022. doi: 10.48550/arXiv.
2203.13513.
[4] A. C. Larsen, N. Blasi, A. Bracco, F. Camera, T. K. Eriksen, A. Görgen, M. Guttormsen, T. W. Hagen, S. Leoni, B. Million, H. T. Nyhus, T. Renstrøm, S. J. Rose, I. E. Ruud, S. Siem, T. Tornyi, G. M. Tveten,
A. V. Voinov, and M. Wiedeking. Evidence for the dipole nature of the low-energy 𝛾 enhancement in 56Fe. Phys. Rev. Lett., 111:242504, Dec 2013. doi: 10.1103/PhysRevLett.111.242504.

Orateur: Yann Beaujeault-Taudiere (IJCLab & LLR)

24 Simulation of the time evolution for one-dimensional wave-functions with quantum computation

Quantum computation is a relatively new field that seeks to harness the power of quantum mechanics to perform calculations that would be impossible with classical computers. The Schrödinger equation lies at the heart of quantum physics. In this poster, we present how to address the one-dimensional time evolution of wave-functions, as governed by the Schrödinger equation, on quantum computation devices. The potential of quantum computation to achieve reliable simulation of the process is demonstrated. In the end, we highlight the challenges of implementing non-unitary boundary conditions.

Orateur: M. Jing Zhang (IJCLab)

mercredi 22 mars

Wednesday 09:15-10:45

Président de session: Gerd Röpke (University of Rostock)

25 Causes and consequences of the emergence of clusters in nuclei

Nuclear clustering refers to the appearance of multi-nucleon localized structures, predominately under the form of α particles, within the interior of a nucleus. The detailed understanding of how quantum correlations among nucleons give rise to nuclear clustering is still lacking, the main difficulty being the need to a priori include four-nucleon correlations.
Using the Energy Density Functional framework, conditions for the emergence of clusters are analyzed and spectroscopic signatures of these structures are discussed.

Orateur: Dr Jean-Paul EBRAN (CEA)

Ebran-QMBC2023.pdf

26 Study of alpha condensate and beyond

After the proposal of the so-called THSR wave function in 2001, it has been developed for 20 years and applied to many nuclear systems. We discuss the achievement in the study of the alpha condensation physics and related issues with the THSR ansatz, which has been done in collaboration with Prof. Peter Schuck.

Orateur: Yasuro Funaki (Kanto Gakuin University)

27 Temperature and Density Conditions for Alpha Clustering in Alpha-Conjugate Nuclei

Orateur: Bernard Borderie (IJCLab)

Borderie_QMBC 2023_indico.pdf

10:45 Coffee break

Wednesday 11:15-12:45

Président de session: Piotr Magierski (Warsaw University of Technology)

28 Pairing in systems with imbalance and BCS-BEC crossover

I will discuss the work initiated in collaboration with Peter Schuck and others on the imbalance of superfluid on the isospin asymmetric nuclear matter. The current understanding of the phase diagram of dilute nuclear matter which includes the phases with broken space symmetries will be discussed covering the temperature-density regime across the BCS-BEC crossover.

Orateur: M. Armen Sedrakian (Frankfurt Institute for Advanced Studies)

Sedrakian-Schuck2023.pdf

29 FFLO correlations in polarized ultracold Fermi gases

Quite generally, an imbalance between the densities of spin-up and spin-down fermions hinders pairing and superfluidity in two-component attractive Fermi gases. The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, in which pairs condense at a finite value of center-of-mass momentum to compensate for the mismatch of the two Fermi surfaces, was proposed many years ago as a possible superfluid phase compatible with a finite polarization. In this talk, I will discuss how significant precursor FFLO fluctuation effects appear already in the normal phase of polarized Fermi gases at finite temperature [1], and how they could be observed experimentally. At zero temperature [2], I will discuss how the quasi-particle parameters of the normal Fermi gas are changed when approaching an FFLO quantum critical point. Within a fully self-consistent t-matrix approach we find that the quasi-particle residues vanish, and the effective masses diverge at the FFLO quantum critical point, with a complete breakdown of the quasi-particle picture that is similar to what is found in heavy-fermion materials at an antiferromagnetic quantum critical point.

References
[1] M. Pini, P. Pieri, and G. Calvanese Strinati, Phys. Rev. Res. 3, 043068 (2021).
[2] M. Pini, P. Pieri, and G. Calvanese Strinati, arXiv:2211.15529.

Orateur: Pierbiagio Pieri (University of Bologna)

Orsay-2023-Pieri.pdf

30 Ferromagnetism in coherently coupled atomic mixtures

In the present talk we review the properties of a spin-1/2 Bose-Einstein condensate (BEC) in presence of both longitudinal and transverse external fields. The system has very peculiar properties not only with respect to the single component BEC, but also with respect to BEC mixtures (where the atom number of each species is conserved). In particular the system can exhibit a para- to ferro-magnetic transition, has a gapped spin spectrum, and can sustain peculiar magnetic defects. Most of the peculiar properties are due to the system being an analog of a (non-dissipative) ferromagnetic material well described by the so-called Landau-Lifshitz equation. We then present some of the system's properties that have very recently been tested experimentally in our lab using $^{23}$Na atomic gases. In particular we measured the excitation spectrum by modulating the trapping potential and generating Faraday patterns, we started characterising the ferromagnetic transition and obtain the first results on magnetic bubble generation from out-of-equilibrium initial states.

Orateur: Dr Alessio Recati (Pitaevskii BEC Center, CNR-INO and Dipartimento di Fisica, Trento, Italy)

Recati.pdf

13:00 Lunch

Wednesday 14:00-15:30

Président de session: Markus Holzmann (LPMMC, Grenoble)

31 Cooperation with Peter Schuck on many-body correlations in nuclei

The first topic describes the Selfconsitent Random Phase Approximation (SCRPA) [1]. The nonlinear SCRPA system of equations is numerically solved for the three level Lipkin model [2]. Goldstone mode and mass parameter in the deformed region are analyzed.
The second topic analyzes simultaneous description of alpha and electromagnetic transitions in $^{212}$Po in terms of the surface alpha clustering [3]. Large dipole electromagnetic transitions from recently discovered unnatural parity states are explained.
[1] P. Schuck, D.S. Delion, J. Dukelsky, M. Jemai, E. Litvinova, G. Ropke, and M. Tohyama, Equation of Motion Method for strongly correlated Fermi systems and Extended RPA approaches, Physics Reports 929 (2021) 1.
[2] D.S. Delion, P.Schuck, and J. Dukelsky, Self Consistent Random Phase Approximation and the restoration of symmetries within the three level Lipkin model, Physical Review C 72 (2005) 064305.
[3] D.S. Delion, R.J. Liotta, P. Schuck. A. Astier, and M.-G. Porquet, Shell model plus cluster description of negative parity states in $^{212}$Po, Physical Review C 85 (2012) 064306.

Orateur: Doru Delion (Horia Hulubei National Institute of Physics and Nuclear Engineering)

Delion-Orsay.pdf

32 dRPA, RPAx, SRPA, SCRPA, rRPA, ...: tell me who you are and I will tell your RPA

The many-body problem, i.e. the solution of the Schrödinger equation in interacting many-body systems, is a formidable problem in condensed matter and nuclear physics, as well as in quantum chemistry, so far insolute despite considerable efforts of a whole generation of physicists and chemists. Every community has attacked the problem from a different point of view; has along the years elaborated varied theories, from wave-function to density-functional, up to Green function based approaches; and has achieved sometimes common findings, but also, and more often, different insights and complementary comprehension. Merging these insights and comprehension was the intuition of Peter Schuck's "RPA multidisciplinary conference" series, which was very stimulating and led to fruitful work in collaboration. Taking inspiration from this work, I will try to present here the largest landscape available, within my comprehension limits, about RPA flavors and many-body approaches encompassing quantum chemistry, condensed-matter and nuclear physics, speaking in a common, as much as possible, Esperanto language.

Orateur: Valerio Olevano (CNRS, Institut NEEL)

olevano.pdf

33 Neutrino-nucleon interactions in dense and hot matter

Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and have among others a strong impact on conditions for heavy element nucleosynthesis and neutron star cooling is dominated by neutrino emission except for very old stars. Many works in the last decades have shown that in dense matter medium effects considerably change the neutrino-matter interaction rates, whereas many astrophysical simulations use analytic approximations which are often far from reproducing more complete calculations. In this talk I will present a scheme which allows to incorporate improved rates into simulations and show as an example some results for core-collapse supernovae and proto-neutron star cooling.

Orateur: Micaela Oertel (LUTH, Observatoire de Paris)

oertel.pdf

15:30 Coffee break

Wednesday 16:00-18:00

Président de session: Jean-Paul Blaizot

34 Nuclear pastas in neutron stars: uncertainties of ETFSI approach

The region between the crust and the core of a neutron star (NS) may consist of a mantle of so-called nuclear pastas. If they exist, these exotic nuclear structures could significantly affect the transport and mechanical properties of dense matter, leaving their imprints on such NS observables as continuous gravitational-wave emission, NS oscillations and their damping, the spin period of x-ray pulsars, and NS cooling.

Such an extremely dense, neutron-rich and inhomogeneous environment represents unique conditions inaccessible to laboratory experiments and challenging for a fully microscopic treatment. Although nuclear pastas have been extensively studied within (semi-) classical methods, limited self-consistent mean-field calculations have been carried out so far due to high computational cost. For this reason, we follow an alternative approach adding proton shell corrections perturbatively via the Strutinsky integral theorem on top of the energy calculated within the 4th-order Extended Thomas-Fermi method [1]. We evaluate the uncertainties of this ETFSI approach using the generalized Skyrme effective interaction BSk24 [2].

First, we show within the ETF approach that the range of densities for which pasta phases are present and the types of pastas can depend on the parameterization of the nucleon profiles employed to speed up the calculations. To improve them, we introduce two new parametrizations for which we have found lower ETF energy solutions, therefore, more accurate nucleon profiles than previously obtained. In the second stage, when adding the SI corrections, we find that the differences in energies corresponding to the adopted parametrizations are (partially) compensated, such that the ETFSI results for the two new profiles are in good agreement. In contrast to the purely ETF calculations, the spaghetti phase is now replaced by spherical clusters, which can also intersperse the lasagna phase and thus shrink the NS mantle. However, results at high densities become sensitive to the imposed boundary conditions to calculate the SI correction. In summary, quantum effects are shown to play an important role for pasta phases, however, fully self-consistent mean-field calculations are still required to identify the configurations that are present near the crust-core boundary.

[1] J. M. Pearson and N. Chamel, Phys. Rev. C 105, 015803 (2022)
[2] S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C 88, 024308 (2013)

Orateur: Nikolai Shchechilin (Universite Libre de Bruxelles)

35 Cluster structures in light nuclei

Orateur: BO ZHOU (Fudan University)

36 Thermodynamics of quark matter with multi-quark clusters in an effective Beth-Uhlenbeck type approach*

We describe multi-quark clusters in quark matter within a generalized Beth-Uhlenbeck approach in a background gluon field that is coupled to the underlying chiral quark dynamics using the Polyakov-gauge and an effective potential for the traced Polyakov-loop. A higher multi-quark cluster of size n is described as a binary composite of smaller subclusters n1 and n2 (n1+n2=n) with a bound state and scattering state spectrum. For the corresponding cluster-cluster phase shifts we use two simple ansätze that capture the Mott dissociation of clusters as a function of temperature and chemical potential. We compare the simple "step-up-step-down" model that ignores continuum correlations with an improved model contains them in a generic form. In order to explain the model, we restrict ourselves here to the cases where n= 1, 2, …, 6.
A striking result is the suppression of the abundance of colored multi-quark clusters at low temperatures by the coupling to the Polyakov loop. This is understood in close analogy to the suppression of quark distributions by the same mechanism and we derive here the corresponding Polyakov-loop generalized distribution functions of n-quark clusters. With the input of the temperature (T) and chemical potential (μ) dependence of the chiral condensate from lattice QCD, we construct the QCD thermodynamics in good agreement with the data. Special emphasis is on the density and entropy density and their ratio in the T- μ plane that show the effect of multi-quark cluster formation and dissociation.

*) This work was initiated by Peter Schuck at the ECT Trento Workshop on “Light Clusters in Nuclei and Nuclear Matter: Nuclear Structure and Decay, Heavy Ion Collisions, and Astrophysics“ (2019)

Orateur: David Blaschke (University of Wroclaw)

Blaschke_QuarkCluster_2023_Orsay.pdf

37 The neutron star crust in the liquid phase

The crust and mantle of neutron stars are exceptional laboratories for the theoretical study of clustering in nuclear matter, and their specific transport properties are known to play an important role in many different observable phenomena, such as the thermal emission of X-ray pulsars, quasi-periodic oscillations, giant flares, pulsar glitches, just to cite a few.
Thanks to density functional theory, the composition and properties of crustal matter is relatively well known at zero temperature. The finite temperature problem has been much less addressed in the literature ; still, neutron stars are born hot, and the crustal ions as well as the non-spherical « pasta » nuclear clusters are in the liquid phase at all temperatures beyond about 10^10 K. In such a configuration, the ions are put into collective motion, the simplifying one-component Wigner-Seitz approximation is not valid, and the collective entropy contributions can strongly affect the thermodynamic and transport properties, as well as the composition of matter. This finite temperature modelling might also have some relevance in the late cooling stage of the neutron stars. Indeed, if the cooling is sufficiently fast, the crust composition might be frozen before the catalyzed configuration is reached. Notably, the cooling history of a neutron star is imprinted in the presence of impurities in the crustal lattice, that are at the origin of the highly resistive behavior of the crustal mantle.

In this contribution, I will present our latest results concerning the thermodynamic and transport properties of the neutron star crust in the liquid phase [1-5]. Concerning static properties, I will discuss the effect of going beyond the one-component plasma approximation and including a renormalization of the ion mass calculated in the hydrodynamic limit. Concerning transport, I will present analytical expressions for the anisotropic collision frequencies in the pasta phase using the Boltzmann equation in the relaxation time approximation, and numerical results for the associated electrical and thermal conductivities.

1. S.Mallik and F.Gulminelli, Statistical treatment of nuclear clusters in the continuum, PRC 103 (2021) 015803
2. T.Carreau, A.F.Fantina and F.Gulminelli, Inner crust of a neutron star at the point of crystallization in a multicomponent approach, A&A 640, A77 (2020).
3. M.R.Pelicer, D.P.Menezes, C.C.Barros Jr., F.Gulminelli, Fluctuations in the nuclear pasta phase, PRC 104 (2021) L022801
4. H.Dinh Thi, A.F.Fantina, F.Gulminelli, The proto-neutron star inner crust in the liquid phase, A&A, submitted
5. M.R.Pelicer, M.Antonelli, D.P.Menezes, F.Gulminelli, Anisotropic electron transport in the nuclear pasta phase, MNRAS, submitted

Orateur: Prof. Francesca Gulminelli (LPC/Ensicaen)

Gulminelli_QMBC.pdf

20:00 Conference Dinner

jeudi 23 mars

Thursday 09:15-10:45

Président de session: Giancarlo Calvanese Strinati (University of Camerino)

38 Josephson effect at finite temperature throughout the BCS-BEC crossover with the inclusion of pairing fluctuations

The BCS-BEC crossover is a research topic common to ultra-cold Fermi gases and the nuclear matter [1]. It represents a useful tool for testing fundamental theories such as the connection between superconductivity and fermionic superfluidity. At zero temperature, a reasonable description of this crossover can be obtained by a mean-field approach. In fact, this approximation provides accurate enough results even at finite temperature provided that the Cooper pair size is much larger than the average inter-particle distance (weak inter-particle interaction). However, a mean-field approach fails when the Cooper pair size is comparable or even smaller than the average-interparticle distance (strong-interparticle interaction). In this case it is necessary to include pairing fluctuations beyond mean-field to obtain reliable results [2]. In the present work, we include pairing fluctuations over and above mean-field in the LPDA (Local Phase Density Approximation) equation, which is a coarse-grained version of the Bogoliubov-de Gennes equations. The LPDA equation is a (highly) non-linear differential equation for the gap parameter, which allows one to study inhomogeneous superfluid systems with a considerably reduced computational time and memory space with respect to the original Bogoliubov-de Gennes equations [3]. In particular, in the present work we address the Josephson effect using a modified (mLPDA) version of the LPDA equation which includes pairing fluctuations on top of the original equation. We show that the outcomes of our numerical simulations for the coupling and temperature dependence on the critical current favorably compare to recent experimental results obtained at LENS with ultra-cold Fermi gases [4,5].
[1] Giancarlo Calvanese Strinati, Pierbiagio Pieri, Gerd Röpke, Peter Schuck, and Michael Urban. The BCS-BEC crossover: From ultra-cold Fermi gases to nuclear systems. Physics Reports, 738:1-76, April 2018.
[2] P. Pieri, L. Pisani, and G. C. Strinati. Bcs-bec crossover at finite temperature in the broken-symmetry phase. Phys. Rev. B, 70:094508, Sep 2004.
[3] S. Simonucci and G. C. Strinati. Equation for the superfluid gap obtained by coarse graining the bogoliubov ̆de gennes equations throughout the bcs-bec crossover. Phys. Rev. B, 89:054511, Feb 2014.
[4] W. J. Kwon, G. Del Pace, R. Panza, M. Inguscio, W. Zwerger, M. Zaccanti, F. Scazza, and G. Roati. Strongly correlated superfluid order parameters from dc josephson supercurrents. Science, 369(6499):84-88, 2020.
[5] G. Del Pace, W. J. Kwon, M. Zaccanti, G. Roati, and F. Scazza. Tunneling transport of unitary fermions across the superfluid transition. Phys. Rev. Lett., 126:055301, Feb 2021.

Orateur: Dr Verdiana Piselli (CNR-INO University of Camerino)

39 Ultracold fermions: from Nuclear to Atomic Physics

Orateur: Christophe Salomon (Ecole Normale supérieure)

Orsay_2023_Salomon_fin.pdf

40 Three-body contact of the resonant Fermi gas

For fermions with two internal states and two-body interactions of large scattering length $a$, we express the number of nearby fermion triplets in terms of a quantity $C_3$ , the three-body contact.

We calculate the three-body contact in a high-temperature regime, similar to a virial expansion. First, at the unitary limit, we use a wave-function approach and we find an analytical formula for the three-body contact in this case. Second, we use a diagrammatic approach in the BEC-BCS crossover. In this approach, the non-trivial scaling of the three-body correlation function is recovered, and we are able to calculate the three-body contact by numerically solving the 3-body problem. At unitarity, we numerically recover the result of the wave-function approach. These calculations could be used as benchmarks for comparisons with experiments.

Orateur: Prof. Xavier Leyronas (LPENS)

Leyronas-QMBC.pdf

10:45 Coffee break

Thursday 11:15-12:45

Président de session: Kutsal Bozkurt

41 Microscopic energy density functional: Success and limitations

Orateur: Dr Marcello Baldo (INFN Catania)

Baldo_Peter_memory.ppt

42 Constructing nuclear functionals for neutron stars and nucleosynthesis applications

Describing all different neutron star layers within a unified framework is a challenge in view of the very wide range of densities encountered. The description of massive pulsars such as J1614-2230 and J0740+662 [1], for example, brings the requirement of a stiff neutron matter equation of state in order to balance the strong gravity field of these compact objects. Additionally, the description of the rapid neutron-capture process (or r-process) nucleosynthesis taking place in neutron-star mergers requires detailed knowledge of nuclear reactions and radioactive decays (hence of the nuclear-structure properties, in particular nuclear masses) for a few thousand exotic neutron-rich nuclei. The challenge for nuclear theory is then the construction of a model that accurately describes: (i) masses of neutron-rich nuclei present in the crust of neutron stars, and (ii) masses of all the neutron-rich potentially produced during the r-process nucleosynthesis, together with (iii) a stiff enough neutron matter equation of state to explain the most massive observed pulsars.
A new family of microscopic nuclear energy density functionals and associated nuclear mass models have been recently proposed [2]. The Brussels-Skyrme-on-a-Grid (BSkG) functionals have the advantage of being based on the concept of utmost symmetry breaking for exotic nuclear configuration, allowing for exotic shapes like for instance triaxial or octupole deformation during the adjustment process. To compensate for the increase in computational cost, machine learning techniques were employed to optimize the parameter adjustment.
We show in this contribution the latest parametrization BSkG3, which greatly improves the infinite nuclear matter properties of BSkG functionals. To do so, we follow the procedure of Ref. [3] and use an extended form of the Skyrme functional to ensure a stiff enough neutron-matter equation of state at high densities. This new functional BSkG3 is consistent with observations of heavy pulsars. Furthermore, we include a pairing interaction designed to match the 1S0 pairing gaps in infinite nuclear matter deduced from ab-initio calculations. The latter is particularly important for a reliable description of superfluids in neutron stars. Both improvements, combined with our state-of-the-art description of atomic nuclei and simultaneous accurate description of many different observables, including in particular nuclear masses for thousands of nuclei, make BSkG3 a tool of choice for applications in nuclear structure and astrophysics.

[1] P. B. Demorest, et. al, Nature 467, 1081 (2010); .M.C. Miller et al., ApJL 918, L28 (2021); T.E. Riley et al., ApJL 918, L27 (2021).
[2] G. Scamps et al., Eur. Phys. J. A 57, 333 (2021); W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022).
[3] N. Chamel, S. Goriely and J. M. Pearson, Phys. Rev. C 80, 065804 (2009).

Orateur: Guilherme Grams (Astronomy and Astrophysics Institute (I.A.A.) - ULB)

Grams_QMBC23.pdf

43 Machine learning light hypernuclei

We employ a feed-forward artificial neural network to extrapolate at large model spaces the results of ab-initio hypernuclear No-Core Shell Model calculations for the $\Lambda$ separation energy $B_\Lambda$ of the lightest hypernuclei, $^3_\Lambda$H, $^4_\Lambda$H and $^4_\Lambda$He, obtained in computationally accessible harmonic oscillator basis spaces using chiral nucleon-nucleon, nucleon-nucleon-nucleon and hyperon-nucleon interactions. The overfitting problem is avoided by enlarging the size of the input dataset and by introducing a Gaussian noise during the training process of the neural network. We find that a network with a single hidden layer of eight neurons is sufficient to extrapolate correctly the value of the $\Lambda$ separation energy to model spaces of size $N_{max}=100$. The results obtained are in agreement with the experimental data in the case of $^3_\Lambda$H and the $0^+$ state of $^4_\Lambda$He, although they are off of the experiment by about $0.3$ MeV for both the $0^+$ and $1^+$states of $^4_\Lambda$H and the $1^+$ state of $^4_\Lambda$He. We find that our results are in excellent agreement with those obtained using other extrapolation schemes of the No-Core Shell Model calculations, showing this that an ANN is a reliable method to extrapolate the results of hypernuclear No-Core Shell Model calculations to large model spaces.

Orateur: Isaac Vidana (INFN)

Vidana_QMBC23.pdf

13:00 Lunch

Thursday 14:00-16:00

Président de session: Peter Ring (Physics Department Technical University Munich)

44 Pairing dynamics in nuclear collisions

I will present the results of nuclear collisions involving medium mass or heavy nuclei, obtained within time-dependent density functional theory (TDDFT) extended to superfluid systems. I will discuss the possible manifestations of pairing dynamics in nuclear collisions, at the vicinity of the Coulomb barrier. These include the mechanism for the increase of the barrier for capture generated by solitonic excitation appearing as a result of pairing phase distortion. Moreover, I will discuss pairing instability occuring in di-nuclear system formed by merging magic nuclei which lead to significant enhancement of pairing correlations.

Orateur: M. Andrzej Makowski (Warsaw University of Technology)

45 Fission properties of r-process nuclei predicted with the BCPM energy density functional

The rapid neutron capture process, or $r$ process, is responsible for the production of about half of the elements heavier than iron found in nature, including the heaviest uranium and thorium. During the $r$ process, several thousands of neutron-rich nuclei are synthesized in few seconds, powering an electromagnetic transient known as kilonova. Since most of such exotic nuclei have never been experimentally observed due to their exceedingly short half-lives, the estimation of abundances and kilonova light curves must rely upon the theoretical predictions of nuclear properties. During this talk, I will present calculations of nuclear properties obtained with the Barcelona-Catania-Paris-Madrid energy density functional (EDF) and their impact on the $r$ process. In particular, I will focus on the nucleosynthesis of translead elements in the merger of two neutron stars, and the role that nuclear masses, beta decays and fission play in shaping the $r$-process abundances and kilonova light curves.

Orateur: Samuel Andrea Giuliani (Universidad Autónoma de Madrid)

Giuliani_OrsayPS_2023.pdf

46 Superfluid dynamics in neutron stars

Produced during gravitational-core collapse supernova explosions with initial temperatures as high as ∼$10^{12}$ K, neutron stars cool down to temperatures $10^9$ K within a few days. The very dense matter in their interior is expected to undergo various quantum phase transitions analogous to those observed in terrestrial laboratories. Similarly to electrons in conventional terrestrial superconductors, free neutrons in the inner crust and the outer core of neutron stars are predicted to form a Bardeen-Cooper-Schrieffer (BCS) condensate of Cooper pairs. Nuclear superfluidity has found support from the rapid decline of luminosity of the Cassiopeia A remnant and has been corroborated by radio-timing observations of frequency glitches in numerous pulsars.

Despite the importance of the superfluid dynamics for interpreting these latter astrophysical phenomena, most microscopic calculations of the nuclear pairing properties have been carried out so far for static situations. We have recently studied the dynamics of hot neutron-proton superfluid mixtures within the self-consistent time-dependent nuclear energy density functional theory [1,2]. In application to neutron stars, we have computed $^{1}S_0$ neutron and proton pairing gaps in the homogeneous core in the presence of arbitrary currents and we have determined the mutual neutron-proton entrainment coupling coefficients [3].

We have also shown within the same framework that there exists a dynamical “gapless” state in which nuclear superfluidity is not destroyed even though the energy spectrum of quasiparticle excitations exhibits no gap. The absence of an energy gap leads to a nucleon specific heat that is very different from that in the classical BCS state (in the absence of superflows). The implications for the cooling of neutron stars will be discussed.

[1] N. Chamel & V. Allard, Phys. Rev. C 100, 065801 (2019).
[2] V. Allard & N. Chamel, Phys. Rev. C 103, 025804 (2021).
[3] V. Allard & N. Chamel, Universe 7, 470 (2021).

Orateur: Valentin Allard (Université Libre de Bruxelles (Institut d'Astronomie et d'Astrophysique))

Allard_QMBC2023.pdf

47 Scalar field, nucleon structure and relativistic chiral theory for nuclear matter

Orateur: Guy Chanfray (IPN Lyon)

Chanfray-QMBC.pptx

16:00 Coffee break