Showing 1–50 of 86 results for author: Mitra, D

Study of ion induced Inner Shell Ionization cross section through electron capture mechanism

Authors: Sumana Ghosh, Debasis Mitra, Soumya Chatterjee

Abstract: Electron Capture (EC) cross-section from K, L and M shells of the target atoms to the vacant K, L and M shells of the projectile ions have been calculated by deriving the accurate momentum transfer to the captured electrons for different charge states. Several correction factors like polarization correction, relativistic effects (R) of the target wave function, Coulomb-deflection factor (C) due to… ▽ More Electron Capture (EC) cross-section from K, L and M shells of the target atoms to the vacant K, L and M shells of the projectile ions have been calculated by deriving the accurate momentum transfer to the captured electrons for different charge states. Several correction factors like polarization correction, relativistic effects (R) of the target wave function, Coulomb-deflection factor (C) due to the effect of the repulsion between the projectile and the target nucleus, correction for projectile energy loss have been introduced. The mean charge state of the projectiles inside the target material has been estimated using suitable empirical models and the fractional charge state distribution has been calculated considering Lorentz distribution. Fractional distribution of charge state of the projectile ions is used to obtain the charge state contributions of the electron capture cross-sections. The effect of Simultaneous Multiple Ionization (SMI) has been considered in the theory of Direct Coulomb Ionization (DCI). The theoretically obtained total cross-sections have been compared with the experimental findings obtained from various literature. The computation scheme has been depicted through sample calculations of ionization cross-sections through electron capture mechanism. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Anisotropic Magnetized Asteroseismic Waves

Authors: B. Tripathi, Dhrubaditya Mitra

Abstract: We solve for waves in a polytropic, stratified medium with a spatially varying background magnetic field that points along a horizontal $x$-direction, and with gravity that is directed along the vertical $z$-direction. Force balance determines the magnitude of the background magnetic field, $B_0^2 \sim z^{n+1}$, where $n$ is the polytropic index. Using numerical and asymptotic methods, we deduce a… ▽ More We solve for waves in a polytropic, stratified medium with a spatially varying background magnetic field that points along a horizontal $x$-direction, and with gravity that is directed along the vertical $z$-direction. Force balance determines the magnitude of the background magnetic field, $B_0^2 \sim z^{n+1}$, where $n$ is the polytropic index. Using numerical and asymptotic methods, we deduce an accurate and explicit dispersion relation for fast pressure-driven waves: $Ω^2 \sim K\left(2m+n\right) \left[1 + (1/M_\mathrm{A})^2 (4-2γ+\cos^2θ-3\cos^4θ)/4 \right]$. Here, $Ω$ is the frequency, $K$ the wavenumber, $θ$ the angle the wave-vector makes with the background magnetic field, $M_\mathrm{A}$ the Alfvénic Mach number, and $m$ an integer representing the eigen state. Applications of our result are in magnetoseismology and nonlinear asteroseismology. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

Anomalous diffusion and effective shear modulus in a semi-solid membrane

Authors: Vikash Pandey, Dhrubaditya Mitra

Abstract: From the perspective of physical properties, the cell membrane is an exotic two-dimensional material that has a dual nature: it exhibits characteristics of fluids, i.e., lipid molecules show lateral diffusion, while also demonstrating properties of solids, evidenced by a non-zero shear modulus. We construct a model for such a $\textit{semi-solid}$ $\textit{membrane}$. Our model is a fluctuating ra… ▽ More From the perspective of physical properties, the cell membrane is an exotic two-dimensional material that has a dual nature: it exhibits characteristics of fluids, i.e., lipid molecules show lateral diffusion, while also demonstrating properties of solids, evidenced by a non-zero shear modulus. We construct a model for such a $\textit{semi-solid}$ $\textit{membrane}$. Our model is a fluctuating randomly triangulated mesh with two different kinds of nodes. The solid nodes never change their neighbors, while the fluid nodes do. As the area fraction occupied by the solid nodes ($Φ$) is increased the motion of fluid nodes transition from diffusion to localization via subdiffusion. Next, the solid nodes are pinned to mimic the pinning of the plasma membrane to the cytoskeleton. For the pinned membrane, there exists a range of $Φ$ over which the model has both a non-zero shear modulus and a non-zero lateral diffusivity. The bending modulus, measured through the spectrum of height fluctuations remains unchanged. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Water-based Quantum Dots Liquid Scintillator for Particle Physics

Authors: M. Zhao, M. Taani, J. Cole, B. Crudele, B. Zou, N. Bhuiyan, E. Chowdhury, Y. Duan, S. Fekri, D. Harvey, D. Mitra, O. Raz, A. Thompson, T. Katori, A. Rakovich

Abstract: Liquid scintillators are typically composed from organic compounds dissolved in organic solvents. However, usage of such material is often restricted due to fire safety and environmental reasons. Because of this, R\&D of water-based liquid scintillators is of extreme relevance; yet, no such scintillators have been made commercially available as yet. Here, we investigate an alternative, water-based… ▽ More Liquid scintillators are typically composed from organic compounds dissolved in organic solvents. However, usage of such material is often restricted due to fire safety and environmental reasons. Because of this, R\&D of water-based liquid scintillators is of extreme relevance; yet, no such scintillators have been made commercially available as yet. Here, we investigate an alternative, water-based quantum dots liquid scintillator. Pre-determined and controllable optical properties of the quantum dots, as well as the existence of large libraries of established protocols for their dispersion in aqueous solutions, make them an attractive option for nuclear and particle physics applications. We characterize the optical properties of water-based quantum dots liquid scintillator and find that most of its optical properties are preserved upon quantum dots' phase transfer into water, through the addition of an oleic acid hydrophilic layer. Using the developed scintillator, the time and charge responses from atmospheric muons are measured, highlighting the practical viability of water-based quantum dots liquid scintillators for nuclear and particle physics, special interest on neutrino physics. △ Less

Submitted 26 June, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 23 pages, 9 figures, 3 tables, v3, accepted by JINST for publication

Journal ref: 2024 JINST 19 P07014

Laser cooling of a fermionic molecule

Authors: Jinyu Dai, Qi Sun, Benjamin C. Riley, Debayan Mitra, Tanya Zelevinsky

Abstract: Only bosonic molecular species have been directly laser cooled to date, primarily due to an abundance of bosonic isotopes in nature and to their simpler hyperfine structure. Fermionic molecules provide new opportunities for ultracold chemistry, quantum simulation, and precision measurements. Here we report direct laser cooling of a fermionic molecular isotopologue, calcium monodeuteride (CaD). Wit… ▽ More Only bosonic molecular species have been directly laser cooled to date, primarily due to an abundance of bosonic isotopes in nature and to their simpler hyperfine structure. Fermionic molecules provide new opportunities for ultracold chemistry, quantum simulation, and precision measurements. Here we report direct laser cooling of a fermionic molecular isotopologue, calcium monodeuteride (CaD). With a nuclear spin I = 1, only 5 hyperfine states need to be addressed for rotational closure in optical cycling. These hyperfine states are unresolved for typical experimental linewidths. We present a method for efficiently producing alkaline-earth metal hydrides and deuterides. We demonstrate rotational closure and show magnetically assisted Sisyphus cooling in one dimension for a beam of CaD molecules. Our results indicate that the experimental complexity for laser cooling CaD is similar to that of calcium monohydride (CaH). Laser cooling of CaD is a promising first step for production of ultracold and trapped atomic deuterium. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 8 pages, 5 figures

Uncovering the multifractality of Lagrangian pair dispersion in shock-dominated turbulence

Authors: Sadhitro De, Dhrubaditya Mitra, Rahul Pandit

Abstract: Lagrangian pair dispersion provides insights into mixing in turbulent flows. By direct numerical simulations (DNS) we show that the statistics of pair dispersion in the randomly forced two-dimensional Burgers equation, which is a typical model of shock-dominated turbulence, is very different from its incompressible counterpart because Lagrangian particles get trapped in shocks. We develop a heuris… ▽ More Lagrangian pair dispersion provides insights into mixing in turbulent flows. By direct numerical simulations (DNS) we show that the statistics of pair dispersion in the randomly forced two-dimensional Burgers equation, which is a typical model of shock-dominated turbulence, is very different from its incompressible counterpart because Lagrangian particles get trapped in shocks. We develop a heuristic theoretical framework that accounts for this -- a generalization of the multifractal model -- whose prediction of the scaling of Lagrangian exit times agrees well with our DNS. △ Less

Submitted 12 November, 2023; originally announced November 2023.

Intermittency in the not-so-smooth elastic turbulence

Authors: Rahul K. Singh, Prasad Perlekar, Dhrubaditya Mitra, Marco E. Rosti

