Discrete Lie Advection of Differential Forms

In this paper, we present a numerical technique for performing Lie advection of arbitrary differential forms. Leveraging advances in high-resolution finite volume methods for scalar hyperbolic conservation laws, we first discretize the interior product (also called contraction) through integrals over Eulerian approximations of extrusions. This, along with Cartan's homotopy formula and a discrete exterior derivative, can then be used to derive a discrete Lie derivative. The usefulness of this operator is demonstrated through the numerical advection of scalar fields and 1-forms on regular grids.Comment: Accepted version; to be published in J. FoC

Coronavirus peplomer charge heterogeneity

Recent advancements in viral hydrodynamics afford the calculation of the transport properties of particle suspensions from first principles, namely, from the detailed particle shapes. For coronavirus suspensions, for example, the shape can be approximated by beading (i) the spherical capsid and (ii) the radially protruding peplomers. The general rigid bead-rod theory allows us to assign Stokesian hydrodynamics to each bead. Thus, viral hydrodynamics yields the suspension rotational diffusivity, but not without first arriving at a configuration for the cationic peplomers. Prior work considered identical peplomers charged identically. However, a recent pioneering experiment uncovers remarkable peplomer size and charge heterogeneities. In this work, we use energy minimization to arrange the spikes, charged heterogeneously to obtain the coronavirus spike configuration required for its viral hydrodynamics. For this, we use the measured charge heterogeneity. We consider 20 000 randomly generated possibilities for cationic peplomers with formal charges ranging from 30 to 55. We find the configurations from energy minimization of all of these possibilities to be nearly spherically symmetric, all slightly oblate, and we report the corresponding breadth of the dimensionless rotational diffusivity, the transport property around which coronavirus cell attachment revolves.journal articl

Coronavirus pleomorphism

The coronavirus is always idealized as a spherical capsid with radially protruding spikes. However, histologically, in the tissues of infected patients, capsids in cross section are elliptical, and only sometimes spherical [Neuman et al., “Supramolecular architecture of severe acute respiratory syndrome coronavirus revealed by electron cryomicroscopy,” J Virol, 80, 7918 (2006)]. This capsid ellipticity implies that coronaviruses are oblate or prolate or both. We call this diversity of shapes, pleomorphism. Recently, the rotational diffusivity of the spherical coronavirus in suspension was calculated, from first principles, using general rigid bead-rod theory [Kanso et al., “Coronavirus rotational diffusivity,” Phys Fluids 32, 113101 (2020)]. We did so by beading the spherical capsid and then also by replacing each of its bulbous spikes with a single bead. In this paper, we use energy minimization for the spreading of the spikes, charged identically, over the oblate or prolate capsids. We use general rigid bead-rod theory to explore the role of such coronavirus cross-sectional ellipticity on its rotational diffusivity, the transport property around which its cell attachment revolves. We learn that coronavirus ellipticity drastically decreases its rotational diffusivity, be it oblate or prolate.journal articl

Teamwork in the viscous oceanic microscale

© The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Kanso, E. A., Lopes, R. M., Strickler, J. R., Dabiri, J. O., & Costello, J. H. Teamwork in the viscous oceanic microscale. Proceedings of the National Academy of Sciences of the United States of America, 118(29), (2021): e2018193118, https://doi.org/10.1073/pnas.2018193118.Nutrient acquisition is crucial for oceanic microbes, and competitive solutions to solve this challenge have evolved among a range of unicellular protists. However, solitary solutions are not the only approach found in natural populations. A diverse array of oceanic protists form temporary or even long-lasting attachments to other protists and marine aggregates. Do these planktonic consortia provide benefits to their members? Here, we use empirical and modeling approaches to evaluate whether the relationship between a large centric diatom, Coscinodiscus wailesii, and a ciliate epibiont, Pseudovorticella coscinodisci, provides nutrient flux benefits to the host diatom. We find that fluid flows generated by ciliary beating can increase nutrient flux to a diatom cell surface four to 10 times that of a still cell without ciliate epibionts. This cosmopolitan species of diatom does not form consortia in all environments but frequently joins such consortia in nutrient-depleted waters. Our results demonstrate that symbiotic consortia provide a cooperative alternative of comparable or greater magnitude to sinking for enhancement of nutrient acquisition in challenging environments.We are grateful to Y. Garcia for help with organism sampling and sorting. E.A.K. is funded by NSF-2100209, NSF RAISE IOS-2034043 and NIH R01 HL 153622-01A1. R.M.L. is a CNPq research fellow (grant # 310642/2017-5). J.H.C. and J.O.D. are funded by Grant NSF-2100705

The motion of the 2D hydrodynamic Chaplygin sleigh in the presence of circulation

We consider the motion of a planar rigid body in a potential flow with circulation and subject to a certain nonholonomic constraint. This model is related to the design of underwater vehicles. The equations of motion admit a reduction to a 2-dimensional nonlinear system, which is integrated explicitly. We show that the reduced system comprises both asymptotic and periodic dynamics separated by a critical value of the energy, and give a complete classification of types of the motion. Then we describe the whole variety of the trajectories of the body on the plane.Comment: 25 pages, 7 figures. This article uses some introductory material from arXiv:1109.321

Measurement of the Bs0→J/ψηB_{s}^{0} \rightarrow J/\psi \eta lifetime

