Exact relaxation in a class of non-equilibrium quantum lattice systems

A reasonable physical intuition in the study of interacting quantum systems says that, independent of the initial state, the system will tend to equilibrate. In this work we study a setting where relaxation to a steady state is exact, namely for the Bose-Hubbard model where the system is quenched from a Mott quantum phase to the strong superfluid regime. We find that the evolving state locally relaxes to a steady state with maximum entropy constrained by second moments, maximizing the entanglement, to a state which is different from the thermal state of the new Hamiltonian. Remarkably, in the infinite system limit this relaxation is true for all large times, and no time average is necessary. For large but finite system size we give a time interval for which the system locally "looks relaxed" up to a prescribed error. Our argument includes a central limit theorem for harmonic systems and exploits the finite speed of sound. Additionally, we show that for all periodic initial configurations, reminiscent of charge density waves, the system relaxes locally. We sketch experimentally accessible signatures in optical lattices as well as implications for the foundations of quantum statistical mechanics.Comment: 8 pages, 3 figures, replaced with final versio

On-lattice agent-based simulation of populations of cells within the open-source chaste framework

Over the years, agent-based models have been developed that combine cell division and reinforced random walks of cells on a regular lattice, reaction-diffusion equations for nutrients and growth factors and ordinary differential equations (ODEs) for the subcellular networks regulating the cell cycle. When linked to a vascular layer, this multiple scale model framework has been applied to tumour growth and therapy. Here we report on the creation of an agent-based multiscale environment amalgamating the characteristics of these models within a Virtual Pysiological Human (VPH) Exemplar Project. This project enables re-use, integration, expansion and sharing of the model and relevant data. The agent-based and reactiondiffusion parts of the multiscale model have been implemented and are available for download as part of the latest public release of Chaste (“Cancer, Heart and Soft Tissue Environment”), (http://www.cs.ox.ac.uk/chaste/) version 3.1, part of the VPH Toolkit (http://toolkit.vph-noe.eu/). The environment functionalities are verified against the original models, in addition to extra validation of all aspects of the code. In this work, we present the details of the implementation of the agent-based environment, including the system description, the conceptual model, the development of the simulation model and the processes of verification and validation of the simulation results. We explore the potential use of the environment by presenting exemplar applications of the “what if” scenarios that can easily be studied in the environment. These examples relate to tumour growth, cellular competition for resources and tumour responses to hypoxia. We conclude our work by summarising the future steps for the expansion of the current system

Solving frustration-free spin systems

We identify a large class of quantum many-body systems that can be solved exactly: natural frustration-free spin-1/2 nearest-neighbor Hamiltonians on arbitrary lattices. We show that the entire ground state manifold of such models can be found exactly by a tensor network of isometries acting on a space locally isomorphic to the symmetric subspace. Thus, for this wide class of models real-space renormalization can be made exact. Our findings also imply that every such frustration-free spin model satisfies an area law for the entanglement entropy of the ground state, establishing a novel large class of models for which an area law is known. Finally, we show that our approach gives rise to an ansatz class useful for the simulation of almost frustration-free models in a simple fashion, outperforming mean field theory.Comment: 5 pages, 1 figur

A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy

Transmission spectroscopy, which consists of measuring the wavelength-dependent absorption of starlight by a planet's atmosphere during a transit, is a powerful probe of atmospheric composition. However, the expected signal is typically orders of magnitude smaller than instrumental systematics, and the results are crucially dependent on the treatment of the latter. In this paper, we propose a new method to infer transit parameters in the presence of systematic noise using Gaussian processes, a technique widely used in the machine learning community for Bayesian regression and classification problems. Our method makes use of auxiliary information about the state of the instrument, but does so in a non-parametric manner, without imposing a specific dependence of the systematics on the instrumental parameters, and naturally allows for the correlated nature of the noise. We give an example application of the method to archival NICMOS transmission spectroscopy of the hot Jupiter HD 189733, which goes some way towards reconciling the controversy surrounding this dataset in the literature. Finally, we provide an appendix giving a general introduction to Gaussian processes for regression, in order to encourage their application to a wider range of problems.Comment: 6 figures, 1 table, accepted for publication in MNRA

The ground state of a class of noncritical 1D quantum spin systems can be approximated efficiently

We study families H_n of 1D quantum spin systems, where n is the number of spins, which have a spectral gap \Delta E between the ground-state and first-excited state energy that scales, asymptotically, as a constant in n. We show that if the ground state |\Omega_m> of the hamiltonian H_m on m spins, where m is an O(1) constant, is locally the same as the ground state |\Omega_n>, for arbitrarily large n, then an arbitrarily good approximation to the ground state of H_n can be stored efficiently for all n. We formulate a conjecture that, if true, would imply our result applies to all noncritical 1D spin systems. We also include an appendix on quasi-adiabatic evolutions.Comment: 9 pages, 1 eps figure, minor change

Evidence for a climate signal in trends of global crop yield variability over the past 50 years

Low variability of crop production from year to year is desirable for many reasons, including reduced income risk and stability of supplies. Therefore, it is important to understand the nature of yield variability, whether it is changing through time, and how it varies between crops and regions. Previous studies have shown that national crop yield variability has changed in the past, with the direction and magnitude dependent on crop type and location. Whilst such studies acknowledge the importance of climate variability in determining yield variability, it has been assumed that its magnitude and its effect on crop production have not changed through time and, hence, that changes to yield variability have been due to non-climatic factors. We address this assumption by jointly examining yield and climate variability for three major crops (rice, wheat and maize) over the past 50 years. National yield time series and growing season temperature and precipitation were de-trended and related using multiple linear regression. Yield variability changed significantly in half of the crop–country combinations examined. For several crop–country combinations, changes in yield variability were related to changes in climate variability

Quasi-Adiabatic Continuation in Gapped Spin and Fermion Systems: Goldstone's Theorem and Flux Periodicity

We apply the technique of quasi-adiabatic continuation to study systems with continuous symmetries. We first derive a general form of Goldstone's theorem applicable to gapped nonrelativistic systems with continuous symmetries. We then show that for a fermionic system with a spin gap, it is possible to insert $\pi$-flux into a cylinder with only exponentially small change in the energy of the system, a scenario which covers several physically interesting cases such as an s-wave superconductor or a resonating valence bond state.Comment: 19 pages, 2 figures, final version in press at JSTA

Topology and Phases in Fermionic Systems

There can exist topological obstructions to continuously deforming a gapped Hamiltonian for free fermions into a trivial form without closing the gap. These topological obstructions are closely related to obstructions to the existence of exponentially localized Wannier functions. We show that by taking two copies of a gapped, free fermionic system with complex conjugate Hamiltonians, it is always possible to overcome these obstructions. This allows us to write the ground state in matrix product form using Grassman-valued bond variables, and show insensitivity of the ground state density matrix to boundary conditions.Comment: 4 pages, see also arxiv:0710.329