Researcher profile

Frédéric Mila

Frédéric Mila contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

19 published item(s)

preprint2024arXiv

Critical line of the triangular Ising antiferromagnet in a field from a $C_3$-symmetric corner transfer matrix algorithm

The corner transfer matrix renormalization group (CTMRG) algorithm has been extensively used to investigate both classical and quantum two-dimensional (2D) lattice models. The convergence of the algorithm can strongly vary from model to model depending on the underlying geometry and symmetries, and the presence of algebraic correlations. An important factor in the convergence of the algorithm is the lattice symmetry, which can be broken due to the necessity of mapping the problem onto the square lattice. We propose a variant of the CTMRG algorithm, designed for models with $C_3$-symmetry, which we apply to the conceptually simple yet numerically challenging problem of the triangular lattice Ising antiferromagnet in a field, at zero and low temperatures. We study how the finite-temperature three-state Potts critical line in this model approaches the ground-state Kosterlitz-Thouless transition driven by a reduced field ($h/T$). In this particular instance, we show that the $C_3$-symmetric CTMRG leads to much more precise results than both existing results from exact diagonalization of transfer matrices and Monte Carlo.

preprint2023arXiv

Critical properties of the quantum Ashkin-Teller chain with chiral perturbations

We investigate the nature of the phase transitions in the quantum Ashkin-Teller chain in the presence of chiral perturbations. We locate the Lifshitz line separating a region of direct chiral transitions from the region where the transition is through an intermediate floating phase. Furthermore, we identify a small region in the vicinity of the four-state Potts point where chiral perturbations are irrelevant and where the transition remains conformal. Implications to Rydberg atoms experiments are briefly discussed.

preprint2022arXiv

Commensurate-incommensurate transition in the chiral Ashkin-Teller model

We investigate the classical chiral Ashkin-Teller model on a square lattice with the corner transfer matrix renormalisation group (CTMRG) algorithm. We show that the melting of the period-4 phase in the presence of a chiral perturbation takes different forms depending on the coefficient of the four-spin term in the Ashkin-Teller model. Close to the clock limit of two decoupled Ising models, the system undergoes a two-step commensurate-incommensurate transition as soon as the chirality is introduced, with an intermediate critical floating phase bounded by a Kosterlitz-Thouless transition at high temperature and a Pokrovsky-Talapov transition at low temperature. By contrast, close to the four-states Potts model, we argue for the existence of a unique commensurate-incommensurate transition that belongs to the chiral universality class, and for the presence of a Lifshitz point where the ordered, disordered and floating phases meet. Finally, we map the whole phase diagram, which turns out to be in qualitative agreement with the 40 year old prediction by Huse and Fisher.

preprint2022arXiv

From SU(2)$_5$ to SU(2)$_3$ Wess-Zumino-Witten transitions in a frustrated spin-5/2 chain

We investigate the properties of a frustrated spin-5/2 chain with next-nearest neighbor two and three-site interactions, with two questions in mind: the nature of the transition into the dimerized phase induced by the three-site interaction, and the possible presence of a critical floating phase at intermediate values of the next-nearest neighbor interaction. We provide strong evidence that the continuous transition into the dimerized phase, which has been found to be generically in the Wess-Zumino-Witten SU(2)$_{2S}$ universality class up to spin $S=2$, is SU(2)$_5$ only at two isolated points of the phase diagram, and that it is SU(2)$_3$ in between, in agreement with the presence of two relevant operators allowed by symmetry for SU(2)$_5$, and with the conservation of the parity of the level index along the renormalization flow between SU(2)$_k$ theories with different values of $k$. We also find that the dimerization induced by the next-nearest neighbor interaction is a three step process, with first a small partially dimerized phase followed by a broad critical floating phase with incommensurate correlations before the fully dimerized phase is reached. Implications for the iron oxide ${\mathrm{Bi}}_{3}{\mathrm{FeMo}}_{2}{\mathrm{O}}_{12}$ are briefly discussed.

preprint2022arXiv

Possibility to detect the bound state of the Heisenberg ferromagnetic chain at intermediate temperature