Abstract: Elastic turbulence is the chaotic fluid motion resulting from elastic instabilities due to the addition of polymers in small concentrations at very small Reynolds ($\mbox{Re}$) numbers. Our direct numerical simulations show that elastic turbulence, though a low $\mbox{Re}$ phenomenon, has more in common with classical, Newtonian turbulence than previously thought. In particular, we find power-law… ▽ More Elastic turbulence is the chaotic fluid motion resulting from elastic instabilities due to the addition of polymers in small concentrations at very small Reynolds ($\mbox{Re}$) numbers. Our direct numerical simulations show that elastic turbulence, though a low $\mbox{Re}$ phenomenon, has more in common with classical, Newtonian turbulence than previously thought. In particular, we find power-law spectra for kinetic energy $E(k) \sim k^{-4}$ and polymeric energy $E_{\rm p}(k) \sim k^{-3/2}$, independent of the Deborah ($\mbox{De}$) number. This is further supported by calculation of scale-by-scale energy budget which shows a balance between the viscous term and the polymeric term in the momentum equation. In real space, as expected, the velocity field is smooth, i.e., the velocity difference across a length scale $r$, $δu \sim r$ but, crucially, with a non-trivial sub-leading contribution $r^{3/2}$ which we extract by using the second difference of velocity. The structure functions of second difference of velocity up to order $6$ show clear evidence of intermittency/multifractality. We provide additional evidence in support of this intermittent nature by calculating moments of rate of dissipation of kinetic energy averaged over a ball of radius $r$, $\varepsilon_{r}$, from which we compute the multifractal spectrum. △ Less

Submitted 23 April, 2024; v1 submitted 14 August, 2023; originally announced August 2023.

Acoustic Kerr Metric in Analogue Gravity

Authors: Deeshani Mitra, Surojit Dalui, Subir Ghosh, Arpan Krishna Mitra

Abstract: The present paper is based on a previous work (involving two of the present authors) where a generalized fluid dynamical model was proposed. The underlying symplectic structure of the Lagrangian discrete degrees of freedom obeyed a Non-Commutative algebra, generated by Berry curvature correction. In an Euler (or Hamiltonian) framework, this is manifested as an extended algebra between the fluid va… ▽ More The present paper is based on a previous work (involving two of the present authors) where a generalized fluid dynamical model was proposed. The underlying symplectic structure of the Lagrangian discrete degrees of freedom obeyed a Non-Commutative algebra, generated by Berry curvature correction. In an Euler (or Hamiltonian) framework, this is manifested as an extended algebra between the fluid variables, leading to the extended fluid model. Here we study the dynamics of sonic fluctuations that live in this effective analogue gravity spacetime. Interestingly enough, the effective metric resembles that of a spinning Black Hole; the spin is induced by the underlying Non-Commutative structure. The effective mass and spin parameters of the Black Hole, in terms of fluid parameters, are also identified. The connection of our model with anomalous Hall systems may lead to observable signatures of the analogue black hole in physical systems. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 7 pages, no figure

Probing the limits of optical cycling in a predissociative diatomic molecule

Authors: Qi Sun, Claire E. Dickerson, Jinyu Dai, Isaac M. Pope, Lan Cheng, Daniel Neuhauser, Anastassia N. Alexandrova, Debayan Mitra, Tanya Zelevinsky

Abstract: Molecular predissociation is the spontaneous, nonradiative bond breaking process that can occur upon excitation. In the context of laser cooling, predissociation is an unwanted consequence of molecular structure that limits the ability to scatter a large number of photons required to reach the ultracold regime. Unlike rovibrational branching, predissociation is irreversible since the fragments fly… ▽ More Molecular predissociation is the spontaneous, nonradiative bond breaking process that can occur upon excitation. In the context of laser cooling, predissociation is an unwanted consequence of molecular structure that limits the ability to scatter a large number of photons required to reach the ultracold regime. Unlike rovibrational branching, predissociation is irreversible since the fragments fly apart with high kinetic energy. Of particular interest is the simple diatomic molecule, CaH, for which the two lowest electronically excited states used in laser cooling lie above the dissociation threshold of the ground potential. In this work, we present measurements and calculations that quantify the predissociation probabilities affecting the cooling cycle. The results allow us to design a laser cooling scheme that will enable the creation of an ultracold and optically trapped cloud of CaH molecules. In addition, we use the results to propose a two-photon pathway to controlled dissociation of the molecules, in order to gain access to their ultracold fragments, including hydrogen. △ Less

Submitted 20 October, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 16 pages, 4 figures

Journal ref: Phys. Rev. Research 5, 043070 (2023)

Flow driven interfacial waves: an inviscid asymptotic study

Authors: A. F. Bonfils, Dhrubaditya Mitra, W. Moon, J. S. Wettlaufer

Abstract: Motivated by wind blowing over water, we use asymptotic methods to study the evolution of short wavelength interfacial waves driven by the combined action of these flows. We solve the Rayleigh equation for the stability of the shear flow, and construct a uniformly valid approximation for the perturbed streamfunction, or eigenfunction. We then expand the real part of the eigenvalue, the phase speed… ▽ More Motivated by wind blowing over water, we use asymptotic methods to study the evolution of short wavelength interfacial waves driven by the combined action of these flows. We solve the Rayleigh equation for the stability of the shear flow, and construct a uniformly valid approximation for the perturbed streamfunction, or eigenfunction. We then expand the real part of the eigenvalue, the phase speed, in a power series of the inverse wavenumber and show that the imaginary part is exponentially small. We give expressions for the growth rates of the Miles (1957) and rippling (e.g., Young & Wolfe 2013) instabilities that are valid for an arbitrary shear flow. The accuracy of the results is demonstrated by a comparison with the exact solution of the eigenvalue problem in the case when both the wind and the current have an exponential profile. △ Less

Submitted 25 October, 2023; v1 submitted 5 November, 2022; originally announced November 2022.

Journal ref: Journal of Fluid Mechanics,976,A19(2023)

Theoretical studies on multiple ionisation and electron capture processes in heavy ion induced M-shell ionisation

Authors: Soumya Chatterjee, Sumana Ghosh, D. Mitra

Abstract: Multiple ionisation and electron capture are found to be vital mechanisms for K and L x-ray emissions along with the direct coulomb ionisation in heavy ion-atom collisions. Naturally, these two mechanisms may also be significant for M x-ray emissions also. However, these mechanisms associating with the M-shell are highly complex and not yet studied convincingly. We, in this work, have discussed ab… ▽ More Multiple ionisation and electron capture are found to be vital mechanisms for K and L x-ray emissions along with the direct coulomb ionisation in heavy ion-atom collisions. Naturally, these two mechanisms may also be significant for M x-ray emissions also. However, these mechanisms associating with the M-shell are highly complex and not yet studied convincingly. We, in this work, have discussed about theoretical techniques to study these mechanisms in more general way, which in turn, almost resolved the wide gap between the theories and experiments regarding M-shell prevalent in literature up to a good extent. We have justified this technique with a few representative cases; for example, the silicon and sulphur projectile ions colliding with the gold and bismuth targets at beam energies ranging from 5 to 10 MeV. Here, we see that the capture contribution is much higher than that coming from the multiple ionisation effect. The combined effects of direct coulomb ionisation associated with simultaneous multiple ionisation and electron capture give a better agreement with the measurements up to a good extent. △ Less

Submitted 3 July, 2024; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: will be published after major modification

Estimate of entropy generation rate can spatiotemporally resolve the active nature of cell flickering

Authors: Sreekanth K Manikandan, Tanmoy Ghosh, Tithi Mandal, Arikta Biswas, Bidisha Sinha, Dhrubaditya Mitra

Abstract: We use the short-time inference scheme (Manikandan, Gupta and Krishnamurthy, Phys. Rev. Lett. 124, 120603, 2020), obtained within the framework of stochastic thermodynamics, to infer a lower-bound to entropy generation rate from flickering data generated by Interference Reflection Microscopy of HeLA cells. We can clearly distinguish active cell membranes from their ATP depleted selves and even spa… ▽ More We use the short-time inference scheme (Manikandan, Gupta and Krishnamurthy, Phys. Rev. Lett. 124, 120603, 2020), obtained within the framework of stochastic thermodynamics, to infer a lower-bound to entropy generation rate from flickering data generated by Interference Reflection Microscopy of HeLA cells. We can clearly distinguish active cell membranes from their ATP depleted selves and even spatio-temporally resolve activity down to the scale of about one $μ$m. Our estimate of activity is model--independent. △ Less

Submitted 26 June, 2024; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: Published version

Report number: NORDITA 2022-034

Journal ref: Phys. Rev. Research 6, 023310 (2024)

Dynamic multiscaling in stochastically forced Burgers turbulence

Authors: Sadhitro De, Dhrubaditya Mitra, Rahul Pandit

