Batas Kecepatan Lieb-Robinson yang Eksplisit

22

Batas Lieb-Robinson menggambarkan bagaimana efek disebarkan melalui sistem karena Hamiltonian lokal. Mereka sering digambarkan dalam bentuk mana dan adalah operator yang dipisahkan oleh jarak pada kisi di mana Hamiltonian memiliki lokal (misalnya tetangga terdekat) interaksi pada kisi itu, dibatasi oleh beberapa kekuatan . The bukti dari Lieb Robinson terikat biasanya menunjukkan adanya kecepatan (yang tergantung pada ). Ini sering sangat berguna untuk membatasi properti dalam sistem ini. Sebagai contoh, ada beberapa hasil yang sangat bagus di sini

|[A,B(t)]|Cevtl,
ABlJvJ tentang berapa lama yang dibutuhkan untuk menghasilkan status GHZ menggunakan Hamiltonian tetangga terdekat.

Masalah yang saya punya adalah bahwa bukti-bukti yang cukup generik yang sulit untuk mendapatkan nilai yang ketat pada apa kecepatan sebenarnya adalah untuk setiap sistem yang diberikan.

Untuk lebih spesifik, bayangkan rantai qubit satu dimensi ditambah dengan Hamiltonian mana untuk semua . Di sini , Y_n dan Z_n mewakili operator Pauli yang diterapkan ke qubit n yang diberikan , dan \ mathbb {I} di tempat lain. Bisakah Anda memberikan batas atas yang baik (yaitu sekencang mungkin) untuk kecepatan Lieb-Robinson v untuk sistem dalam Persamaan. (1)

(1)H=n=1NBn2Zn+n=1N1Jn2(XnXn+1+YnYn+1),
JnJnXnYnZnnIv

Pertanyaan ini dapat ditanyakan dengan dua asumsi berbeda:

  • The Jn dan Bn semua tetap dalam waktu
  • The Jn dan Bn diperbolehkan untuk bervariasi dalam waktu.

Yang pertama adalah asumsi yang lebih kuat yang dapat membuat bukti lebih mudah, sedangkan yang terakhir biasanya dimasukkan dalam pernyataan batas Lieb-Robinson.


Motivasi

Komputasi kuantum, dan informasi kuantum yang lebih umum, muncul untuk membuat keadaan kuantum yang menarik. Melalui karya-karya seperti ini , kita melihat bahwa informasi membutuhkan sejumlah waktu untuk menyebar dari satu tempat ke tempat lain dalam sistem kuantum yang mengalami evolusi karena orang Hamilton seperti Persamaan. (1), dan status kuantum, seperti status GHZ, atau status dengan urutan topologi, membutuhkan waktu tertentu untuk diproduksi. Apa yang ditunjukkan oleh hasil saat ini adalah relasi penskalaan, misalnya waktu yang diperlukan adalah Ω(N) .

Jadi, katakanlah saya datang dengan skema yang tidak mentransfer informasi, atau menghasilkan GHZ negara dll dengan cara yang skala linear di . Seberapa baik skema itu sebenarnya? Jika saya memiliki kecepatan eksplisit, saya dapat melihat seberapa dekat koefisien penskalaan dalam skema saya dibandingkan dengan batas bawah.N

Jika saya berpikir bahwa suatu hari yang ingin saya lihat adalah protokol yang diterapkan di lab, maka saya sangat peduli untuk mengoptimalkan koefisien penskalaan ini, bukan hanya fungsi penskalaan yang luas, karena semakin cepat saya dapat mengimplementasikan protokol, semakin sedikit kesempatan di sana. adalah untuk kebisingan datang dan mengacaukan semuanya.


Informasi lebih lanjut

