Search

When the intensity of turbulence is increased (by increasing the Reynolds number, e.g. by reducing the viscosity of the fluid), the rate of the dissipation of kinetic energy decreases but does not tend asymptotically to zero: it levels off to a non-zero constant as smaller and smaller vortical flow structures are generated. This fundamental property, called the dissipation anomaly, is sometimes referred to as the zeroth law of turbulence. The question of what happens in the limit of vanishing viscosity (purely hypothetical in classical fluids) acquires a particular physical significance in the context of liquid helium, a quantum fluid which becomes effectively inviscid at low temperatures achievable in the laboratory. By performing numerical simulations and identifying the superfluid Reynolds number, here we show evidence for a superfluid analog to the classical dissipation anomaly. Our numerics indeed show that as the superfluid Reynolds number increases, smaller and smaller structures are generated on the quantized vortex lines on which the superfluid vorticity is confined, balancing the effect of weaker and weaker dissipation., 14 pages , 9 figures. Accepted in Phys Rev Fluids

We probe the collective dynamics of a quantum degenerate Bose-Bose mixture of Cs and $^{174}$Yb with attractive interspecies interactions. Specifically, we excite vertical center of mass oscillations of the Cs condensate, and observe significant damping for the Cs dipole mode, due to the rapid transfer of energy to the larger Yb component, and the ensuing acoustic dissipation. Numerical simulations based on coupled Gross-Pitaevskii equations provide excellent agreement, and additionally reveal the possibility of late-time revivals (beating) which are found to be highly sensitive to the Cs and Yb atom number combinations. By further tuning the interaction strength of Cs using a broad Feshbach resonance, we explore the stability of the degenerate mixture, and observe collapse of the Cs condensate mediated by the attractive Cs-Yb interaction when $a_{\mathrm{Cs}}, 13 pages, 7 figures (inc. supplemental material), supplemental movie

Human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) have promising clinical applications which often rely on clonally-homogeneous cell populations. To achieve this, it is important to ensure that each colony originates from a single founding cell and to avoid subsequent merging of colonies during their growth. Clonal homogeneity can be obtained with low seeding densities; however, this leads to low yield and viability. It is therefore important to quantitatively assess how seeding density affects clonality loss so that experimental protocols can be optimised to meet the required standards. Here we develop a quantitative framework for modelling the growth of hESC colonies from a given seeding density based on stochastic exponential growth. This allows us to identify the timescales for colony merges and over which colony size no longer predicts the number of founding cells. We demonstrate the success of our model by applying it to our own experiments of hESC colony growth; while this is based on a particular experimental set-up, the model can be applied more generally to other cell lines and experimental conditions to predict these important timescales.

Reconnections of coherent filamentary structures play a key role in the dynamics of fluids, redistributing energy and helicity among the length scales, triggering dissipative effects and inducing fine-scale mixing. Unlike ordinary (classical) fluids where vorticity is a continuous field, in superfluid helium and in atomic Bose-Einstein condensates (BECs) vorticity takes the form of isolated quantised vortex lines, which are conceptually easier to study. New experimental techniques now allow visualisation of individual vortex reconnections in helium and condensates. It has long being suspected that reconnections obey universal laws, particularly a universal scaling with time of the minimum distance between vortices $\delta$. Here we perform a comprehensive analysis of this scaling across a range of scenarios relevant to superfluid helium and trapped condensates, combining our own numerical simulations with the previous results in the literature. We reveal that the scaling exhibit two distinct fundamental regimes: a $\delta \sim t^{1/2}$ scaling arising from the mutual interaction of the reconnecting strands and a $\delta \sim t$ scaling when extrinsic factors drive the individual vortices., Comment: 13 pages, 8 Figures main manuscript; 9 pages, 4 Figures Supporting Information. Accepted in PNAS

The aim of this chapter is to convey the importance and usefulness of mathematical modeling as a tool to achieve a deeper understanding of stem cell biology. We introduce key mathematical concepts (random walk theory, differential equations, and agent-based modeling), which form the basis of current descriptions of induced pluripotent stem cells. We hope to encourage a meaningful dialogue between biologists and mathematicians and highlight the value of such an interdisciplinary approach.

