@article{SeroussiNowickiPayneetal.2020, author = {Seroussi, Helene and Nowicki, Sophie and Payne, Antony J. and Goelzer, Heiko and Lipscomb, William H. and Abe-Ouchi, Ayako and Agosta, Cecile and Albrecht, Torsten and Asay-Davis, Xylar and Barthel, Alice and Calov, Reinhard and Cullather, Richard and Dumas, Christophe and Galton-Fenzi, Benjamin K. and Gladstone, Rupert and Golledge, Nicholas R. and Gregory, Jonathan M. and Greve, Ralf and Hattermann, Tore and Hoffman, Matthew J. and Humbert, Angelika and Huybrechts, Philippe and Jourdain, Nicolas C. and Kleiner, Thomas and Larour, Eric and Leguy, Gunter R. and Lowry, Daniel P. and Little, Chistopher M. and Morlighem, Mathieu and Pattyn, Frank and Pelle, Tyler and Price, Stephen F. and Quiquet, Aurelien and Reese, Ronja and Schlegel, Nicole-Jeanne and Shepherd, Andrew and Simon, Erika and Smith, Robin S. and Straneo, Fiammetta and Sun, Sainan and Trusel, Luke D. and Van Breedam, Jonas and van de Wal, Roderik S. W. and Winkelmann, Ricarda and Zhao, Chen and Zhang, Tong and Zwinger, Thomas}, title = {ISMIP6 Antarctica}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3033-2020}, pages = {3033 -- 3070}, year = {2020}, abstract = {Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015-2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between 7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between 6 :1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.}, language = {en} } @article{BarthelPinedaEisert2009, author = {Barthel, Thomas and Pineda, Carlos and Eisert, Jens}, title = {Contraction of fermionic operator circuits and the simulation of strongly correlated fermions}, issn = {1050-2947}, doi = {10.1103/Physreva.80.042333}, year = {2009}, abstract = {A fermionic operator circuit is a product of fermionic operators of usually different and partially overlapping support. Further elements of fermionic operator circuits (FOCs) are partial traces and partial projections. The presented framework allows for the introduction of fermionic versions of known qudit operator circuits (QUOC), important for the simulation of strongly correlated d-dimensional systems: the multiscale entanglement renormalization ansaumltze (MERA), tree tensor networks (TTN), projected entangled pair states (PEPS), or their infinite-size versions (iPEPS etc.). After the definition of a FOC, we present a method to contract it with the same computation and memory requirements as a corresponding QUOC, for which all fermionic operators are replaced by qudit operators of identical dimension. A given scheme for contracting the QUOC relates to an analogous scheme for the corresponding fermionic circuit, where additional marginal computational costs arise only from reordering of modes for operators occurring in intermediate stages of the contraction. Our result hence generalizes efficient schemes for the simulation of d- dimensional spin systems, as MERA, TTN, or PEPS to the fermionic case.}, language = {en} } @article{BarthelSchollwockWhite2009, author = {Barthel, Thomas and Schollwock, Ulrich and White, Steven R.}, title = {Spectral functions in one-dimensional quantum systems at finite temperature using the density matrix renormalization group}, issn = {1098-0121}, doi = {10.1103/Physrevb.79.245101}, year = {2009}, abstract = {We present time-dependent density matrix renormalization group simulations (t-DMRG) at finite temperatures. It is demonstrated how a combination of finite-temperature t-DMRG and time-series prediction allows for an easy and very accurate calculation of spectral functions in one-dimensional quantum systems, irrespective of their statistics for arbitrary temperatures. This is illustrated with spin structure factors of XX and XXX spin-1/2 chains. For the XX model we can compare against an exact solution, and for the XXX model (Heisenberg antiferromagnet) against a Bethe ansatz solution and quantum Monte Carlo data.}, language = {en} } @article{BarthelKlieschEisert2010, author = {Barthel, Thomas and Kliesch, Martin and Eisert, Jens}, title = {Real-space renormalization yields finite correlations}, issn = {0031-9007}, doi = {10.1103/Physrevlett.105.010502}, year = {2010}, abstract = {Real-space renormalization approaches for quantum lattice systems generate certain hierarchical classes of states that are subsumed by the multiscale entanglement renormalization Ansatz (MERA). It is shown that, with the exception of one spatial dimension, MERA states are actually states with finite correlations, i.e., projected entangled pair states (PEPS) with a bond dimension independent of the system size. Hence, real-space renormalization generates states which can be encoded with local effective degrees of freedom, and MERA states form an efficiently contractible class of PEPS that obey the area law for the entanglement entropy. It is further pointed out that there exist other efficiently contractible schemes violating the area law.}, language = {en} } @article{KlieschBarthelGogolinetal.2011, author = {Kliesch, Martin and Barthel, Thomas and Gogolin, C. and Kastoryano, M. and Eisert, J.}, title = {Dissipative quantum church-turing theorem}, series = {Physical review letters}, volume = {107}, journal = {Physical review letters}, number = {12}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.107.120501}, pages = {5}, year = {2011}, abstract = {We show that the time evolution of an open quantum system, described by a possibly time dependent Liouvillian, can be simulated by a unitary quantum circuit of a size scaling polynomially in the simulation time and the size of the system. An immediate consequence is that dissipative quantum computing is no more powerful than the unitary circuit model. Our result can be seen as a dissipative Church-Turing theorem, since it implies that under natural assumptions, such as weak coupling to an environment, the dynamics of an open quantum system can be simulated efficiently on a quantum computer. Formally, we introduce a Trotter decomposition for Liouvillian dynamics and give explicit error bounds. This constitutes a practical tool for numerical simulations, e.g., using matrix-product operators. We also demonstrate that most quantum states cannot be prepared efficiently.}, language = {en} } @article{BarthelKliesch2012, author = {Barthel, Thomas and Kliesch, Martin}, title = {Quasilocality and efficient simulation of Markovian Quantum Dynamics}, series = {Physical review letters}, volume = {108}, journal = {Physical review letters}, number = {23}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.108.230504}, pages = {5}, year = {2012}, abstract = {We consider open many-body systems governed by a time-dependent quantum master equation with short-range interactions. With a generalized Lieb-Robinson bound, we show that the evolution in this very generic framework is quasilocal; i.e., the evolution of observables can be approximated by implementing the dynamics only in a vicinity of the observables' support. The precision increases exponentially with the diameter of the considered subsystem. Hence, time evolution can be simulated on classical computers with a cost that is independent of the system size. Providing error bounds for Trotter decompositions, we conclude that the simulation on a quantum computer is additionally efficient in time. For experiments and simulations in the Schrodinger picture, our result can be used to rigorously bound finite-size effects.}, language = {en} }