BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Weizhu Bao

Weizhu Bao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAcond-mat.mtrl-scicond-mat.quant-gasNumerical Analysisquant-phmath-phmath.MPphysics.comp-phcond-mat.softmath.APphysics.flu-dyn

Catalog footprint

What is connected

25works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Diffuse-interface approach to competition between viscous flow and diffusion in pinch-off dynamics

The pinch-off dynamics of a liquid thread has been studied through numerical simulations and theoretical analysis. Occurring at small length scales, the pinch-off dynamics admits similarity solutions that can be classified into the Stokes regime and the diffusion-dominated regime, with the latter being recently experimentally observed in aqueous two-phase systems [Phys. Rev. Lett. 123, 134501 (2019)]. Derived by applying Onsager's variational principle, the Cahn-Hilliard-Navier-Stokes model is employed as a minimal model capable of describing the interfacial motion driven by not only advection but also diffusion. By analyzing the free energy dissipation mechanisms in the model, a characteristic length scale is introduced to measure the competition between diffusion and viscous flow in interfacial motion. This length scale is typically of nanometer scale for systems far from the critical point, but can approach micrometer scale for aqueous two-phase systems close to the critical point. The Cahn-Hilliard-Navier-Stokes model is solved by using an accurate and efficient spectral method in a cylindrical domain with axisymmetry. Ample numerical examples are presented to show the pinch-off processes in the Stokes regime and the diffusion-dominated regime respectively. In particular, the crossover between these two regimes is investigated numerically and analytically to reveal how the scaling behaviors of similarity solutions are to be qualitatively changed as the characteristic length scale is inevitably accessed by the pinching neck of the interface. Discussions are also provided for numerical examples that are performed for the breakup of long liquid filaments and show qualitatively different phenomena in different scaling regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved uniform error bounds of the time-splitting methods for the long-time (nonlinear) Schrödinger equation

We establish improved uniform error bounds for the time-splitting methods for the long-time dynamics of the Schrödinger equation with small potential and the nonlinear Schrödinger equation (NLSE) with weak nonlinearity. For the Schrödinger equation with small potential characterized by a dimensionless parameter $\varepsilon \in (0, 1]$ representing the amplitude of the potential, we employ the unitary flow property of the (second-order) time-splitting Fourier pseudospectral (TSFP) method in $L^2$-norm to prove a uniform error bound at $C(T)(h^m +τ^2)$ up to the long time $T_\varepsilon= T/\varepsilon$ for any $T>0$ and uniformly for $0<\varepsilon\le1$, while $h$ is the mesh size, $τ$ is the time step, $m \ge 2$ depends on the regularity of the exact solution, and $C(T) =C_0+C_1T$ grows at most linearly with respect to $T$ with $C_0$ and $C_1$ two positive constants independent of $T$, $\varepsilon$, $h$ and $τ$. Then by introducing a new technique of {\sl regularity compensation oscillation} (RCO) in which the high frequency modes are controlled by regularity and the low frequency modes are analyzed by phase cancellation and energy method, an improved uniform error bound at $O(h^{m-1} + \varepsilon τ^2)$ is established in $H^1$-norm for the long-time dynamics up to the time at $O(1/\varepsilon)$ of the Schrödinger equation with $O(\varepsilon)$-potential with $m \geq 3$, which is uniformly for $\varepsilon\in(0,1]$. Moreover, the RCO technique is extended to prove an improved uniform error bound at $O(h^{m-1} + \varepsilon^2τ^2)$ in $H^1$-norm for the long-time dynamics up to the time at $O(1/\varepsilon^2)$ of the cubic NLSE with $O(\varepsilon^2)$-nonlinearity strength, uniformly for $\varepsilon \in (0, 1]$. Extensions to the first-order and fourth-order time-splitting methods are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved uniform error bounds on time-splitting methods for the long-time dynamics of the weakly nonlinear Dirac equation