Ada beberapa fitur bagus dari Hamiltonian ini yang saya anggap membuat perhitungan lebih mudah. Secara khusus, Hamiltonian memiliki struktur ruang bagian berdasarkan jumlah 1s dalam standar dasar (dikatakan mempertahankan eksitasi) dan, lebih baik lagi, transformasi Jordan-Wigner menunjukkan bahwa semua sifat dari ruang bagian eksitasi yang lebih tinggi dapat diturunkan. dari subruang 1-eksitasi.N×Nh2N×2NHh = N n = 1 B n | n n | + N - 1 Σ n = 1 J n ( | n n + 1 | + | n + 1 n | ) . v = 2 J 2 J, di mana Ada beberapa bukti bahwa kecepatan Lieb-Robinson adalah , seperti di sini dan di sini , tetapi ini semua menggunakan rantai yang mendekati gabungan, yang memiliki kecepatan grup (dan saya berasumsi bahwa kecepatan grup terkait erat dengan Kecepatan Lieb-Robinson). Itu tidak membuktikan bahwa semua pilihan kekuatan kopling yang memungkinkan memiliki kecepatan yang sangat dibatasi.

h=n=1NBn|nn|+n=1N1Jn(|nn+1|+|n+1n|).
v=2J2J

Saya dapat menambahkan sedikit lebih jauh pada motivasi. Pertimbangkan evolusi waktu dari satu eksitasi tunggal yang dimulai pada satu ujung rantai, , dan berapa amplitudo untuk tiba di ujung rantai yang lain , waktu yang singkat nanti. Untuk urutan pertama dalam , ini adalah Anda dapat melihat fungsionalitas eksponensial yang Anda harapkan berada di luar 'kerucut cahaya' yang ditentukan oleh sistem Lieb-Robinson, tetapi yang lebih penting, jika Anda ingin memaksimalkan amplitudo itu, Anda akan mengatur semua|1|Nδtδt

N|eihδt|1=δtN1(N1)!n=1N1Jn+O(δtN).
Jn=J. Jadi, dalam waktu singkat, sistem yang digabungkan secara seragam mengarah ke transfer paling cepat. Mencoba untuk mendorong ini lebih jauh, Anda dapat bertanya, sebagai sedikit kebohongan, kapan bisa Mengambil batas besar , dan menggunakan rumus Stirling pada faktorial mengarah ke yang menunjukkan kecepatan maksimum sekitar . Tutup, tetapi tidak terlalu ketat (karena persyaratan pesanan tinggi tidak dapat diabaikan)!
tN1(N1)!n=1N1Jn1
N
etJN11,
eJ

DaftWullie
sumber
Sudahkah Anda menghitung LR-terikat terbaik dari bukti untuk model itu? Bagaimana perbandingannya dengan kecepatan yang Anda kutip?
Norbert Schuch
1
Ok, saya akui itu adalah pertanyaan komputasi kuantum, setidaknya cara saya menafsirkannya sekarang: "Apa pilihan dan (tergantung pada beberapa kendala) yang menghasilkan kecepatan maksimum untuk transfer informasi / keadaan / .... " --- Apakah ini interpretasi yang benar? B nJnBn
Norbert Schuch
@NorbertSchuch Tidak cukup. Saya ingin dapat mengatakan, "Saya telah datang dengan satu set sambungan yang mencapai protokol dengan penskalaan tertentu. Protokol itu diketahui dibatasi oleh batas Lieb-Robinson. Seberapa dekat saya dengan menjenuhkan kendala itu?" sebagai ukuran seberapa cepat protokol saya.
DaftWullie
@DaftWullie Jadi - apakah Anda mempertanyakan: "Seberapa dekat saya dengan optimal", atau "Seberapa dekat saya dengan semacam ikatan (mengambil yang seketat mungkin)"?
Norbert Schuch
1
@ user1271772 Itu benar. B(t)=eiHtB(0)eiHt
DaftWullie

Jawaban:

4

Biarkan saya menjawab pertanyaan umum bagaimana cara mendapatkan kecepatan Lieb-Robinson (LR) yang cukup ketat ketika Anda menghadapi model kisi yang berinteraksi secara lokal secara generik, dan kemudian saya akan kembali ke model 1D XY dalam pertanyaan Anda, yang sangat khusus untuk dipecahkan secara tepat.


Metode umum