Abstract: We carry out a detailed study of dynamic multiscaling in the turbulent nonequilibrium, but statistically steady, state of the stochastically forced one-dimensional Burgers equation. We introduce the concept of $\textit{interval collapse times}$ $τ_{\rm col}$, the time taken for an interval of length $\ell$, demarcated by a pair of Lagrangian tracers, to collapse at a shock. By calculating the dyna… ▽ More We carry out a detailed study of dynamic multiscaling in the turbulent nonequilibrium, but statistically steady, state of the stochastically forced one-dimensional Burgers equation. We introduce the concept of $\textit{interval collapse times}$ $τ_{\rm col}$, the time taken for an interval of length $\ell$, demarcated by a pair of Lagrangian tracers, to collapse at a shock. By calculating the dynamic scaling exponent of the order-$p$ moment of $τ_{\rm col}$, we show that (a) there is $\textit{not one but an infinity of characteristic time scales}$ and (b) the probability distribution function of $τ_{\rm col}$ is non-Gaussian and has a power-law tail. Our study is based on (a) a theoretical framework that allows us to obtain dynamic-multiscaling exponents analytically, (b) extensive direct numerical simulations, and (c) a careful comparison of the results of (a) and (b). We discuss possible generalizations of our work to dimensions $d >1 $, for the stochastically forced Burgers equation, and to other compressible flows that exhibit turbulence with shocks. △ Less

Submitted 7 June, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Quantum control of molecules for fundamental physics

Authors: D. Mitra, K. H. Leung, T. Zelevinsky

Abstract: The extraordinary success in laser cooling, trapping, and coherent manipulation of atoms has energized the efforts in extending this exquisite control to molecules. Not only are molecules ubiquitous in nature, but the control of their quantum states offers unparalleled access to fundamental constants and possible physics beyond the Standard Model. Quantum state manipulation of molecules can enable… ▽ More The extraordinary success in laser cooling, trapping, and coherent manipulation of atoms has energized the efforts in extending this exquisite control to molecules. Not only are molecules ubiquitous in nature, but the control of their quantum states offers unparalleled access to fundamental constants and possible physics beyond the Standard Model. Quantum state manipulation of molecules can enable high-precision measurements including tests of fundamental symmetries and searches for new particles and fields. At the same time, their complex internal structure presents experimental challenges to overcome in order to gain sensitivity to new physics. In this Perspective, we review recent developments in this thriving new field. Moreover, throughout the text we discuss many current and future research directions that have the potential to place molecules at the forefront of fundamental science. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Comments: Perspective article

Journal ref: Phys. Rev. A 105, 040101 (2022)

Kolmogorov Turbulence Coexists with Pseudo-Turbulence in Buoyancy-Driven Bubbly Flows

Authors: Vikash Pandey, Dhrubaditya Mitra, Prasad Perlekar

Abstract: We investigate spectral properties of buoyancy driven bubbly flows. Using high-resolution numerical simulations and phenomenology of homogeneous turbulence, we identify the relevant energy transfer mechanisms. We find: (a) At high enough Galilei number (ratio of the buoyancy to viscous forces) the kinetic energy spectrum shows the Kolmogorov scaling with a power law exponent $-5/3$ for the range o… ▽ More We investigate spectral properties of buoyancy driven bubbly flows. Using high-resolution numerical simulations and phenomenology of homogeneous turbulence, we identify the relevant energy transfer mechanisms. We find: (a) At high enough Galilei number (ratio of the buoyancy to viscous forces) the kinetic energy spectrum shows the Kolmogorov scaling with a power law exponent $-5/3$ for the range of scales between the bubble diameter and the dissipation scale ($η$). (b) For scales smaller than $η$, the physics of pseudo-turbulence is recovered. △ Less

Submitted 23 August, 2023; v1 submitted 9 April, 2022; originally announced April 2022.

MeMC: A package for monte-carlo simulations of spherical shells

Authors: Vipin Agrawal, Vikash Pandey, Hanna Kylhammar, Apurba Dev, Dhrubaditya Mitra

Abstract: The MeMC is an open-source software package for monte-carlo simulation of elastic shells. It is designed as a tool to interpret the force-distance data generated by indentation of biological nano-vesicles by atomic force microscopes. The code is written in c++ and python. The code is customizable -- new modules can be added in a straightforward manner. The MeMC is an open-source software package for monte-carlo simulation of elastic shells. It is designed as a tool to interpret the force-distance data generated by indentation of biological nano-vesicles by atomic force microscopes. The code is written in c++ and python. The code is customizable -- new modules can be added in a straightforward manner. △ Less

Submitted 25 April, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 6 pages,4 figures

Direct laser cooling of calcium monohydride molecules

Authors: S. F. Vázquez-Carson, Q. Sun, J. Dai, D. Mitra, T. Zelevinsky

Abstract: We demonstrate optical cycling and sub-Doppler laser cooling of a cryogenic buffer-gas beam of calcium monohydride (CaH) molecules. We measure vibrational branching ratios for laser cooling transitions for both excited electronic states A and B. We measure further that repeated photon scattering via the $A\leftarrow X$ transition is achievable at a rate of $\sim$ $1.6\times10^6$ photons/s and demo… ▽ More We demonstrate optical cycling and sub-Doppler laser cooling of a cryogenic buffer-gas beam of calcium monohydride (CaH) molecules. We measure vibrational branching ratios for laser cooling transitions for both excited electronic states A and B. We measure further that repeated photon scattering via the $A\leftarrow X$ transition is achievable at a rate of $\sim$ $1.6\times10^6$ photons/s and demonstrate the interaction-time limited scattering of $\sim$ $200$ photons by repumping the largest vibrational decay channel. We also demonstrate the ability to sub-Doppler cool a molecular beam of CaH through the magnetically assisted Sisyphus effect. Using a standing wave of light, we lower the molecular beam's transverse temperature from $12.2(1.2)$ mK to $5.7(1.1)$ mK. We compare these results to sub-Doppler forces modeled using optical Bloch equations and Monte Carlo simulations of the molecular beam trajectories. This work establishes a clear pathway for creating a magneto-optical trap (MOT) of CaH molecules. Such a MOT could serve as a starting point for production of ultracold hydrogen gas via dissociation of a trapped CaH cloud. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 11 pages, 4 figures, 3 tables

Functionalizing Aromatic Compounds with Optical Cycling Centers

Authors: Guo-Zhu Zhu, Debayan Mitra, Benjamin L. Augenbraun, Claire E. Dickerson, Michael J. Frim, Guanming Lao, Zack D. Lasner, Anastassia N. Alexandrova, Wesley C. Campbell, Justin R. Caram, John M. Doyle, Eric R. Hudson

Abstract: Molecular design principles provide guidelines for augmenting a molecule with a smaller group of atoms to realize a desired property or function. We demonstrate that these concepts can be used to create an optical cycling center that can be attached to a number of aromatic ligands, allowing the scattering of many photons from the resulting molecules without changing the molecular vibrational state… ▽ More Molecular design principles provide guidelines for augmenting a molecule with a smaller group of atoms to realize a desired property or function. We demonstrate that these concepts can be used to create an optical cycling center that can be attached to a number of aromatic ligands, allowing the scattering of many photons from the resulting molecules without changing the molecular vibrational states. We provide further design principles that indicate the ability to expand this work. This represents a significant step towards a quantum functional group, which may serve as a generic qubit moiety that can be attached to a wide range of molecular structures and surfaces. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Comments: 6 pages, 3 figures

Journal ref: Nat.Chem.14, 995-999 (2022)

Pathway Towards Optical Cycling and Laser Cooling of Functionalized Arenes

Authors: Debayan Mitra, Zack D. Lasner, Guo-Zhu Zhu, Claire E. Dickerson, Benjamin L. Augenbraun, Austin D. Bailey, Anastassia N. Alexandrova, Wesley C. Campbell, Justin R. Caram, Eric R. Hudson, John M. Doyle

Abstract: Rapid and repeated photon cycling has enabled precision metrology and the development of quantum information systems using a variety of atoms and simple molecules. Extending optical cycling to structurally complex molecules would provide new capabilities in these areas, as well as in ultracold chemistry. Increased molecular complexity, however, makes realizing closed optical transitions more diffi… ▽ More Rapid and repeated photon cycling has enabled precision metrology and the development of quantum information systems using a variety of atoms and simple molecules. Extending optical cycling to structurally complex molecules would provide new capabilities in these areas, as well as in ultracold chemistry. Increased molecular complexity, however, makes realizing closed optical transitions more difficult. Building on the already established strong optical cycling of diatomic, linear triatomic, and symmetric top molecules, recent theoretical and experimental work has indicated that cycling will be extendable to phenol containing molecules, as well as other asymmetric species. The paradigm for these systems is the use of an optical cycling center bonded to a molecular ligand. Theory has suggested that cycling may be extended to even larger ligands, like naphthalene, pyrene and coronene. Here, we study the optical excitation and vibrational branching of the molecules CaO-2-naphthyl, SrO-2-naphthyl and CaO-1-naphthyl and find only weak decay to excited vibrational states, indicating a promising path to full quantum control and laser cooling of large arene-based molecules. △ Less

Submitted 9 March, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 9 pages, 3 figures, 2 tables

Magneto-Optical Trapping and Sub-Doppler Cooling of a Polyatomic Molecule

Authors: Nathaniel B. Vilas, Christian Hallas, Loïc Anderegg, Paige Robichaud, Andrew Winnicki, Debayan Mitra, John M. Doyle

