Researcher profile

A. M. Somoza

A. M. Somoza contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
7topics
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

15 published item(s)

preprint2015arXiv

Deconfined Quantum Criticality, Scaling Violations, and Classical Loop Models

Numerical studies of the Néel to valence-bond solid phase transition in 2D quantum antiferromagnets give strong evidence for the remarkable scenario of deconfined criticality, but display strong violations of finite-size scaling that are not yet understood. We show how to realise the universal physics of the Neel-VBS transition in a 3D classical loop model (this includes the interference effect that suppresses Néel hedgehogs). We use this model to simulate unprecedentedly large systems (of size $L\leq 512$). Our results are compatible with a direct continuous transition at which both order parameters are critical, and we do not see conventional signs of first-order behaviour. However, we find that the scaling violations are stronger than previously realised and are incompatible with conventional finite-size scaling over the size range studied, even if allowance is made for a weakly/marginally irrelevant scaling variable. In particular, different determinations of the anomalous dimensions $η_\text{VBS}$ and $η_\text{Néel}$ yield very different results. The assumption of conventional finite-size scaling gives estimates which drift to negative values at large $L$, in violation of unitarity bounds. In contrast, the behaviour of correlators on scales much smaller than $L$ is consistent with large positive anomalous dimensions. Barring an unexpected reversal in behaviour at still larger sizes, this implies that the transition, if continuous, must show unconventional finite-size scaling, e.g. from a dangerously irrelevant scaling variable. Another possibility is an anomalously weak first-order transition. By analysing the renormalisation group flows for the non-compact $CP^{n-1}$ model ($n$-component Abelian Higgs model) between two and four dimensions, we give the simplest scenario by which an anomalously weak first-order transition can arise without fine-tuning of the Hamiltonian.

preprint2015arXiv

Emergent $SO(5)$ Symmetry at the Néel to Valence-Bond-Solid Transition