Improved uniform error bounds on time-splitting methods are rigorously proven for the long-time dynamics of the weakly nonlinear Dirac equation (NLDE), where the nonlinearity strength is characterized by a dimensionless parameter $\varepsilon \in (0, 1]$ . We adopt a second order Strang splitting method to discretize the NLDE in time and combine the Fourier pseudospectral method in space for the full-discretization. By employing the {\sl regularity compensation oscillation} (RCO) technique where the high frequency modes are controlled by the regularity of the exact solution and the low frequency modes are analyzed by phase cancellation and energy method, we establish improved uniform error bounds at $O(\varepsilon^2τ^2)$ and $O(h^{m-1}+ \varepsilon^2τ^2)$ for the second-order Strang splitting semi-discretizaion and full-discretization up to the long-time $T_{\varepsilon} = T/\varepsilon^2$ with $T>0$ fixed, respectively. Furthermore, the numerical scheme and error estimates are extended to an oscillatory NLDE which propagates waves with $O(\varepsilon^2)$ wavelength in time. Finally, numerical examples verifying our analytical results are given.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Volume-preserving parametric finite element methods for axisymmetric geometric evolution equations

We propose and analyze volume-preserving parametric finite element methods for surface diffusion, conserved mean curvature flow and an intermediate evolution law in an axisymmetric setting. The weak formulations are presented in terms of the generating curves of the axisymmetric surfaces. The proposed numerical methods are based on piecewise linear parametric finite elements. The constructed fully practical schemes satisfy the conservation of the enclosed volume. In addition, we prove the unconditional stability and consider the distribution of vertices for the discretized schemes. The introduced methods are implicit and the resulting nonlinear systems of equations can be solved very efficiently and accurately via the Newton's iterative method. Numerical results are presented to show the accuracy and efficiency of the introduced schemes for computing the considered axisymmetric geometric flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An energy-stable parametric finite element method for simulating solid-state dewetting

We propose an energy-stable parametric finite element method (ES-PFEM) for simulating solid-state dewetting of thin films in two dimensions via a sharp-interface model, which is governed by surface diffusion and contact line (point) migration together with proper boundary conditions. By reformulating the relaxed contact angle condition into a Robin-type boundary condition and then treating it as a natural boundary condition, we obtain a new variational formulation for the problem, in which the interface curve and its contact points are evolved simultaneously. Then, the variational problem is discretized in space by using piecewise linear elements. A full discretization is presented by adopting the backward Euler method in time, and the well-posedness and energy dissipation of the full discretization are established. The numerical method is semi-implicit (i.e., a linear system to be solved at each time step and thus efficient), unconditionally energy-stable with respect to the time step, and second-order in space measured by a manifold distance between two curves. In addition, it demonstrates equal mesh distribution when the solution reaches its equilibrium, i.e., long-time dynamics. Numerical results are reported to show accuracy and efficiency as well as some good properties of the proposed numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Power-law scaling for solid-state dewetting of thin films: an Onsager variational approach

We examine the kinetics of surface diffusion-controlled, solid-state dewetting by consideration of the retraction of the contact in a semi-infinite solid thin film on a flat rigid substrate. The analysis is performed within the framework of the Onsager variational principle applied to surface diffusion-controlled morphology evolution. Based on this approach, we derive a simple, reduced-order model to quantitatively analyse the power-law scaling of the dewetting process. Using asymptotic analysis and numerical simulations for the reduced-order model, we find that the retraction distance grows as the $2/5$ power of time and the height of the ridge, adjacent to the contact, grows as the $1/5$ power of time for late time. While the asymptotic analysis focuses on late time and a relatively simple geometric model, the Onsager approach is applicable to all times and descriptions of the morphology of arbitrary complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sharp-interface model for simulating solid-state dewetting in three dimensions

The problem of simulating solid-state dewetting of thin films in three dimensions (3D) by using a sharp-interface approach is considered in this paper. Based on the thermodynamic variation, a speed method is used for calculating the first variation to the total surface energy functional. The speed method shares more advantages than the traditional use of parameterized curves (or surfaces), e.g., it is more intrinsic and its variational structure (related with Cahn-Hoffman $\boldsymbolξ$-vector) is clearer and more direct. By making use of the first variation, necessary conditions for the equilibrium shape of the solid-state dewetting problem is given, and a kinetic sharp-interface model which includes the surface energy anisotropy is also proposed. This sharp-interface model describes the interface evolution in 3D which occurs through surface diffusion and contact line migration. By solving the proposed model, we perform lots of numerical simulations to investigate the evolution of patterned films, e.g., the evolution of a short cuboid and pinch-off of a long cuboid. Numerical simulations in 3D demonstrate the accuracy and efficacy of the sharp-interface approach to capture many of the complexities observed in solid-state dewetting experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A parametric finite element method for solid-state dewetting problems in three dimensions