Abstract: We report magneto-optical trapping (MOT) of a polyatomic molecule, calcium monohydroxide (CaOH). The MOT contains $2.0(5)\times 10^4$ CaOH molecules at a peak density of $3.0(8)\times10^{6}$ cm$^{-3}$. CaOH molecules are further sub-Doppler laser cooled in an optical molasses, to a temperature of 110(4) $μ$K. The temperatures and densities achieved here make CaOH a viable candidate for a wide vari… ▽ More We report magneto-optical trapping (MOT) of a polyatomic molecule, calcium monohydroxide (CaOH). The MOT contains $2.0(5)\times 10^4$ CaOH molecules at a peak density of $3.0(8)\times10^{6}$ cm$^{-3}$. CaOH molecules are further sub-Doppler laser cooled in an optical molasses, to a temperature of 110(4) $μ$K. The temperatures and densities achieved here make CaOH a viable candidate for a wide variety of quantum science applications, including the creation of optical tweezer arrays of CaOH molecules. This work also suggests that laser cooling and magneto-optical trapping of many other polyatomic species will be both feasible and practical. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 6 pages, 4 figures

Large is different: non-monotonic behaviour of elastic range scaling in polymeric turbulence at large Reynolds and Deborah numbers

Authors: Marco E Rosti, Prasad Perlekar, Dhrubaditya Mitra

Abstract: We use direct numerical simulations to study homogeneous, and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wavenumbers $k$, the kinetic energy spectrum shows Kolmogorov--like behavior which crosses over at a larger $k$ to a novel, elastic scaling regime, $E(k) \sim k^{-ξ}$, with $ξ\approx 2.3$. We study the contribution of the p… ▽ More We use direct numerical simulations to study homogeneous, and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wavenumbers $k$, the kinetic energy spectrum shows Kolmogorov--like behavior which crosses over at a larger $k$ to a novel, elastic scaling regime, $E(k) \sim k^{-ξ}$, with $ξ\approx 2.3$. We study the contribution of the polymers to the flux of kinetic energy through scales, and find that it can be decomposed into two parts: one increase in effective viscous dissipation, and a purely elastic contribution that dominates over the nonlinear flux in the range of $k$ over which the elastic scaling is observed. The multiscale balance between the two fluxes determines the crossover wavenumber which depends non-monotically on the Deborah number. Consistently, structure functions also show two scaling ranges, with intermittency present in both of them in equal measure. △ Less

Submitted 12 February, 2023; v1 submitted 22 November, 2021; originally announced November 2021.

Rate of formation of caustics in heavy particles advected by turbulence

Authors: Akshay Bhatnagar, Vikash Pandey, Prasad Perlekar, Dhrubaditya Mitra

Abstract: The rate of collision and the relative velocities of the colliding particles in turbulent flows is a crucial part of several natural phenomena, e.g., rain formation in warm clouds and planetesimal formation in a protoplanetary disks. The particles are often modeled as passive, but heavy and inertial. Within this model, large relative velocities emerge due to formation of singularities (caustics) o… ▽ More The rate of collision and the relative velocities of the colliding particles in turbulent flows is a crucial part of several natural phenomena, e.g., rain formation in warm clouds and planetesimal formation in a protoplanetary disks. The particles are often modeled as passive, but heavy and inertial. Within this model, large relative velocities emerge due to formation of singularities (caustics) of in the gradient matrix of the velocities of the particles. Using extensive direct numerical simulations of heavy particles in both two (direct and inverse cascade) and three dimensional turbulent flows we calculate the rate of formation of caustics, $J$ as a function of the Stokes number (${\rm St}$).The best approximation to our data is $J \sim \exp(-C/{\rm St})$, in the limit ${\rm St} \to 0 $ where $C$ is a non-universal constant. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: Consistent with the version accepted for publication in Philosophical Transactions of the Royal Society A

Zeptosecond dynamics in atoms: fact or fiction?

Authors: T. Nandi, Prashant Sharma, Soumya Chatterjee, D. Mitra, Adya P Mishra, Y. Azuma, F. Koike, A. S. Kheifets

Abstract: Photon exchange due to nuclear bremsstrahlung during nuclear collisions can cause Coulomb excitation in the projectile and the target nuclei. The corresponding process originated in nuclear timescales can also be observed in atomic phenomenon experimentally if it delayed by at least with an attosecond or longer timescales. We have found that this happens due to a mechanism involving the Eisenbud-W… ▽ More Photon exchange due to nuclear bremsstrahlung during nuclear collisions can cause Coulomb excitation in the projectile and the target nuclei. The corresponding process originated in nuclear timescales can also be observed in atomic phenomenon experimentally if it delayed by at least with an attosecond or longer timescales. We have found that this happens due to a mechanism involving the Eisenbud-Wigner-Smith time delay process. We have estimated photoionization time delays in atomic collisions utilizing the nonrelativistic version of random phase approximation with exchange and Hartree-Fock methods. We present three representative processes in which we can observe the phenomena in attosecond timescales even though they originate from excitations in the zeptosecond timescales. Thus the work represents an investigation of parallels between two neighboring areas of physics. Furthermore the present work suggests new possibilities for atomic physics research near the Coulomb barrier energy, where the laser is replaced by nuclear bremsstrahlung. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 6 pages, 4 figures

Asymptotic interpretation of the Miles mechanism of wind-wave instability

Authors: A. F. Bonfils, Dhrubaditya Mitra, W. Moon, J. S. Wettlaufer

Abstract: When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because o… ▽ More When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air-water density ratio, $ρ_{\rm{a}}/ρ_{\rm{w}}\equivε<<1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $ε$, the eigenvalue has a finite imaginary part, which implies growth. Miles (1957) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically. △ Less

Submitted 25 June, 2022; v1 submitted 14 July, 2021; originally announced July 2021.

Journal ref: Journal of Fluid Mechanics,944,A8(2021)

Understanding the L-shell ionization mechanism through osmium atoms bombarded by 4-6 MeV/u fluorine ions

Authors: Soumya Chatterjee, Sunil Kumar, Sarvesh Kumar, M Oswal, Biraja Mohanty, D. Mehta, D. Mitra, A. M. P. Mendez, D. M. Mitnik, C. C. Montanari, L. Sarkadi, T. Nandi

Abstract: The L-subshell ionization mechanism is studied in an ultra-thin osmium target bombarded by 4-6 MeV/u fluorine ions. Multiple ionization effects in the collisions are considered through the change of fluorescence and Coster-Kronig yields while determining L-subshell ionization cross sections from L-line x-ray production cross sections. The L-subshell ionization, as well as L-shell x-ray production… ▽ More The L-subshell ionization mechanism is studied in an ultra-thin osmium target bombarded by 4-6 MeV/u fluorine ions. Multiple ionization effects in the collisions are considered through the change of fluorescence and Coster-Kronig yields while determining L-subshell ionization cross sections from L-line x-ray production cross sections. The L-subshell ionization, as well as L-shell x-ray production cross sections so obtained, are compared with various theoretical approximations. The Coulomb direct ionization contributions is studied by (i) the relativistic semi-classical approximations (RSCA), (ii) the shellwise local plasma approximation (SLPA), and (iii) the ECUSAR theory, along with the inclusion of the vacancy sharing among the subshells by the coupled-states model (CSM) and the electron capture (EC) by a standard formalism. We find that the ECUSAR-CSM-EC describes the measured excitation function curves the best. However, the theoretical calculations are still about a factor of two smaller than the measured values. Such differences are resolved by re-evaluating the fluorescence and the Coster-Kronig yields. This work demonstrates that, in the present energy range, the heavy-ion induced inner-shell ionization of heavy atoms can be understood by combining the basic mechanisms of the direct Coulomb ionization, the electron capture, the multiple ionization, and the vacancy sharing among subshells, together with optimized atomic parameters. △ Less

Submitted 13 June, 2021; originally announced June 2021.

Comments: 11 pages, 10 figures

K-shell ionization cross sections of Cu, Zn and Ge by 3-5 MeV/U Si-ion bombardment

Authors: Shashank Singh, Mumtaz Oswal, Sunil Kumar, K. P. Singh, D. Mitra, T. Nandi

Abstract: The K x-ray spectra of different targets (Cu, Zn, and Ge) induced by 3 to 5 MeV/u Si projectile ions have been measured to determine the K-shell ionization cross-section. A significant difference is observed between the measurements and theoretical estimates, where the theoretical ones are about 28-35% of the experimental results. Such difference is reduced to a good extent 51-56% if multiple ioni… ▽ More The K x-ray spectra of different targets (Cu, Zn, and Ge) induced by 3 to 5 MeV/u Si projectile ions have been measured to determine the K-shell ionization cross-section. A significant difference is observed between the measurements and theoretical estimates, where the theoretical ones are about 28-35% of the experimental results. Such difference is reduced to a good extent 51-56% if multiple ionization effects are taken into account. Remaining discrepancy may be attributed to the electron capture contribution. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 7 pages, 4 figures. arXiv admin note: substantial text overlap with arXiv:2103.08299

