bit circus, copenhagen, 23-24 of august 2012
DESCRIPTION
On the dynamics of the Fermi-Bose model Magnus Ögren Nano-Science Center, Copenhagen University. DTU-Mathematics, Technical University of Denmark. In collaboration with Marcus Carlsson Center for Mathematical Sciences, Lund University. BIT Circus, Copenhagen, 23-24 of August 2012. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/1.jpg)
On the dynamics of the Fermi-Bose model
Magnus Ögren Nano-Science Center, Copenhagen University.
DTU-Mathematics, Technical University of Denmark.
In collaboration with Marcus CarlssonCenter for Mathematical Sciences, Lund University.
BIT Circus, Copenhagen, 23-24 of August 2012
![Page 2: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/2.jpg)
![Page 3: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/3.jpg)
Background: What is the problem we would like to say something about using numerical calculations?Quantum dynamics of molecular BEC dissociation!
Formulation: How can we write up the dynamicalevolution of the system in differential equations?Linear ODEs for operators, evolve a complex matrix!
Improvements: What have we done to be able to treatlarge (i.e. realistic) arbitrary shaped 3D systems?Symmetries for block-matrices, D-block-Hankel matrix!
Some numerical results!
Outline of my talk:
![Page 4: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/4.jpg)
i) Conceptual:Molecular dissociation as a fermionic analog of optical parametric down-conversion, a good candidate for developing the paradigm of fermionic quantum atom optics
in fundamental physics and a test bench for simulations.
ii) Pragmatic:Can we explain the experimentally observed pair-correlations. (Molecules made up of fermions have longer lifetime.)
Motivation to study dissociation into fermions:
dimers fermions
![Page 5: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/5.jpg)
Fermi-Bose Hamiltonian and applications
![Page 6: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/6.jpg)
Exact simulation of molecular dissociation: MÖ, KK, JC, E.P.L. 2010.
We have earlier applied the Gaussian phase-space representation to stochastically model a 1D uniform molecular BEC dissociating into fermionic atoms.
![Page 7: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/7.jpg)
Implement a molecular-field approximation
Linear operator equations!
![Page 8: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/8.jpg)
Fourier transformation to momentum-space
Represent the BEC geometry with a D-dimensional Fourier series.
![Page 9: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/9.jpg)
Linear ODEs for momentum-space operators
Fourier coefficients are delta spikes for uniform systems.
![Page 10: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/10.jpg)
Uniform (even and real) condensate wavefunction
Connects to alternative formulation “PMFT”, but this require two indices per unknown for non-uniform systems.
![Page 11: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/11.jpg)
Uniform (even and real) condensate wavefunction
Valuable with analytic solutions for software tests!
![Page 12: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/12.jpg)
General formulation for a complex BEC wavefunction
![Page 13: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/13.jpg)
General formulation for a complex BEC wavefunction
![Page 14: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/14.jpg)
Theory 1: Define all necessary physical observables in terms of pairs of raws of the matrixexponential.Numerics 1: Use efficient software (expokit) for the calculation of only these raws from a sparse (truncated) system matrix.
Theory 2: Prove block-matrix symmetries.Numerics 2: Find block-matrix symmetries and implement them in the corresponding algorithms.
Theory 3: Define a D-block-Hankel matrix structure.Numerics 3: Implement algorithm for multiplication between a D-block-Hankel matrix and a vector and incorporate them into efficient matrixexponentiation software (expokit).
Major 3 steps towards a realistic 3D simulation:
![Page 15: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/15.jpg)
Theory 1: Define all necessary physical observables in terms of pairs of raws of the matrixexponential.
Numerics 1: Use efficient software (expokit) for the calculation of theses raws from a sparse (truncation) system matrix.
1.st step towards a realistic 3D simulation:
![Page 16: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/16.jpg)
What do we need to calculate?
![Page 17: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/17.jpg)
Physical observables are formed by pairs of raws
![Page 18: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/18.jpg)
Any observable is available (Wick approximated)
![Page 19: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/19.jpg)
Theory 2: Prove block-matrix symmetries.
Numerics 2: Find block-matrix symmetries and implement them in the corresponding algorithms.
2.nd step towards a realistic 3D simulation:
![Page 20: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/20.jpg)
General formulation for a complex BEC wavefunction
![Page 21: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/21.jpg)
Real and even BEC wavefunctions
![Page 22: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/22.jpg)
From symmetries in the system matrix to the observables
![Page 23: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/23.jpg)
From symmetries in the system matrix to the propagator
![Page 24: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/24.jpg)
From symmetries in the system matrix to the propagator
![Page 25: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/25.jpg)
Real BEC wavefunction
![Page 26: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/26.jpg)
Real and even BEC wavefunction (common in exp.)
![Page 27: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/27.jpg)
Theory 3: Define a D-block-Hankel matrix structure.
Numerics 3: Implement algorithm for multiplication between a D-block-Hankel matrix and a vector and incorporate them into efficient matrixexponentiation software (expokit).
3.rd step towards a realistic 3D simulation:
![Page 28: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/28.jpg)
General formulation for a complex BEC wavefunction
![Page 29: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/29.jpg)
Visualization of a D-block-Hankel matrix (D=3, K=30)
![Page 30: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/30.jpg)
Visualization of a D-block-Hankel matrix (D=3, K=30)
![Page 31: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/31.jpg)
Visualization of a D-block-Hankel matrix (D=3, K=30)
![Page 32: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/32.jpg)
Numerical results for fermionic atom-atom correlations
![Page 33: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/33.jpg)
Numerical evaluations of analytic asymptotes
![Page 34: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/34.jpg)
Collinear (CL) correlations, molecular dissociation
(b) Collinear (CL) correlations due to particle statistics, (like Hanbury Brown and Twiss for photons).
We have derived an analytical asymptote (dashed lines), strictly valid for short times (t/t0<<1). But useful even for t/t0~1 as here. Solid lines are from a numerical calculation at t/t0=0.5.
+-
![Page 35: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/35.jpg)
Observations from the field of ultra-cold atoms:
T. Jeltes et al., Nature445 (2007) 402.
See also: M. Henny et al., Science 284, 296 (1999). For ‘anti-bunching of electrons’ in a solid state device.
(CL) gj,j(2)(k,k’,t), j=1,2
Bosons
Fermions
![Page 36: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/36.jpg)
First 3D calculation for general BEC wavefunction
![Page 37: BIT Circus, Copenhagen, 23-24 of August 2012](https://reader036.vdocument.in/reader036/viewer/2022062409/568145fd550346895db30bb8/html5/thumbnails/37.jpg)
Related work:
On the dynamics of the Fermi-Bose model.M. Ögren and M. Carlsson, To be submitted to J. Phys. A: Math. Gen. 2012.
Stochastic simulations of fermionic dynamics with phase-space representations. M. Ögren, K. V. Kheruntsyan and J. F. Corney, Comp. Phys. Comm. 182 1999 (2011).
First-principles quantum dynamics for fermions: application to molecular dissociation.M. Ögren, K. V. Kheruntsyan and J. F. Corney, Europhys. Lett. 92, (2010) 36003.
Role of spatial inhomogeneity in dissociation of trapped molecular condensates. M. Ögren and K. V. Kheruntsyan, Phys. Rev. A 82, 013641 (2010). Directional effects due to quantum statistics in dissociation of elongated molecular condensates. M. Ögren, C. M. Savage and K. V. Kheruntsyan, Phys. Rev. A 79, 043624 (2009).
Atom-atom correlations from condensate collisions.M. Ögren and K. V. Kheruntsyan, Phys. Rev. A 79, 021606(R) (2009).
Atom-atom correlations and relative number squeezing in dissociation of spatially inhomogeneous molecular condensates.M. Ögren and K. V. Kheruntsyan, Phys. Rev. A 78, 011602(R) (2008).