We propose a parametric finite element method (PFEM) for efficiently solving the morphological evolution of solid-state dewetting of thin films on a flat rigid substrate in three dimensions (3D). The interface evolution of the dewetting problem in 3D is described by a sharp-interface model, which includes surface diffusion coupled with contact line migration. A variational formulation of the sharp-interface model is presented, and a PFEM is proposed for spatial discretization. For temporal discretization, at each time step, we first update the position of the contact line according to the relaxed contact angle condition; then, by using the position of the new contact line as the boundary condition, we solve a linear algebra system resulted from the discretization of PFEM to obtain the new interface surface for the next step. The well-posedness of the solution of the PFEM is also established. Extensive numerical results are reported to demonstrate the accuracy and efficiency of the proposed PFEM and to show the complexities of the dewetting morphology evolution observed in solid-state dewetting experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dimension reduction for dipolar Bose-Einstein condensates in the strong interaction regime

We study dimension reduction for the three-dimensional Gross-Pitaevskii equation with a long-range and anisotropic dipole-dipole interaction modeling dipolar Bose-Einstein condensation in a strong interaction regime. The cases of disk shaped condensates (confinement from dimension three to dimension two) and cigar shaped condensates (confinement to dimension one) are analyzed. In both cases, the analysis combines averaging tools and semiclassical techniques. Asymptotic models are derived, with rates of convergence in terms of two small dimensionless parameters characterizing the strength of the confinement and the strength of the interaction between atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fractional quantum mechanics in polariton condensates with velocity dependent mass

We introduce and analyze a novel mean-field model for polariton condensates with velocity dependence of the effective polariton mass due the photon and exciton components. The effective mass depends on the in-plane wave vector k, which at the inflection point of the lower polariton energy branch becomes infinite and above negative. The polariton condensate modes of the new mean-field theory are now sensitive to mass variations and for certain points of the energy dispersion the polariton condensate mode represents fractional quantum mechanics. The impact of the generalized kinetic energy term is elucidated by numerical studies in 1D and 2D showing significant differences for large velocities. Analytical expressions for plane wave solutions as well as a linear waves analysis show the significance of this new model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A uniformly accurate multiscale time integrator pseudospectral method for the Klein-Gordon equation in the nonrelativistic limit regime

We propose and analyze a multiscale time integrator Fourier pseudospectral (MTI-FP) method for solving the Klein-Gordon (KG) equation with a dimensionless parameter $0<\varepsilon\leq1$ which is inversely proportional to the speed of light. In the nonrelativistic limit regime, i.e. $0<\varepsilon\ll1$, the solution to the KG equation propagates waves with amplitude at $O(1)$ and wavelength at $O(\varepsilon^2)$ in time and $O(1)$ in space, which causes significantly numerical burdens due to the high oscillation in time. The MTI-FP method is designed by adapting a multiscale decomposition by frequency (MDF) to the solution at each time step and applying an exponential wave integrator to the nonlinear Schrödinger equation with wave operator under well-prepared initial data for $\varepsilon^2$-frequency and $O(1)$-amplitude waves and a KG-type equation with small initial data for the reminder waves in the MDF. We rigorously establish two independent error bounds in $H^2$-norm to the MTI-FP method at $O(h^{m_0}+τ^2+\varepsilon^2)$ and $O(h^{m_0}+τ^2/\varepsilon^2)$ with $h$ mesh size, $τ$ time step and $m_0\ge2$ an integer depending on the regularity of the solution, which immediately imply that the MTI-FP converges uniformly and optimally in space with exponential convergence rate if the solution is smooth, and uniformly in time with linear convergence rate at $O(τ)$ for all $\varepsilon\in(0,1]$ and optimally with quadratic convergence rate at $O(τ^2)$ in the regimes when either $\varepsilon=O(1)$ or $0<\varepsilon\le τ$. Numerical results are reported to confirm the error bounds and demonstrate the efficiency and accuracy of the MTI-FP method for the KG equation, especially in the nonrelativistic limit regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computing the ground state and dynamics of the nonlinear Schrödinger equation with nonlocal interactions via the nonuniform FFT