Fast and High-Yield Loading of a D$_2$ MOT of Potassium from a Cryogenic Buffer Gas Beam

Authors: Zack Lasner, Debayan Mitra, Maryam Hiradfar, Benjamin Augenbraun, Lawrence Cheuk, Eunice Lee, Sridhar Prabhu, John Doyle

Abstract: We demonstrate the direct loading of a D$_2$ MOT of potassium-39 atoms from a cryogenic buffer gas beam source. We load $10^8$ atoms in a 10 ms pulse, with no degradation in performance up to a 10 Hz repetition rate. Observed densities reach $\sim10^{11}$ atoms/cm$^3$ in a single pulse, achieved with no sub-Doppler cooling or transverse compression. This system produces an ideal starting point for… ▽ More We demonstrate the direct loading of a D$_2$ MOT of potassium-39 atoms from a cryogenic buffer gas beam source. We load $10^8$ atoms in a 10 ms pulse, with no degradation in performance up to a 10 Hz repetition rate. Observed densities reach $\sim10^{11}$ atoms/cm$^3$ in a single pulse, achieved with no sub-Doppler cooling or transverse compression. This system produces an ideal starting point for ultracold atom experiments where high experimental repetition rates are desirable and initial high densities are critical. Extension to other atomic species (e.g., refractory metals) that present technical challenges to high-yield production using oven-based sources is straightforward. △ Less

Submitted 8 June, 2021; v1 submitted 31 May, 2021; originally announced May 2021.

Comments: 7 pages, 4 figures

Turbulence modulation in buoyancy-driven bubbly flows

Authors: Vikash Pandey, Dhrubaditya Mitra, Prasad Perlekar

Abstract: We present a Direct Numerical Simulation (DNS) study of buoyancy-driven bubbly flows in the presence of large scale driving that generates turbulence. On increasing the turbulence intensity: (a) the bubble trajectories become more curved, and (b) the average rise velocity of the bubbles decreases. We find that the energy spectrum of the flow shows a pseudo-turbulence scaling for length scales smal… ▽ More We present a Direct Numerical Simulation (DNS) study of buoyancy-driven bubbly flows in the presence of large scale driving that generates turbulence. On increasing the turbulence intensity: (a) the bubble trajectories become more curved, and (b) the average rise velocity of the bubbles decreases. We find that the energy spectrum of the flow shows a pseudo-turbulence scaling for length scales smaller than the bubble diameter and a Kolmogorov scaling for scales larger than the bubble diameter. We conduct a scale-by-scale energy budget analysis to understand the scaling behaviour observed in the spectrum. Although our bubbles are weakly buoyant, the statistical properties of our DNS are consistent with the experiments that investigate turbulence modulation by air bubbles in water. △ Less

Submitted 21 November, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Journal ref: J. Fluid Mech. 932 (2022) A19

Exploring a novel model for projectile charge state distribution inside a solid-target

Authors: Soumya Chatterjee, Prashant Sharma, D. Mitra, T. Nandi

Abstract: For the first time, we report a theoretical methodology to predict charge state distribution of projectile ions inside a solid-target. The method utilizes either a simple Fermi gas model or an ab initio theoretical method and a certain parameterization of width for the Lorentzian charge state distributions. Results obtained from the two approaches are comparable, but the former has a certain edge… ▽ More For the first time, we report a theoretical methodology to predict charge state distribution of projectile ions inside a solid-target. The method utilizes either a simple Fermi gas model or an ab initio theoretical method and a certain parameterization of width for the Lorentzian charge state distributions. Results obtained from the two approaches are comparable, but the former has a certain edge over the latter. The projectile charge state distribution inside a solid-target plays a significant role in estimating electron capture cross-sections and then to describe the observed K-shell ionization dynamics. The electron capture process plays a certain role in L-shell ionization dynamics too, but in a test case of Si on Au target the subshell charge sharing contributes a more vital role than the electron capture. Thus, we have validated the present model as a reliable as well as useful for many solid-target based applications viz. tumour therapy, biophysics, accelerators, material science etc. △ Less

Submitted 5 April, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 7 pages, 3 figures. arXiv admin note: substantial text overlap with arXiv:2103.08299

Significance of high charge state of projectile ions inside the target and its role on electron capture leading to target ionization phenomenon

Authors: Soumya Chatterjee, Prashant Sharma, Montanary Claudia, D. Mitra, T. Nandi

Abstract: The K x-ray spectra of different targets (Cu, Zn, and Ge) induced by 3 to 5 MeV/u Si projectile ions have been measured to determine the K-shell ionization cross-section. A significant difference is observed between the measurements and theoretical estimates, with the latter being about 50% below the experimental results. This underestimation} is attributed to the charge-exchange from target K-she… ▽ More The K x-ray spectra of different targets (Cu, Zn, and Ge) induced by 3 to 5 MeV/u Si projectile ions have been measured to determine the K-shell ionization cross-section. A significant difference is observed between the measurements and theoretical estimates, with the latter being about 50% below the experimental results. This underestimation} is attributed to the charge-exchange from target K-shell to projectile K- and L-shells. Such observation can only be possible if the projectile ions attain up to H- and He-like charge states. Corresponding projectile charge state fractions have been evaluated from the Lorentzian charge state distribution, where mean charge state is taken from the Fermi gas model [Phys. Rev. Lett. 30, 358 (1973)] and width from the Novikov and Teplova approach [Phys. Lett. A378, 1286-1289 (2014)]. The sum of the direct ionization cross-section and K-K + K-L capture cross-sections gives a good agreement with the measured cross-sections. Furthermore, we have validated this methodology with available data for Si-ion on Ti target. Such results may be useful in many solid target based applications. △ Less

Submitted 17 June, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Comments: 13 pages, 12 figures

Journal ref: Phys. Rev. A 104, 022810 (2021)

Evaluation of accurate uncertainty of measurement in L subshell ionization cross-section

Authors: Shashank Singh, Soumya Chatterjee, Debasis Mitra, Tapan Nandi

Abstract: To have a better understanding of a physical process, a comparison of ex-perimental data with theoretical values is mandatory. The comparison ismeaningful if the uncertainty in the experiment is accounted well. However,it is seldom seen, especially for a complex phenomenon. We take a test casethrough L subshell ionization of atoms by particle impact. Experimentally,x-ray production cross-sections… ▽ More To have a better understanding of a physical process, a comparison of ex-perimental data with theoretical values is mandatory. The comparison ismeaningful if the uncertainty in the experiment is accounted well. However,it is seldom seen, especially for a complex phenomenon. We take a test casethrough L subshell ionization of atoms by particle impact. Experimentally,x-ray production cross-sections are measured, but ionization cross-sectionsare calculated theoretically. Furthermore, the uncertainty of the x-ray pro-duction cross-section is mainly statistics and detector-efficiency driven. Butionization cross-section involves many other factors because of the relation-ship between the production and ionization cross section, having wide uncer-tainty spectrum. Consequently, determining the measurement uncertainty inL subshell ionization cross-section is always difficult. We have studied thisissue in the simplest way, where the rule of weighted propagation of rela-tive uncertainty is utilised. We notice that larger uncertainties are involvedin atomic parameters relevant to L1(2s1/2) subshell than those associatedwith the other two L2(2p1/2) and L3(2p3/2) subshells. Hence, comparisonbetween theory and experiment would give higher emphasis onL2andL3subshell ionization cross sections. We believe this work aware us that theappropriate uncertainty evaluation is extremely important for providing theright judgment on the data. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 31 pages, 2 figures

Paths to caustic formation in turbulent aerosols

Authors: Jan Meibohm, Vikash Pandey, Akshay Bhatnagar, Kristian Gustavsson, Dhrubaditya Mitra, Prasad Perlekar, B. Mehlig

Abstract: The dynamics of small, yet heavy, identical particles in turbulence exhibits singularities, called caustics, that lead to large fluctuations in the spatial particle-number density, and in collision velocities. For large particle, inertia the fluid velocity at the particle position is essentially a white-noise signal and caustic formation is analogous to Kramers escape. Here we show that caustic fo… ▽ More The dynamics of small, yet heavy, identical particles in turbulence exhibits singularities, called caustics, that lead to large fluctuations in the spatial particle-number density, and in collision velocities. For large particle, inertia the fluid velocity at the particle position is essentially a white-noise signal and caustic formation is analogous to Kramers escape. Here we show that caustic formation at small particle inertia is different. Caustics tend to form in the vicinity of particle trajectories that experience a specific history of fluid-velocity gradients, characterised by low vorticity and a violent strain exceeding a large threshold. We develop a theory that explains our findings in terms of an optimal path to caustic formation that is approached in the small inertia limit. △ Less

Submitted 15 June, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 6 pages, 3 figures, supplementary material

Journal ref: Phys. Rev. Fluids 6 (2021) L062302

Planetesimals on eccentric orbits erode rapidly