Motivated by the lack of direct evidence with inelastic neutron scattering of the well documented bound state of Heisenberg ferromagnets, we use the time-dependent Thermal Density Matrix Renormalization Group algorithm to study the temperature dependence of the dynamical spin structure factor of Heisenberg ferromagnetic spin chains. For spin-1/2, we show that the bound state appears as a well defined excitation with significant spectral weight in the temperature range $J/12 \lesssim T \lesssim J/3$, pointing to the possibility of detecting it with inelastic neutron scattering near $k=π$ provided the temperature is neither too low nor too high - at low temperature, the spectral weight only grows as $T^{3/2}$, and at high temperature the bound state peak merges with the two-magnon continuum. For spin-1, the situation is more subtle because the bound state with two neighboring spin flips competes with an anti-bound state with two spin-flips on the same site. As a consequence, the relative spectral weight of the bound state is smaller than for spin-1/2, and a weak resonance due to the anti-bound state appears in the continuum. A clearer signature of the bound state (resp. anti-bound state) can be obtained if a negative (resp. positive) biquadratic interaction is present.

preprint2022arXiv

Spinon confinement and deconfinement in spin-1 chains

Motivated by the rich phase diagrams recently unveiled in frustrated spin-1 chains with competing next-nearest neighbour and four-spin interactions, we investigate the nature of the elementary excitations of spin-1 chains in the vicinity of phase transitions using the matrix-product state (MPS) representation of elementary excitations. First, we show that both spinons and magnons naturally arise in SU(2) invariant spin chains when describing ground states and elementary excitations using MPS. Then we investigate the nature of the elementary excitations across the first-oder transition between the Haldane phase and the topologically trivial next-nearest neighbor Haldane phase. We show explicitly that spinons deconfine at the transition, and we calculate the dispersion of these deconfined spinons with MPS. We also show that, immediately away from the transition line, spinons confine, forming dispersive spinon/anti-spinon bound states on both sides of the transition. Finally, we show that deconfined spinons also appear at the transition between the Haldane phase and the spontaneously dimerized phase when this transition is first order.

preprint2022arXiv

Thermal Ising transition in the spin-1/2 J1-J2 Heisenberg model

Using an SU(2) invariant finite-temperature tensor network algorithm, we provide strong numerical evidence in favor of an Ising transition in the collinear phase of the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. In units of $J_2$, the critical temperature reaches a maximal value of $T_c/J_2\simeq 0.18$ around $J_2/J_1\simeq 1.0$. It is strongly suppressed upon approaching the zero-temperature boundary of the collinear phase $J_2/J_1\simeq 0.6$, and it vanishes as $1/\log(J_2/J_1)$ in the large $J_2/J_1$ limit, as predicted by Chandra, Coleman and Larkin [Phys. Rev. Lett. 64, 88, 1990]. Enforcing the SU(2) symmetry is crucial to avoid the artifact of finite-temperature SU(2) symmetry breaking of U(1) algorithms, opening new perspectives in the investigation of the thermal properties of quantum Heisenberg antiferromagnets.

preprint2021arXiv

Phase diagram of the Shastry-Sutherland Compound SrCu2(BO3)2 under extreme combined conditions of field and pressure

Motivated by the intriguing properties of the Shastry-Sutherland compound SrCu2(BO3)2 under pressure, with a still debated intermediate plaquette phase appearing at around 20 kbar and a possible deconfined critical point at higher pressure upon entering the antiferromagnetic phase, we have investigated its high-field properties in this pressure range using tunnel diode oscillator (TDO) measurements. The two main new phases revealed by these measurements are fully consistent with those identified by infinite Projected Entangled Pair states (iPEPS) calculations of the Shastry-Sutherland model, a 1/5 plateau and a 10 x 2 supersolid. Remarkably, these phases are descendants of the full-plaquette phase, the prominent candidate for the intermediate phase of SrCu2(BO3)2. The emerging picture for SrCu2(BO3)2 is shown to be that of a system dominated by a tendency to an orthorhombic distortion at intermediate pressure, an important constraint on any realistic description of the transition into the antiferromagnetic phase.

preprint2021arXiv

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Quantum Monte Carlo (QMC) simulations constitute nowadays one of the most powerful methods to study strongly correlated quantum systems, provided that no "sign problem" arises. However, many systems of interest, including highly frustrated magnets, suffer from an average sign that is close to zero in standard QMC simulations. Nevertheless, a possible sign problem depends on the simulation basis, and here we demonstrate how a suitable choice of cluster basis can be used to eliminate or at least reduce the sign problem in highly frustrated magnets that were so far inaccessible to efficient QMC simulations. We focus in particular on the application of a two-spin (dimer)-based QMC method to the thermodynamics of the spin-1/2 Shastry-Sutherland model for SrCu$_2$(BO$_3$)$_2$.

preprint2020arXiv

Floating, critical and dimerized phases in a frustrated spin-3/2 chain

