@inproceedings{7479,
  abstract     = {Multi-exit architectures, in which a stack of processing layers is interleaved with early output layers, allow the processing of a test example to stop early and thus save computation time and/or energy.  In this work, we propose a new training procedure for multi-exit architectures based on the principle of knowledge distillation. The method encourage searly exits to mimic later, more accurate exits, by matching their output probabilities.
Experiments  on  CIFAR100  and  ImageNet  show  that distillation-based training significantly improves the accuracy of early exits while maintaining state-of-the-art accuracy  for  late  ones.   The  method  is  particularly  beneficial when  training  data  is  limited  and  it  allows  a  straightforward extension to semi-supervised learning,i.e. making use of unlabeled data at training time. Moreover, it takes only afew lines to implement and incurs almost no computational overhead at training time, and none at all at test time.},
  author       = {Bui Thi Mai, Phuong and Lampert, Christoph},
  booktitle    = {IEEE International Conference on Computer Vision},
  isbn         = {9781728148038},
  issn         = {15505499},
  location     = {Seoul, Korea},
  pages        = {1355--1364},
  publisher    = {IEEE},
  title        = {{Distillation-based training for multi-exit architectures}},
  doi          = {10.1109/ICCV.2019.00144},
  volume       = {2019-October},
  year         = {2019},
}

@inbook{7513,
  abstract     = {Social insects (i.e., ants, termites and the social bees and wasps) protect their colonies from disease using a combination of individual immunity and collectively performed defenses, termed social immunity. The first line of social immune defense is sanitary care, which is performed by colony members to protect their pathogen-exposed nestmates from developing an infection. If sanitary care fails and an infection becomes established, a second line of social immune defense is deployed to stop disease transmission within the colony and to protect the valuable queens, which together with the males are the reproductive individuals of the colony. Insect colonies are separated into these reproductive individuals and the sterile worker force, forming a superorganismal reproductive unit reminiscent of the differentiated germline and soma in a multicellular organism. Ultimately, the social immune response preserves the germline of the superorganism insect colony and increases overall fitness of the colony in case of disease. },
  author       = {Cremer, Sylvia and Kutzer, Megan},
  booktitle    = {Encyclopedia of Animal Behavior},
  editor       = {Choe, Jae},
  isbn         = {9780128132517},
  pages        = {747--755},
  publisher    = {Elsevier},
  title        = {{Social immunity}},
  doi          = {10.1016/B978-0-12-809633-8.90721-0},
  year         = {2019},
}

@unpublished{7524,
  abstract     = {We prove a lower bound for the free energy (per unit volume) of the two-dimensional Bose gas in the thermodynamic limit. We show that the free energy at density $\rho$ and inverse temperature $\beta$ differs from the one of the non-interacting system by the correction term $4 \pi \rho^2 |\ln a^2 \rho|^{-1} (2 - [1 - \beta_{\mathrm{c}}/\beta]_+^2)$. Here $a$ is the scattering length of the interaction potential, $[\cdot]_+ = \max\{ 0, \cdot \}$ and $\beta_{\mathrm{c}}$ is the inverse Berezinskii--Kosterlitz--Thouless critical temperature for superfluidity. The result is valid in the dilute limit
$a^2\rho \ll 1$ and if $\beta \rho \gtrsim 1$.},
  author       = {Deuchert, Andreas and Mayer, Simon and Seiringer, Robert},
  booktitle    = {arXiv:1910.03372},
  pages        = {61},
  publisher    = {ArXiv},
  title        = {{The free energy of the two-dimensional dilute Bose gas. I. Lower bound}},
  year         = {2019},
}

@inproceedings{7542,
  abstract     = {We present a novel class of convolutional neural networks (CNNs) for set functions,i.e., data indexed with the powerset of a finite set. The convolutions are derivedas linear, shift-equivariant functions for various notions of shifts on set functions.The framework is fundamentally different from graph convolutions based on theLaplacian, as it provides not one but several basic shifts, one for each element inthe ground set. Prototypical experiments with several set function classificationtasks on synthetic datasets and on datasets derived from real-world hypergraphsdemonstrate the potential of our new powerset CNNs.},
  author       = {Wendler, Chris and Alistarh, Dan-Adrian and Püschel, Markus},
  issn         = {1049-5258},
  location     = {Vancouver, Canada},
  pages        = {927--938},
  publisher    = {Neural Information Processing Systems Foundation},
  title        = {{Powerset convolutional neural networks}},
  volume       = {32},
  year         = {2019},
}