Metode untuk mendapatkan terikat paling ketat saat ini (untuk model interaksi jarak pendek generik) diperkenalkan dalam Ref1 = arXiv: 1908.03997 . Gagasan dasarnya adalah norma komutator waktu yang tidak merataantara operator lokal sewenang-wenang dapat dibatasi oleh solusi untuk satu set persamaan diferensial linear orde pertama yang hidup pada grafik komutatifitas model. Grafik komutatif, seperti yang diperkenalkan pada Sec.II A dari Ref1, dapat dengan mudah diambil dari model Hamiltonian , dan dirancang untuk mencerminkan hubungan pergantian antara berbagai operator lokal yang disajikan dalam[AX(t),BY(0)]H H ω max ( i κH^H^. Dalam sistem invarian terjemahan, rangkaian persamaan diferensial ini dapat dengan mudah dipecahkan oleh transformasi Fourier, dan batas atas kecepatan LR dapat dihitung dari frekuensi eigenfrequency terbesar menggunakan Persamaan (31) dari Ref1 . Berikut ini saya akan menerapkan metode ini untuk model 1D XY sebagai contoh pedagogis. Untuk kesederhanaan, saya akan fokus pada case invarian-independen dan terjemahan , (batas yang dihasilkan tidak tergantung pada tanda-tandaωmax(iκ)| B n | = B > 0 | J n | = J > 0 B n , J n | J n ( t ) | J , | B n (|Bn|=B>0|Jn|=J>0Bn,Jn). Untuk terjemahan yang tidak berubah-ubah, bergantung pada waktu, Anda dapat menyelesaikan persamaan diferensial secara numerik (yang merupakan tugas komputasi yang mudah untuk sistem ribuan situs), atau Anda dapat menggunakan keseluruhan batas atas dan lanjutkan untuk menggunakan metode di bawah ini (tapi ini sedikit kompromi sesaknya dibandingkan dengan metode numerik).|Jn(t)|J, |Bn(t)|B

  1. Pertama-tama kita menggambar grafik komutatifitas, seperti di bawah ini. Setiap operator di Hamiltonian ~ ( , , ) diwakili oleh sebuah vertex, dan kami menautkan dua simpul jika dan hanya jika operator terkait tidak bepergian ( atau, dalam kasus saat ini, anti-perjalanan). XnXn+1YnYn+1Znmasukkan deskripsi gambar di sini

  2. Kemudian tuliskan persamaan diferensial Persamaan (10) dari Ref1 :

    γ¯˙α,n=J[γ¯α,n1(t)+γ¯α,n+1(t)]+B[γ¯3,n(t)+γ¯3,n+1(t)],  α=1,2,γ¯˙3,n=Jα=1,2[γ¯α,n1(t)+γ¯α,n(t)].

  3. Fourier mentransformasikan persamaan di atas, kita memiliki Frekuensi eigen adalah . Kecepatan LR diberikan oleh Pers. (31) dari Ref1 : mana

    ddt(γ¯1,kγ¯2,kγ¯3,k)=(2Jcosk0B(1+eik)02JcoskB(1+eik)J(1+eik)J(1+eik)0)(γ¯1,kγ¯2,kγ¯3,k).
    2Jcosk,Jcosk±(Jcosk)2+2BJ(1+cosk)
    vLRminκ>0ωmax(iκ)κ=ZB2JJ,
    Zyminκ>0coshκ+cosh2κ+4y(1+coshκ)κ.

Catatan: Ini terikat berbeda ketika , sedangkan kecepatan propagasi informasi fisik tetap terbatas. Kita bisa menyingkirkan masalah ini dengan menggunakan metode di Sec. VI dari Ref1 . Hasilnya adalah dalam batas ini, di mana didefinisikan sebagai solusi untuk persamaan .B/JvLR4X0JXyxarcsinh(x)=x2+1+y


Batas kecepatan untuk beberapa model klasik

Metode di atas sepenuhnya umum. Jika Anda tertarik lebih dalam, saya mencantumkan batas kecepatan untuk beberapa model klasik pada tabel berikut, yang diperoleh dengan cara yang sama. Perhatikan bahwa kecepatan LR dibatasi oleh yang terkecil dari semua ekspresi yang terdaftar (jadi di wilayah parameter yang berbeda, ekspresi yang berbeda harus digunakan). Fungsi didefinisikan sebagai root terbesarSemua parameter diasumsikan positif (ambil saja nilai absolut untuk kasus negatif).vLRF(Jx,Jy,Jz)x3(JxJy+JxJz+JyJz)x2JxJyJz=0.

ModelvLRd-dimensional TFIM2X0dJh=3.02dJhH^=JmnXmXn+hnZn4Xd1ddJ8.93dJ4X0dh=6.04dhd-dimensional Fermi-Hubbard2X3U4dJdJH^=Jmn,s=↑,(am,san,s+H.c.)8Xd1ddJ17.9dJ  +UnananananZU/JJ (d=1)1D Heisenberg XYZ4X0F(Jx,Jy,Jz)H^=n(JxXnXn+1+JyYnYn+1+JzZnZn+1)34.6max{Jx,Jy}

Adapun seberapa baik batas-batas ini, saya belum menyelidiki secara umum, tetapi untuk TFIM 1D pada titik kritis , solusi tepat memberikan , sedangkan batas di atas memberikan . Demikian pula, pada titik FH dan titik Heisenberg XYZ, batas di atas semuanya lebih besar dari solusi tepat dengan faktor . [Sebenarnya pada titik-titik khusus ini dua yang terakhir setara dengan rantai TFIM yang terpisah, seperti yang dapat dinilai secara langsung dari grafik komutatifitas mereka.]J=hvLR=2J2X0J3.02JU=0Jx=Jy,Jz=0X01.50888


Terikat ketat untuk 1D XY dengan memetakan ke fermion gratis

Sekarang mari kita bicara lebih banyak tentang model 1D XY. Seperti yang Anda perhatikan, ini dapat dipecahkan dengan memetakan fermion gratis: Untuk umum Anda perlu menyelesaikan masalah free-fermion secara numerik, tetapi izinkan saya menyebutkan dua kasus khusus yang secara analitik dapat ditelusuri.

H^=nBn(anan1/2)+nJn(anan+1+H.c.).
Bn(t),Jn(t)

  1. Bn(t)=B,Jn(t)=J adalah tetap dan terjemahan tidak berubah. Maka solusi yang tepat adalah mana adalah fungsi Bessel dari urutan. Jadi kecepatan LR adalah .

    an(t)=12πππa~kei2Jtcoskeikxdk=mJ|nm|(2Jt)am(0),
    J|nm|(2Jt)|nm|vLRXY=2J

  2. Bn,Jn diperbaiki dalam waktu tetapi benar-benar acak (gangguan padam). Kemudian karena lokalisasi banyak-tubuh (atau lokalisasi Anderson dalam gambar fermion), informasi tidak menyebar di sistem ini, jadi . Lebih tepatnya, di arXiv: quant-ph / 0703209 , batasan berikut ini terbukti untuk kasus tidak teratur: dengan kerucut cahaya logaritmik yang melambat .v LR = 0 [ A X (vLR=0

    [AX(t),BY(0)]const. t edXY/ξ,
    dXY=ξlnt

Lagrenge
sumber
Haruskah saya menyimpulkan dari apa yang Anda katakan bahwa untuk setiap model (termasuk yang tanpa terjemahan invarian) dengan , bahwa kecepatannya adalah ? XY|Jn|JvLRXY2J
DaftWullie
@ DaftWullie Tidak, Anda hanya dapat menggunakan batas atas keseluruhan untuk parameter dalam metode umum, karena metode umum selalu memberikan batas yang benar-benar tidak berkurang dalam nilai absolut dari koefisien apa pun. terikat diperoleh dari solusi eksak fermion bebas, di mana Anda tidak dapat menggunakan batas atas keseluruhan untuk parameter, dan harus menyelesaikan kasus per kasus. Jika adalah terjemahan invarian, maka Anda dapat mengatur dalam metode umum karena istilah bepergian dengan , dan dapatkan . B n ( t ) B = 0 B H v LR2 X 0 J = 3,02 J2JBn(t)B=0BH^vLR2X0J=3.02J
Lagrenge
@DaftWullie Dear DaftWullie, jika Anda berpikir ada sesuatu yang hilang dalam jawaban saya, atau poin apa pun masih belum jelas, tolong beri tahu saya.
Lagrenge
jawabannya terlihat berpotensi bermanfaat. Saya belum punya waktu untuk melihat kertas Anda (mungkin beberapa minggu). Dengan asumsi saya mengerti semuanya OK, itu intinya saya akan menerima jawaban Anda.
DaftWullie