Source author record

Jordi Boronat

Jordi Boronat 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

18works
10topics
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

18 published item(s)

preprint2023arXiv

Absence of off-diagonal long-range order in hcp $^{\bf 4}$He dislocation cores

The mass transport properties along dislocation cores in hcp $^4$He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation core structures and their off-diagonal long-range order properties. It is found that the Bose-Einstein condensate fraction of such defective $^4$He systems is practically null ($\le 10^{-6}$), just as in the bulk defect-free crystal. These results provide compelling evidence for the absence of intrinsic superfluidity in dislocation cores in hcp $^4$He and challenge the superfluid dislocation-network interpretation of the mass-flux-experiment observations, calling for further experimental investigation.

preprint2022arXiv

Hole-induced anomaly in the thermodynamic behavior of a one-dimensional Bose gas

We reveal an intriguing anomaly in the temperature dependence of the specific heat of a one-dimensional Bose gas. The observed peak holds for arbitrary interaction and remembers a superfluid-to-normal phase transition in higher dimensions, but phase transitions are not allowed in one dimension. The presence of the anomaly signals a region of unpopulated states which behaves as an energy gap and is located below the hole branch in the excitation spectrum. The anomaly temperature is found to be of the same order of the energy of the maximum of the hole branch. We rely on the Bethe Ansatz to obtain the specific heat exactly and provide interpretations of the analytically tractable limits. The dynamic structure factor is computed with the Path Integral Monte Carlo method for the first time. We notice that at temperatures similar to the anomaly threshold, the energy of the thermal fluctuations become comparable with the maximal hole energy, leading to a qualitative change in the structure of excitations. This excitation pattern experiences the breakdown of the quasi-particle description for any value of the interaction strength at the anomaly, similarly to any superfluid phase transition at the critical temperature. We provide indications for future observations and how the hole anomaly can be employed for in-situ thermometry, identifying different collisional regimes and understanding other anomalies in atomic, solid-state, electronic, spin-chain and ladder systems.

preprint2022arXiv

Van der Waals five-body size-energy universality

A universal relationship between scaled size and scaled energy is explored in five-body self-bound quantum systems. The ground-state binding energy and structure properties are obtained by means of the diffusion Monte Carlo method. We use pure estimators to eliminate any residual bias in the estimation of the cluster size. Strengthening the inter-particle interaction, we extend the exploration from the halo region to classical systems. Universal scaled size-scaled energy line, which does not depend on the short-range potential details and binding strength, is found for homogeneous pentamers with interaction potentials decaying at long range predominantly as $r^{-6}$. For mixed pentamers, we discuss under which conditions the universal line can approximately describe the size-energy ratio. Our data is compatible with generalized Tjon lines, which assume a linear dependence between the binding energy of the pentamers and the one of tetramers, when both are divided by the trimer energies.

preprint2020arXiv

Finite-range effects in ultradilute quantum drops

In the first experimental realization of dilute Bose-Bose liquid drops using two hyperfine states of $^{39}$K some discrepancies between theory and experiment were observed. The standard analysis of the data using the Lee-Huang-Yang beyond mean-field theory predicted critical numbers which were significantly off the experimental measurements. Also, the radial size of the drops in the experiment proved to be larger than expected from this theory. Using a new functional, which is based on quantum Monte Carlo results of the bulk phase incorporating finite-range effects, we can explain the origin of the discrepancies in the critical number. This result proves the necessity of including finite-range corrections to deal with the observed properties in this setup. The controversy on the radial size is reasoned in terms of the departure from the optimal concentration ratio between the two species of the mixture.

preprint2016arXiv

Elusive structure of helium trimers

