Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Collective excitations of trapped Bose-Einstein condensates in the energy and time domains

A time-dependent method for calculating the collective excitation frequencies and densities of a trapped, inhomogeneous Bose-Einstein condensate with circulation is presented. The results are compared with time-independent solutions of the Bogoliubov-de Gennes equations. The method is based on time-dependent linear-response theory combined with spectral analysis of moments of the excitation modes of interest. The technique is straightforward to apply, extremely efficient in our implementation with parallel fast Fourier transform methods, and produces highly accurate results. For high dimensionality or low symmetry the time-dependent approach is a more practical computational scheme and produces accurate and reliable data. The method is suitable for general trap geometries, condensate flows and condensates permeated with defects and vortex structures.

Vortex line and ring dynamics in trapped Bose-Einstein condensates

Vortex dynamics in inhomogeneous Bose-Einstein condensates are studied numerically in two and three dimensions. We simulate the precession of a single vortex around the center of a trapped condensate, and use the Magnus force to estimate the precession frequency. Vortex ring dynamics in a spherical trap are also simulated, and we discover that a ring undergoes oscillatory motion around a circle of maximum energy. The position of this locus is calculated as a function of the number of condensed atoms. In the presence of dissipation, the amplitude of the oscillation will increase, eventually resulting in self-annihilation of the ring.

Vortex shedding and drag in dilute Bose-Einstein condensates

Above a critical velocity, the dominant mechanism of energy transfer between a moving object and a dilute Bose-Einstein condensate is vortex formation. In this paper, we discuss the critical velocity for vortex formation and the link between vortex shedding and drag in both homogeneous and inhomogeneous condensates. We find that at supersonic velocities sound radiation also contributes significantly to the drag force.

The distorted-wave impulse approximation for electron capture into excited states

Total cross sections for electron capture are calculated for collisions of fast protons and a-particles with atomic hydrogen. The distorted-wave impulse approximation is applied over the energy range 10-1500 keV/u. State-selective results are given for the 1s, 2s and 2p levels. Both the post and prior forms of the model are calculated and compared with results from other theories and experimental measurements. In general the model performs very well in comparison with experiment over this energy range though discrepancies arise at lower energies.

Dissipation and vortex creation in Bose-Einstein condensed gases

We solve the Gross-Pitaevskii equation to study energy transfer from an oscillating

Wave mixing of hybrid Bogoliubov modes in a Bose-Einstein condensate

Mode-mixing of coherent excitations of a trapped Bose-Einstein condensate is modeled using the Bogoliubov approximation. Calculations are presented for second-harmonic generation between the two lowest-lying even-parity m=0 modes in an oblate spheroidal trap. Hybridization of the modes of the breather (l=0) and surface (l=4) states leads to the formation of a Bogoliubov dark state near phase-matching resonance so that a single mode is coherently populated. Efficient harmonic generation requires a strong coupling rate, sharply-defined and well-separated frequency spectrum, and good phase matching. We find that in all three respects the quantal results are significantly different from hydrodynamic predictions. Typically the second-harmonic conversion rate is half that given by an equivalent hydrodynamic estimate.

Phylogeographic analysis of North American populations of the parasitic herbaceous plant Monotropa hypopitys reveals a complex history of range expansion from multiple late glacial refugia

Aim We carried out a phylogeographic study across the range of the herbaceous plant species Monotropa hypopitys L. in North America to determine whether its current disjunct distribution is due to recolonization from separate eastern and western refugia after the Last Glacial Maximum (LGM). Location North America: Pacific Northwest and north-eastern USA/south-eastern Canada. Methods Palaeodistribution modelling was carried out to determine suitable climatic regions for M. hypopitys at the LGM. We analysed between 155 and 176 individuals from 39 locations spanning the species' entire range in North America. Sequence data were obtained for the chloroplast rps2 gene (n=168) and for the nuclear ITS region (n=158). Individuals were also genotyped for eight microsatellite loci (n=176). Interpolation of diversity values was used to visualize the range-wide distribution of genetic diversity for each of the three marker classes. Minimum spanning networks were constructed showing the relationships between the rps2 and ITS haplotypes, and the geographical distributions of these haplotypes were plotted. The numbers of genetic clusters based on the microsatellite data were estimated using Bayesian clustering approaches. Results The palaeodistribution modelling indicated suitable climate envelopes for M. hypopitys at the LGM in both the Pacific Northwest and south-eastern USA. High levels of genetic diversity and endemic haplotypes were found in Oregon, the Alexander Archipelago, Wisconsin, and in the south-eastern part of the species' distribution range. Main conclusions Our results suggest a complex recolonization history for M. hypopitys in North America, involving persistence in separate eastern and western refugia. A generally high degree of congruence between the different marker classes analysed indicated the presence of multiple refugia, with at least two refugia in each area. In the west, putative refugia were identified in Oregon and the Alexander Archipelago, whereas eastern refugia may have been located in the southern part of the species' current distribution, as well as in the 'Driftless Area'. These findings are in contrast to a previous study on the related species Orthilia secunda, which has a similar disjunct distribution to M. hypopitys, but which appears to have recolonized solely from western refugia. © 2011 Blackwell Publishing Ltd.