OCT4 is one of the transcription factors required to maintain an undifferentiated state in human embryonic stem cells (hESCs). Thus, it is crucial to understand how OCT4 transcription is regulated both at the single-cell and colony level. Here we analyse the changes of OCT4-mCherry intensity expression in hESCs in the presence and absence of the BMP4 morphogenetic protein.We show that OCT4 expression is dynamic, reaching a maximum response 10 h after BMP4 treatment. We obtain the stationary probability distributions that govern the hESCs transitions amongst the different cell states in the presence/absence of BMP4 and establish the times at which the hESCs, that lead to differentiated and pluripotent cells, cluster in the colony. Furthermore, by quantifying the similarities between the OCT4 expression amongst neighbouring hESCs, we show that hESCs express, on average, similar values in their local neighbourhood within the first two days of the experiment and before BMP4 treatment. These results are relevant for the development of mathematical and computational models of adherent hESC colonies.

We study numerically the formation of a vortex lattice inside a rotating bucket containing superfluid helium, paying attention to an important feature which is practically unavoidable in all experiments: the microscopic roughness of the bucket's surface. We model this using the Gross-Pitaevskii for a weakly-interacting Bose gas, a model which is idealised when applied to superfluid helium but captures the key physics of the vortex dynamics which we are interested in. We find that the vortex lattice arises from the interaction and reconnections of nucleated U-shaped vortex lines, which merge and align along the axis of rotation. We quantify the effects which the surface roughness and remanent vortex lines play in this process., 11 pages, 12 figures

The improved in vitro regulation of human embryonic stem cell (hESC) pluripotency and differentiation trajectories is required for their promising clinical applications. The temporal and spatial quantification of the molecular interactions controlling pluripotency is also necessary for the development of successful mathematical and computational models. Here we use time-lapse experimental data of OCT4-mCherry fluorescence intensity to quantify the temporal and spatial dynamics of the pluripotency transcription factor OCT4 in a growing hESC colony in the presence and absence of BMP4. We characterise the internal self-regulation of OCT4 using the Hurst exponent and autocorrelation analysis, quantify the intra-cellular fluctuations and consider the diffusive nature of OCT4 evolution for individual cells and pairs of their descendants. We find that OCT4 abundance in the daughter cells fluctuates sub-diffusively, showing anti-persistent self-regulation. We obtain the stationary probability distributions governing hESC transitions amongst the different cell states and establish the times at which pro-fate cells (which later give rise to pluripotent or differentiated cells) cluster in the colony. By quantifying the similarities between the OCT4 expression amongst neighbouring cells, we show that hESCs express similar OCT4 to cells within their local neighbourhood within the first two days of the experiment and before BMP4 treatment. Our framework allows us to quantify the relevant properties of proliferating hESC colonies and the procedure is widely applicable to other transcription factors and cell populations.

The leapfrogging of coaxial vortex rings is a famous effect which has been noticed since the times of Helmholtz. Recent advances in ultra-cold atomic gases show that the effect can now be studied in quantum fluids. The strong confinement which characterizes these systems motivates the study of leapfrogging of vortices within narrow channels. Using the two-dimensional point vortex model, we show that in the constrained geometry of a two-dimensional channel the dynamics is richer than in an unbounded domain: alongsize the known regimes of standard leapfrogging and the absence of it, we identify new regimes of backward leapfrogging and periodic orbits. Moreover, by solving the Gross-Pitaevskii equation for a Bose-Einstein condensate, we show that all four regimes exist for quantum vortices too. Finally, we discuss the differences between classical and quantum vortex leapfrogging which appear when the quantum healing length becomes significant compared to the vortex separation or the channel size, and when, due to high velocity, compressibility effects in the condensate becomes significant., Comment: 18 pages, 7 figures

Ultracold dipolar droplets have been realized in a series of ground-breaking experiments, where the stability of the droplet state is attributed to beyond-mean-field effects in the form of the celebrated Lee-Huang-Yang (LHY) correction. We scrutinize the dipolar droplet states in a one-dimensional context using a combination of analytical and numerical approaches, and identify experimentally viable parameters for accessing our findings for future experiments. In particular we identify regimes of stability in the restricted geometry, finding multiple roton instabilities as well as regions supporting quasi-one-dimensional droplet states. By applying an interaction quench to the droplet, a modulational instability is induced and multiple droplets are produced, along with bright solitons and atomic radiation. We also assess the droplets robustness to collisions, revealing population transfer and droplet fission., Comment: Expanded and revised version, J. Phys. Commun. in press