Over the years many He-He interaction potentials have been developed, some very sophisticated, including various corrections beyond Born-Oppenheimer approximation. Most of them were used to predict properties of helium dimers and trimers, examples of exotic quantum states, whose experimental study proved to be very challenging. Recently, detailed structural properties of helium trimers were measured for the first time, allowing a comparison with theoretical predictions and possibly enabling the evaluation of different interaction potentials. The comparisons already made included adjusting the maxima of both theoretical and experimental correlation functions to one, so the overall agreement between theory and experiment appeared satisfactory. However, no attempt was made to evaluate the quality of the interaction potentials used in the calculations. In this work, we calculate the experimentally measured correlation functions using both new and old potentials, compare them with experimental data and rank the potentials. We use diffusion Monte Carlo simulations at $T=0\ \mathrm{K}$, which give within statistical noise exact results of the ground state properties. All models predict both trimers $^4\mathrm{He}_{3}$ and $^4\mathrm{He}_2{}^3\mathrm{He}$ to be in a quantum halo state.

preprint2015arXiv

First-principles modeling of three-body interactions in highly compressed solid helium

We present a new set of three-body interaction models based on the Bruch-McGee (BM) potential that are suitable for the study of the energy, structural and elastic properties of solid 4He at high pressure. Our ab initio three-body potentials are obtained from the fit to total energies and atomic forces computed with the van der Waals density functional theory method due to Grimme, and represent an improvement with respect to previously reported three-body interaction models. In particular, we show that some of the introduced BM parametrizations reproduce closely the experimental equation of state and bulk modulus of solid helium up to a pressure of ~ 60 GPa, when used in combination with standard pairwise interaction models in diffusion Monte Carlo simulations. Importantly, we find that recent predictions reporting a surprisingly small variation of the kinetic energy and Lindeman ratio on quantum crystals under increasing pressure are likely to be artifacts produced by the use of incomplete interaction models. Also, we show that the experimental variation of the shear modulus, C44, at P < 25 GPa can be quantitatively described with the new set of three-body BM potentials. At higher pressures, however, the agreement between our C44 results and experiments deteriorates and thus we argue that higher order many-body terms in the expansion of the atomic interactions probably are necessary in order to better describe elasticity in very dense solid 4He.

preprint2015arXiv

Quantum Monte Carlo estimation of complex-time correlations for the study of the ground-state dynamic structure function

We present a method based on the Path Integral Monte Carlo formalism for the calculation of ground-state time correlation functions in quantum systems. The key point of the method is the consideration of time as a complex variable whose phase $δ$ acts as an adjustable parameter. By using high-order approximations for the quantum propagator, it is possible to obtain Monte Carlo data all the way from purely imaginary time to $δ$ values near the limit of real time. As a consequence, it is possible to infer accurately the spectral functions using simple inversion algorithms. We test this approach in the calculation of the dynamic structure function $S(q,ω)$ of two one-dimensional model systems, harmonic and quartic oscillators, for which $S(q,ω)$ can be exactly calculated. We notice a clear improvement in the calculation of the dynamic response with respect to the common approach based on the inverse Laplace transform of the imaginary-time correlation function.

preprint2014arXiv

Atomic monolayer deposition on the surface of nanotube mechanical resonators

We studied monolayers of noble gas atoms (Xe, Kr, Ar, and Ne) deposited on individual ultra- clean suspended nanotubes. For this, we recorded the resonance frequency of the mechanical motion of the nanotube, since it provides a direct measure of the coverage. The latter is the number of adsorbed atoms divided by the number of the carbon atoms of the suspended nanotube. Monolayers formed when the temperature was lowered in a constant pressure of noble gas atoms. The coverage of Xe monolayers remained constant at 1/6 over a large temperature range. This finding reveals that Xe monolayers are solid phases with a triangular atomic arrangement, and are commensurate with the underlying carbon nanotube. By comparing our measurements to theoretical calculations, we identify the phases of Ar and Ne monolayers as fluids, and we tentatively describe Kr monolayers as solid phases. These results underscore that mechanical resonators made from single nanotubes are excellent probes for surface science.

preprint2014arXiv

First-Principles Modeling of Quantum Nuclear Effects and Atomic Interactions in Solid 4He at High Pressure