Authors: Lukas Cedenblad, Noemi Schaffer, Anders Johansen, B. Mehlig, Dhrubaditya Mitra

Abstract: We investigate the possibility of erosion of planetesimals in a protoplanetary disk. We use theory and direct numerical simulations (Lattice Boltzmann Method) to calculate the erosion of large -- much larger than the mean-free-path of gas molecules -- bodies of different shapes in flows. We find that erosion follows a universal power-law in time, at intermediate times, independent of the Reynolds… ▽ More We investigate the possibility of erosion of planetesimals in a protoplanetary disk. We use theory and direct numerical simulations (Lattice Boltzmann Method) to calculate the erosion of large -- much larger than the mean-free-path of gas molecules -- bodies of different shapes in flows. We find that erosion follows a universal power-law in time, at intermediate times, independent of the Reynolds number of the flow and the initial shape of the body. Consequently, we estimate that planetesimals in eccentric orbits, of even very small eccentricity, rapidly (in about hundred years) erodes away if the semi-major axis of their orbit lies in the inner disk -- less than about $10$ au. Even planetesimals in circular orbits erode away in approximately ten thousand years if the semi-major axis of their orbits are $\lessapprox 0.6$au. △ Less

Submitted 30 September, 2021; v1 submitted 29 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astrophysical J

Report number: NORDITA 2021-084

Journal ref: The Astrophysical Journal 921 (2021) 123

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

Authors: A. Brandenburg, A. Johansen, P. A. Bourdin, W. Dobler, W. Lyra, M. Rheinhardt, S. Bingert, N. E. L. Haugen, A. Mee, F. Gent, N. Babkovskaia, C. -C. Yang, T. Heinemann, B. Dintrans, D. Mitra, S. Candelaresi, J. Warnecke, P. J. Käpylä, A. Schreiber, P. Chatterjee, M. J. Käpylä, X. -Y. Li, J. Krüger, J. R. Aarnes, G. R. Sarson , et al. (12 additional authors not shown)

Abstract: The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the… ▽ More The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 7 pages, submitted to the Journal for Open Source Software (JOSS)

Report number: NORDITA-2020-087

Journal ref: Journal of Open Source Software 6, 2807 (2021)

Establishing a nearly closed cycling transition in a polyatomic molecule

Authors: Louis Baum, Nathaniel B. Vilas, Christian Hallas, Benjamin L. Augenbraun, Shivam Rava, Debayan Mitra, John M. Doyle

Abstract: We study optical cycling in the polar free radical calcium monohydroxide (CaOH) and establish an experimental path towards scattering $\sim$$10^4$ photons. We report rovibronic branching ratio measurements with precision at the $\sim10^{-4}$ level and observe weak symmetry-forbidden decays to bending modes with non-zero vibrational angular momentum. Calculations are in excellent agreement with the… ▽ More We study optical cycling in the polar free radical calcium monohydroxide (CaOH) and establish an experimental path towards scattering $\sim$$10^4$ photons. We report rovibronic branching ratio measurements with precision at the $\sim10^{-4}$ level and observe weak symmetry-forbidden decays to bending modes with non-zero vibrational angular momentum. Calculations are in excellent agreement with these measurements and predict additional decay pathways. Additionally, we perform high-resolution spectroscopy of the $\widetilde{\text{X}}\,^2Σ^+(12^00)$ and $\widetilde{\text{X}}\,^2Σ^+(12^20)$ hybrid vibrational states of CaOH. These advances establish a path towards radiative slowing, 3D magneto-optical trapping, and sub-Doppler cooling of CaOH. △ Less

Submitted 15 June, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. A 103, 043111 (2021)

Direct Laser Cooling of a Symmetric Top Molecule

Authors: Debayan Mitra, Nathaniel B. Vilas, Christian Hallas, Loïc Anderegg, Benjamin L. Augenbraun, Louis Baum, Calder Miller, Shivam Raval, John M. Doyle

Abstract: We report direct laser cooling of a symmetric top molecule, reducing the transverse temperature of a beam of calcium monomethoxide (CaOCH$_3$) to $1.8\pm0.7$ mK while addressing two distinct nuclear spin isomers. These results open a path to efficient production of ultracold chiral molecules and conclusively demonstrate that by using proper rovibronic optical transitions, both photon cycling and l… ▽ More We report direct laser cooling of a symmetric top molecule, reducing the transverse temperature of a beam of calcium monomethoxide (CaOCH$_3$) to $1.8\pm0.7$ mK while addressing two distinct nuclear spin isomers. These results open a path to efficient production of ultracold chiral molecules and conclusively demonstrate that by using proper rovibronic optical transitions, both photon cycling and laser cooling of complex molecules can be as efficient as for much simpler linear species. △ Less

Submitted 9 November, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

1D Magneto-Optical Trap of Polyatomic Molecules

Authors: Louis Baum, Nathaniel B. Vilas, Christian Hallas, Benjamin L. Augenbraun, Shivam Raval, Debayan Mitra, John M. Doyle

Abstract: We demonstrate a 1D magneto-optical trap of the polar free radical calcium monohydroxide (CaOH). A quasi-closed cycling transition is established to scatter $\sim 10^3$ photons per molecule, predominantly limited by interaction time. This enables radiative laser cooling of CaOH while compressing the molecular beam, leading to a significant increase in on-axis beam brightness and reduction in tempe… ▽ More We demonstrate a 1D magneto-optical trap of the polar free radical calcium monohydroxide (CaOH). A quasi-closed cycling transition is established to scatter $\sim 10^3$ photons per molecule, predominantly limited by interaction time. This enables radiative laser cooling of CaOH while compressing the molecular beam, leading to a significant increase in on-axis beam brightness and reduction in temperature from 8.4 mK to 1.4 mK. △ Less

Submitted 29 January, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. Lett. 124, 133201 (2020)

The spreading of viruses by airborne aerosols: lessons from a first-passage-time problem for tracers in turbulent flows

Authors: Akhilesh Kumar Verma, Akshay Bhatnagar, Dhrubaditya Mitra, Rahul Pandit

Abstract: We study the spreading of viruses, such as SARS-CoV-2, by airborne aerosols, via a new first-passage-time problem for Lagrangian tracers that are advected by a turbulent flow: By direct numerical simulations of the three-dimensional (3D) incompressible, Navier-Stokes equation, we obtain the time $t_R$ at which a tracer, initially at the origin of a sphere of radius $R$, crosses the surface of the… ▽ More We study the spreading of viruses, such as SARS-CoV-2, by airborne aerosols, via a new first-passage-time problem for Lagrangian tracers that are advected by a turbulent flow: By direct numerical simulations of the three-dimensional (3D) incompressible, Navier-Stokes equation, we obtain the time $t_R$ at which a tracer, initially at the origin of a sphere of radius $R$, crosses the surface of the sphere \textit{for the first time}. We obtain the probability distribution function $\mathcal{P}(R,t_R)$ and show that it displays two qualitatively different behaviors: (a) for $R \ll L_{\rm I}$, $\mathcal{P}(R,t_R)$ has a power-law tail $\sim t_R^{-α}$, with the exponent $α= 4$ and $L_{\rm I}$ the integral scale of the turbulent flow; (b) for $l_{\rm I} \lesssim R $, the tail of $\mathcal{P}(R,t_R)$ decays exponentially. We develop models that allow us to obtain these asymptotic behaviors analytically. We show how to use $\mathcal{P}(R,t_R)$ to develop social-distancing guidelines for the mitigation of the spreading of airborne aerosols with viruses such as SARS-CoV-2. △ Less

Submitted 18 June, 2020; v1 submitted 5 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. Research 2, 033239 (2020)

Study of bremsstrahlung radiation coming from a transmission type x-ray generator and its application on EDXRF technique

Authors: Soumya Chatterjee, Sumana Ghosh, D. Mitra

Abstract: We have demonstrated the use of the bremsstrahlung radiation in energy dispersive x-ray fluorescence technique as a tool to perform non-destructive elemental analysis of solid samples employed in inter-disciplinary science research. The bremsstrahlung radiation can be taken from a small, portable, transmission type x-ray generator. As our knowledge goes, till date, this is the first attempt to gen… ▽ More We have demonstrated the use of the bremsstrahlung radiation in energy dispersive x-ray fluorescence technique as a tool to perform non-destructive elemental analysis of solid samples employed in inter-disciplinary science research. The bremsstrahlung radiation can be taken from a small, portable, transmission type x-ray generator. As our knowledge goes, till date, this is the first attempt to generate bremsstrahlung spectrum theoretically from a portable transmission type x-ray generator and compared it with the actual observation. Theoretically generated bremsstrahlung spectra are also found to be in good agreement with the experimentally observed spectra obtained with various operating anode voltages of the x-ray generator. A computer program has been developed utilizing few atomic parameters to obtain the elemental concentrations in the sample by a single run using the whole bremsstrahlung. The knowledge of incoming x-ray flux, geometry of experimental arrangements is not required in this technique. To validate the technique, we have taken two samples: an Indian one rupee coin of the year 2000 and a NIST made brass sample, whose compositions are well known. We used these samples to expose to the bremsstrahlung radiation, produced by the operating anode voltages 20, 25, and 30 KV. The relative concentrations of different elements are determined, which is in good agreement with the earlier results. △ Less