The behaviour of a harmonically trapped dipolar Bose-Einstein condensate with its dipole moments rotating at angular frequencies lower than the transverse harmonic trapping frequency is explored in the co-rotating frame. We obtain semi-analytical solutions for the stationary states in the Thomas-Fermi limit of the corresponding dipolar Gross-Pitaevskii equation and utilise linear stability analysis to elucidate a phase diagram for the dynamical stability of these stationary solutions with respect to collective modes. These results are verified via direct numerical simulations of the dipolar Gross-Pitaevskii equation, which demonstrate that dynamical instabilities of the co-rotating stationary solutions lead to the seeding of vortices that eventually relax into a triangular lattice configuration. Our results illustrate that rotation of the dipole polarization represents a new route to vortex formation in dipolar Bose-Einstein condensates., 11 pages, 5 figures; updated with responses to referee reports and minor typographical corrections

The possibility of effectively inverting the sign of the dipole-dipole interaction, by fast rotation of the dipole polarization, is examined within a harmonically trapped dipolar Bose-Einstein condensate. Our analysis is based on the stationary states in the Thomas-Fermi limit, in the corotating frame, as well as direct numerical simulations in the Thomas-Fermi regime, explicitly accounting for the rotating polarization. The condensate is found to be inherently unstable due to the dynamical instability of collective modes. This ultimately prevents the realization of robust and long-lived rotationally tuned states. Our findings have major implications for experimentally accessing this regime., Comment: 9 pages with 5 figures

Summary 1. Infectious disease outbreaks in plants threaten ecosystems, agricultural crops and food trade. Currently, several fungal diseases are affecting forests worldwide, posing a major risk to tree species, habitats and consequently ecosystem decay. Prediction and control of disease spread are difficult, mainly due to the complexity of the interaction between individual components involved. 2. In this work, we introduce a lattice-based epidemic model coupled with a stochastic process that mimics, in a very simplified way, the interaction between the hosts and pathogen. We studied the disease spread by measuring the propagation velocity of the pathogen on the susceptible hosts. Quantitative results indicate the occurrence of a critical transition between two stable phases: local confinement and an extended epiphytotic outbreak that depends on the density of the susceptible individuals. 3. Quantitative predictions of epiphytotics are performed using the framework early-warning indicators for impending regime shifts, widely applied on dynamical systems. These signals forecast successfully the outcome of the critical shift between the two stable phases before the system enters the epiphytotic regime. 4. Synthesis: Our study demonstrates that early-warning indicators could be useful for the prediction of forest disease epidemics through mathematical and computational models suited to more specific pathogen-host-environmental interactions., Comment: 11 pages

Human pluripotent stem cells hold great promise for developments in regenerative medicine and drug design. The mathematical modelling of stem cells and their properties is necessary to understand and quantify key behaviours and develop non-invasive prognostic modelling tools to assist in the optimisation of laboratory experiments. Here, the recent advances in the mathematical modelling of hPSCs are discussed, including cell kinematics, cell proliferation and colony formation, and pluripotency and differentiation.

Anaerobic digestion enables the water industry to treat wastewater as a resource for generating energy and recovering valuable by-products. The complexity of the anaerobic digestion process has motivated the development of complex models. However, this complexity makes it intractable to pin-point stability and emergent behaviour. Here, the widely used Anaerobic Digestion Model No. 1 (ADM1) has been reduced to its very backbone, a syntrophic two-tiered microbial ‘food chain’ and a slightly more complex three-tiered microbial ‘food web’, with their stability analysed as a function of the inflowing substrate concentration and dilution rate. Parameterised for phenol and chlorophenol degradation, steady-states were always stable and non-oscillatory. Low input concentrations of chlorophenol were sufficient to maintain chlorophenol- and phenol-degrading populations but resulted in poor conversion and a hydrogen flux that was too low to sustain hydrogenotrophic methanogens. The addition of hydrogen and phenol boosted the populations of all three organisms, resulting in the counterintuitive phenomena that (i) the phenol degraders were stimulated by adding hydrogen, even though hydrogen inhibits phenol degradation, and (ii) the dechlorinators indirectly benefitted from measures that stimulated their hydrogenotrophic competitors; both phenomena hint at emergent behaviour., Highlights • A chlorophenol-mineralising food web is mathematically analysed for local stability. • Adding hydrogen allows for full chlorophenol-mineralisation at high dilution rates. • At low chlorophenol, adding phenol allows for mutualism rather than competition to occur. • The dechlorinator may partially replace the methanogen as the syntroph.