@article{7547,
  abstract     = {The BH3-only family of proteins is key for initiating apoptosis in a variety of contexts, and may also contribute to non-apoptotic cellular processes. Historically, the nematode Caenorhabditis elegans has provided a powerful system for studying and identifying conserved regulators of BH3-only proteins. In C. elegans, the BH3-only protein egl-1 is expressed during development to cell-autonomously trigger most developmental cell deaths. Here we provide evidence that egl-1 is also transcribed after development in the sensory neuron pair URX without inducing apoptosis. We used genetic screening and epistasis analysis to determine that its transcription is regulated in URX by neuronal activity and/or in parallel by orthologs of Protein Kinase G and the Salt-Inducible Kinase family. Because several BH3-only family proteins are also expressed in the adult nervous system of mammals, we suggest that studying egl-1 expression in URX may shed light on mechanisms that regulate conserved family members in higher organisms.},
  author       = {Cohn, Jesse and Dwivedi, Vivek and Valperga, Giulio and Zarate, Nicole and de Bono, Mario and Horvitz, H. Robert and Pierce, Jonathan T.},
  issn         = {2160-1836},
  journal      = {G3: Genes, Genomes, Genetics},
  number       = {11},
  pages        = {3703--3714},
  publisher    = {Genetics Society of America},
  title        = {{Activity-dependent regulation of the proapoptotic BH3-only gene egl-1 in a living neuron pair in Caenorhabditis elegans}},
  doi          = {10.1534/g3.119.400654},
  volume       = {9},
  year         = {2019},
}

@article{7548,
  abstract     = {Although the aggregation of the amyloid-β peptide (Aβ) into amyloid fibrils is a well-established hallmark of Alzheimer’s disease, the complex mechanisms linking this process to neurodegeneration are still incompletely understood. The nematode worm C. elegans is a valuable model organism through which to study these mechanisms because of its simple nervous system and its relatively short lifespan. Standard Aβ-based C. elegans models of Alzheimer’s disease are designed to study the toxic effects of the overexpression of Aβ in the muscle or nervous systems. However, the wide variety of effects associated with the tissue-level overexpression of Aβ makes it difficult to single out and study specific cellular mechanisms related to the onset of Alzheimer’s disease. Here, to better understand how to investigate the early events affecting neuronal signalling, we created a C. elegans model expressing Aβ42, the 42-residue form of Aβ, from a single-copy gene insertion in just one pair of glutamatergic sensory neurons, the BAG neurons. In behavioural assays, we found that the Aβ42-expressing animals displayed a subtle modulation of the response to CO2, compared to controls. Ca2+ imaging revealed that the BAG neurons in young Aβ42-expressing nematodes were activated more strongly than in control animals, and that neuronal activation remained intact until old age. Taken together, our results suggest that Aβ42-expression in this very subtle model of AD is sufficient to modulate the behavioural response but not strong enough to generate significant neurotoxicity, suggesting that slightly more aggressive perturbations will enable effectively studies of the links between the modulation of a physiological response and its associated neurotoxicity.},
  author       = {Sinnige, Tessa and Ciryam, Prashanth and Casford, Samuel and Dobson, Christopher M. and de Bono, Mario and Vendruscolo, Michele},
  issn         = {1932-6203},
  journal      = {PLOS ONE},
  number       = {5},
  publisher    = {Public Library of Science},
  title        = {{Expression of the amyloid-β peptide in a single pair of C. elegans sensory neurons modulates the associated behavioural response}},
  doi          = {10.1371/journal.pone.0217746},
  volume       = {14},
  year         = {2019},
}

@article{7550,
  abstract     = {We consider an optimal control problem for an abstract nonlinear dissipative evolution equation. The differential constraint is penalized by augmenting the target functional by a nonnegative global-in-time functional which is null-minimized in the evolution equation is satisfied. Different variational settings are presented, leading to the convergence of the penalization method for gradient flows, noncyclic and semimonotone flows, doubly nonlinear evolutions, and GENERIC systems. },
  author       = {Portinale, Lorenzo and Stefanelli, Ulisse},
  issn         = {1343-4373},
  journal      = {Advances in Mathematical Sciences and Applications},
  number       = {2},
  pages        = {425--447},
  publisher    = {Gakko Tosho},
  title        = {{Penalization via global functionals of optimal-control problems for dissipative evolution}},
  volume       = {28},
  year         = {2019},
}

