Added paper reference, changed date

paper.bib

@@ -53,6 +53,21 @@ @article{Oriordan2016
   url = {https://link.aps.org/doi/10.1103/PhysRevA.93.023609}
 }
 
+@article{Oriordan2016b,
+  title = {Topological defect dynamics of vortex lattices in Bose-Einstein condensates},
+  author = {O'Riordan, L. J. and Busch, Th.},
+  journal = {Phys. Rev. A},
+  volume = {94},
+  issue = {5},
+  pages = {053603},
+  numpages = {8},
+  year = {2016},
+  month = {Nov},
+  publisher = {American Physical Society},
+  doi = {10.1103/PhysRevA.94.053603},
+  url = {https://link.aps.org/doi/10.1103/PhysRevA.94.053603}
+}
+
 @article{Oriordan2013,
   title = {Coherent transport by adiabatic passage on atom chips},
   author = {Morgan, T. and O'Riordan, L. J. and Crowley, N. and O'Sullivan, B. and Busch, Th.},

paper.md

@@ -17,7 +17,7 @@ authors:
 affiliations:
  - name: Okinawa Institute of Science and Technology Graduate University, Onna-son, Okinawa 904-0495, Japan.
    index: 1
-date: 01 August 2033
+date: 21 September 2018
 bibliography: paper.bib
 ---
 
@@ -26,25 +26,26 @@ bibliography: paper.bib
 Bose--Eintein Condensates (BECs) are superfluid systems consisting of bosonic atoms that have been cooled and condensed into a single, macroscopic ground state.
 These systems can be created in an expermental laboratory and allow for the the exploration of many interesting physical phenomenon, such as superfluid turbulence.
 Simulations of BECs allow for new discoveries that directly mimic what can be seen in experiments and are thus highly valuable for fundamental research.
-In practice, almost all dynamics of BEC systems can be found by solving the non-linear Schrödinger equation known as the Gross-Pitaevskii Equation (GPE):
+In practice, almost all dynamics of BEC systems can be found by solving the non-linear Schrödinger equation known as the Gross--Pitaevskii Equation (GPE):
 
 $$
 \frac{\partial\Psi(\mathbf{r},t)}{\partial t} = \left( -\frac{\hbar^2}{2m} \frac{\partial}{\partial\mathbf{r}^2} + V(\mathbf{r}) + g|\Psi(\mathbf{r},t)|^2\right)\Psi(\mathbf{r},t)
 $$
 
-Where $\Psi(\mathbf{r},t)$ is the many-body wavefunction of the quantum system, $m$ is the atomic mass, $V(\mathbf{r})$ is a potential to trap the atomic system, $g = \frac{4\pi\hbar^2a_s}{m}$ is a coupling factor, and $a_s$ is the scattering length of the atomic species.
+Where $\Psi(\mathbf{r},t)$ is the many-body wavefunction of the quantum system, $m$ is the atomic mass, $V(\mathbf{r})$ is a potential to trap the atomic system, $g = \frac{4\pi\hbar^2a\_s}{m}$ is a coupling factor, and $a\_s$ is the scattering length of the atomic species.
 Here, the GPE is shown in 1 dimension, but it can easily be extended to 2 or 3 dimensions, if necessary.
 Though there are many methods to solve the GPE, one of the most straightforward is the split-operator method, which has previously been accelerated with GPU devices [@Ruf2009]; however, there are no general software packages available using this method.
 Even so, there are several software packages designed to simulate BECs with other methods, including GPELab [@Antoine2014] and the Massively Parallel Trotter-Suzuki Solver [@Wittek2013].
 
-GPUE is a GPU-based Gross-Pitaevskii Equation solver via the split-operator method for superfluid simulations of Bose-Einstein Condensates with an emphasis on vortex dynamics in 2 and 3 dimensions.
+GPUE is a GPU-based Gross-Pitaevskii Equation solver via the split-operator method for superfluid simulations of Bose--Einstein Condensates with an emphasis on vortex dynamics in 2 and 3 dimensions.
 For this purpose, GPUE provides a number of unique features:
 1. Dynamic field generation for trapping potentials and other variables on the GPU device
 2. Vortex tracking in 2D and vortex highlighting in 3D
 3. Configurable gauge fields for the generation of artificial magnetic fields and corresponding vortex distributions
+4. Vortex manipulation via direct control of the wavefunction phase
 
-All of these features are essential for usability and performance and have all been adequately described in the documentation [@documentation].
-GPUE provides a fast, robust, and accessible method to simulate superfluid physics for fundamental research in the area and has been used to simulate large vortex lattices in two dimensions [@Oriordan2016] and spatial adiabatic passage in atom chips [@Oriordan2013], along with ongoing studies on vortex turbulence in 2 dimensions and vortex structures in 3 dimensions.
+All of these features are essential for usability and performance, and have all been adequately described in the documentation [@documentation].
+GPUE provides a fast, robust, and accessible method to simulate superfluid physics for fundamental research in the area and has been used to simulate and manipulate large vortex lattices in two dimensions [@Oriordan2016, @Oriordan2016b] and spatial adiabatic passage in atom chips [@Oriordan2013], along with ongoing studies on vortex turbulence in two dimensions and vortex structures in three dimensions.
 
 # Acknowledgements
 This work has been supported by the Okinawa Institute of Science and Technology Graduate University and by JSPS KAKENHI Grant Number JP17J01488.