We show numerically that the `deconfined' quantum critical point between the Néel antiferromagnet and the columnar valence-bond-solid, for a square lattice of spin-1/2s, has an emergent $SO(5)$ symmetry. This symmetry allows the Néel vector and the valence-bond-solid order parameter to be rotated into each other. It is a remarkable 2+1-dimensional analogue of the $SO(4)= [SU(2)\times SU(2)]/Z_2$ symmetry that appears in the scaling limit for the spin-1/2 Heisenberg chain. The emergent $SO(5)$ is strong evidence that the phase transition in the 2+1D system is truly continuous, despite the violations of finite-size scaling observed previously in this problem. It also implies surprising relations between correlation functions at the transition. The symmetry enhancement is expected to apply generally to the critical two-component Abelian Higgs model (non-compact $CP^1$ model). The result indicates that in three dimensions there is an $SO(5)$-symmetric conformal field theory which has no relevant singlet operators, so is radically different to conventional Wilson-Fisher-type conformal field theories.

preprint2015arXiv

Unbinding transition in semi-infinite two-dimensional localized systems

We consider a two-dimensional strongly localized system defined in a half-space and whose transfer integral in the edge can be different than in the bulk. We predict an unbinding transition, as the edge transfer integral is varied, from a phase where conduction paths are distributed across the bulk to a bound phase where propagation is mainly along the edge. At criticality the logarithm of the conductance follows the $F_1$ Tracy-Widom distribution. We verify numerically these predictions for both the Anderson and the Nguyen, Spivak and Shklovskii models. We also check that for a half-space, i.e., when the edge transfer integral is equal to the bulk transfer integral, the distribution of the conductance is the $F_4$ Tracy-Widom distribution. These findings are strong indications that random signs directed polymer models and their quantum extensions belong to the Kardar-Parisi- Zhang universality class. We have analyzed finite-size corrections at criticality and for a half-plane.

preprint2013arXiv

Length Distributions in Loop Soups

Statistical lattice ensembles of loops in three or more dimensions typically have phases in which the longest loops fill a finite fraction of the system. In such phases it is natural to ask about the distribution of loop lengths. We show how to calculate moments of these distributions using $CP^{n-1}$ or $RP^{n-1}$ and O(n) $σ$ models together with replica techniques. The resulting joint length distribution for macroscopic loops is Poisson-Dirichlet with a parameter $θ$ fixed by the loop fugacity and by symmetries of the ensemble. We also discuss features of the length distribution for shorter loops, and use numerical simulations to test and illustrate our conclusions.

preprint2013arXiv

Loop models with crossings

The universal behaviour of two-dimensional loop models can change dramatically when loops are allowed to cross. We study models with crossings both analytically and with extensive Monte Carlo simulations. Our main focus (the 'completely packed loop model with crossings') is a simple generalisation of well-known models which shows an interesting phase diagram with continuous phase transitions of a new kind. These separate the unusual 'Goldstone' phase observed previously from phases with short loops. Using mappings to Z_2 lattice gauge theory, we show that the continuum description of the model is a replica limit of the sigma model on real projective space (RP^{n-1}). This field theory sustains Z_2 point defects which proliferate at the transition. In addition to studying the new critical points, we characterise the universal properties of the Goldstone phase in detail, comparing renormalisation group (RG) calculations with numerical data on systems of linear size up to L=10^6 at loop fugacity n=1. (Very large sizes are necessary because of the logarithmic form of correlation functions and other observables.) The model is relevant to polymers on the verge of collapse, and a particular point in parameter space maps to self-avoiding trails at their Θ-point; we use the RG treatment of a perturbed sigma model to resolve some perplexing features in the previous literature on trails. Finally, one of the phase transitions considered here is a close analogue of those in disordered electronic systems --- specifically, Anderson metal-insulator transitions --- and provides a simpler context in which to study the properties of these poorly-understood (central-charge-zero) critical points.

preprint2013arXiv

Phase transitions in 3D loop models and the $CP^{n-1}$ $σ$ model

We consider the statistical mechanics of a class of models involving close-packed loops with fugacity $n$ on three-dimensional lattices. The models exhibit phases of two types as a coupling constant is varied: in one, all loops are finite, and in the other, some loops are infinitely extended. We show that the loop models are discretisations of $CP^{n-1}$ $σ$ models. The finite and infinite loop phases represent, respectively, disordered and ordered phases of the $σ$ model, and we discuss the relationship between loop properties and $σ$ model correlators. On large scales, loops are Brownian in an ordered phase and have a non-trivial fractal dimension at a critical point. We simulate the models, finding continuous transitions between the two phases for $n=1,2,3$ and first order transitions for $n\geq 4$. We also give a renormalisation group treatment of the $CP^{n-1}$ model that shows how a continuous transition can survive for values of $n$ larger than (but close to) two, despite the presence of a cubic invariant in the Landau-Ginzburg description. The results we obtain are of broader relevance to a variety of problems, including SU(n) quantum magnets in (2+1) dimensions, Anderson localisation in symmetry class C, and the statistics of random curves in three dimensions.

preprint2012arXiv

Quantum Coulomb gap in low dimensions

We study the single-particle density of states of one-dimensional and two-dimensional quantum disordered systems with long-range interactions. We consider a $1/\sqrt{r}$ interaction in one dimension and a Coulomb interaction in two dimensions, which produce linear gaps in the density of states in both cases. We focus on the strong localization regime where the localization length is small but non-zero. We use an exact diagonalization technique for small system sizes and a perturbation approach for larger sizes. We find that, with both methods, the inclusion of a finite hopping contribution does not change the linear character of the gap, but reduces its slope, widening the gap.

preprint2012arXiv

Re-entrance and entanglement in the one-dimensional Bose-Hubbard model

Re-entrance is a novel feature where the phase boundaries of a system exhibit a succession of transitions between two phases A and B, like A-B-A-B, when just one parameter is varied monotonically. This type of re-entrance is displayed by the 1D Bose Hubbard model between its Mott insulator (MI) and superfluid phase as the hopping amplitude is increased from zero. Here we analyse this counter-intuitive phenomenon directly in the thermodynamic limit by utilizing the infinite time-evolving block decimation algorithm to variationally minimize an infinite matrix product state (MPS) parameterized by a matrix size chi. Exploiting the direct restriction on the half-chain entanglement imposed by fixing chi, we determined that re-entrance in the MI lobes only emerges in this approximate when chi >= 8. This entanglement threshold is found to be coincident with the ability an infinite MPS to be simultaneously particle-number symmetric and capture the kinetic energy carried by particle-hole excitations above the MI. Focussing on the tip of the MI lobe we then applied, for the first time, a general finite-entanglement scaling analysis of the infinite order Kosterlitz-Thouless critical point located there. By analysing chi's up to a very moderate chi = 70 we obtained an estimate of the KT transition as t_KT = 0.30 +/- 0.01, demonstrating the how a finite-entanglement approach can provide not only qualitative insight but also quantitatively accurate predictions.

preprint2011arXiv

3D loop models and the CP^{n-1} sigma model

Many statistical mechanics problems can be framed in terms of random curves; we consider a class of three-dimensional loop models that are prototypes for such ensembles. The models show transitions between phases with infinite loops and short-loop phases. We map them to $CP^{n-1}$ sigma models, where $n$ is the loop fugacity. Using Monte Carlo simulations, we find continuous transitions for $n=1,2,3$, and first order transitions for $n\geq 5$. The results are relevant to line defects in random media, as well as to Anderson localization and $(2+1)$-dimensional quantum magnets.

preprint2011arXiv

Effects of many-electron jumps in relaxation and conductivity of Coulomb glasses

A numerical study of the energy relaxation and conductivity of the Coulomb glass is presented. The role of many-electron transitions is studied by two complementary methods: a kinetic Monte Carlo algorithm and a master equation in configuration space method. A calculation of the transition rate for two-electron transitions is presented, and the proper extension of this to multi-electron transitions is discussed. It is shown that two-electron transitions are important in bypassing energy barriers which effectively block sequential one-electron transitions. The effect of two-electron transitions is also discussed.

preprint2011arXiv

Nonlinear conductivity of two-dimensional Coulomb glasses

We have studied the nonlinear conductivity of two-dimensional Coulomb glasses. We have used a Monte Carlo algorithm to simulate the dynamic of the system under an applied electric field $E$. We found that in the nonlinear regime the site occupancy in the Coulomb gap follows a Fermi-Dirac distribution with an effective temperature $T_{\rm eff}$, higher than the phonon bath temperature $T$. The value of the effective temperature is compatible with that obtained for slow modes from the generalized fluctuation-dissipation theorem. The nonlinear conductivity for a given electric field and $T$ is fairly similar to the linear conductivity at the corresponding $T_{\rm eff}$. We found that the dissipated power and the effective temperature are related by an expression of the form $(T_{\rm eff}^α-T^α)T_{\rm eff}^{β-α}$.

preprint2011arXiv

Phase diagram of weak-magnetic-field quantum Hall transition quantified from classical percolation

We consider magnetotransport in high-mobility 2D electron gas in a non-quantizing magnetic field. We employ a weakly chiral network model to test numerically the prediction of the scaling theory that the transition from an Anderson to a quantum Hall insulator takes place when the Drude value of the non-diagonal conductivity is equal to 1/2. The weaker is the magnetic field the harder it is to locate a delocalization transition using quantum simulations. The main idea of the present study is that the position of the transition does not change when a strong local inhomogeneity is introduced. Since the strong inhomogeneity suppresses interference, transport reduces to classical percolation. We show that the corresponding percolation problem is bond percolation over two sublattices coupled to each other by random bonds. Simulation of this percolation allows to access the domain of very weak magnetic fields. Simulation results confirm the criterion σ_{xy}=1/2 for values σ_{xx}\sim 10, where they agree with earlier quantum simulation results. However for larger σ_{xx} we find that the transition boundary is described by σ_{xy} σ_{xx}^k with k= 0.5, i.e., the transition takes place at higher magnetic fields. The strong inhomogeneity limit of magnetotransport in the presence of a random magnetic field, pertinent to composite fermions, corresponds to a different percolation problem. In this limit we find for the delocalization transition boundary σ_{xy} σ_{xx}^{0.6}.

preprint2011arXiv

Spin quantum Hall effect and plateau transitions in multilayer network models

We study the spin quantum Hall effect and transitions between Hall plateaus in quasi two-dimensional network models consisting of several coupled layers. Systems exhibiting the spin quantum Hall effect belong to class C in the symmetry classification for Anderson localisation, and for network models in this class there is an established mapping between the quantum problem and a classical one involving random walks. This mapping permits numerical studies of plateau transitions in much larger samples than for other symmetry classes, and we use it to examine localisation in systems consisting of $n$ weakly coupled layers. Standard scaling ideas lead one to expect $n$ distinct plateau transitions, but in the case of the unitary symmetry class this conclusion has been questioned. Focussing on a two-layer model, we demonstrate that there are two separate plateau transitions, with the same critical properties as in a single-layer model, even for very weak interlayer coupling.

preprint2009arXiv

Crossover from diffusive to strongly localized regime in two-dimensional systems

We have studied the conductance distribution function of two-dimensional disordered noninteracting systems in the crossover regime between the diffusive and the localized phases. The distribution is entirely determined by the mean conductance, g, in agreement with the strong version of the single-parameter scaling hypothesis. The distribution seems to change drastically at a critical value very close to one. For conductances larger than this critical value, the distribution is roughly Gaussian while for smaller values it resembles a log-normal distribution. The two distributions match at the critical point with an often appreciable change in behavior. This matching implies a jump in the first derivative of the distribution which does not seem to disappear as system size increases. We have also studied 1/g corrections to the skewness to quantify the deviation of the distribution from a Gaussian function in the diffusive regime.

preprint2009arXiv

Random Walks and Anderson Localisation in a Three-Dimensional Class C Network Model

We study the disorder-induced localisation transition in a three-dimensional network model that belongs to symmetry class C. The model represents quasiparticle dynamics in a gapless spin-singlet superconductor without time-reversal invariance. It is a special feature of network models with this symmetry that the conductance and density of states can be expressed as averages in a classical system of dense, interacting random walks. Using this mapping, we present a more precise numerical study of critical behaviour at an Anderson transition than has been possible previously in any context.