We determine analytically the phase diagram of the toric code model in a parallel magnetic field which displays three distinct regions. Our study relies on two high-order perturbative expansions in the strong- and weak-field limit, as well as a large-spin analysis. Calculations in the topological phase establish a quasiparticle picture for the anyonic excitations. We obtain two second-order transition lines that merge with a first-order line giving rise to a multicritical point as recently suggested by numerical simulations. We compute the values of the corresponding critical fields and exponents that drive the closure of the gap. We also give the one-particle dispersions of the anyonic quasiparticles inside the topological phase.
The Polchinski version of the exact renormalisation group equations is applied to multicritical fixed points, which are present for dimensions between two and four, for scalar theories using both the local potential approximation and its extension, the derivative expansion. The results are compared with the epsilon expansion by showing that the non linear differential equations may be linearised at each multicritical point and the epsilon expansion treated as a perturbative expansion. The results for critical exponents are compared with corresponding epsilon expansion results from standard perturbation theory. The results provide a test for the validity of the local potential approximation and also the derivative expansion. An alternative truncation of the exact RG equation leads to equations which are similar to those found in the derivative expansion but which gives correct results for critical exponents to order $epsilon$ and also for the field anomalous dimension to order $epsilon^2$. An exact marginal operator for the full RG equations is also constructed.
On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.
Applicability of classical Lifshitz-Slyozov theory of Ostwald ripening is analyzed and found limited by relatively large cluster sizes due to restrictions imposed by theoretical assumptions. An assumption about the steady state ripening regime poses an upper limit, while another, implicit assumption of continuous description poses a cluster size-dependent lower limit on the supersaturation level. These two limits mismatch for the clusters under certain size in the nanometer scale making the theory inapplicable. We present a more generic, molecular theory of Ostwald ripening, which reproduces classical Lifshitz-Slyozov and Wagner theories in appropriate extreme cases. This theory has a wider applicability than classical theories, especially at lower supersaturation levels, and is more suitable for nanoscale systems.
A non-Hermitian complex symmetric 2x2 matrix toy model is used to study projective Hilbert space structures in the vicinity of exceptional points (EPs). The bi-orthogonal eigenvectors of a diagonalizable matrix are Puiseux-expanded in terms of the root vectors at the EP. It is shown that the apparent contradiction between the two incompatible normalization conditions with finite and singular behavior in the EP-limit can be resolved by projectively extending the original Hilbert space. The complementary normalization conditions correspond then to two different affine charts of this enlarged projective Hilbert space. Geometric phase and phase jump behavior are analyzed and the usefulness of the phase rigidity as measure for the distance to EP configurations is demonstrated. Finally, EP-related aspects of PT-symmetrically extended Quantum Mechanics are discussed and a conjecture concerning the quantum brachistochrone problem is formulated.
We present a formal demonstration that light can simultaneously exhibit a superfluid behavior and spatial long-range order when propagating in a photonic crystal with self-focussing nonlinearity. In this way, light presents the distinguishing features of matter in a supersolid phase. We show that this supersolid phase provides the stability conditions for nonlinear Bloch waves and, at the same time, permits the existence of topological solitons or defects for the envelope of these waves. We use a condensed matter analysis instead of a standard nonlinear optics approach and provide numerical evidence of these theoretical findings.
Investigations of the Fermi surface via the electron momentum distribution reconstructed from either angular correlation of annihilation radiation (or Compton scattering) experimental spectra are presented. The basis of these experiments and mathematical methods applied in reconstructing three-dimensional densities from line (or plane) projections measured in these experiments are described. The review of papers where such techniques have been applied to study the Fermi surface of metallic materials with showing their main results is also done.
A significant obstacle for practical quantum computation is the loss of physical qubits in quantum computers, a decoherence mechanism most notably in optical systems. Here we experimentally demonstrate, both in the quantum circuit model and in the one-way quantum computer model, the smallest non-trivial quantum codes to tackle this problem. In the experiment, we encode single-qubit input states into highly-entangled multiparticle codewords, and we test their ability to protect encoded quantum information from detected one-qubit loss error. Our results prove the in-principle feasibility of overcoming the qubit loss error by quantum codes.
Scatterings of electrons at quasiparticles or photons are very important for many topics in solid state physics, e.g., spintronics, magnonics or photonics, and therefore a correct numerical treatment of these scatterings is very important. For a quantum-mechanical description of these scatterings Fermis golden rule is used in order to calculate the transition rate from an initial state to a final state in a first-order time-dependent perturbation theory. One can calculate the total transition rate from all initial states to all final states with Boltzmann rate equations involving Brillouin zone integrations. The numerical treatment of these integrations on a finite grid is often done via a replacement of the Dirac delta distribution by a Gaussian. The Dirac delta distribution appears in Fermis golden rule where it describes the energy conservation among the interacting particles. Since the Dirac delta distribution is a not a function it is not clear from a mathematical point of view that this procedure is justified. We show with physical and mathematical arguments that this numerical procedure is in general correct, and we comment on critical points.
In a recent paper [Phys.~Rev.~A {bf 89}, 022503 (2014)], the average density approximation (ADA) was implemented to develop a parameter-free, nonlocal kinetic energy functional to be used in the orbital-free density-functional theory of an inhomogenous, two-dimensional (2D), Fermi gas. In this work, we provide a detailed comparison of self-consistent calculations within the ADA with the exact results of the Kohn-Sham density-functional theory, and the elementary Thomas-Fermi (TF) approximation. We demonstrate that the ADA for the 2D kinetic energy functional works very well under a wide variety of confinement potentials, even for relatively small particle numbers. Remarkably, the TF approximation for the kinetic energy functional, {em without any gradient corrections}, also yields good agreement with the exact kinetic energy for all confining potentials considered, although at the expense of the spatial and kinetic energy densities exhibiting poor point-wise agreement, particularly near the TF radius. Our findings illustrate that the ADA kinetic energy functional yields accurate results for {em both} the local and global equilibrium properties of an inhomogeneous 2D Fermi gas, without the need for any fitting parameters.