We present efficient and accurate numerical methods for computing the ground state and dynamics of the nonlinear Schrödinger equation (NLSE) with nonlocal interactions based on a fast and accurate evaluation of the long-range interactions via the nonuniform fast Fourier transform (NUFFT). We begin with a review of the fast and accurate NUFFT based method in \cite{JGB} for nonlocal interactions where the singularity of the Fourier symbol of the interaction kernel at the origin can be canceled by switching to spherical or polar coordinates. We then extend the method to compute other nonlocal interactions whose Fourier symbols have stronger singularity at the origin that cannot be canceled by the coordinate transform. Many of these interactions do not decay at infinity in the physical space, which adds another layer of complexity since it is more difficult to impose the correct artificial boundary conditions for the truncated bounded computational domain. The performance of our method against other existing methods is illustrated numerically, with particular attention on the effect of the size of the computational domain in the physical space. Finally, to study the ground state and dynamics of the NLSE, we propose efficient and accurate numerical methods by combining the NUFFT method for potential evaluation with the normalized gradient flow using backward Euler Fourier pseudospectral discretization and time-splitting Fourier pseudospectral method, respectively. Extensive numerical comparisons are carried out between these methods and other existing methods for computing the ground state and dynamics of the NLSE with various nonlocal interactions. Numerical results show that our scheme performs much better than those existing methods in terms of both accuracy and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dimension reduction for anisotropic Bose-Einstein condensates in the strong interaction regime

We study the problem of dimension reduction for the three dimensional Gross-Pitaevskii equation (GPE) describing a Bose-Einstein condensate confined in a strongly anisotropic harmonic trap. Since the gas is assumed to be in a strong interaction regime, we have to analyze two combined singular limits: a semi-classical limit in the transport direction and the strong partial confinement limit in the transversal direction. We prove that both limits commute together and we provide convergence rates. The by-products of this work are approximated models in reduced dimension for the GPE, with a priori estimates of the approximation errors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast and Accurate Evaluation of Nonlocal Coulomb and Dipole-Dipole Interactions via the Nonuniform FFT

We present a fast and accurate algorithm for the evaluation of nonlocal (long-range) Coulomb and dipole-dipole interactions in free space. The governing potential is simply the convolution of an interaction kernel $U(\bx)$ and a density function $ρ(\bx)=|ψ(\bx)|^2$, for some complex-valued wave function $ψ(\bx)$, permitting the formal use of Fourier methods. These are hampered by the fact that the Fourier transform of the interaction kernel $\widehat{U}(\bk)$ has a singularity at the origin $\bk={\bf 0}$ in Fourier (phase) space. Thus, accuracy is lost when using a uniform Cartesian grid in $\bk$ which would otherwise permit the use of the FFT for evaluating the convolution. Here, we make use of a high-order discretization of the Fourier integral, accelerated by the nonuniform fast Fourier transform (NUFFT). By adopting spherical and polar phase-space discretizations in three and two dimensions, respectively, the singularity in $\hat{U}(\bk)$ at the origin is canceled, so that only a modest number of degrees of freedom are required to evaluate the Fourier integral, assuming that the density function $ρ(\bx)$ is smooth and decays sufficiently fast as $\bx \rightarrow \infty$. More precisely, the calculation requires $O(N\log N)$ operations, where $N$ is the total number of discretization points in the computational domain. Numerical examples are presented to demonstrate the performance of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ground states and dynamics of spin-orbit-coupled Bose-Einstein condensates

We study analytically and asymptotically as well as numerically ground states and dynamics of two-component spin-orbit-coupled Bose-Einstein condensates (BECs) modeled by the coupled Gross-Pitaevskii equations (CGPEs). In fact, due to the appearance of the spin-orbit (SO) coupling in the two-component BEC with a Raman coupling, the ground state structures and dynamical properties become very rich and complicated. For the ground states, we establish the existence and non-existence results under different parameter regimes, and obtain their limiting behaviors and/or structures with different combinations of the SO and Raman coupling strengths. For the dynamics, we show that the motion of the center-of-mass is either non-periodic or with different frequency to the trapping frequency when the external trapping potential is taken as harmonic and the initial data is chosen as a stationary state (e.g. ground state) with a shift, which is completely different from the case of a two-component BEC without the SO coupling, and obtain the semiclassical limit of the CGPEs in the linear case via the Wigner transform method. Efficient and accurate numerical methods are proposed for computing the ground states and dynamics, especially for the case of box potentials. Numerical results are reported to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomenon in the SO-coupled BECs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sharp interface model for solid-state dewetting problems with weakly anisotropic surface energies

