@article{14650, abstract = {We study the out-of-equilibrium quantum dynamics of dipolar polarons, i.e., impurities immersed in a dipolar Bose-Einstein condensate, after a quench of the impurity-boson interaction. We show that the dipolar nature of the condensate and of the impurity results in anisotropic relaxation dynamics, in particular, anisotropic dressing of the polaron. More relevantly for cold-atom setups, quench dynamics is strongly affected by the interplay between dipolar anisotropy and trap geometry. Our findings pave the way for simulating impurities in anisotropic media utilizing experiments with dipolar mixtures.}, author = {Volosniev, Artem and Bighin, Giacomo and Santos, Luis and Peña Ardila, Luisllu A.}, issn = {2542-4653}, journal = {SciPost Physics}, keywords = {General Physics and Astronomy}, number = {6}, publisher = {SciPost Foundation}, title = {{Non-equilibrium dynamics of dipolar polarons}}, doi = {10.21468/scipostphys.15.6.232}, volume = {15}, year = {2023}, } @article{13966, abstract = {We present a low-scaling diagrammatic Monte Carlo approach to molecular correlation energies. Using combinatorial graph theory to encode many-body Hugenholtz diagrams, we sample the Møller-Plesset (MPn) perturbation series, obtaining accurate correlation energies up to n=5, with quadratic scaling in the number of basis functions. Our technique reduces the computational complexity of the molecular many-fermion correlation problem, opening up the possibility of low-scaling, accurate stochastic computations for a wide class of many-body systems described by Hugenholtz diagrams.}, author = {Bighin, Giacomo and Ho, Quoc P and Lemeshko, Mikhail and Tscherbul, T. V.}, issn = {2469-9969}, journal = {Physical Review B}, number = {4}, publisher = {American Physical Society}, title = {{Diagrammatic Monte Carlo for electronic correlation in molecules: High-order many-body perturbation theory with low scaling}}, doi = {10.1103/PhysRevB.108.045115}, volume = {108}, year = {2023}, } @article{11592, abstract = {We compare recent experimental results [Science 375, 528 (2022)] of the superfluid unitary Fermi gas near the critical temperature with a thermodynamic model based on the elementary excitations of the system. We find good agreement between experimental data and our theory for several quantities such as first sound, second sound, and superfluid fraction. We also show that mode mixing between first and second sound occurs. Finally, we characterize the response amplitude to a density perturbation: Close to the critical temperature both first and second sound can be excited through a density perturbation, whereas at lower temperatures only the first sound mode exhibits a significant response.}, author = {Bighin, Giacomo and Cappellaro, Alberto and Salasnich, L.}, issn = {2469-9934}, journal = {Physical Review A}, number = {6}, publisher = {American Physical Society}, title = {{Unitary Fermi superfluid near the critical temperature: Thermodynamics and sound modes from elementary excitations}}, doi = {10.1103/PhysRevA.105.063329}, volume = {105}, year = {2022}, } @article{11997, abstract = {We study the fate of an impurity in an ultracold heteronuclear Bose mixture, focusing on the experimentally relevant case of a ⁴¹K - ⁸⁷Rb mixture, with the impurity in a ⁴¹K hyperfine state. Our paper provides a comprehensive description of an impurity in a BEC mixture with contact interactions across its phase diagram. We present results for the miscible and immiscible regimes, as well as for the impurity in a self-bound quantum droplet. Here, varying the interactions, we find exotic states where the impurity localizes either at the center or at the surface of the droplet. }, author = {Bighin, Giacomo and Burchianti, A. and Minardi, F. and Macrì, T.}, issn = {2469-9934}, journal = {Physical Review A}, number = {2}, publisher = {American Physical Society}, title = {{Impurity in a heteronuclear two-component Bose mixture}}, doi = {10.1103/PhysRevA.106.023301}, volume = {106}, year = {2022}, } @article{11998, abstract = {Recently it became possible to study highly excited rotational states of molecules in superfluid helium through nonadiabatic alignment experiments (Cherepanov et al 2021 Phys. Rev. A 104 L061303). This calls for theoretical approaches that go beyond explaining renormalized values of molecular spectroscopic constants, which suffices when only the lowest few rotational states are involved. As the first step in this direction, here we present a basic quantum mechanical model describing highly excited rotational states of molecules in superfluid helium nanodroplets. We show that a linear molecule immersed in a superfluid can be seen as an effective symmetric top, similar to the rotational structure of radicals, such as OH or NO, but with the angular momentum of the superfluid playing the role of the electronic angular momentum in free molecules. The simple theory sheds light onto what happens when the rotational angular momentum of the molecule increases beyond the lowest excited states accessible by infrared spectroscopy. In addition, the model allows to estimate the effective rotational and centrifugal distortion constants for a broad range of species and to explain the crossover between light and heavy molecules in superfluid 4He in terms of the many-body wavefunction structure. Some of the above mentioned insights can be acquired by analyzing a simple 2 × 2 matrix.}, author = {Cherepanov, Igor and Bighin, Giacomo and Schouder, Constant A. and Chatterley, Adam S. and Stapelfeldt, Henrik and Lemeshko, Mikhail}, issn = {1367-2630}, journal = {New Journal of Physics}, number = {7}, publisher = {IOP}, title = {{A simple model for high rotational excitations of molecules in a superfluid}}, doi = {10.1088/1367-2630/ac8113}, volume = {24}, year = {2022}, } @article{12154, abstract = {We review our theoretical results of the sound propagation in two-dimensional (2D) systems of ultracold fermionic and bosonic atoms. In the superfluid phase, characterized by the spontaneous symmetry breaking of the U(1) symmetry, there is the coexistence of first and second sound. In the case of weakly-interacting repulsive bosons, we model the recent measurements of the sound velocities of 39K atoms in 2D obtained in the weakly-interacting regime and around the Berezinskii–Kosterlitz–Thouless (BKT) superfluid-to-normal transition temperature. In particular, we perform a quite accurate computation of the superfluid density and show that it is reasonably consistent with the experimental results. For superfluid attractive fermions, we calculate the first and second sound velocities across the whole BCS-BEC crossover. In the low-temperature regime, we reproduce the recent measurements of first-sound speed with 6Li atoms. We also predict that there is mixing between sound modes only in the finite-temperature BEC regime.}, author = {Salasnich, Luca and Cappellaro, Alberto and Furutani, Koichiro and Tononi, Andrea and Bighin, Giacomo}, issn = {2073-8994}, journal = {Symmetry}, keywords = {Physics and Astronomy (miscellaneous), General Mathematics, Chemistry (miscellaneous), Computer Science (miscellaneous)}, number = {10}, publisher = {MDPI}, title = {{First and second sound in two-dimensional bosonic and fermionic superfluids}}, doi = {10.3390/sym14102182}, volume = {14}, year = {2022}, } @article{9606, abstract = {Sound propagation is a macroscopic manifestation of the interplay between the equilibrium thermodynamics and the dynamical transport properties of fluids. Here, for a two-dimensional system of ultracold fermions, we calculate the first and second sound velocities across the whole BCS-BEC crossover, and we analyze the system response to an external perturbation. In the low-temperature regime we reproduce the recent measurements [Phys. Rev. Lett. 124, 240403 (2020)] of the first sound velocity, which, due to the decoupling of density and entropy fluctuations, is the sole mode excited by a density probe. Conversely, a heat perturbation excites only the second sound, which, being sensitive to the superfluid depletion, vanishes in the deep BCS regime and jumps discontinuously to zero at the Berezinskii-Kosterlitz-Thouless superfluid transition. A mixing between the modes occurs only in the finite-temperature BEC regime, where our theory converges to the purely bosonic results.}, author = {Tononi, A. and Cappellaro, Alberto and Bighin, Giacomo and Salasnich, L.}, issn = {24699934}, journal = {Physical Review A}, number = {6}, publisher = {American Physical Society}, title = {{Propagation of first and second sound in a two-dimensional Fermi superfluid}}, doi = {10.1103/PhysRevA.103.L061303}, volume = {103}, year = {2021}, } @article{10631, abstract = {We combine experimental and theoretical approaches to explore excited rotational states of molecules embedded in helium nanodroplets using CS2 and I2 as examples. Laser-induced nonadiabatic molecular alignment is employed to measure spectral lines for rotational states extending beyond those initially populated at the 0.37 K droplet temperature. We construct a simple quantum-mechanical model, based on a linear rotor coupled to a single-mode bosonic bath, to determine the rotational energy structure in its entirety. The calculated and measured spectral lines are in good agreement. We show that the effect of the surrounding superfluid on molecular rotation can be rationalized by a single quantity, the angular momentum, transferred from the molecule to the droplet.}, author = {Cherepanov, Igor and Bighin, Giacomo and Schouder, Constant A. and Chatterley, Adam S. and Albrechtsen, Simon H. and Muñoz, Alberto Viñas and Christiansen, Lars and Stapelfeldt, Henrik and Lemeshko, Mikhail}, issn = {2469-9934}, journal = {Physical Review A}, number = {6}, publisher = {American Physical Society}, title = {{Excited rotational states of molecules in a superfluid}}, doi = {10.1103/PhysRevA.104.L061303}, volume = {104}, year = {2021}, } @article{8170, abstract = {Alignment of OCS, CS2, and I2 molecules embedded in helium nanodroplets is measured as a function of time following rotational excitation by a nonresonant, comparatively weak ps laser pulse. The distinct peaks in the power spectra, obtained by Fourier analysis, are used to determine the rotational, B, and centrifugal distortion, D, constants. For OCS, B and D match the values known from IR spectroscopy. For CS2 and I2, they are the first experimental results reported. The alignment dynamics calculated from the gas-phase rotational Schrödinger equation, using the experimental in-droplet B and D values, agree in detail with the measurement for all three molecules. The rotational spectroscopy technique for molecules in helium droplets introduced here should apply to a range of molecules and complexes.}, author = {Chatterley, Adam S. and Christiansen, Lars and Schouder, Constant A. and Jørgensen, Anders V. and Shepperson, Benjamin and Cherepanov, Igor and Bighin, Giacomo and Zillich, Robert E. and Lemeshko, Mikhail and Stapelfeldt, Henrik}, issn = {10797114}, journal = {Physical Review Letters}, number = {1}, publisher = {American Physical Society}, title = {{Rotational coherence spectroscopy of molecules in Helium nanodroplets: Reconciling the time and the frequency domains}}, doi = {10.1103/PhysRevLett.125.013001}, volume = {125}, year = {2020}, } @article{8587, abstract = {Inspired by the possibility to experimentally manipulate and enhance chemical reactivity in helium nanodroplets, we investigate the effective interaction and the resulting correlations between two diatomic molecules immersed in a bath of bosons. By analogy with the bipolaron, we introduce the biangulon quasiparticle describing two rotating molecules that align with respect to each other due to the effective attractive interaction mediated by the excitations of the bath. We study this system in different parameter regimes and apply several theoretical approaches to describe its properties. Using a Born–Oppenheimer approximation, we investigate the dependence of the effective intermolecular interaction on the rotational state of the two molecules. In the strong-coupling regime, a product-state ansatz shows that the molecules tend to have a strong alignment in the ground state. To investigate the system in the weak-coupling regime, we apply a one-phonon excitation variational ansatz, which allows us to access the energy spectrum. In comparison to the angulon quasiparticle, the biangulon shows shifted angulon instabilities and an additional spectral instability, where resonant angular momentum transfer between the molecules and the bath takes place. These features are proposed as an experimentally observable signature for the formation of the biangulon quasiparticle. Finally, by using products of single angulon and bare impurity wave functions as basis states, we introduce a diagonalization scheme that allows us to describe the transition from two separated angulons to a biangulon as a function of the distance between the two molecules.}, author = {Li, Xiang and Yakaboylu, Enderalp and Bighin, Giacomo and Schmidt, Richard and Lemeshko, Mikhail and Deuchert, Andreas}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {16}, publisher = {AIP Publishing}, title = {{Intermolecular forces and correlations mediated by a phonon bath}}, doi = {10.1063/1.5144759}, volume = {152}, year = {2020}, } @article{8644, abstract = {Determining the phase diagram of systems consisting of smaller subsystems 'connected' via a tunable coupling is a challenging task relevant for a variety of physical settings. A general question is whether new phases, not present in the uncoupled limit, may arise. We use machine learning and a suitable quasidistance between different points of the phase diagram to study layered spin models, in which the spin variables constituting each of the uncoupled systems (to which we refer as layers) are coupled to each other via an interlayer coupling. In such systems, in general, composite order parameters involving spins of different layers may emerge as a consequence of the interlayer coupling. We focus on the layered Ising and Ashkin–Teller models as a paradigmatic case study, determining their phase diagram via the application of a machine learning algorithm to the Monte Carlo data. Remarkably our technique is able to correctly characterize all the system phases also in the case of hidden order parameters, i.e. order parameters whose expression in terms of the microscopic configurations would require additional preprocessing of the data fed to the algorithm. We correctly retrieve the three known phases of the Ashkin–Teller model with ferromagnetic couplings, including the phase described by a composite order parameter. For the bilayer and trilayer Ising models the phases we find are only the ferromagnetic and the paramagnetic ones. Within the approach we introduce, owing to the construction of convolutional neural networks, naturally suitable for layered image-like data with arbitrary number of layers, no preprocessing of the Monte Carlo data is needed, also with regard to its spatial structure. The physical meaning of our results is discussed and compared with analytical data, where available. Yet, the method can be used without any a priori knowledge of the phases one seeks to find and can be applied to other models and structures.}, author = {Rzadkowski, Wojciech and Defenu, N and Chiacchiera, S and Trombettoni, A and Bighin, Giacomo}, issn = {13672630}, journal = {New Journal of Physics}, number = {9}, publisher = {IOP Publishing}, title = {{Detecting composite orders in layered models via machine learning}}, doi = {10.1088/1367-2630/abae44}, volume = {22}, year = {2020}, } @article{6940, abstract = {We study the effect of a linear tunneling coupling between two-dimensional systems, each separately exhibiting the topological Berezinskii-Kosterlitz-Thouless (BKT) transition. In the uncoupled limit, there are two phases: one where the one-body correlation functions are algebraically decaying and the other with exponential decay. When the linear coupling is turned on, a third BKT-paired phase emerges, in which one-body correlations are exponentially decaying, while two-body correlation functions exhibit power-law decay. We perform numerical simulations in the paradigmatic case of two coupled XY models at finite temperature, finding evidences that for any finite value of the interlayer coupling, the BKT-paired phase is present. We provide a picture of the phase diagram using a renormalization group approach.}, author = {Bighin, Giacomo and Defenu, Nicolò and Nándori, István and Salasnich, Luca and Trombettoni, Andrea}, issn = {1079-7114}, journal = {Physical Review Letters}, number = {10}, publisher = {American Physical Society}, title = {{Berezinskii-Kosterlitz-Thouless paired phase in coupled XY models}}, doi = {10.1103/physrevlett.123.100601}, volume = {123}, year = {2019}, } @article{5886, abstract = {Problems involving quantum impurities, in which one or a few particles are interacting with a macroscopic environment, represent a pervasive paradigm, spanning across atomic, molecular, and condensed-matter physics. In this paper we introduce new variational approaches to quantum impurities and apply them to the Fröhlich polaron–a quasiparticle formed out of an electron (or other point-like impurity) in a polar medium, and to the angulon–a quasiparticle formed out of a rotating molecule in a bosonic bath. We benchmark these approaches against established theories, evaluating their accuracy as a function of the impurity-bath coupling.}, author = {Li, Xiang and Bighin, Giacomo and Yakaboylu, Enderalp and Lemeshko, Mikhail}, issn = {00268976}, journal = {Molecular Physics}, publisher = {Taylor and Francis}, title = {{Variational approaches to quantum impurities: from the Fröhlich polaron to the angulon}}, doi = {10.1080/00268976.2019.1567852}, year = {2019}, } @article{6339, abstract = {We introduce a diagrammatic Monte Carlo approach to angular momentum properties of quantum many-particle systems possessing a macroscopic number of degrees of freedom. The treatment is based on a diagrammatic expansion that merges the usual Feynman diagrams with the angular momentum diagrams known from atomic and nuclear structure theory, thereby incorporating the non-Abelian algebra inherent to quantum rotations. Our approach is applicable at arbitrary coupling, is free of systematic errors and of finite-size effects, and naturally provides access to the impurity Green function. We exemplify the technique by obtaining an all-coupling solution of the angulon model; however, the method is quite general and can be applied to a broad variety of systems in which particles exchange quantum angular momentum with their many-body environment.}, author = {Bighin, Giacomo and Tscherbul, Timur and Lemeshko, Mikhail}, journal = {Physical Review Letters}, number = {16}, publisher = {American Physical Society}, title = {{Diagrammatic Monte Carlo approach to angular momentum in quantum many-particle systems}}, doi = {10.1103/physrevlett.121.165301}, volume = {121}, year = {2018}, } @article{417, abstract = {We introduce a Diagrammatic Monte Carlo (DiagMC) approach to complex molecular impurities with rotational degrees of freedom interacting with a many-particle environment. The treatment is based on the diagrammatic expansion that merges the usual Feynman diagrams with the angular momentum diagrams known from atomic and nuclear structure theory, thereby incorporating the non-Abelian algebra inherent to quantum rotations. Our approach works at arbitrary coupling, is free of systematic errors and of finite size effects, and naturally provides access to the impurity Green function. We exemplify the technique by obtaining an all-coupling solution of the angulon model, however, the method is quite general and can be applied to a broad variety of quantum impurities possessing angular momentum degrees of freedom. }, author = {Bighin, Giacomo and Tscherbul, Timur and Lemeshko, Mikhail}, journal = {Physical Review Letters}, number = {16}, publisher = {American Physical Society}, title = {{Diagrammatic Monte Carlo approach to rotating molecular impurities}}, doi = {10.1103/PhysRevLett.121.165301}, volume = {121}, year = {2018}, } @article{420, abstract = {We analyze the theoretical derivation of the beyond-mean-field equation of state for two-dimensional gas of dilute, ultracold alkali-metal atoms in the Bardeen–Cooper–Schrieffer (BCS) to Bose–Einstein condensate (BEC) crossover. We show that at zero temperature our theory — considering Gaussian fluctuations on top of the mean-field equation of state — is in very good agreement with experimental data. Subsequently, we investigate the superfluid density at finite temperature and its renormalization due to the proliferation of vortex–antivortex pairs. By doing so, we determine the Berezinskii–Kosterlitz–Thouless (BKT) critical temperature — at which the renormalized superfluid density jumps to zero — as a function of the inter-atomic potential strength. We find that the Nelson–Kosterlitz criterion overestimates the BKT temperature with respect to the renormalization group equations, this effect being particularly relevant in the intermediate regime of the crossover.}, author = {Bighin, Giacomo and Salasnich, Luca}, journal = {International Journal of Modern Physics B}, number = {17}, pages = {1840022}, publisher = {World Scientific Publishing}, title = {{Renormalization of the superfluid density in the two-dimensional BCS-BEC crossover}}, doi = {10.1142/S0217979218400222}, volume = {32}, year = {2018}, } @article{1015, abstract = {Vortices are commonly observed in the context of classical hydrodynamics: from whirlpools after stirring the coffee in a cup to a violent atmospheric phenomenon such as a tornado, all classical vortices are characterized by an arbitrary circulation value of the local velocity field. On the other hand the appearance of vortices with quantized circulation represents one of the fundamental signatures of macroscopic quantum phenomena. In two-dimensional superfluids quantized vortices play a key role in determining finite-temperature properties, as the superfluid phase and the normal state are separated by a vortex unbinding transition, the Berezinskii-Kosterlitz-Thouless transition. Very recent experiments with two-dimensional superfluid fermions motivate the present work: we present theoretical results based on the renormalization group showing that the universal jump of the superfluid density and the critical temperature crucially depend on the interaction strength, providing a strong benchmark for forthcoming investigations.}, author = {Bighin, Giacomo and Salasnich, Luca}, issn = {20452322}, journal = {Scientific Reports}, publisher = {Nature Publishing Group}, title = {{Vortices and antivortices in two-dimensional ultracold Fermi gases}}, doi = {10.1038/srep45702}, volume = {7}, year = {2017}, } @article{995, abstract = {Recently it was shown that an impurity exchanging orbital angular momentum with a surrounding bath can be described in terms of the angulon quasiparticle [Phys. Rev. Lett. 118, 095301 (2017)]. The angulon consists of a quantum rotor dressed by a many-particle field of boson excitations, and can be formed out of, for example, a molecule or a nonspherical atom in superfluid helium, or out of an electron coupled to lattice phonons or a Bose condensate. Here we develop an approach to the angulon based on the path-integral formalism, which sets the ground for a systematic, perturbative treatment of the angulon problem. The resulting perturbation series can be interpreted in terms of Feynman diagrams, from which, in turn, one can derive a set of diagrammatic rules. These rules extend the machinery of the graphical theory of angular momentum - well known from theoretical atomic spectroscopy - to the case where an environment with an infinite number of degrees of freedom is present. In particular, we show that each diagram can be interpreted as a 'skeleton', which enforces angular momentum conservation, dressed by an additional many-body contribution. This connection between the angulon theory and the graphical theory of angular momentum is particularly important as it allows to systematically and substantially simplify the analytical representation of each diagram. In order to exemplify the technique, we calculate the 1- and 2-loop contributions to the angulon self-energy, the spectral function, and the quasiparticle weight. The diagrammatic theory we develop paves the way to investigate next-to-leading order quantities in a more compact way compared to the variational approaches.}, author = {Bighin, Giacomo and Lemeshko, Mikhail}, issn = {24699950}, journal = {Physical Review B - Condensed Matter and Materials Physics}, number = {8}, publisher = {American Physical Society}, title = {{Diagrammatic approach to orbital quantum impurities interacting with a many-particle environment}}, doi = {10.1103/PhysRevB.96.085410}, volume = {96}, year = {2017}, }