Submitted 3 March, 2023; v1 submitted 23 November, 2019; originally announced November 2019.

Comments: 6 pages, 3 figures

Distinguishing between paediatric brain tumour types using multi-parametric magnetic resonance imaging and machine learning: a multi-site study

Authors: James T. Grist, Stephanie Withey, Lesley MacPherson, Adam Oates, Stephen Powell, Jan Novak, Laurence Abernethy, Barry Pizer, Richard Grundy, Simon Bailey, Dipayan Mitra, Theodoros N. Arvanitis, Dorothee P. Auer, Shivaram Avula, Andrew C Peet

Abstract: The imaging and subsequent accurate diagnosis of paediatric brain tumours presents a radiological challenge, with magnetic resonance imaging playing a key role in providing tumour specific imaging information. Diffusion weighted and perfusion imaging are commonly used to aid the non invasive diagnosis of paediatric brain tumours, but are usually evaluated by expert qualitative review. Quantitative… ▽ More The imaging and subsequent accurate diagnosis of paediatric brain tumours presents a radiological challenge, with magnetic resonance imaging playing a key role in providing tumour specific imaging information. Diffusion weighted and perfusion imaging are commonly used to aid the non invasive diagnosis of paediatric brain tumours, but are usually evaluated by expert qualitative review. Quantitative studies are mainly single centre and single modality. The aim of this work was to combine multi centre diffusion and perfusion imaging, with machine learning, to develop machine learning based classifiers to discriminate between three common paediatric tumour types. The results show that diffusion and perfusion weighted imaging of both the tumour and whole brain provide significant features which differ between tumour types, and that combining these features gives the optimal machine learning classifier with greater than 80 percent predictive precision. This work represents a step forward to aid in the non invasive diagnosis of paediatric brain tumours, using advanced clinical imaging. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Comments: 22 pages, 7 figures and tables

Clustering and energy spectra in two-dimensional dusty gas turbulence

Authors: Vikash Pandey, Dhrubaditya Mitra, Prasad Perlekar

Abstract: We present Direct Numerical Simulation (DNS) of heavy inertial particles (dust) immersed in two-dimensional turbulent flow (gas). The dust are modeled as mono-dispersed heavy particles capable of modifying the flow through two-way coupling. By varying the Stokes number (St) and the mass-loading parameter $(φ_{\rm m})$, we study the clustering phenomenon and the gas phase kinetic energy spectra. We… ▽ More We present Direct Numerical Simulation (DNS) of heavy inertial particles (dust) immersed in two-dimensional turbulent flow (gas). The dust are modeled as mono-dispersed heavy particles capable of modifying the flow through two-way coupling. By varying the Stokes number (St) and the mass-loading parameter $(φ_{\rm m})$, we study the clustering phenomenon and the gas phase kinetic energy spectra. We find that the dust-dust correlation dimension $(d_2)$ also depends on $φ_{\rm m}$. In particular, clustering decreases as mass-loading $(φ_{\rm m})$, is increased. In the kinetic energy spectra of gas we show: (i) emergence of a new scaling regime, (ii) the scaling exponent in this regime is not unique but rather a function of both St and $(φ_{\rm m} )$. Using a scale-by-scale enstrophy budget analysis we show in the new scaling regime, viscous dissipation due to the gas balances back-reaction from the dust. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. E 100, 013114 (2019)

The breakdown of Darcy's law in a soft porous material

Authors: Marco E. Rosti, Satyajit Pramanik, Luca Brandt, Dhrubaditya Mitra

Abstract: We perform direct numerical simulations of the flow through a model of a deformable porous medium. Our model is a two-dimensional hexagonal lattice, with defects, of soft elastic cylindrical pillars, with elastic shear modulus $G$, immersed in a liquid. We use a two-phase approach: the liquid phase is a viscous fluid and the solid phase is modeled as an incompressible viscoelastic material, whose… ▽ More We perform direct numerical simulations of the flow through a model of a deformable porous medium. Our model is a two-dimensional hexagonal lattice, with defects, of soft elastic cylindrical pillars, with elastic shear modulus $G$, immersed in a liquid. We use a two-phase approach: the liquid phase is a viscous fluid and the solid phase is modeled as an incompressible viscoelastic material, whose complete nonlinear structural response is considered. We observe that the Darcy flux ($q$) is a nonlinear function -- steeper than linear -- of the pressure-difference ($ΔP$) across the medium. Furthermore, the flux is larger for a softer medium (smaller $G$). We construct a theory of this super-linear behavior by modelling the channels between the solid cylinders as elastic channels whose walls are made of a material with a linear constitutive relation but can undergo large deformation. Our theory further predicts that the flow permeability is a universal function of $ΔP/G$, which is confirmed by the present simulations. △ Less

Submitted 9 December, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: 6 pages, 3 figures, Some minor changes (including the title) from the previous submission

Journal ref: Soft Matter, 2020, 16, 939-944

Relative velocities in bi-disperse turbulent aerosols: simulations and theory

Authors: Akshay Bhatnagar, K. Gustavsson, B. Mehlig, Dhrubaditya Mitra