@unpublished{7552,
  abstract     = {There is increasing evidence that protein binding to specific sites along DNA can activate the reading out of genetic information without coming into direct physical contact with the gene. There also is evidence that these distant but interacting sites are embedded in a liquid droplet of proteins which condenses out of the surrounding solution. We argue that droplet-mediated interactions can account for crucial features of gene regulation only if the droplet is poised at a non-generic point in its phase diagram. We explore a minimal model that embodies this idea, show that this model has a natural mechanism for self-tuning, and suggest direct experimental tests. },
  author       = {Bialek, William and Gregor, Thomas and Tkačik, Gašper},
  booktitle    = {arXiv:1912.08579},
  pages        = {5},
  publisher    = {ArXiv},
  title        = {{Action at a distance in transcriptional regulation}},
  year         = {2019},
}

@inproceedings{7576,
  abstract     = {We present the results of a friendly competition for formal verification of continuous and hybrid systems with nonlinear continuous dynamics. The friendly competition took place as part of the workshop Applied Verification for Continuous and Hybrid Systems (ARCH) in 2019. In this year, 6 tools Ariadne, CORA, DynIbex, Flow*, Isabelle/HOL, and JuliaReach (in alphabetic order) participated. They are applied to solve reachability analysis problems on four benchmark problems, one of them with hybrid dynamics. We do not rank the tools based on the results, but show the current status and discover the potential advantages of different tools.},
  author       = {Immler, Fabian and Althoff, Matthias and Benet, Luis and Chapoutot, Alexandre and Chen, Xin and Forets, Marcelo and Geretti, Luca and Kochdumper, Niklas and Sanders, David P. and Schilling, Christian},
  booktitle    = {EPiC Series in Computing},
  issn         = {23987340},
  location     = {Montreal, Canada},
  pages        = {41--61},
  publisher    = {EasyChair Publications},
  title        = {{ARCH-COMP19 Category Report: Continuous and hybrid systems with nonlinear dynamics}},
  doi          = {10.29007/m75b},
  volume       = {61},
  year         = {2019},
}

@inproceedings{7606,
  abstract     = {We derive a tight lower bound on equivocation (conditional entropy), or equivalently a tight upper bound on mutual information between a signal variable and channel outputs. The bound is in terms of the joint distribution of the signals and maximum a posteriori decodes (most probable signals given channel output). As part of our derivation, we describe the key properties of the distribution of signals, channel outputs and decodes, that minimizes equivocation and maximizes mutual information. This work addresses a problem in data analysis, where mutual information between signals and decodes is sometimes used to lower bound the mutual information between signals and channel outputs. Our result provides a corresponding upper bound.},
  author       = {Hledik, Michal and Sokolowski, Thomas R and Tkačik, Gašper},
  booktitle    = {IEEE Information Theory Workshop, ITW 2019},
  isbn         = {9781538669006},
  location     = {Visby, Sweden},
  publisher    = {IEEE},
  title        = {{A tight upper bound on mutual information}},
  doi          = {10.1109/ITW44776.2019.8989292},
  year         = {2019},
}

@unpublished{7627,
  abstract     = {Electrodepositing  insulating  and  insoluble  Li2O2  is  the  key  process  during  discharge  of  aprotic  Li-O2 batteries  and  determines  rate,  capacity,  and  reversibility.  Current  understanding  states  that  the partition  between  surface  adsorbed  and  solvated  LiO2  governs  whether  Li2O2  grows  as  surface  film, leading to low capacity even at low rates, or in solution, leading to particles and high capacities. Here we show that Li2O2 forms to the widest extent as particles via solution mediated LiO2 disproportionation. We  describe  a  unified  Li2O2  growth  model  that  conclusively  explains  capacity  limitations  across  the whole range of electrolytes. Deciding for particle morphology, achievable rate and capacities are species mobilities,   electrode   specific   surface   area   (determining  true   areal   rate)   and   the  concentration distribution of associated LiO2 in solution. Provided that species mobilities and surface are high, high, capacities are possible even with low-donor-number electrolytes, previously considered prototypical for low   capacity   via   surface   growth.   The   tools   for   these   insights   are   microscopy,   hydrodynamic voltammetry, a numerical reaction model, and in situ small/wide angle X-ray scattering (SAXS/WAXS). Combined with sophisticated data analysis, SAXS allows retrieving rich quantitative information from complex multi-phase systems. On a wider perspective, this SAXS method is a powerful in situ metrology with  atomic  to  sub-micron  resolution  to  study  mechanisms  in  complex  electrochemical  systems  and beyond. },
  author       = {Prehal, Christian and Samojlov, Aleksej and Nachtnebel, Manfred and Kriechbaum, Manfred and Amenitsch, Heinz and Freunberger, Stefan Alexander},
  pages        = {50},
  publisher    = {ChemRxiv},
  title        = {{A revised O2 reduction model in Li-O2 batteries as revealed by in situ small angle X-ray scattering}},
  year         = {2019},
}