We study spontaneous dimerization and emergent criticality in a spin-3/2 chain with antiferromagnetic nearest-neighbor $J_1$, next-nearest-neighbor $J_2$ and three-site $J_3$ interactions. In the absence of three-site interaction $J_3$, we provide evidence that the model undergoes a remarkable sequence of three phase transitions as a function of $J_2/J_1$, going successively through a critical commesurate phase, a partially dimerized gapped phase, a critical floating phase with quasi-long-range incommensurate order, to end up in a fully dimerized phase at very large $J_2/J_1$. In the field theory language, this implies that the coupling constant of the marginal operator responsible for dimerization changes sign three times. For large enough $J_3$, the fully dimerized phase is stabilized for all $J_2$, and the phase transitions between the critical phases and this phase are both Wess-Zumino-Witten (WZW) SU(2)$_3$ along part of the boundary and turn first order at some point due to the presence of a marginal operator in the WZW SU(2)$_3$ model. By contrast, the transition between the two dimerized phase is always first order, and the phase transitions between the partially dimerized phase and the critical phases are Kosterlitz-Thouless. Finally, we discuss the intriguing spin-1/2 edge states that emerge in the partially dimerized phase for even chains. Unlike their counterparts in the spin-1 chain, they are not confined and disappear upon increasing $J_2$ in favour of a reorganization of the dimerization pattern.

preprint2020arXiv

Haldane Gap of the Three-Box Symmetric $\mathrm{SU}(3)$ Chain

Motivated by the recent generalization of the Haldane conjecture to $\mathrm{SU}(3)$ chains [M. Lajkó et al., Nucl. Phys. B924, 508 (2017)] according to which a Haldane gap should be present for symmetric representations if the number of boxes in the Young diagram is a multiple of three, we develop a density matrix renormalization group algorithm based on standard Young tableaus to study the model with three boxes directly in the representations of the global $\mathrm{SU}(3)$ symmetry. We show that there is a finite gap between the singlet and the symmetric $[3\,0\,0]$ sector $Δ_{[3\,0\,0]}/J = 0.040\pm0.006$ where $J$ is the antiferromagnetic Heisenberg coupling, and we argue on the basis of the structure of the low energy states that this is sufficient to conclude that the spectrum is gapped.

preprint2020arXiv

Kibble-Zurek exponent and chiral transition of the period-4 phase of Rydberg chains

Chains of Rydberg atoms have emerged as an amazing playground to study quantum physics in 1D. Playing with inter-atomic distances and laser detuning, one can in particular explore the commensurate-incommensurate transition out of charge-density waves through the Kibble-Zurek mechanism, and the possible presence of a chiral transition with dynamical exponent $z>1$. Here we address this problem theoretically with effective blockade models where the short-distance repulsions are replaced by a constraint of no double occupancy. For the period-4 phase, we show there is an Ashkin-Teller transition point with exponent $ν=0.78$ surrounded by a direct chiral transition with a dynamical exponent $z=1.14$ and a Kibble-Zurek exponent $μ=0.4$. For Rydberg atoms with a van der Waals potential, we suggest that the experimental value $μ=0.25$ is due to a chiral transition with $z\simeq 1.9$ and $ν\simeq 0.47$ surrounding an Ashkin-Teller transition close to the 4-state Potts universality.

preprint2020arXiv

Magnetic field induced bound states in spin-$\frac{1}{2}$ ladders

Motivated by the intriguing mode splittings in a magnetic field recently observed with inelastic neutron scattering in the spin ladder compound (C$_5$H$_{12}$N)$_2$CuBr$_4$ (BPCB), we investigate the nature of the spin ladder excitations using DMRG and analytical arguments. Starting from the fully frustrated ladder, for which we derive the low-energy spectrum, we show that bound states are generically present close to $k=0$ in the dynamical structure factor of spin ladders above $H_{c 1}$, and that they are characterized by a field-independent binding energy and an intensity that grows with $H-H_{c 1}$. These predictions are shown to explain quantitatively the split modes observed in BPCB.

preprint2020arXiv

Reduction of the sign problem near $T=0$ in quantum Monte Carlo simulations

Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian---obtained by making all off-diagonal elements negative in a given basis---have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.

preprint2020arXiv

Solving frustrated Ising models using tensor networks

Motivated by the recent success of tensor networks to calculate the residual entropy of spin ice and kagome Ising models, we develop a general framework to study frustrated Ising models in terms of infinite tensor networks %, i.e. tensor networks that can be contracted using standard algorithms for infinite systems. This is achieved by reformulating the problem as local rules for configurations on overlapping clusters chosen in such a way that they relieve the frustration, i.e. that the energy can be minimized independently on each cluster. We show that optimizing the choice of clusters, including the weight on shared bonds, is crucial for the contractibility of the tensor networks, and we derive some basic rules and a linear program to implement them. We illustrate the power of the method by computing the residual entropy of a frustrated Ising spin system on the kagome lattice with next-next-nearest neighbour interactions, vastly outperforming Monte Carlo methods in speed and accuracy. The extension to finite-temperature is briefly discussed.