We demonstrate an experimentally feasible method for generating the classical Kelvin-Helmholtz instability in a single component atomic Bose-Einstein condensate. By progressively reducing a potential barrier between two counter-flowing channels we seed a line of quantised vortices, which precede to form progressively larger clusters, mimicking the classical roll-up behaviour of the Kelvin-Helmholtz instability. This cluster formation leads to an effective superfluid shear layer, formed through the collective motion of many quantised vortices. From this we demonstrate a straightforward method to measure the effective viscosity of a turbulent quantum fluid in a system with a moderate number of vortices, within the range of current experimental capabilities.

We determine the evolution of a cluster of quantum vortices initially placed at the centre of a larger vortex-free region. We find that the cluster spreads out spatially. This spreading motion consists of two effects: the rapid evaporation of vortex dipoles from the cluster and the slow expansion of the cluster itself. The latter is akin to a diffusion process controlled by the quantum of circulation. Numerical simulations of the Gross-Pitaevskii equation show that this phenomenon is qualitatively unaffected by the presence of sound waves, vortex annihilations, and boundaries, and it should be possible to create it in the laboratory., Comment: Revised version accepted for publication in PRA; supplementary material not available in this version

Experiments and numerical simulations show that quantum turbulence exists in two distinct limiting regimes: Kolmogorov turbulence (which shares with classical turbulence the important property of a cascade of kinetic energy from large eddies to small eddies) and Vinen turbulence (which is more similar to a random flow). In this work, we define a mesoscale helicity for the superfluid, which, tested in numerical experiments, distinguishes the two turbulent regimes, quantifying the amount of nonlocal vortex interactions and the orientation of the vortex lines., Comment: To be published in Phys. Rev. B (2021)

We numerically study the spatial spreading of quantized vortex lines in low-temperature liquid helium. The vortex lines, initially concentrated in a small region, diffuse into the surrounding vortex-free helium, a situation which is typical of many experiments. We find that this spreading, which occurs in the absence of viscosity, emerges from the interactions between the vortex lines and can be interpreted as a diffusion process with effective coefficient equal to approximately $0.5\ensuremath{\kappa}$ where $\ensuremath{\kappa}$ is the quantum of circulation.

We study the elementary characteristics of turbulence in a quantum ferrofluid through the context of a dipolar Bose gas condensing from a highly non-equilibrium thermal state. Our simulations reveal that the dipolar interactions drive the emergence of polarized turbulence and density corrugations. The superfluid vortex lines and density fluctuations adopt a columnar or stratified configuration, depending on the sign of the dipolar interactions, with the vortices tending to form in the low density regions to minimize kinetic energy. When the interactions are dominantly dipolar, the decay of vortex line length is enhanced, closely following a $t^{-3/2}$ behaviour. This system poses exciting prospects for realizing stratified quantum turbulence and new levels of generating and controlling turbulence using magnetic fields., to be published in Phys. Rev. Lett

Motivated by recent experiments, we study theoretically the dynamics of vortices in the crossover from two dimensions to one in atomic condensates in elongated traps. We explore the transition from the dynamics of a vortex to that of a dark soliton as the one-dimensional limit is approached, mapping this transition out as a function of the key system parameters. Moreover, we probe this transition dynamically through the hysteresis under time-dependent deformation of the trap at the dimensionality crossover. When the solitonic regime is probed during the hysteresis, significant angular momentum is lost from the system, but remarkably, the vortex can reemerge.