We present a first-principles computational study of solid 4He at T=0K and pressures up to 160GPa. Our computational strategy consists in using van der Waals density functional theory (DFT-vdW) to describe the electronic degrees of freedom in this material, and the diffusion Monte Carlo (DMC) method to solve the Schrodinger equation describing the behavior of the quantum nuclei. For this, we construct an analytical interaction function based on the pairwise Aziz potential that closely matches the volume variation of the cohesive energy calculated with DFT-vdW in dense helium. Interestingly, we find that the kinetic energy of solid 4He does not increase appreciably with compression for P > 85GPa. Also, we show that the Lindemann ratio in dense solid 4He amounts to 0.10 almost independently of pressure. The reliability of customary quasi-harmonic DFT (QH DFT) approaches in the description of quantum nuclear effects in solids is also studied. We find that QH DFT simulations, although provide a reasonable equation of state in agreement with experiments, are not able to reproduce correctly these critical effects in compressed 4He. In particular, we disclose huge discrepancies of at least ~50% in the calculated 4He kinetic energies using both the QH DFT and present DFT-DMC methods.

preprint2014arXiv

The Limit of Mechanical Stability in Quantum Crystals: A Diffusion Monte Carlo Study of Solid 4He

We present a first-principles study of the energy and elastic properties of solid helium at pressures below the range in which is energetically stable. We find that the limit of mechanical stability in hcp 4He is $P_{s}$ = -33.82 bar, which lies significantly below the spinodal pressure found in the liquid phase (i.e., -9.6 bar). Furthermore, we show that the pressure variation of the transverse and longitudinal sound velocities close to $P_{s}$ do not follow a power law of the form $\propto \left( P - P_{s} \right)^γ$, in contrast to what is observed on the fluid.

preprint2013arXiv

Elastic Constants of Incommensurate Solid 4He

We study the elastic properties of incommensurate solid 4He in the limit of zero temperature. Specifically, we calculate the pressure dependence of the five elastic constants (C11, C12, C13, C33, and C44), longitudinal and transversal speeds of sound, and the T = 0 Debye temperature of incommensurate and commensurate hcp 4He using the diffusion Monte Carlo method. Our results show that under compression the commensurate crystal is globally stiffer than the incommensurate, however at pressures close to melting (i.e. P ~ 25 bars) some of the elastic constants accounting for strain deformations of the hcp basal plane (C12 and C13) are slightly larger in the incommensurate solid. Also, we find that upon the introduction of tiny concentrations of point defects the shear modulus of 4He (C44) undergoes a small reduction.

preprint2013arXiv

Possible superfluidity of molecular hydrogen in a two-dimensional crystal phase of sodium

We theoretically investigate the ground-state properties of a molecular para-hydrogen (p-H2) film in which crystallization is energetically frustrated by embedding sodium (Na) atoms periodically distributed in a triangular lattice. In order to fully deal with the quantum nature of p-H2 molecules, we employ the diffusion Monte Carlo method and realistic semi-empirical pairwise potentials describing the interactions between H2-H2 and Na-H2 species. In particular, we calculate the energetic, structural and superfluid properties of two-dimensional Na-H2 systems within a narrow density interval around equilibrium at zero temperature. In contrast to previous computational studies considering other alkali metal species such as rubidium and potassium, we find that the p-H2 ground-state is a liquid with a significantly large superfluid fraction of ~30%. The appearance of p-H2 superfluid response is due to the fact that the interactions between Na atoms and H2 molecules are less attractive than between H2 molecules. This induces a considerable reduction of the hydrogen density which favours the stabilization of the liquid phase.

preprint2012arXiv

Excitations and Stripe Phase Formation in a 2D Dipolar Bose Gas with Tilted Polarization

