TY - JOUR

T1 - Reduction-based strategy for optimal control of Bose-Einstein condensates

AU - Adriazola, J.

AU - Goodman, R. H.

N1 - Publisher Copyright:
© 2022 American Physical Society.

PY - 2022/2/1

Y1 - 2022/2/1

N2 - Applications of Bose-Einstein condensates (BEC) often require that the condensate be prepared in a specific complex state. Optimal control is a reliable framework to prepare such a state while avoiding undesirable excitations, and, when applied to the time-dependent Gross-Pitaevskii equation (GPE) model of BEC in multiple space dimensions, results in a large computational problem. We propose a control method based on first reducing the problem, using a Galerkin expansion, from a partial differential equation to a low-dimensional Hamiltonian ordinary differential equation system. We then apply a two-stage hybrid control strategy. At the first stage, we approximate the control using a second Galerkin-like method known as the chopped random basis to derive a finite-dimensional nonlinear programing problem, which we solve with a differential evolution algorithm. This search method then yields a candidate local minimum which we further refine using a variant of gradient descent. This hybrid strategy allows us to greatly reduce excitations both in the reduced model and the full GPE system.

AB - Applications of Bose-Einstein condensates (BEC) often require that the condensate be prepared in a specific complex state. Optimal control is a reliable framework to prepare such a state while avoiding undesirable excitations, and, when applied to the time-dependent Gross-Pitaevskii equation (GPE) model of BEC in multiple space dimensions, results in a large computational problem. We propose a control method based on first reducing the problem, using a Galerkin expansion, from a partial differential equation to a low-dimensional Hamiltonian ordinary differential equation system. We then apply a two-stage hybrid control strategy. At the first stage, we approximate the control using a second Galerkin-like method known as the chopped random basis to derive a finite-dimensional nonlinear programing problem, which we solve with a differential evolution algorithm. This search method then yields a candidate local minimum which we further refine using a variant of gradient descent. This hybrid strategy allows us to greatly reduce excitations both in the reduced model and the full GPE system.

UR - http://www.scopus.com/inward/record.url?scp=85125593164&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85125593164&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.105.025311

DO - 10.1103/PhysRevE.105.025311

M3 - Article

C2 - 35291192

AN - SCOPUS:85125593164

VL - 105

JO - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

JF - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

SN - 1063-651X

IS - 2

M1 - 025311

ER -