Using a data set corresponding to an integrated luminosity of $3 fb^{-1}$, collected by the LHCb experiment in $pp$ collisions at centre-of-mass energies of 7 and 8 TeV, the effective lifetime in the $B^0_s \rightarrow J/\psi \eta$ decay mode, $\tau_{\textrm{eff}}$, is measured to be $\tau_{\textrm{eff}} = 1.479 \pm 0.034~\textrm{(stat)} \pm 0.011 ~\textrm{(syst)}$ ps. Assuming $CP$ conservation, $\tau_{\textrm{eff}}$ corresponds to the lifetime of the light $B_s^0$ mass eigenstate. This is the first measurement of the effective lifetime in this decay mode.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2016-017.htm

Observation of two new Ξb−\Xi_b^- baryon resonances

Two structures are observed close to the kinematic threshold in the $\Xi_b^0 \pi^-$ mass spectrum in a sample of proton-proton collision data, corresponding to an integrated luminosity of 3.0 fb$^{-1}$ recorded by the LHCb experiment. In the quark model, two baryonic resonances with quark content $bds$ are expected in this mass region: the spin-parity $J^P = \frac{1}{2}^+$ and $J^P=\frac{3}{2}^+$ states, denoted $\Xi_b^{\prime -}$ and $\Xi_b^{*-}$. Interpreting the structures as these resonances, we measure the mass differences and the width of the heavier state to be $m(\Xi_b^{\prime -}) - m(\Xi_b^0) - m(\pi^{-}) = 3.653 \pm 0.018 \pm 0.006$ MeV$/c^2$, $m(\Xi_b^{*-}) - m(\Xi_b^0) - m(\pi^{-}) = 23.96 \pm 0.12 \pm 0.06$ MeV$/c^2$, $\Gamma(\Xi_b^{*-}) = 1.65 \pm 0.31 \pm 0.10$ MeV, where the first and second uncertainties are statistical and systematic, respectively. The width of the lighter state is consistent with zero, and we place an upper limit of $\Gamma(\Xi_b^{\prime -}) < 0.08$ MeV at 95% confidence level. Relative production rates of these states are also reported.Comment: 17 pages, 2 figure

Measurement of the mass and lifetime of the Ωb−\Omega_b^- baryon

A proton-proton collision data sample, corresponding to an integrated luminosity of 3 fb$^{-1}$ collected by LHCb at $\sqrt{s}=7$ and 8 TeV, is used to reconstruct $63\pm9$ $\Omega_b^-\to\Omega_c^0\pi^-$, $\Omega_c^0\to pK^-K^-\pi^+$ decays. Using the $\Xi_b^-\to\Xi_c^0\pi^-$, $\Xi_c^0\to pK^-K^-\pi^+$ decay mode for calibration, the lifetime ratio and absolute lifetime of the $\Omega_b^-$ baryon are measured to be \begin{align*} \frac{\tau_{\Omega_b^-}}{\tau_{\Xi_b^-}} &= 1.11\pm0.16\pm0.03, \\ \tau_{\Omega_b^-} &= 1.78\pm0.26\pm0.05\pm0.06~{\rm ps}, \end{align*} where the uncertainties are statistical, systematic and from the calibration mode (for $\tau_{\Omega_b^-}$ only). A measurement is also made of the mass difference, $m_{\Omega_b^-}-m_{\Xi_b^-}$, and the corresponding $\Omega_b^-$ mass, which yields \begin{align*} m_{\Omega_b^-}-m_{\Xi_b^-} &= 247.4\pm3.2\pm0.5~{\rm MeV}/c^2, \\ m_{\Omega_b^-} &= 6045.1\pm3.2\pm 0.5\pm0.6~{\rm MeV}/c^2. \end{align*} These results are consistent with previous measurements.Comment: 11 pages, 5 figures, All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2016-008.htm

Measurement of the relative rate of prompt χc0, χc1 and χc2 production at √s=7TeV

Prompt production of charmonium χc0, χc1 and χc2 mesons is studied using proton-proton collisions at the LHC at a centre-of-mass energy of √s=7TeV. The χc mesons are identified through their decay to J/ψγ, with J/ψ→μ+mu− using photons that converted in the detector. A data sample, corresponding to an integrated luminosity of 1.0fb−1 collected by the LHCb detector, is used to measure the relative prompt production rate of χc1 and χc2 in the rapidity range 2.0<y<4.5 as a function of the J/ψ transverse momentum from 3 to 20 GeV/c. First evidence for χc0 meson production at a hadron collider is also presented

Model-independent evidence for J/ψpJ/\psi p contributions to Λb0→J/ψpK−\Lambda_b^0\to J/\psi p K^- decays

The data sample of $\Lambda_b^0\to J/\psi p K^-$ decays acquired with the LHCb detector from 7 and 8~TeV $pp$ collisions, corresponding to an integrated luminosity of 3 fb$^{-1}$, is inspected for the presence of $J/\psi p$ or $J/\psi K^-$ contributions with minimal assumptions about $K^- p$ contributions. It is demonstrated at more than 9 standard deviations that $\Lambda_b^0\to J/\psi p K^-$ decays cannot be described with $K^- p$ contributions alone, and that $J/\psi p$ contributions play a dominant role in this incompatibility. These model-independent results support the previously obtained model-dependent evidence for $P_c^+\to J/\psi p$ charmonium-pentaquark states in the same data sample.Comment: 21 pages, 12 figures (including the supplemental section added at the end