We present calculations of the ground state and excitations of an anisotropic dipolar Bose gas in two dimensions, realized by a non-perpendicular polarization with respect to the system plane. For sufficiently high density an increase of the polarization angle leads to a density instability of the gas phase in the direction where the anisotropic interaction is strongest. Using a dynamic many-body theory, we calculate the dynamic structure function in the gas phase which shows the anisotropic dispersion of the excitations. We find that the energy of roton excitations in the strongly interacting direction decreases with increasing polarization angle and almost vanishes close to the instability. Exact path integral ground state Monte Carlo simulations show that this instability is indeed a quantum phase transition to a stripe phase, characterized by long-range order in the strongly interacting direction.

preprint2012arXiv

Onset temperature of Bose-Einstein condensation in incommensurate solid 4He

The temperature dependence of the one-body density matrix in 4He crystals presenting vacancies is computed with Path Integral Monte Carlo. The main purpose of this study is to estimate the onset temperature T_0 of Bose-Einstein condensation in these systems. We see that T_0 depends on the vacancy concentration X_v of the simulated system, but not following the law $T_0 \sim X_v^{2/3}$ obtained assuming non-interacting vacancies. For the lowest X_v we have studied, that is X_v = 1/256, we get T_0 = (0.15 \pm 0.05) K, close to the temperatures at which a finite fraction of non-classical rotational inertia is experimentally observed. Below T_0, vacancies do not act as classical point defects becoming completely delocalized entities.

preprint2011arXiv

Condensate fraction in liquid 4He at zero temperature

We present results of the one-body density matrix (OBDM) and the condensate fraction n_0 of liquid 4He calculated at zero temperature by means of the Path Integral Ground State Monte Carlo method. This technique allows to generate a highly accurate approximation for the ground state wave function Psi_0 in a totally model-independent way, that depends only on the Hamiltonian of the system and on the symmetry properties of Psi_0. With this unbiased estimation of the OBDM, we obtain precise results for the condensate fraction n_0 and the kinetic energy K of the system. The dependence of n_0 with the pressure shows an excellent agreement of our results with recent experimental measurements. Above the melting pressure, overpressurized liquid 4He shows a small condensate fraction that has dropped to 0.8% at the highest pressure of p = 87 bar.

preprint2011arXiv

Microscopic approach to the bcc phase of solid 4He

The coexistence between the hcp and bcc phases of solid 4He at fixed pressure, and at three different temperatures, is studied by means of the path integral Monte Carlo method. Microscopic results for the main energetic and structure properties of the two phases are reported. The differences between hcp and bcc 4He are shown to be small with the obvious exception of the static structure factor. When crossing the phase transition line, the most significant changes are observed in the kinetic energy per particle and in the Lindemann ratio, both pointing to a less correlated quantum solid in the bcc side. No off-diagonal long-range order is observed in the the hcp and bcc perfect lattices.

preprint2010arXiv

Path Integral Monte Carlo calculation of momentum distribution in solid \^4\He

We perform calculations of the momentum distribution n(k) in solid \^4\He by means of path integral Monte Carlo methods. We see that, in perfect crystal, n(k) does not depend on temperature T and that is different from the classical Gaussian shape of Maxwell-Boltzmann distribution, even though these discrepancies decrease when the density of the system increases. In crystal presenting vacancies, we see that for T \ge 0.75K, n(k) presents the same behavior as in the perfect crystal, but, at lower T, it presents a peak when k\to0.

preprint2009arXiv

Thermal and quantum fluctuations in chains of ultracold polar molecules

Ultracold polar molecules, in highly anisotropic traps and interacting via a repulsive dipolar potential, may form one-dimensional chains at high densities. According to classical theory, at low temperatures there exists a critical value of the density at which a second order phase transition from a linear to a zigzag chain occurs. We study the effect of thermal and quantum fluctuations on these self-organized structures using classical and quantum Monte Carlo methods, by means of which we evaluate the pair correlation function and the static structure factor. Depending on the parameters, these functions exhibit properties typical of a crystalline or of a liquid system. We compare the thermal and the quantum results, identifying analogies and differences. Finally, we discuss experimental parameter regimes where the effects of quantum fluctuations on the linear - zigzag transition can be observed.