We propose a sharp interface model for simulating solid-state dewetting where the surface energy is (weakly) anisotropic. The morphology evolution of thin films is governed by surface diffusion and contact line migration. The mathematical model is based on an energy variational approach. Anisotropic surface energies lead to multiple solutions of the contact angle equation at contact points. Introduction of a finite contact point mobility is both physically based and leads to robust, unambiguous determination of the contact angles. We implement the mathematical model in an explicit finite difference scheme with cubic spline interpolation for evolving marker points. Following validation of the mathematical and numerical approaches, we simulate the evolution of thin film islands, semi-infinite films, and films with holes as a function of film dimensions, Young's angle $θ_i$, anisotropy strength and crystal symmetry, and film crystal orientation relative to the substrate normal. We find that the contact point retraction rate can be well described by a power-law, $l \sim t^n$. Our results demonstrate that the exponent $n$ is not universal -- it is sensitive to the Young's angle $θ_i$ (and insensitive to anisotropy). In addition to classical wetting (where holes in a film heal) and dewetting (where holes in a film grow), we observe cases where a hole through the film heals but leave a finite size hole/bubble between the continuous film and substrate or where the hole heals leaving a continuous film that is not bonded to the substrate. Surface energy anisotropy (i) increases the instability that leads to island break-up into multiple islands, (ii) enhances hole healing, and (iii) leads to finite island size even under some conditions where the Young's angle $θ_i$ suggests that the film wets the substrate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A simple and efficient numerical method for computing the dynamics of rotating Bose-Einstein condensates via a rotating Lagrangian coordinate

We propose a simple, efficient and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with/without a long-range dipole-dipole interaction. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range dipole-dipole interaction, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential and review some dynamical laws related to the 2D and 3D GPE. By introducing a rotating Lagrangian coordinate system, the original GPEs are re-formulated to GPEs without the angular momentum rotation which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable and very efficient in computation. It is spectral order accurate in space and second-order accurate in time, and conserves the mass in the discrete level. Extensive numerical results are reported to demonstrate the efficiency and accuracy of the new numerical method. Finally, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation and center-of-mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without/with the long-range dipole-dipole interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computational methods for the dynamics of the nonlinear Schrodinger/Gross-Pitaevskii equations

In this paper, we begin with the nonlinear Schrodinger/Gross-Pitaevskii equation (NLSE/GPE) for modeling Bose-Einstein condensation (BEC) and nonlinear optics as well as other applications, and discuss their dynamical properties ranging from time reversible, time transverse invariant, mass and energy conservation, dispersion relation to soliton solutions. Then, we review and compare different numerical methods for solving the NLSE/GPE including finite difference time domain methods and time-splitting spectral method, and discuss different absorbing boundary conditions. In addition, these numerical methods are extended to the NLSE/GPE with damping terms and/or an angular momentum rotation term as well as coupled NLSEs/GPEs. Finally, applications to simulate a quantized vortex lattice dynamics in a rotating BEC are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Uniform and optimal error estimates of an exponential wave integrator sine pseudospectral method for the nonlinear Schrodinger equation with wave operator

We propose an exponential wave integrator sine pseudospectral (EWI-SP) method for the nonlinear Schrödinger equation (NLS) with wave operator (NLSW), and carry out rigorous error analysis. The NLSW is NLS perturbed by the wave operator with strength described by a dimensionless parameter $\varepsilon\in(0,1]$. As $\varepsilon\to0^+$, the NLSW converges to the NLS and for the small perturbation, i.e. $0<\varepsilon\ll1$, the solution of the NLSW differs from that of the NLS with a function oscillating in time with $O(\varepsilon^2)$-wavelength at $O(\varepsilon^2)$ and $O(\varepsilon^4)$ amplitudes for ill-prepared and well-prepared initial data, respectively. This rapid oscillation in time brings significant difficulties in designing and analyzing numerical methods with error bounds uniformly in $\varepsilon$. In this work, we show that the proposed EWI-SP possesses the optimal uniform error bounds at $O(τ^2)$ and $O(τ)$ in $τ$ (time step) for well-prepared initial data and ill-prepared initial data, respectively, and spectral accuracy in $h$ (mesh size) for the both cases, in the $L^2$ and semi-$H^1$ norms. This result significantly improves the error bounds of the finite difference methods for the NLSW. Our approach involves a careful study of the error propagation, cut-off of the nonlinearity and the energy method. Numerical examples are provided to confirm our theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effective dipole-dipole interactions in multilayered dipolar Bose-Einstein condensates

