Finally, we investigate the fidelity susceptibility as an indicator of the topological phase transitions, and find a simple scaling law as a function of the number of bands of our multiband model which is found to be the same for both bulk and boundary fidelity susceptibilities. Expand abstract.

We investigate the Loschmidt amplitude and dynamical quantum phase transitions in multiband one dimensional topological insulators. For this purpose we introduce a new solvable multiband model based on the Su-Schrieffer-Heeger model, generalized to unit cells containing many atoms but with the same symmetry properties. Such models have a richer structure of dynamical quantum phase transitions than the simple two-band topological insulator models typically considered previously, with both quasiperiodic and aperiodic dynamical quantum phase transitions present. Moreover the aperiodic transitions can still occur for quenches within a single topological phase. We also investigate the boundary contributions from the presence of the topologically protected edge states of this model. Plateaus in the boundary return rate are related to the topology of the time evolving Hamiltonian, and hence to a dynamical bulk-boundary correspondence. We go on to consider the dynamics of the entanglement entropy generated after a quench, and its potential relation to the critical times of the dynamical quantum phase transitions. Finally, we investigate the fidelity susceptibility as an indicator of the topological phase transitions, and find a simple scaling law as a function of the number of bands of our multiband model which is found to be the same for both bulk and boundary fidelity susceptibilities.

Recently quantum simulators have been constructed to investigate experimentally the most prominent theoretical four-point many-body system described by the Hubbard model. Expand abstract.

Recently quantum simulators have been constructed to investigate experimentally the most prominent theoretical four-point many-body system described by the Hubbard model. By varying the coupling strength of the four-point interaction in relation to the kinetic term, one can analyze the phase structure of the model. This intriguing fact leads us to investigate whether similar Hamiltonians with four-point interactions can also be studied as a function of their four-point coupling strength. In this paper, we reexamine the Nambu-Jona-Lasinio model, regarding it generally beyond the context of quantum chromodynamics. Essentially, it is a model in which particle-antiparticle pairing leads to a BCS-like condensate, with the result that chiral symmetry is broken dynamically in the strong-coupling regime. To study the behavior of the system, it is necessary to move from this regime to a hypothetical regime of weak coupling, altering the coupling strength of the interaction arbitrarily. In order to do this, the gap equation must be regarded as complex and its Riemann surface structure must be known. We do this and obtain a continuous quantum phase transition characterized by the development of a complex order parameter (the dynamically generated mass) from the second sheet of the Riemann surface, as we move into the weak-coupling regime. The power-law behavior of the order parameter in the vicinity of the phase transition point is demonstrated to be independent of the choice of the regularization scheme with the critical exponent as $\beta \approx 0.55$. At the same time, the isovector pseudoscalar modes retain their feature as Goldstone modes and still have zero mass, while the isoscalar scalar meson follows the behavior of the order parameter and gains a width. Energetically, this mode is not favored over the normal, uncondensed mode but would have to be accessed through an excitation process.

We present results on phase transitions of local and global survival in a two-species model on Gilbert graphs. Expand abstract.

We present results on phase transitions of local and global survival in a two-species model on Gilbert graphs. At initial time there is an infection at the origin that propagates on the Gilbert graph according to a continuous-time nearest-neighbor interacting particle system. The Gilbert graph consists of susceptible nodes and nodes of a second type, which we call white knights. The infection can spread on susceptible nodes without restriction. If the infection reaches a white knight, this white knight starts to spread on the set of infected nodes according to the same mechanism, with a potentially different rate, giving rise to a competition of chase and escape. We show well-definedness of the model, isolate regimes of global survival and extinction of the infection and present estimates on local survival. The proofs rest on comparisons to the process on trees, percolation arguments and finite-degree approximations of the underlying random graphs.

Layered perovskite ruthenium oxides exhibit a striking series of metal-insulator and magnetic-nonmagnetic phase transitions easily tuned by temperature, pressure, epitaxy, and nonlinear drive. Expand abstract.