@inproceedings{7639,
  abstract     = {Deep neural networks (DNNs) have become increasingly important due to their excellent empirical performance on a wide range of problems. However, regularization is generally achieved by indirect means, largely due to the complex set of functions defined by a network and the difficulty in measuring function complexity. There exists no method in the literature for additive regularization based on a norm of the function, as is classically considered in statistical learning theory. In this work, we study the tractability of function norms for deep neural networks with ReLU activations. We provide, to the best of our knowledge, the first proof in the literature of the NP-hardness of computing function norms of DNNs of 3 or more layers. We also highlight a fundamental difference between shallow and deep networks. In the light on these results, we propose a new regularization strategy based on approximate function norms, and show its efficiency on a segmentation task with a DNN.},
  author       = {Rannen-Triki, Amal and Berman, Maxim and Kolmogorov, Vladimir and Blaschko, Matthew B.},
  booktitle    = {Proceedings of the 2019 International Conference on Computer Vision Workshop},
  isbn         = {9781728150239},
  location     = {Seoul, South Korea},
  publisher    = {IEEE},
  title        = {{Function norms for neural networks}},
  doi          = {10.1109/ICCVW.2019.00097},
  year         = {2019},
}

@inproceedings{7640,
  abstract     = {We propose a new model for detecting visual relationships, such as "person riding motorcycle" or "bottle on table". This task is an important step towards comprehensive structured mage understanding, going beyond detecting individual objects. Our main novelty is a Box Attention mechanism that allows to model pairwise interactions between objects using standard object detection pipelines. The resulting model is conceptually clean, expressive and relies on well-justified training and prediction procedures. Moreover, unlike previously proposed approaches, our model does not introduce any additional complex components or hyperparameters on top of those already required by the underlying detection model. We conduct an experimental evaluation on two datasets, V-COCO and Open Images, demonstrating strong quantitative and qualitative results.},
  author       = {Kolesnikov, Alexander and Kuznetsova, Alina and Lampert, Christoph and Ferrari, Vittorio},
  booktitle    = {Proceedings of the 2019 International Conference on Computer Vision Workshop},
  isbn         = {9781728150239},
  location     = {Seoul, South Korea},
  publisher    = {IEEE},
  title        = {{Detecting visual relationships using box attention}},
  doi          = {10.1109/ICCVW.2019.00217},
  year         = {2019},
}

@article{7710,
  abstract     = {The number of human genomes being genotyped or sequenced increases exponentially and efficient haplotype estimation methods able to handle this amount of data are now required. Here we present a method, SHAPEIT4, which substantially improves upon other methods to process large genotype and high coverage sequencing datasets. It notably exhibits sub-linear running times with sample size, provides highly accurate haplotypes and allows integrating external phasing information such as large reference panels of haplotypes, collections of pre-phased variants and long sequencing reads. We provide SHAPEIT4 in an open source format and demonstrate its performance in terms of accuracy and running times on two gold standard datasets: the UK Biobank data and the Genome In A Bottle.},
  author       = {Delaneau, Olivier and Zagury, Jean-François and Robinson, Matthew Richard and Marchini, Jonathan L. and Dermitzakis, Emmanouil T.},
  issn         = {2041-1723},
  journal      = {Nature Communications},
  publisher    = {Springer Nature},
  title        = {{Accurate, scalable and integrative haplotype estimation}},
  doi          = {10.1038/s41467-019-13225-y},
  volume       = {10},
  year         = {2019},
}