We propose a two-dimensional model for a multilayer stack of dipolar Bose-Einstein condensates formed by a strong optical lattice. We derive effective intra- and interlayer dipole-dipole interaction potentials and provide simple analytical approximations for a given number of lattice sites at arbitrary polarization. We find that the interlayer dipole-dipole interaction changes the transverse aspect ratio of the ground state in the central layers depending on its polarization and the number of lattice sites. The changing aspect ratio should be observable in time of flight images. Furthermore, we show that the interlayer dipole-dipole interaction reduces the excitation energy of local perturbations, affecting the development of a roton minimum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gross-Pitaevskii-Poisson equations for dipolar Bose-Einstein condensate with anisotropic confinement

Ground states and dynamical properties of dipolar Bose-Einstein condensate are analyzed based on the Gross-Pitaevskii-Poisson system (GPPS) and its dimension reduction models under anisotropic confining potential. We begin with the three-dimensional (3D) Gross-Pitaevskii-Poisson system and review its quasi-2D approximate equations when the trap is strongly confined in $z$-direction and quasi-1D approximate equations when the trap is strongly confined in $x$-, $y$-directions. In fact, in the quasi-2D equations, a fractional Poisson equation with the operator $(-Δ)^{1/2}$ is involved which brings significant difficulties into the analysis. Existence and uniqueness as well as nonexistence of the ground state under different parameter regimes are established for the quasi-2D and quasi-1D equations. Well-posedness of the Cauchy problem for both equations and finite time blowup in 2D are analyzed. Finally, we rigorously prove the convergence and linear convergence rate between the solutions of the 3D GPPS and its quasi-2D and quasi-1D approximate equations in weak interaction regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scattering and bound states in two-dimensional anisotropic potentials

We propose a framework for calculating scattering and bound state properties in anisotropic two-dimensional potentials. Using our method, we derive systematic approximations of partial wave phase shifts and binding energies. Moreover, the method is suitable for efficient numerical computations. We calculate the s-wave phase shift and binding energy of polar molecules in two layers polarized by an external field along an arbitrary direction. We find that scattering depends strongly on their polarization direction and that absolute interlayer binding energies are larger than thermal energies at typical ultracold temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mean-field regime of trapped dipolar Bose-Einstein condensates in one and two dimensions

We derive rigorous one- and two-dimensional mean-field equations for cigar- and pancake-shaped dipolar Bose-Einstein condensates with arbitrary polarization angle. We show how the dipolar interaction modifies the contact interaction of the strongly confined atoms. In addition, our equations introduce a nonlocal potential, which is anisotropic for pancake-shaped condensates. We propose to observe this anisotropy via measurement of the condensate aspect ratio. We also derive analytically approximate density profiles from our equations. Both the numerical solutions of our reduced mean-field equations and the analytical density profiles agree well with numerical solutions of the full Gross-Pitaevskii equation while being more efficient to compute.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Computing Ground States of Spin-1 Bose-Einstein Condensates by the Normalized Gradient Flow

In this paper, we propose an efficient and accurate numerical method for computing the ground state of spin-1 Bose-Einstein condensates (BEC) by using the normalized gradient flow or imaginary time method. The key idea is to find a third projection or normalization condition based on the relation between the chemical potentials so that the three projection parameters used in the projection step of the normalized gradient flow are uniquely determined by this condition as well as the other two physical conditions given by the conservation of total mass and total magnetization. This allows us to successfully extend the most popular and powerful normalized gradient flow or imaginary time method for computing the ground state of single component BEC to compute the ground state of spin-1 BEC. An efficient and accurate discretization scheme, the backward-forward Euler sine-pseudospectral method (BFSP), is proposed to discretize the normalized gradient flow. Extensive numerical results on ground states of spin-1 BEC with ferromagnetic/antiferromagnetic interaction and harmonic/optical lattice potential in one/three dimensions are reported to demonstrate the efficiency of our new numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Convergence rate of dimension reduction in Bose-Einstein condensates

In this paper, we study dimension reduction of the three-dimensional (3D) Gross-Pitaevskii equation (GPE) modelling Bose-Einstein condensation under different limiting interaction and trapping frequencies parameter regimes. Convergence rates for the dimension reduction of 3D ground state and dynamics of the GPE in the case of disk-shaped condensation and cigar-shaped condensation are reported based on our asymptotic and numerical results. In addition, the parameter regimes in which the 3D GPE cannot be reduced to lower dimensions are identified.

0 citations0 reviewsNo code linkedOpen access