Layered perovskite ruthenium oxides exhibit a striking series of metal-insulator and magnetic-nonmagnetic phase transitions easily tuned by temperature, pressure, epitaxy, and nonlinear drive. In this work, we combine results from two complementary state of the art many-body methods, Auxiliary Field Quantum Monte Carlo and Dynamical Mean Field Theory, to determine the low-temperature phase diagram of Ca$_2$RuO$_4$. Both methods predict a low temperature, pressure-driven metal-insulator transition accompanied by a ferromagnetic-antiferromagnetic transition. The properties of the ferromagnetic state vary non-monotonically with pressure and are dominated by the ruthenium $d_{xy}$ orbital, while the properties of the antiferromagnetic state are dominated by the $d_{xz}$ and $d_{yz}$ orbitals. Differences of detail in the predictions of the two methods are analyzed. This work is theoretically important as it presents the first application of the Auxiliary Field Quantum Monte Carlo method to an orbitally-degenerate system with both Mott and Hunds physics, and provides an important comparison of the Dynamical Mean Field and Auxiliary Field Quantum Monte Carlo methods.

The fully-connected Ising $p$-spin model has for $p >2$ a discontinuous phase transition from the paramagnetic phase to a stable state with one-step replica symmetry breaking (1RSB). Expand abstract.

The fully-connected Ising $p$-spin model has for $p >2$ a discontinuous phase transition from the paramagnetic phase to a stable state with one-step replica symmetry breaking (1RSB). However, simulations in three dimension do not look like these mean-field results and have features more like those which would arise with full replica symmetry breaking (FRSB). To help understand how this might come about we have studied in the fully connected $p$-spin model the state of two-step replica symmetry breaking (2RSB). It has a free energy degenerate with that of 1RSB, but the weight of the additional peak in $P(q)$ vanishes. We expect that the state with full replica symmetry breaking (FRSB) is also degenerate with that of 1RSB. We suggest that finite size effects will give a non-vanishing weight to the FRSB features, as also will fluctuations about the mean-field solution. Our conclusion is that outside the fully connected model in the thermodynamic limit, FRSB is to be expected rather than 1RSB.

Loops are one of the fundamental structures that trace the geometry of the magnetic field in the solar atmosphere. Expand abstract.

Loops are one of the fundamental structures that trace the geometry of the magnetic field in the solar atmosphere. Their evolution and dynamics provide a crucial proxy for studying how the magnetized structures are formed and heated in the solar atmosphere. Here, we report on spectroscopic observations of a set of transition region loops taken by the Interface Region Imaging Spectrograph (IRIS) at Si IV 1394 \AA\ with a sit-and-stare mode. The loops are corresponding to the flux emergence at its very late phase when the emerged magentic features in the photosphere have fully developed. We find the transition region loops are still expanding and moving upward with a velocity of a few kilometers per second ($\lesssim$10 km/s) at this stage. The expansion of the loops leads to interactions between themselves and the ambient field, which can drive magnetic reconnection evidenced by multiple intense brightenings, including transition region explosive events and IRIS bombs in the footpoint region associated with the moving polarity. A set of quasi-periodic brightenings with a period of about 130 s is found at the loop apex, from which the Si IV 1394 \AA\ profiles are significantly non-Gaussian with enhancements at both blue and red wings at Doppler velocities of about 50 km/s. We suggest that the transition region loops in the very late phase of flux emergence can be powered by heating events generated by the interactions between the expanding loops and the ambient fields and also by (quasi-)periodic processes, such as oscillation-modulated braiding reconnection.

Many extensions of the Standard Model include a new $U(1)$ gauge group that is broken spontaneously at a scale much above TeV. Expand abstract.

Many extensions of the Standard Model include a new $U(1)$ gauge group that is broken spontaneously at a scale much above TeV. If a $U(1)$-breaking phase transition occurs at nucleation temperature of $O(100)$-$O(1000)$~TeV, it can generate stochastic gravitational waves in $O(10)$-$O(100)$ Hz range if $\beta_{\rm n}/H_{\rm n}=1000$, which can be detected by ground-based detectors. Meanwhile, supersymmetry (SUSY) may play a crucial role in the dynamics of such high-scale $U(1)$ gauge symmetry breaking, because SUSY breaking scale is expected to be at TeV to solve the hierarchy problem. In this paper, we study the phase transition of $U(1)$ gauge symmetry breaking in a SUSY model in the SUSY limit. We consider a particular example, the minimal SUSY $U(1)_{B-L}$ model. We derive the finite temperature effective potential of the model in the SUSY limit, study a $U(1)_{B-L}$-breaking phase transition, and estimate gravitational waves generated from it.