We study the sensitivity of coupled condensate formation dynamics on the history of initial stochastic domain formation in the context of instantaneously quenched elongated harmonically-trapped immiscible two-component atomic Bose gases. The spontaneous generation of defects in the fastest condensing component, and subsequent coarse-graining dynamics, can lead to a deep oscillating microtrap into which the other component condenses, thereby establishing a long-lived composite defect in the form of a dark-bright solitary wave. We numerically map out diverse key aspects of these competing growth dynamics, focussing on the role of shot-to-shot fluctuations and global parameter changes (initial state choices, quench parameters and condensate growth rates). We conclude that phase-separated structures observable on experimental timescales are likely to be metastable states whose form is influenced by the stability and dynamics of the spontaneously-emerging dark-bright solitary wave., Comment: 5 pages + Refs, 5 figures. Revised and expanded from previous version

As well as being free from viscosity, the Bose–Einstein condensate has another striking property—it is constrained to circulate only through the presence of whirlpools of fixed size and quantized circulation. In contrast, in conventional fluids, the eddies can have arbitrary size and circulation. Here we establish the form of these quantum vortices, their key properties, and how they are formed and modelled.

Bose and Einstein’s prediction of Bose–Einstein condensation came out of their theory for how quantum particles in a gas behaved, and was built on the pioneering statistical approach of Boltzmann for classical particles. Here we follow Boltzmann, Bose and Einstein’s footsteps, leading to the derivation of Bose–Einstein condensation for an ideal gas and its key properties.

During the past year of the MCTDHB project at the HLRS, we continued to strive and conquest further applications, developments, and expansion of the MultiConfigurational Time-Dependent Hartree for Bosons (MCTDHB) method in the context of ultracold atomic systems. We also announce the MCTDH-X package, the Multiconfigurational Time-Dependent Hartree for Indistinguishable Particles X package, which is able to treat identical bosons and fermions, with or without spin/internal degrees of freedom, alike. Here we report on a plethora of results and versatile applications which include: (i) single-shot imaging of fluctuating vortices in a fragmented Bose-Einstein condensate (BEC); (ii) the many-body tunneling and fragmetnation of vortices in 2D trapped BECs; (iii) the transition from vortices to solitonic vortices in 2D trapped BECs; (iv) the variance of a many-particle system being very sensitive to correlations even in the infinite-particle limit; (v) the consequences of the latter on the out-of-equilibrium uncertainty product of an evolving BEC; (vi) the mechanism of tunneling to open space of a few interacting polarized fermions; and (vii) composite fragmentation of multi-components BECs (i.e., with internal degrees of freedom). These are all exciting results made throughout the allocation of computer time by the HLRS to the MCTDHB project. Finally, further perspectives and future research plans are briefly discussed.

In the previous chapter we considered the shape of steady state condensates, either homogeneous or confined by trapping potentials. We have seen that the condensate described by the GPE is a special kind of fluid, similar to the idealized Euler fluid without viscosity that appears in classical fluid dynamics textbooks. Not surprisingly for a fluid, the dynamics of the condensate exhibit a variety of interesting time-dependent phenomena, from sound waves and shape oscillations, to solitons and vortices.

The Gross-Pitaevskii equation (GPE) is a successful and well-established model for describing an atomic Bose-Einstein condensate. Here we introduce this model, along with its assumptions. Throughout the rest of this chapter we explore its properties and key time-independent solutions.

We map the ultrasonic (8 MHz) speed and attenuation of edible-grade gelatin in water, exploring the key dependencies on temperature, concentration and time. The ultrasonic signatures of the sol-gel transition, confirmed by rheological measurements, and incomplete gel formation at low concentrations, enable a phase diagram of the system to be constructed. Sensitivity is also demonstrated to the kinetics of gel formation and melting, and associated hysteresis effects upon cyclic temperature sweeps. Furthermore, simple acoustic models of the sol and gel state enable estimation of the speed of sound and compressibility of gelatin. Our results demonstrate the potential of ultrasonic measurements to characterise the structure and visco-elasticity of gelatin hydrogels., 15 pages, 8 figures

A theoretical study of vortex-lattice formation in atomic Bose-Einstein condensates confined by a rotating elliptical trap is presented. For the conventional case of purely s-wave interatomic interactions, this is done through a consideration of both hydrodynamic equations and time-dependent simulations of the Gross-Pitaevskii equation. We discriminate three distinct, experimentally testable regimes of instability: ripple, interbranch, and catastrophic. Additionally, we generalize the classical hydrodynamical approach to include long-range dipolar interactions, showing how the static solutions and their stability in the rotating frame are significantly altered. This enables us to examine the routes towards unstable dynamics, which, in analogy to conventional condensates, may lead to vortex-lattice formation.