preprint2020arXiv

Spin-$\frac{1}{2}$ kagome Heisenberg antiferromagnet with strong breathing anisotropy

We study the zero-temperature phase diagram of the spin-$\frac{1}{2}$ Heisenberg model with breathing anisotropy (i.e., with different coupling strength on the upward and downward triangles) on the kagome lattice. Our study relies on large scale tensor network simulations based on infinite projected entangled-pair state and infinite projected entangled-simplex state methods adapted to the kagome lattice. Our energy analysis suggests that the U(1) algebraic quantum spin-liquid (QSL) ground-state of the isotropic Heisenberg model is stable up to very large breathing anisotropy until it breaks down to a critical lattice-nematic phase that breaks rotational symmetry in real space through a first-order quantum phase transition. Our results also provide further insight into the recent experiment on vanadium oxyfluoride compounds which has been shown to be relevant platforms for realizing QSL in the presence of breathing anisotropy.

preprint2020arXiv

SrCu$_2$(BO$_3$)$_2$ under pressure: a first-principle study

Using density-functional theory (DFT) band-structure calculations, we study the crystal structure, the lattice dynamics, and the magnetic interactions in the Shastry-Sutherland magnet SrCu$_2$(BO$_3$)$_2$ under pressure, concentrating on experimentally relevant pressures up to 4 GPa. All tendencies inside the candidate intermediate phases are thoroughly worked out, including specific predictions for some Raman active phonon modes that could be used to pin down the nature of the intermediate phase.

preprint2019arXiv

Generalization of the Haldane conjecture to SU($n$) chains

Recently, SU(3) chains in the symmetric and self-conjugate representations have been studied using field theory techniques. For certain representations, namely rank-$p$ symmetric ones with $p$ not a multiple of 3, it was argued that the ground state exhibits gapless excitations. For the remaining representations considered, a finite energy gap exists above the ground state. In this paper, we extend these results to SU($n$) chains in the symmetric representation. For a rank-$p$ symmetric representation with $n$ and $p$ coprime, we predict gapless excitations above the ground state. If $p$ is a multiple of $n$, we predict a unique ground state with a finite energy gap. Finally, if $p$ and $n$ have a greatest common divisor $1<q<n$, we predict a ground state degeneracy of $n/q$, with a finite energy gap. To arrive at these results, we derive a non-Lorentz invariant flag manifold sigma model description of the SU($n$) chains, and use the renormalization group to show that Lorentz invariance is restored at low energies. We then make use of recently developed anomaly matching conditions for these Lorentz-invariant models. We also review the Lieb-Shultz-Mattis-Affleck theorem, and extend it to SU($n$) models with longer range interactions.

preprint2019arXiv

Perturbative approach to quantum tunneling and quantum interferences in spin clusters

Collective tunneling is a ubiquitous phenomenon in finite-size spin clusters that shows up in systems as diverse as molecular magnets or spin clusters adsorbed at surfaces. The problem we explore is to understand how small flipping terms can cooperate to flip a large spin to the opposite direction or a cluster of interacting elementary Ising spins into the time-reversed state. These high order processes will involve at least two channels, a single spin-flip channel due to a transverse field and a two-spin flip channel due to exchange or due to single-ion anisotropies. In the present paper, we show that high-order perturbation theory can be formulated and evaluated with the help of simple recurrence relations, leading to a compact theory of tunnelling in macroscopic spins, in one-dimensional clusters, as well as in small higher-dimensional clusters. This is demonstrated explicitly in the case of the Ising model with a transverse field and transverse exchange, and in the case of macroscopic spins with uniaxial anisotropy. Our approach provides a transparent theory of level crossings, where the tunneling between time reversed configurations vanishes as a function of the external field. The crossings result from destructive quantum interferences between competing flipping channels. Our theory consistently predicts N crossings in chains of N Ising spins, 2S crossings in single spins of magnitude S, and yields explicit analytical formulae for the level crossings of open chains and macroscopic spins. Disorder can be easily implemented in this perturbative formalism. Leading disorder effects can be treated analytically for spin rings. At the smallest transverse field crossing, the suppression of tunneling is most robust with respect to disorder. We briefly discuss the implications of our findings for the use of realistic spin clusters on surfaces to store information.