Abstract: We perform direct numerical simulations of a bi-disperse suspension of heavy spherical particles in forced, homogeneous, and isotropic three-dimensional turbulence. We compute the joint distribution of relative particle distances and longitudinal relative velocities between particles of different sizes, and compare the results with recent theoretical predictions [Meibohm et al. Phys. Rev. E 96 (20… ▽ More We perform direct numerical simulations of a bi-disperse suspension of heavy spherical particles in forced, homogeneous, and isotropic three-dimensional turbulence. We compute the joint distribution of relative particle distances and longitudinal relative velocities between particles of different sizes, and compare the results with recent theoretical predictions [Meibohm et al. Phys. Rev. E 96 (2017) 061102] for the shape of this distribution. We also compute the moments of relative velocities as a function of particle separation, and compare with the theoretical predictions. We observe good agreement. △ Less

Submitted 20 December, 2018; v1 submitted 27 September, 2018; originally announced September 2018.

Comments: 12 pages, 3 figures

Journal ref: Phys. Rev. E 98, 063107 (2018)

Kazantsev dynamo in turbulent compressible flows

Authors: Marco Martins Afonso, Dhrubaditya Mitra, Dario Vincenzi

Abstract: We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling e… ▽ More We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling exponents for the solenoidal and potential parts differ, in particular if they correspond to typical Kolmogorov and Burgers values, we again find that an increase in compressibility slows down the growth rate but does not turn it off. The slow down is, however, weaker and the critical magnetic Reynolds number is lower than when both the solenoidal and potential components display the Kolmogorov scaling. Intriguingly, we find that there exist cases, when the potential part is smoother than the solenoidal part, for which an increase in compressibility increases the growth rate. We also find that the critical value of the scaling exponent above which a dynamo is seen is unity irrespective of the compressibility. Finally, we realize that the dimension $d = 3$ is special, since for all other values of $d$ the critical exponent is higher and depends on the compressibility. △ Less

Submitted 1 March, 2019; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: 12 pages, 6 figures

Journal ref: Proc. R. Soc. A 475, 20180591 (2019)

Computation of a Theoretical Membrane Phase Diagram, and the Role of Phase in Lipid Raft-Mediated Protein Organization

Authors: Eshan D. Mitra, Samuel C. Whitehead, David Holowka, Barbara Baird, James P. Sethna

Abstract: Lipid phase heterogeneity in the plasma membrane is thought to be crucial for many aspects of cell signaling, but the physical basis of participating membrane domains such as "lipid rafts" remains controversial. Here we consider a lattice model yielding a phase diagram that includes several states proposed to be relevant for the cell membrane, including microemulsion - which can be related to memb… ▽ More Lipid phase heterogeneity in the plasma membrane is thought to be crucial for many aspects of cell signaling, but the physical basis of participating membrane domains such as "lipid rafts" remains controversial. Here we consider a lattice model yielding a phase diagram that includes several states proposed to be relevant for the cell membrane, including microemulsion - which can be related to membrane curvature - and Ising critical behavior. Using a neural network-based machine learning approach, we compute the full phase diagram of this lattice model. We analyze selected regions of this phase diagram in the context of a signaling initiation event in mast cells: recruitment of the membrane-anchored tyrosine kinase Lyn to a cluster of transmembrane of IgE-FcεRI receptors. We find that model membrane systems in microemulsion and Ising critical states can mediate roughly equal levels of kinase recruitment (binding energy ~ -0.6 kBT), whereas a membrane near a tricritical point can mediate much stronger kinase recruitment (-1.7 kBT). By comparing several models for lipid heterogeneity within a single theoretical framework, this work points to testable differences between existing models. We also suggest the tricritical point as a new possibility for the basis of membrane domains that facilitate preferential partitioning of signaling components. △ Less

Submitted 3 February, 2018; originally announced February 2018.

Comments: 33 pages, 7 figures, 16 supplementary pages, 10 supplementary figures

Heavy inertial particles in turbulent flows gain energy slowly but lose it rapidly

Authors: Akshay Bhatnagar, Anupam Gupta, Dhrubaditya Mitra, Rahul Pandit

Abstract: We present an extensive numerical study of the time irreversibility of the dynamics of heavy inertial particles in three-dimensional, statistically homogeneous and isotropic turbulent flows. We show that the probability density function (PDF) of the increment, $W(τ)$, of a particle's energy over a time-scale $τ$ is non-Gaussian, and skewed towards negative values. This implies that, on average, pa… ▽ More We present an extensive numerical study of the time irreversibility of the dynamics of heavy inertial particles in three-dimensional, statistically homogeneous and isotropic turbulent flows. We show that the probability density function (PDF) of the increment, $W(τ)$, of a particle's energy over a time-scale $τ$ is non-Gaussian, and skewed towards negative values. This implies that, on average, particles gain energy over a period of time that is longer than the duration over which they lose energy. We call this $\textit{slow gain}$ and $\textit{fast loss}$. We find that the third moment of $W(τ)$ scales as $τ^3$, for small values of $τ$. We show that the PDF of power-input $p$ is negatively skewed too; we use this skewness ${\rm Ir}$ as a measure of the time-irreversibility and we demonstrate that it increases sharply with the Stokes number ${\rm St}$, for small ${\rm St}$; this increase slows down at ${\rm St} \simeq 1$. Furthermore, we obtain the PDFs of $t^+$ and $t^-$, the times over which $p$ has, respectively, positive or negative signs, i.e., the particle gains or loses energy. We obtain from these PDFs a direct and natural quantification of the the slow-gain and fast-loss of the particles, because these PDFs possess exponential tails, whence we infer the characteristic loss and gain times $t_{\rm loss}$ and $t_{\rm gain}$, respectively; and we obtain $t_{\rm loss} < t_{\rm gain}$, for all the cases we have considered. Finally, we show that the slow-gain in energy of the particles is equally likely in vortical or strain-dominated regions of the flow; in contrast, the fast-loss of energy occurs with greater probability in the latter than in the former. △ Less

Submitted 19 November, 2017; originally announced November 2017.

Report number: NORDITA 2017-117

Journal ref: Phys. Rev. E 97, 033102 (2018)

Topology of two-dimensional turbulent flows of dust and gas

Authors: Dhrubaditya Mitra, Prasad Perlekar

Abstract: We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, ${\rm St} < 1$, using both Lagrangian and Eulerian approach (with a shock-capturing scheme). We find that: The dust-density field in our Eulerian simulations have the same correlation dimension $d_2$ as obtained fr… ▽ More We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, ${\rm St} < 1$, using both Lagrangian and Eulerian approach (with a shock-capturing scheme). We find that: The dust-density field in our Eulerian simulations have the same correlation dimension $d_2$ as obtained from the clustering of particles in the Lagrangian simulations for ${\rm St} < 1$; The cumulative probability distribution function of the dust-density coarse-grained over a scale $r$ in the inertial range has a left-tail with a power-law fall-off indicating presence of voids; The energy spectrum of the dust-velocity has a power-law range with an exponent that is same as the gas-velocity spectrum except at very high Fourier modes; The compressibility of the dust-velocity field is proportional to ${\rm St}^2$. We quantify the topological properties of the dust-velocity and the gas-velocity through their gradient matrices, called $\mathcal{A}$ and $\mathcal{B}$, respectively. The topological properties of $\mathcal{B}$ are the same in Eulerian and Lagrangian frames only if the Eulerian data are weighed by the dust-density -- a correspondence that we use to study Lagrangian properties of $\mathcal{A}$. In the Lagrangian frame, the mean value of the trace of $\mathcal{A} \sim - \exp(-C/{\rm St}$, with a constant $C\approx 0.1$. The topology of the dust-velocity fields shows that as ${\rm St} increases the contribution to negative divergence comes mostly from saddles and the contribution to positive divergence comes from both vortices and saddles. Compared to the Eulerian case, the density-weighed Eulerian case has less inward spirals and more converging saddles. Outward spirals are the least probable topological structures in both cases. △ Less

Submitted 9 November, 2017; originally announced November 2017.

Comments: abstract condensed to fit arxiv format

Report number: NORDITA 2017-115

Journal ref: Phys. Rev. Fluids 3, 044303 (2018)

Probing quench dynamics across a quantum phase transition into a 2D Ising antiferromagnet

Authors: Elmer Guardado-Sanchez, Peter T. Brown, Debayan Mitra, Trithep Devakul, David A. Huse, Peter Schauss, Waseem S. Bakr

Abstract: Simulating the real-time evolution of quantum spin systems far out of equilibrium poses a major theoretical challenge, especially in more than one dimension. We experimentally explore the dynamics of a two-dimensional Ising spin system with transverse and longitudinal fields as we quench it across a quantum phase transition from a paramagnet to an antiferromagnet. We realize the system with a near… ▽ More Simulating the real-time evolution of quantum spin systems far out of equilibrium poses a major theoretical challenge, especially in more than one dimension. We experimentally explore the dynamics of a two-dimensional Ising spin system with transverse and longitudinal fields as we quench it across a quantum phase transition from a paramagnet to an antiferromagnet. We realize the system with a near unit-occupancy atomic array of over 200 atoms obtained by loading a spin-polarized band insulator of fermionic lithium into an optical lattice and induce short-range interactions by direct excitation to a low-lying Rydberg state. Using site-resolved microscopy, we probe the correlations in the system after a sudden quench from the paramagnetic state and compare our measurements to exact calculations in the regime where it is possible. We achieve many-body states with longer-range antiferromagnetic correlations by implementing a near-adiabatic quench and study the buildup of correlations as we cross the quantum phase transition at different rates. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Journal ref: Phys. Rev. X 8, 021069 (2018)

Statistics of the relative velocity of particles in turbulent flows : monodisperse particles

Authors: Akshay Bhatnagar, K. Gustavsson, Dhrubaditya Mitra

Abstract: We use direct numerical simulations to calculate the joint probability density function of the relative distance $R$ and relative radial velocity component $V_R$ for a pair of heavy inertial particles suspended in homogeneous and isotropic turbulent flows. At small scales the distribution is scale invariant, with a scaling exponent that is related to the particle-particle correlation dimension in… ▽ More We use direct numerical simulations to calculate the joint probability density function of the relative distance $R$ and relative radial velocity component $V_R$ for a pair of heavy inertial particles suspended in homogeneous and isotropic turbulent flows. At small scales the distribution is scale invariant, with a scaling exponent that is related to the particle-particle correlation dimension in phase space, $D_2$. It was argued [1, 2] that the scale invariant part of the distribution has two asymptotic regimes: (1) $|V_R| \ll R$ where the distribution depends solely on $R$; and (2) $|V_R| \gg R$ where the distribution is a function of $|V_R|$ alone. The probability distributions in these two regimes are matched along a straight line $|V_R| = z^\ast R$. Our simulations confirm that this is indeed correct. We further obtain $D_2$ and $z^\ast$ as a function of the Stokes number, ${\rm St}$. The former depends non-monotonically on ${\rm St}$ with a minimum at about ${\rm St} \approx 0.7$ and the latter has only a weak dependence on ${\rm St}$. △ Less

Submitted 13 October, 2017; originally announced October 2017.

Report number: NORDITA 2017-108

Journal ref: Phys. Rev. E 97, 023105 (2018)

Rheology of suspensions of viscoelastic spheres: deformability as an effective volume fraction

Authors: Marco Edoardo Rosti, Luca Brandt, Dhrubaditya Mitra

Abstract: We study suspensions of deformable (viscoelastic) spheres in a Newtonian solvent in plane Couette geometry, by means of direct numerical simulations. We find that in the limit of vanishing inertia the effective viscosity $μ$ of the suspension increases as the volume-fraction occupied by the spheres $Φ$ increases and decreases as the elastic modulus of the spheres $G$ decreases; the function… ▽ More We study suspensions of deformable (viscoelastic) spheres in a Newtonian solvent in plane Couette geometry, by means of direct numerical simulations. We find that in the limit of vanishing inertia the effective viscosity $μ$ of the suspension increases as the volume-fraction occupied by the spheres $Φ$ increases and decreases as the elastic modulus of the spheres $G$ decreases; the function $μ(Φ,G)$ collapses to an universal function, $μ(Φ_e)$, with a reduced effective volume fraction $Φ_e(Φ,G)$. Remarkably, the function $μ(Φ_e)$ is the well-known Eilers fit that describes the rheology of suspension of rigid spheres at all $Φ$. Our results suggest new ways to interpret macro-rheology of blood. △ Less

Submitted 10 January, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Journal ref: Phys. Rev. Fluids 3, 012301 (2018)