We theoretically examine the rotation of an atomic Bose-Einstein condensate in an elliptical trap, both in the absence and presence of a quantized vortex. Two methods of introducing the rotating potential are considered - adiabatically increasing the rotation frequency at fixed ellipticity, and adiabatically increasing the trap ellipticity at fixed rotation frequency. Extensive simulations of the Gross-Pitaevskii equation are employed to map out the points where the condensate becomes unstable and ultimately forms a vortex lattice. We highlight the key features of having a quantized vortex in the initial condensate. In particular, we find that the presence of the vortex causes the instabilities to shift to lower or higher rotation frequencies, depending on the direction of the vortex relative to the trap rotation., Comment: 15 pages, 8 figures

Water-in-oil emulsions (WOE) are two-phase colloidal systems formed during crude oil production and spills. The high viscosity and stability of WOEs imply challenges during their clean-up and removal. Such emulsions are difficult to disaggregate due to a combination of chemical and physical factors. Ultrasound spectrometry can be used to characterise the WOE physical properties, providing access to, e.g. the droplet size distribution (DSD), density and viscosity. The DSD has been identified as a significant property impacting emulsion stability. This study focuses on the data post-acquisition stage modelling the droplet size growth as a stochastic process. Geometric Brownian motion (GBM) and Ito stochastic differential equations (SDEs) are used. Bayesian inference is introduced as a tool aiding in conditions of poor sample quality. The obtained model could predict emulsion separation indicated by a sufficiently large mean and standard deviation of the droplet growth process. It could be used for emulsions of different chemical compositions, including with added dispersants, allowing to characterise their impact on the WOE stability over time.

Medical ultrasound scanners are typically calibrated to the soft tissue average of 1540 m s$^{-1}$. In regions of different sound speed, for example, organs and tumours, the $B$-scan image then becomes a distortion of the true tissue cross-section, due to the misrepresentation of length and refraction. To quantify this distortion we develop a general geometric ray model for an object with an atypical speed of sound embedded in an ambient medium. We analyse the ensuing area distortion for circular and elliptical objects, mapping it out as a function of the key parameters, including the speed of sound mismatch, the object size and its elongation. We find that the area distortion can become significant, even for small-scale speed of sound mismatches. Our findings are verified by ultrasound imaging of a test object., Reduced in length, modified title

The motion of a dark soliton is investigated in a one-dimensional dilute Bose-Einstein condensate confined in a harmonic trap and an optical lattice. The harmonic trap induces a dynamical instability of the soliton, culminating in sound emission. The presence of the optical lattice enhances the instability, and in addition, dephases the emitted sound waves, thus preventing stabilisation of the soliton by sound reabsorption. This instability can be probed experimentally by monitoring the soliton oscillations under various lattice configurations, which can be realised by changing the intensity and angle between the laser beams that form the lattice. For short enough times, such that the emitted sound does not reinteract with the soliton, the power emitted by the soliton is found to be proportional to the square of the local soliton acceleration, which is in turn proportional to the deformation of the soliton profile., 12 pages, 6 figures. Submitted to J. Phys. B special issue for "Theory of Quantum Gases and Quantum Coherence" workshop (Levico, Italy, June 2003)

A dark soliton becomes unstable when it is incident on a background density gradient, and the induced instability results in the emission of sound. Detailed quantitative studies of sound emission are performed for various potentials, such as steps, linear ramps and gaussian traps. The amount of sound emission is found to be a significant fraction of the soliton energy for typical potentials. Continuous emission of sound is found to lead to an apparent deformation of the soliton profile. The power emitted by the soliton is shown to be parametrised by the square of the displacement of the centre of mass of the soliton from its density minimum, thus highlighting the significance of the inhomogeneity-induced soliton deformation., 22 pages, 11 figures, J. Phys. B accepted