Extending our models to allow for non-equilibrium terms yields new stable phases which qualitatively agrees with observations in vivo. Expand abstract.

Epigenetics is a driving force of important and ubiquitous phenomena in nature such as cell differentiation or even metamorphosis. Oppositely to its widespread role, understanding the biophysical principles that allow epigenetics to control and rewire gene regulatory networks remains an open challenge. In this work we study the effects of epigenetic modifications on the spatial folding of chromosomes -- and hence on the expression of the underlying genes -- by mapping the problem to a class of models known as magnetic polymers. In this work we show that a first order phase transition underlies the simultaneous spreading of certain epigenetic marks and the conformational collapse of a chromosome. Further, we describe Brownian Dynamics simulations of the model in which the topology of the polymer and thermal fluctuations are fully taken into account and that confirm our mean field predictions. Extending our models to allow for non-equilibrium terms yields new stable phases which qualitatively agrees with observations in vivo. Our results show that statistical mechanics techniques applied to models of magnetic polymers can be successfully exploited to rationalize the outcomes of experiments designed to probe the interplay between a dynamic epigenetic landscape and chromatin organization.

Interestingly, we find that there is no finite-temperature phase transition between these two regions though rather a smooth crossover. Expand abstract.

Using extensive Monte Carlo simulations, we investigate the surface adsorption of self-avoiding trails on the triangular lattice with two- and three-body on-site monomer-monomer interactions. In the parameter space of two-body, three-body, and surface interaction strengths, the phase diagram displays four phases: swollen (coil), globule, crystal, and adsorbed. For small values of the surface interaction, we confirm the presence of swollen, globule, and crystal bulk phases. For sufficiently large values of the surface interaction, the system is in an adsorbed state, and the adsorption transition can be continuous or discontinuous, depending on the bulk phase. As such, the phase diagram contains a rich phase structure with transition surfaces that meet in multicritical lines joining in a single special multicritical point. The adsorbed phase displays two distinct regions with different characteristics, dominated by either single or double layer adsorbed ground states. Interestingly, we find that there is no finite-temperature phase transition between these two regions though rather a smooth crossover.

Finally, we measure the associated magnetic--monopole current induced around the chromo-flux tube and observe that the confinement/deconfinement phase transition is associated with the appearance/disappearance of the induced magnetic--monopole current, respectively. Expand abstract.

We investigate the confinement-deconfinement phase transition at finite temperature of the SU(3) Yang-Mills(YM) theory on the lattice from a viewpoint of the dual superconductor picture based on the novel reformulation of the YM theory. In particular, we compare the conventional Abelian dual superconductor picture with the non-Abelian dual superconductor picture proposed in our previous works as the mechanism of quark confinement in the SU(3) YM theory. For the SU(3) YM theory, the reformulation allows two possible options called maximal and minimal. The maximal option corresponds to the manifestly gauge-invariant extension of the Abelian projection scheme, while the minimal option is really new to give the non-Abelian dual superconductor picture. Keeping these differences in mind, we present the numerical evidences that the confinement/deconfinement phase transition is caused by appearance/disappearance of the dual Meissner effects. First, we measure the Polyakov loop average at various temperatures to determine the critical temperature separating the low-temperature confined phase and the high-temperature deconfined phase. Second, we measure the static quark-antiquark potential. Third, we measure the chromoelectric and chromomagnetic flux created by a pair of quark and antiquark at temperatures below and above the critical temperature. We observe no more squeezing of the chromoelectric flux tube above the critical temperature. Finally, we measure the associated magnetic--monopole current induced around the chromo-flux tube and observe that the confinement/deconfinement phase transition is associated with the appearance/disappearance of the induced magnetic--monopole current, respectively. We confirm that these results are also obtained by the restricted field alone in both options, indicating the restricted field dominance in quark confinement at finite temperature.