@article{7711,
  abstract     = {The nature and extent of mitochondrial DNA variation in a population and how it affects traits is poorly understood. Here we resequence the mitochondrial genomes of 169 Drosophila Genetic Reference Panel lines, identifying 231 variants that stratify along 12 mitochondrial haplotypes. We identify 1,845 cases of mitonuclear allelic imbalances, thus implying that mitochondrial haplotypes are reflected in the nuclear genome. However, no major fitness effects are associated with mitonuclear imbalance, suggesting that such imbalances reflect population structure at the mitochondrial level rather than genomic incompatibilities. Although mitochondrial haplotypes have no direct impact on mitochondrial respiration, some haplotypes are associated with stress- and metabolism-related phenotypes, including food intake in males. Finally, through reciprocal swapping of mitochondrial genomes, we demonstrate that a mitochondrial haplotype associated with high food intake can rescue a low food intake phenotype. Together, our findings provide new insight into population structure at the mitochondrial level and point to the importance of incorporating mitochondrial haplotypes in genotype–phenotype relationship studies.},
  author       = {Bevers, Roel P. J. and Litovchenko, Maria and Kapopoulou, Adamandia and Braman, Virginie S. and Robinson, Matthew Richard and Auwerx, Johan and Hollis, Brian and Deplancke, Bart},
  issn         = {2522-5812},
  journal      = {Nature Metabolism},
  number       = {12},
  pages        = {1226--1242},
  publisher    = {Springer Nature},
  title        = {{Mitochondrial haplotypes affect metabolic phenotypes in the Drosophila Genetic Reference Panel}},
  doi          = {10.1038/s42255-019-0147-3},
  volume       = {1},
  year         = {2019},
}

@unpublished{7782,
  abstract     = {As genome-wide association studies (GWAS) increased in size, numerous gene-environment interactions (GxE) have been discovered, many of which however explore only one environment at a time and may suffer from statistical artefacts leading to biased interaction estimates. Here we propose a maximum likelihood method to estimate the contribution of GxE to complex traits taking into account all interacting environmental variables at the same time, without the need to measure any. This is possible because GxE induces fluctuations in the conditional trait variance, the extent of which depends on the strength of GxE. The approach can be applied to continuous outcomes and for single SNPs or genetic risk scores (GRS). Extensive simulations demonstrated that our method yields unbiased interaction estimates and excellent confidence interval coverage. We also offer a strategy to distinguish specific GxE from general heteroscedasticity (scale effects). Applying our method to 32 complex traits in the UK Biobank reveals that for body mass index (BMI) the GRSxE explains an additional 1.9% variance on top of the 5.2% GRS contribution. However, this interaction is not specific to the GRS and holds for any variable similarly correlated with BMI. On the contrary, the GRSxE interaction effect for leg impedance Embedded Image is significantly (P < 10−56) larger than it would be expected for a similarly correlated variable Embedded Image. We showed that our method could robustly detect the global contribution of GxE to complex traits, which turned out to be substantial for certain obesity measures.},
  author       = {Sulc, Jonathan and Mounier, Ninon and Günther, Felix and Winkler, Thomas and Wood, Andrew R. and Frayling, Timothy M. and Heid, Iris M. and Robinson, Matthew Richard and Kutalik, Zoltán},
  booktitle    = {bioRxiv},
  pages        = {20},
  publisher    = {Cold Spring Harbor Laboratory},
  title        = {{Maximum likelihood method quantifies the overall contribution of gene-environment interaction to continuous traits: An application to complex traits in the UK Biobank}},
  year         = {2019},
}

@unpublished{7950,
  abstract     = {The input to the token swapping problem is a graph with vertices v1, v2, . . . , vn, and n tokens with labels 1,2, . . . , n, one on each vertex.  The goal is to get token i to vertex vi for all i= 1, . . . , n using a minimum number of swaps, where a swap exchanges the tokens on the endpoints of an edge.Token swapping on a tree, also known as “sorting with a transposition tree,” is not known to be in P nor NP-complete.  We present some partial results:
1.  An optimum swap sequence may need to perform a swap on a leaf vertex that has the correct token (a “happy leaf”), disproving a conjecture of Vaughan.
2.  Any algorithm that fixes happy leaves—as all known approximation algorithms for the problem do—has approximation factor at least 4/3.  Furthermore, the two best-known 2-approximation algorithms have approximation factor exactly 2.
3.  A generalized problem—weighted coloured token swapping—is NP-complete on trees, but solvable in polynomial time on paths and stars.  In this version, tokens and  vertices  have  colours,  and  colours  have  weights.   The  goal  is  to  get  every token to a vertex of the same colour, and the cost of a swap is the sum of the weights of the two tokens involved.},
  author       = {Biniaz, Ahmad and Jain, Kshitij and Lubiw, Anna and Masárová, Zuzana and Miltzow, Tillmann and Mondal, Debajyoti and Naredla, Anurag Murty and Tkadlec, Josef and Turcotte, Alexi},
  booktitle    = {arXiv},
  title        = {{Token swapping on trees}},
  year         = {2019},
}