We reinvestigate numerically the classic problem of two-dimensional superfluid flow past an obstacle. Taking the obstacle to be elongated (perpendicular to the flow), rather than the usual circular form, is shown to promote the nucleation of quantized vortices, enhance their subsequent interactions, and lead to wakes which bear striking similarity to their classical (viscous) counterparts. Then, focussing on the recent experiment of Kwon et al. (arXiv:1403.4658) in a trapped condensate, we show that an elliptical obstacle leads to a cleaner and more efficient means to generate two-dimensional quantum turbulence., 7 pages, 6 figures. Conference proceedings for the International Laser Physics Workshop (Sofia, July 2014)

In a recent experiment, Kwon et. al (arXiv:1403.4658 [cond-mat.quant-gas]) generated a disordered state of quantum vortices by translating an oblate Bose-Einstein condensate past a laser-induced obstacle and studying the subsequent decay of vortex number. Using mean-field simulations of the Gross-Pitaevskii equation, we shed light on the various stages of the observed dynamics. We find that the flow of the superfluid past the obstacle leads initially to the formation of a classical-like wake, which later becomes disordered. Following removal of the obstacle, the vortex number decays due to vortices annihilating and reaching the boundary. Our results are in excellent agreement with the experimental observations. Furthermore, we probe thermal effects through phenomenological dissipation., 6 pages, 6 figures

We consider an equilibrium single-species homogeneous Bose gas from which a proportion of the atoms are instantaneously and coherently transferred to a second species, thereby forming a binary Bose gas in a non-equilibrium initial state. We study the ensuing evolution towards a new equilibrium, mapping the dynamics and final equilibrium state out as a function of the population transfer and the interspecies interactions by means of classical field methods. While in certain regimes, the condensate fractions are largely unaffected by the population transfer process, in others, particularly for immiscible interactions, one or both condensate fractions are vastly reduced to a new equilibrium value., 4 figures, 9 pages

The study of vortex reconnections is an essential ingredient of understanding superfluid turbulence, a phenomenon recently also reported in trapped atomic Bose-Einstein condensates. In this work we show that, despite the established dependence of vortex motion on temperature in such systems, vortex reconnections are actually temperature independent on the typical length/time scales of atomic condensates. Our work is based on a dissipative Gross-Pitaevskii equation for the condensate, coupled to a semiclassical Boltzmann equation for the thermal cloud (the Zaremba-Nikuni-Griffin formalism). Comparison to vortex reconnections in homogeneous condensates further show reconnections to be insensitive to the inhomogeneity in the background density., 6 pages, 4 figures

The generation of isotropic vortex configurations in trapped atomic Bose-Einstein condensates offers a platform to elucidate quantum turbulence on mesoscopic scales. We demonstrate that a laser-induced obstacle moving in a figure-eight path within the condensate provides a simple and effective means to generate an isotropic three-dimensional vortex tangle due to its minimal net transfer of angular momentum to the condensate. Our characterisation of vortex structures and their isotropy is based on projected vortex lengths and velocity statistics obtained numerically via the Gross-Pitaevskii equation. Our methodology provides a possible experimental route for generating and characterising vortex tangles and quantum turbulence in atomic Bose-Einstein condensates., Comment: 5 pages, 4 figures

We show that an elliptical obstacle moving through a Bose-Einstein condensate generates wakes of quantum vortices which resemble those of classical viscous flow past a cylinder or sphere. The role of ellipticity is to facilitate the interaction of the vortices nucleated by the obstacle. Initial steady symmetric wakes lose their symmetry and form clusters of like-signed vortices, in analogy to the classical B\'enard-von K\'arm\'an vortex street. Our findings, demonstrated numerically in both two and three dimensions, confirm the intuition that a sufficiently large number of quanta of circulation reproduce classical physics., Comment: 10 pages, 5 figures

We perform a theoretical study into how dipole-dipole interactions modify the properties of superfluid vortices within the context of a two-dimensional atomic Bose gas of co-oriented dipoles. The reduced density at a vortex acts like a giant anti-dipole, changing the density profile and generating an effective dipolar potential centred at the vortex core whose most slowly decaying terms go as $1/\rho^2$ and $\ln(\rho)/\rho^3$. These effects modify the vortex-vortex interaction which, in particular, becomes anisotropic for dipoles polarized in the plane. Striking modifications to vortex-vortex dynamics are demonstrated, i.e. anisotropic co-rotation dynamics and the suppression of vortex annihilation., Comment: PRL accepted, 6 pages, 5 figures