@article{8,
  abstract     = {Despite their different origins, Drosophila glia and hemocytes are related cell populations that provide an immune function. Drosophila hemocytes patrol the body cavity and act as macrophages outside the nervous system whereas glia originate from the neuroepithelium and provide the scavenger population of the nervous system. Drosophila glia are hence the functional orthologs of vertebrate microglia, even though the latter are cells of immune origin that subsequently move into the brain during development. Interestingly, the Drosophila immune cells within (glia) and outside the nervous system (hemocytes) require the same transcription factor Glide/Gcm for their development. This raises the issue of how do glia specifically differentiate in the nervous system and hemocytes in the procephalic mesoderm. The Repo homeodomain transcription factor and pan-glial direct target of Glide/Gcm is known to ensure glial terminal differentiation. Here we show that Repo also takes center stage in the process that discriminates between glia and hemocytes. First, Repo expression is repressed in the hemocyte anlagen by mesoderm-specific factors. Second, Repo ectopic activation in the procephalic mesoderm is sufficient to repress the expression of hemocyte-specific genes. Third, the lack of Repo triggers the expression of hemocyte markers in glia. Thus, a complex network of tissue-specific cues biases the potential of Glide/Gcm. These data allow us to revise the concept of fate determinants and help us understand the bases of cell specification. Both sexes were analyzed.SIGNIFICANCE STATEMENTDistinct cell types often require the same pioneer transcription factor, raising the issue of how does one factor trigger different fates. In Drosophila, glia and hemocytes provide a scavenger activity within and outside the nervous system, respectively. While they both require the Glide/Gcm transcription factor, glia originate from the ectoderm, hemocytes from the mesoderm. Here we show that tissue-specific factors inhibit the gliogenic potential of Glide/Gcm in the mesoderm by repressing the expression of the homeodomain protein Repo, a major glial-specific target of Glide/Gcm. Repo expression in turn inhibits the expression of hemocyte-specific genes in the nervous system. These cell-specific networks secure the establishment of the glial fate only in the nervous system and allow cell diversification.},
  author       = {Trébuchet, Guillaume and Cattenoz, Pierre B and Zsámboki, János and Mazaud, David and Siekhaus, Daria E and Fanto, Manolis and Giangrande, Angela},
  journal      = {Journal of Neuroscience},
  number       = {2},
  pages        = {238--255},
  publisher    = {Society for Neuroscience},
  title        = {{The Repo homeodomain transcription factor suppresses hematopoiesis in Drosophila and preserves the glial fate}},
  doi          = {10.1523/JNEUROSCI.1059-18.2018},
  volume       = {39},
  year         = {2019},
}

@article{80,
  abstract     = {We consider an interacting, dilute Bose gas trapped in a harmonic potential at a positive temperature. The system is analyzed in a combination of a thermodynamic and a Gross–Pitaevskii (GP) limit where the trap frequency ω, the temperature T, and the particle number N are related by N∼ (T/ ω) 3→ ∞ while the scattering length is so small that the interaction energy per particle around the center of the trap is of the same order of magnitude as the spectral gap in the trap. We prove that the difference between the canonical free energy of the interacting gas and the one of the noninteracting system can be obtained by minimizing the GP energy functional. We also prove Bose–Einstein condensation in the following sense: The one-particle density matrix of any approximate minimizer of the canonical free energy functional is to leading order given by that of the noninteracting gas but with the free condensate wavefunction replaced by the GP minimizer.},
  author       = {Deuchert, Andreas and Seiringer, Robert and Yngvason, Jakob},
  journal      = {Communications in Mathematical Physics},
  number       = {2},
  pages        = {723--776},
  publisher    = {Springer},
  title        = {{Bose–Einstein condensation in a dilute, trapped gas at positive temperature}},
  doi          = {10.1007/s00220-018-3239-0},
  volume       = {368},
  year         = {2019},
}

@article{8013,
  author       = {Currin, Christopher B. and Khoza, Phumlani N. and Antrobus, Alexander D. and Latham, Peter E. and Vogels, Tim P and Raimondo, Joseph V.},
  issn         = {1553-7358},
  journal      = {PLOS Computational Biology},
  number       = {7},
  publisher    = {Public Library of Science},
  title        = {{Think: Theory for Africa}},
  doi          = {10.1371/journal.pcbi.1007049},
  volume       = {15},
  year         = {2019},
}

