Showing 1–50 of 496 results for author: Price, J

Observational Signatures of Circumbinary Discs II: Kinematic Signatures in Velocity Residuals

Authors: Josh Calcino, Brodie Norfolk, Daniel J. Price, Thomas Hilder, Jessica Speedie, Christophe Pinte, Himanshi Garg, Richard Teague, Cassandra Hall, Jochen Stadler

Abstract: Kinematic studies of protoplanetary discs are a valuable method for uncovering hidden companions. In the first paper of this series, we presented five morphological and kinematic criteria that aid in asserting the binary nature of a protoplanetary disc. In this work we study the kinematic signatures of circumbinary discs in the residuals of their velocity maps. We show that Doppler-flips, spiral a… ▽ More Kinematic studies of protoplanetary discs are a valuable method for uncovering hidden companions. In the first paper of this series, we presented five morphological and kinematic criteria that aid in asserting the binary nature of a protoplanetary disc. In this work we study the kinematic signatures of circumbinary discs in the residuals of their velocity maps. We show that Doppler-flips, spiral arms, eccentric gas motion, fast flows inside of the cavity, and vortex-like kinematic signatures are commonly observed. Unlike in the planetary mass companion case, Doppler-flips in circumbinary discs are not necessarily centred on a companion, and can extend towards the cavity edge. We then compare the kinematic signatures in our simulations with observations and see similarities to the Doppler-flip signal in HD 100546 and the vortex-like kinematic signatures in HD 142527. Our analysis also reveals kinematic evidence for binarity in several protoplantary disks typically regarded as circumstellar rather than circumbinary, including AB Aurigae and HD 100546. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 15 pages, 12 Figures. Submitted to MNRAS, revised version. Comments welcomed and appreciated!

Chemical tracers of a highly eccentric AGB-main sequence star binary

Authors: T. Danilovich, J. Malfait, M. Van de Sande, M. Montargès, P. Kervella, F. De Ceuster, A. Coenegrachts, T. J. Millar, A. M. S. Richards, L. Decin, C. A. Gottlieb, C. Pinte, E. De Beck, D. J. Price, K. T. Wong, J. Bolte, K. M. Menten, A. Baudry, A. de Koter, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, E. Lagadec , et al. (10 additional authors not shown)

Abstract: Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chem… ▽ More Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chemical and dynamical analysis to reveal a highly eccentric and long-period orbit for W Aquilae, a binary system containing an AGB star and a main sequence companion. Our results are based on anisotropic SiN emission, the first detections of NS and SiC towards an S-type star, and density structures observed in the CO emission. These features are all interpreted as having formed during periastron interactions. Our astrochemistry-based method can yield stringent constraints on the orbital parameters of long-period binaries containing AGB stars, and will be applicable to other systems. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Comments: This version of the article has been accepted for publication, after peer review (when applicable) and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at: https://doi.org/10.1038/s41550-023-02154-y

Journal ref: Nature Astronomy, 8, pages 308-327 (2024)

On shock capturing in smoothed particle hydrodynamics

Authors: Daniel J. Price

Abstract: For the past 20 years, our approach to shock capturing in smoothed particle hydrodynamics (SPH) has been to use artificial viscosity and conductivity terms supplemented by switches to control excess dissipation away from shocks (Monaghan 1997; Morris & Monaghan 1997). This approach has been demonstrated to be superior to approximate Riemann solvers in a recent comparison (Puri & Ramachandran 2014)… ▽ More For the past 20 years, our approach to shock capturing in smoothed particle hydrodynamics (SPH) has been to use artificial viscosity and conductivity terms supplemented by switches to control excess dissipation away from shocks (Monaghan 1997; Morris & Monaghan 1997). This approach has been demonstrated to be superior to approximate Riemann solvers in a recent comparison (Puri & Ramachandran 2014). The Cullen & Dehnen (2010) switch is regarded as the state of the art. But are we missing something? I will present a novel approach to shock capturing in SPH that utilises the philosophy of approximate Riemann solvers but provides a direct improvement on the ability to reduce excess dissipation away from shocks while preserving the fidelity of the shock itself. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 6 pages, 12 figures, proceedings of 2019 international SPHERIC workshop, Exeter UK. Despite what I said in these proceedings, I later concluded that switches are indeed better, because one can ensure positive definite entropy increase, which is not the case with reconstruction (see Price & Laibe 2020; arXiv:2005.06562)

Dust formation during the interaction of binary stars by common envelope

Authors: Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Miguel González-Bolívar, Mike Y. M. Lau, Chunliang Mu, Ryosuke Hirai, Taïssa Danilovich, Mansi Kasliwal

Abstract: We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as $\sim$1-3 yrs after the onset of the CE, and are smaller than grains formed later.… ▽ More We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as $\sim$1-3 yrs after the onset of the CE, and are smaller than grains formed later. As the simulations progress, a high-opacity dusty shell forms, resulting in the CE photosphere being up to an order of magnitude larger than it would be without the inclusion of dust. At the end of the simulations, the total dust yield is $0.0082~M_{\odot}$ ($0.022~M_{\odot}$) for a CE with a $1.7~M_{\odot}$ ($3.7~M_{\odot}$) AGB star. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 6 pages, 6 figures, contributed talk to the Proceedings IAU Symposium No. 384

Compaction during fragmentation and bouncing produces realistic dust grain porosities in protoplanetary discs

Authors: Stéphane Michoulier, Jean-François Gonzalez, Daniel J. Price

Abstract: Context: In protoplanetary discs, micron-sized dust grows to form millimetre- to centimetre-sized pebbles but encounters several barriers during its evolution. Collisional fragmentation and radial drift impede further dust growth to planetesimal size. Fluffy grains have been hypothesised to solve these problems. While porosity leads to faster grain growth, the implied porosity values obtained from… ▽ More Context: In protoplanetary discs, micron-sized dust grows to form millimetre- to centimetre-sized pebbles but encounters several barriers during its evolution. Collisional fragmentation and radial drift impede further dust growth to planetesimal size. Fluffy grains have been hypothesised to solve these problems. While porosity leads to faster grain growth, the implied porosity values obtained from previous simulations were larger than suggested by observations. Aims: In this paper, we study the influence of porosity on dust evolution taking into account growth, bouncing, fragmentation, compaction, rotational disruption and snow lines, in order to understand their impact on dust evolution. Methods: We develop a module for porosity evolution for the 3D Smoothed Particle Hydrodynamics (SPH) code Phantom that accounts for dust growth and fragmentation. This mono-disperse model is integrated into both a 1D code and the 3D code to capture the overall evolution of dust and gas. Results: We show that porosity helps dust growth and leads to the formation of larger solids than when considering compact grains, as predicted by previous work. Our simulations taking into account compaction during fragmentation show that large millimetre grains are still formed, but are 10 to 100 times more compact. Thus, mm sizes with typical filling factors of ~0.1 match the values measured on comets or via polarimetric observations of protoplanetary discs. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 24 pages, 21 figures, accepted for publication in A&A

Journal ref: A&A 688, A31 (2024)

First Measurement of Deeply Virtual Compton Scattering on the Neutron with Detection of the Active Neutron

Authors: CLAS Collaboration, A. Hobart, S. Niccolai, M. Čuić, K. Kumerički, P. Achenbach, J. S. Alvarado, W. R. Armstrong, H. Atac, H. Avakian, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, S. Boiarinov, M. Bondi, W. A. Booth, F. Bossù, K. -Th. Brinkmann, W. J. Briscoe , et al. (124 additional authors not shown)

Abstract: Measuring Deeply Virtual Compton Scattering on the neutron is one of the necessary steps to understand the structure of the nucleon in terms of Generalized Parton Distributions (GPDs). Neutron targets play a complementary role to transversely polarized proton targets in the determination of the GPD $E$. This poorly known and poorly constrained GPD is essential to obtain the contribution of the qua… ▽ More Measuring Deeply Virtual Compton Scattering on the neutron is one of the necessary steps to understand the structure of the nucleon in terms of Generalized Parton Distributions (GPDs). Neutron targets play a complementary role to transversely polarized proton targets in the determination of the GPD $E$. This poorly known and poorly constrained GPD is essential to obtain the contribution of the quarks' angular momentum to the spin of the nucleon. DVCS on the neutron was measured for the first time selecting the exclusive final state by detecting the neutron, using the Jefferson Lab longitudinally polarized electron beam, with energies up to 10.6 GeV, and the CLAS12 detector. The extracted beam-spin asymmetries, combined with DVCS observables measured on the proton, allow a clean quark-flavor separation of the imaginary parts of the GPDs $H$ and $E$. △ Less

Submitted 25 June, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: 7 pages, 6 figures

Report number: JLAB-PHY-24-4089

V892 Tau: A tidally perturbed circumbinary disc in a triple stellar system

Authors: Antoine Alaguero, Nicolás Cuello, François Ménard, Simone Ceppi, Álvaro Ribas, Rebecca Nealon, Miguel Vioque, Andrés Izquierdo, James Miley, Enrique Macías, Daniel J. Price

Abstract: V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data… ▽ More V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data with V892 Tau as a triple system and perform a grid of hydrodynamical simulations testing several orbits of the companion. The simulation outputs are then post-processed to build synthetic maps that we compare to the observations. The 12CO emission of the disc shows clear non-Keplerian features such as spiral arms. When comparing the data with our synthetic observations, we interpret these features as ongoing interactions with the companion. Our simulations indicate that an eccentricity of 0.5 of the companion is needed to reproduce the observed disc extent and that a mutual inclination of approximately 60° with the inner binary reproduces the measured disc tilt. In order to explain most of the features of the circumbinary disc, we propose that V892 Tau NE follows a misaligned eccentric orbit, with an eccentricity between 0.2 and 0.5 and a mutual inclination between 30° and 60°. Such a misaligned companion suggests the disc is oscillating and precessing with time, stabilising in an intermediate plane with a non-zero mutual inclination with the inner binary. Given that orbital configuration, we show that the stability of future planets is compromised in the second half of the disc once the gas has dissipated. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A, 29 pages, 15 figures

Journal ref: A&A 687, A311 (2024)

Decretion disc size in Be/X-ray binaries depends upon the disc aspect ratio

Authors: Rebecca G. Martin, Stephen H. Lubow, Philip J. Armitage, Daniel J. Price

Abstract: With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation… ▽ More With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation. For low disc aspect ratio, $H/R\lesssim 0.1$, the disc is efficiently tidally truncated by the neutron star. Most material that escapes the Roche lobe of the Be star is accreted by the neutron star through tidal streams. For larger disc aspect ratio, the outflow rate through the Be star disc is higher, tidal truncation becomes inefficient, the disc fills the Roche lobe and extends to the orbit of the companion. Some material escapes the binary as a gas stream that begins near the L2 point. While the accretion rate onto the neutron star is higher, the fraction of the outflow that is accreted by the neutron star is smaller. Low density Be star discs are expected to be approximately isothermal, such that $H/R$ increases with radius. Tidal truncation is therefore weaker for larger separation binaries, and lower mass primaries. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS

Cross-modal Diffusion Modelling for Super-resolved Spatial Transcriptomics

Authors: Xiaofei Wang, Xingxu Huang, Stephen J. Price, Chao Li

Abstract: The recent advancement of spatial transcriptomics (ST) allows to characterize spatial gene expression within tissue for discovery research. However, current ST platforms suffer from low resolution, hindering in-depth understanding of spatial gene expression. Super-resolution approaches promise to enhance ST maps by integrating histology images with gene expressions of profiled tissue spots. Howeve… ▽ More The recent advancement of spatial transcriptomics (ST) allows to characterize spatial gene expression within tissue for discovery research. However, current ST platforms suffer from low resolution, hindering in-depth understanding of spatial gene expression. Super-resolution approaches promise to enhance ST maps by integrating histology images with gene expressions of profiled tissue spots. However, current super-resolution methods are limited by restoration uncertainty and mode collapse. Although diffusion models have shown promise in capturing complex interactions between multi-modal conditions, it remains a challenge to integrate histology images and gene expression for super-resolved ST maps. This paper proposes a cross-modal conditional diffusion model for super-resolving ST maps with the guidance of histology images. Specifically, we design a multi-modal disentangling network with cross-modal adaptive modulation to utilize complementary information from histology images and spatial gene expression. Moreover, we propose a dynamic cross-attention modelling strategy to extract hierarchical cell-to-tissue information from histology images. Lastly, we propose a co-expression-based gene-correlation graph network to model the co-expression relationship of multiple genes. Experiments show that our method outperforms other state-of-the-art methods in ST super-resolution on three public datasets. △ Less

Submitted 27 May, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Eddington envelopes: The fate of stars on parabolic orbits tidally disrupted by supermassive black holes

Authors: Daniel J. Price, David Liptai, Ilya Mandel, Joanna Shepherd, Giuseppe Lodato, Yuri Levin

Abstract: Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the retur… ▽ More Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the returning debris stream. These reveal the formation of an expanding asymmetric bubble of material extending to hundreds of astronomical units -- an outflowing Eddington envelope with an optically thick inner region. Such envelopes have been hypothesised as the reprocessing layer needed to explain optical/UV emission in tidal disruption events, but never produced self-consistently in a simulation. Our model broadly matches the observed light curves with low temperatures, faint luminosities, and line widths of 10,000--20,000 km/s. △ Less

Submitted 9 July, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 12 pages + appendices, 15 figures, accepted to ApJL. Data+movies on zenodo at https://doi.org/10.5281/zenodo.11438154 . Movies on YouTube... watch this one first: https://www.youtube.com/watch?v=KlffYCtYxFU then: https://www.youtube.com/watch?v=bPurbeQNgvI&t=19s and finally watch the expanding "bubble" with rho and T at photosphere: https://www.youtube.com/watch?v=D4dyZeUbfvM

Synthetic Data for Robust Stroke Segmentation

Authors: Liam Chalcroft, Ioannis Pappas, Cathy J. Price, John Ashburner

Abstract: Deep learning-based semantic segmentation in neuroimaging currently requires high-resolution scans and extensive annotated datasets, posing significant barriers to clinical applicability. We present a novel synthetic framework for the task of lesion segmentation, extending the capabilities of the established SynthSeg approach to accommodate large heterogeneous pathologies with lesion-specific augm… ▽ More Deep learning-based semantic segmentation in neuroimaging currently requires high-resolution scans and extensive annotated datasets, posing significant barriers to clinical applicability. We present a novel synthetic framework for the task of lesion segmentation, extending the capabilities of the established SynthSeg approach to accommodate large heterogeneous pathologies with lesion-specific augmentation strategies. Our method trains deep learning models, demonstrated here with the UNet architecture, using label maps derived from healthy and stroke datasets, facilitating the segmentation of both healthy tissue and pathological lesions without sequence-specific training data. Evaluated against in-domain and out-of-domain (OOD) datasets, our framework demonstrates robust performance, rivaling current methods within the training domain and significantly outperforming them on OOD data. This contribution holds promise for advancing medical imaging analysis in clinical settings, especially for stroke pathology, by enabling reliable segmentation across varied imaging sequences with reduced dependency on large annotated corpora. Code and weights available at https://github.com/liamchalcroft/SynthStroke. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Partial tidal disruption events: The elixir of life

Authors: Megha Sharma, Daniel J. Price, Alexander Heger

Abstract: In our Galactic Center, about 10,000 to 100,000 stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the… ▽ More In our Galactic Center, about 10,000 to 100,000 stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the tidal disruption of 1, 3, and 10 solar mass stars at zero-age (ZAMS), middle-age (MAMS), and terminal-age main-sequence (TAMS). We map the disruption remnants into Kepler in order to understand their post-distribution evolution. We find distinct characteristics in the remnants, including increased radius, rapid core rotation, and differential rotation in the envelope. The remnants undergo composition mixing that affects their stellar evolution. Whereas the remnants formed by disruption of ZAMS models evolve similarly to unperturbed models of the same mass, for MAMS and TAMS stars, the remnants have higher luminosity and effective temperature. Potential observational signatures include peculiarities in nitrogen and carbon abundances, higher luminosity, rapid rotation, faster evolution, and unique tracks in the Hertzsprung-Russell diagram △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 22 pages, 27 figures, submitted to MNRAS, comments welcome

Detailed Report on the Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm

Authors: D. P. Aguillard, T. Albahri, D. Allspach, A. Anisenkov, K. Badgley, S. Baeßler, I. Bailey, L. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, E. Barzi, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, S. Braun, M. Bressler, G. Cantatore, R. M. Carey, B. C. K. Casey , et al. (168 additional authors not shown)

Abstract: We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference b… ▽ More We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference between the muon spin precession frequency and its cyclotron frequency. This difference is normalized to the strength of the magnetic field, measured using Nuclear Magnetic Resonance (NMR). The ratio is then corrected for small contributions from beam motion, beam dispersion, and transient magnetic fields. We measure $a_μ= 116 592 057 (25) \times 10^{-11}$ (0.21 ppm). This is the world's most precise measurement of this quantity and represents a factor of $2.2$ improvement over our previous result based on the 2018 dataset. In combination, the two datasets yield $a_μ(\text{FNAL}) = 116 592 055 (24) \times 10^{-11}$ (0.20 ppm). Combining this with the measurements from Brookhaven National Laboratory for both positive and negative muons, the new world average is $a_μ$(exp) $ = 116 592 059 (22) \times 10^{-11}$ (0.19 ppm). △ Less

Submitted 22 May, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 48 pages, 29 figures; 4 pages of Supplement Material; version accepted for publication in Physical Review D

Report number: FERMILAB-PUB-24-0084-AD-CSAID-PPD

Adaptive hybrid density functionals

Authors: Danish Khan, Alastair James Arthur Price, Maximilian L. Ach, O. Anatole von Lilienfeld

Abstract: Exact exchange contributions are known to crucially affect electronic states, which in turn govern covalent bond formation and breaking in chemical species. Empirically averaging the exact exchange admixture over compositional degrees of freedom, hybrid density functional approximations have been widely successful, yet have fallen short to reach high level quantum chemistry accuracy, primarily due… ▽ More Exact exchange contributions are known to crucially affect electronic states, which in turn govern covalent bond formation and breaking in chemical species. Empirically averaging the exact exchange admixture over compositional degrees of freedom, hybrid density functional approximations have been widely successful, yet have fallen short to reach high level quantum chemistry accuracy, primarily due to delocalization errors. We propose to `adaptify` hybrid functionals by generating optimal admixture ratios of exact exchange on the fly, i.e. specifically for any chemical compound, using extremely data efficient quantum machine learning models that carry negligible overhead. The adaptive Perdew-Burke-Ernzerhof based hybrid density functional (aPBE0) is shown to yield atomization energies with sufficient accuracy to effectively cure the infamous spin gap problem in open shell systems, such as carbenes. aPBE0 further improves energetics, electron densities, and HOMO-LUMO gaps in organic molecules drawn from the QM9 and QM7b data set. Obtained with aPBE0 in a large basis, we present a revision of the entire QM9 data set (revQM9) with an estimated quality vastly superior to the original containing on average, stronger covalent binding, larger band-gaps, more localized electron densities, and larger dipole-moments. While aPBE0 is applicable in the equilibrium regime, outstanding limitations include covalent bond dissociation when going beyond the Coulson-Fisher point. △ Less

Submitted 10 April, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

The orbit of HD 142527 B is too compact to explain many of the disc features

Authors: M. Nowak, S. Rowther, S. Lacour, F. Meru, R. Nealon, D. J. Price

Abstract: HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13… ▽ More HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13 $\pm$ 0.03 $M_\odot{}$), which is thought to be responsible for most of the structures observed in the surrounding disc. We gathered VLTI/GRAVITY observations of HD 142527, either from our own programmes or from the ESO archive. We used this inhomogeneous set of data to extract a total of seven high-precision measurements of the relative astrometry between HD 142527 A and B, spread from mid-2017 to early 2021. Combined with what is available in the literature, we now have 9 yr of astrometric monitoring on HD 142527. We used orbit fitting tools to determine the orbital parameters of HD 142527 B, and used them as inputs for a 3D hydrodynamical model of the disc to determine whether or not the binary is able to create the structures observed in the disc. Our VLTI/GRAVITY astrometry gives excellent constraints on the orbit of HD 142527 B. We show that the secondary is following an orbit of semi-major axis a = 10.80 $\pm$ 0.22 au, with moderate eccentricity (e = 0.47 $\pm$ 0.01). With such a compact orbit, we show that HD 142527 B can only generate a gap and spiral arms of $\sim$30 au in the disc, which is much smaller than what is revealed by observations. Even from a theoretical standpoint, the observed cavity size of $\sim$100 au far exceeds even the most generous predictions for a companion like HD 142527 B on such a compact orbit. Thus, we conclude that the low-mass companion cannot be solely responsible for the observed morphology of the disc surrounding the system. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 11 pages, 6 figures. Accepted for publication in A&A

Estimating the construct validity of Principal Components Analysis

Authors: Thomas M. H. Hope, Cathy J. Price, Ajay Halai, Carola Salvi, Jenny Crinion, Merel Keijsers, Christoph Sperber, Howard Bowman

Abstract: In many scientific disciplines, the features of interest cannot be observed directly, so must instead be inferred from observed behaviour. Latent variable analyses are increasingly employed to systematise these inferences, and Principal Components Analysis (PCA) is perhaps the simplest and most popular of these methods. Here, we examine how the assumptions that we are prepared to entertain, about… ▽ More In many scientific disciplines, the features of interest cannot be observed directly, so must instead be inferred from observed behaviour. Latent variable analyses are increasingly employed to systematise these inferences, and Principal Components Analysis (PCA) is perhaps the simplest and most popular of these methods. Here, we examine how the assumptions that we are prepared to entertain, about the latent variable system, mediate the likelihood that PCA-derived components will capture the true sources of variance underlying data. As expected, we find that this likelihood is excellent in the best case, and robust to empirically reasonable levels of measurement noise, but best-case performance is also: (a) not robust to violations of the method's more prominent assumptions, of linearity and orthogonality; and also (b) requires that other subtler assumptions be made, such as that the latent variables should have varying importance, and that weights relating latent variables to observed data have zero mean. Neither variance explained, nor replication in independent samples, could reliably predict which (if any) PCA-derived components will capture true sources of variance in data. We conclude by describing a procedure to fit these inferences more directly to empirical data, and use it to find that components derived via PCA from two different empirical neuropsychological datasets, are less likely to have meaningful referents in the brain than we hoped. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 26 pages, 3 figures, 3 tables

Probing initial distributions of orbital eccentricity and disc misalignment via polar discs

Authors: Simone Ceppi, Nicolás Cuello, Giuseppe Lodato, Cristiano Longarini, Daniel J. Price, Daniel Elsender, Matthew R. Bate

Abstract: In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the di… ▽ More In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the distributions of initial eccentricity and mutual inclination in multiple stellar systems at birth. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: accepted for publication in A&A

Dust formation in common envelope binary interactions -- II: 3D simulations with self-consistent dust formation

Authors: Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Miguel González-Bolívar, Mike Y. M. Lau, Chunliang Mu, Ryosuke Hirai, Taïssa Danilovich, Mansi M. Kasliwal

Abstract: We performed numerical simulations of the common envelope (CE) interaction between two thermally-pulsing asymptotic giant branch (AGB) stars of 1.7 $M_\odot$ and 3.7 $M_\odot$, and their 0.6 $M_\odot$ compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust in the envelope is calculated explicitly, using a ca… ▽ More We performed numerical simulations of the common envelope (CE) interaction between two thermally-pulsing asymptotic giant branch (AGB) stars of 1.7 $M_\odot$ and 3.7 $M_\odot$, and their 0.6 $M_\odot$ compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust in the envelope is calculated explicitly, using a carbon dust nucleation network with a gas phase C/O number ratio of 2.5. By the end of the simulations, the total dust yield are $\sim8.2\times10^{-3}~M_\odot$ and $\sim2.2\times10^{-2}~M_\odot$ for the CE with a 1.7 $M_\odot$ and a 3.7 $M_\odot$ AGB star, respectively, close to the theoretical limit. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution. The first dust grains appear as early as $\sim1-3$ yrs after the onset of the CE rapidly forming an optically thick shell at $\sim10-20$ au, growing in thickness and radius to values of $\sim400-500$ au by $\sim40$ yrs. These large objects have approximate temperatures of 400 K. While dust yields are commensurate with those of single AGB stars of comparable mass, the dust in CE ejections forms over decades as opposed to tens of thousands of years. It is likely that these rapidly evolving IR objects correspond to the post-optically-luminous tail of the lightcurve of some luminous red novae. The simulated characteristics of dusty CEs also lend further support to the idea that extreme carbon stars and the so called ``water fountains" may be objects observed in the immediate aftermath of a CE event. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: 17 pages, 17 figures, 1 table

Optical Appearance of Eccentric Tidal Disruption Events

Authors: Fangyi, Hu, Daniel J. Price, Ilya Mandel

Abstract: Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($β= 1$ or $5$) using… ▽ More Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($β= 1$ or $5$) using general relativistic smoothed particle hydrodynamics. We follow the simulation for up to a year post-disruption. We show that accretion disks in eccentric TDEs are masked by unbound material outflowing at $\sim10,000~$km/s. Assuming electron scattering opacity, this material would be visible as a $\sim100~$au photosphere at $\sim10^4~$K, in line with observations of candidate TDEs. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 11 pages, 5 figures, submitted to ApJL

Journal ref: ApJL 963 L27 (2024)

The Automatic Identification and Tracking of Coronal Flux Ropes -- Part II: New Mathematical Morphology-based Flux Rope Extraction Method and Deflection Analysis

Authors: Andreas Wagner, Slava Bourgeois, Emilia K. J. Kilpua, Ranadeep Sarkar, Daniel J. Price, Anshu Kumari, Jens Pomoell, Stefaan Poedts, Teresa Barata, Robertus Erdélyi, Orlando Oliveira, Ricardo Gafeira

Abstract: We present a magnetic flux rope (FR) extraction tool for solar coronal magnetic field modelling data, which builds upon the methodology from Wagner et al. (2023). We apply the scheme to magnetic field simulations of active regions AR12473 and AR11176. We compare the method to its predecessor and study the 3D movement of the newly extracted FRs up to heights of 200 and 300 Mm, respectively. The ext… ▽ More We present a magnetic flux rope (FR) extraction tool for solar coronal magnetic field modelling data, which builds upon the methodology from Wagner et al. (2023). We apply the scheme to magnetic field simulations of active regions AR12473 and AR11176. We compare the method to its predecessor and study the 3D movement of the newly extracted FRs up to heights of 200 and 300 Mm, respectively. The extraction method is based on the twist parameter and a variety of mathematical morphology algorithms, including the opening transform and the morphological gradient. We highlight the differences between the methods by investigating the circularity of the FRs in the plane we extract from. The simulations for the active regions are carried out with a time-dependent data-driven magnetofrictional model (TMFM; Pomoell et al. (2019)). We investigate the FR trajectories by tracking their apex throughout the full simulation time span. We demonstrate that this upgraded methodology provides the user with more tools and less a-priori assumptions about the FR shape that, in turn, leads to a more accurate set of field lines. The propagation analysis yields that the erupting FR from AR12473 showcases stronger dynamics than the AR11176 FR and a significant deflection during its ascent through the domain. The AR11176 FR appears more stable, though there still is a notable deflection. This confirms that at these low coronal heights, FRs do undergo significant changes in the direction of their propagation even for less dynamic cases. The modelling results are also verified with observations, with AR12473 being indeed dynamic and eruptive, while AR11176 only features an eruption outside of our simulation time window. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Predicting recovery following stroke: deep learning, multimodal data and feature selection using explainable AI

Authors: Adam White, Margarita Saranti, Artur d'Avila Garcez, Thomas M. H. Hope, Cathy J. Price, Howard Bowman

Abstract: Machine learning offers great potential for automated prediction of post-stroke symptoms and their response to rehabilitation. Major challenges for this endeavour include the very high dimensionality of neuroimaging data, the relatively small size of the datasets available for learning, and how to effectively combine neuroimaging and tabular data (e.g. demographic information and clinical characte… ▽ More Machine learning offers great potential for automated prediction of post-stroke symptoms and their response to rehabilitation. Major challenges for this endeavour include the very high dimensionality of neuroimaging data, the relatively small size of the datasets available for learning, and how to effectively combine neuroimaging and tabular data (e.g. demographic information and clinical characteristics). This paper evaluates several solutions based on two strategies. The first is to use 2D images that summarise MRI scans. The second is to select key features that improve classification accuracy. Additionally, we introduce the novel approach of training a convolutional neural network (CNN) on images that combine regions-of-interest extracted from MRIs, with symbolic representations of tabular data. We evaluate a series of CNN architectures (both 2D and a 3D) that are trained on different representations of MRI and tabular data, to predict whether a composite measure of post-stroke spoken picture description ability is in the aphasic or non-aphasic range. MRI and tabular data were acquired from 758 English speaking stroke survivors who participated in the PLORAS study. The classification accuracy for a baseline logistic regression was 0.678 for lesion size alone, rising to 0.757 and 0.813 when initial symptom severity and recovery time were successively added. The highest classification accuracy 0.854 was observed when 8 regions-of-interest was extracted from each MRI scan and combined with lesion size, initial severity and recovery time in a 2D Residual Neural Network.Our findings demonstrate how imaging and tabular data can be combined for high post-stroke classification accuracy, even when the dataset is small in machine learning terms. We conclude by proposing how the current models could be improved to achieve even higher levels of accuracy using images from hospital scanners. △ Less

Submitted 29 October, 2023; originally announced October 2023.

Non-ergodicity in reinforcement learning: robustness via ergodicity transformations

Authors: Dominik Baumann, Erfaun Noorani, James Price, Ole Peters, Colm Connaughton, Thomas B. Schön

Abstract: Envisioned application areas for reinforcement learning (RL) include autonomous driving, precision agriculture, and finance, which all require RL agents to make decisions in the real world. A significant challenge hindering the adoption of RL methods in these domains is the non-robustness of conventional algorithms. In this paper, we argue that a fundamental issue contributing to this lack of robu… ▽ More Envisioned application areas for reinforcement learning (RL) include autonomous driving, precision agriculture, and finance, which all require RL agents to make decisions in the real world. A significant challenge hindering the adoption of RL methods in these domains is the non-robustness of conventional algorithms. In this paper, we argue that a fundamental issue contributing to this lack of robustness lies in the focus on the expected value of the return as the sole ``correct'' optimization objective. The expected value is the average over the statistical ensemble of infinitely many trajectories. For non-ergodic returns, this average differs from the average over a single but infinitely long trajectory. Consequently, optimizing the expected value can lead to policies that yield exceptionally high returns with probability zero but almost surely result in catastrophic outcomes. This problem can be circumvented by transforming the time series of collected returns into one with ergodic increments. This transformation enables learning robust policies by optimizing the long-term return for individual agents rather than the average across infinitely many trajectories. We propose an algorithm for learning ergodicity transformations from data and demonstrate its effectiveness in an instructive, non-ergodic environment and on standard RL benchmarks. △ Less

Submitted 10 April, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Reducing Training Data Needs with Minimal Multilevel Machine Learning (M3L)

Authors: Stefan Heinen, Danish Khan, Guido Falk von Rudorff, Konstantin Karandashev, Daniel Jose Arismendi Arrieta, Alastair J. A. Price, Surajit Nandi, Arghya Bhowmik, Kersti Hermansson, O. Anatole von Lilienfeld

Abstract: For many machine learning applications in science, data acquisition, not training, is the bottleneck even when avoiding experiments and relying on computation and simulation. Correspondingly, and in order to reduce cost and carbon footprint, training data efficiency is key. We introduce minimal multilevel machine learning (M3L) which optimizes training data set sizes using a loss function at multi… ▽ More For many machine learning applications in science, data acquisition, not training, is the bottleneck even when avoiding experiments and relying on computation and simulation. Correspondingly, and in order to reduce cost and carbon footprint, training data efficiency is key. We introduce minimal multilevel machine learning (M3L) which optimizes training data set sizes using a loss function at multiple levels of reference data in order to minimize a combination of prediction error with overall training data acquisition costs (as measured by computational wall-times). Numerical evidence has been obtained for calculated atomization energies and electron affinities of thousands of organic molecules at various levels of theory including HF, MP2, DLPNO-CCSD(T), DFHFCABS, PNOMP2F12, and PNOCCSD(T)F12, and treating tens with basis sets TZ, cc-pVTZ, and AVTZ-F12. Our M3L benchmarks for reaching chemical accuracy in distinct chemical compound sub-spaces indicate substantial computational cost reductions by factors of $\sim$ 1.01, 1.1, 3.8, 13.8 and 25.8 when compared to heuristic sub-optimal multilevel machine learning (M2L) for the data sets QM7b, QM9$^\mathrm{LCCSD(T)}$, EGP, QM9$^\mathrm{CCSD(T)}_\mathrm{AE}$, and QM9$^\mathrm{CCSD(T)}_\mathrm{EA}$, respectively. Furthermore, we use M2L to investigate the performance for 76 density functionals when used within multilevel learning and building on the following levels drawn from the hierarchy of Jacobs Ladder:~LDA, GGA, mGGA, and hybrid functionals. Within M2L and the molecules considered, mGGAs do not provide any noticeable advantage over GGAs. Among the functionals considered and in combination with LDA, the three on average top performing GGA and Hybrid levels for atomization energies on QM9 using M3L correspond respectively to PW91, KT2, B97D, and $τ$-HCTH, B3LYP$\ast$(VWN5), TPSSH. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Large-kernel Attention for Efficient and Robust Brain Lesion Segmentation

Authors: Liam Chalcroft, Ruben Lourenço Pereira, Mikael Brudfors, Andrew S. Kayser, Mark D'Esposito, Cathy J. Price, Ioannis Pappas, John Ashburner

Abstract: Vision transformers are effective deep learning models for vision tasks, including medical image segmentation. However, they lack efficiency and translational invariance, unlike convolutional neural networks (CNNs). To model long-range interactions in 3D brain lesion segmentation, we propose an all-convolutional transformer block variant of the U-Net architecture. We demonstrate that our model pro… ▽ More Vision transformers are effective deep learning models for vision tasks, including medical image segmentation. However, they lack efficiency and translational invariance, unlike convolutional neural networks (CNNs). To model long-range interactions in 3D brain lesion segmentation, we propose an all-convolutional transformer block variant of the U-Net architecture. We demonstrate that our model provides the greatest compromise in three factors: performance competitive with the state-of-the-art; parameter efficiency of a CNN; and the favourable inductive biases of a transformer. Our public implementation is available at https://github.com/liamchalcroft/MDUNet . △ Less

Submitted 14 August, 2023; originally announced August 2023.

Distinguishing Risk Preferences using Repeated Gambles

Authors: James Price, Colm Connaughton

Abstract: Sequences of repeated gambles provide an experimental tool to characterize the risk preferences of humans or artificial decision-making agents. The difficulty of this inference depends on factors including the details of the gambles offered and the number of iterations of the game played. In this paper we explore in detail the practical challenges of inferring risk preferences from the observed ch… ▽ More Sequences of repeated gambles provide an experimental tool to characterize the risk preferences of humans or artificial decision-making agents. The difficulty of this inference depends on factors including the details of the gambles offered and the number of iterations of the game played. In this paper we explore in detail the practical challenges of inferring risk preferences from the observed choices of artificial agents who are presented with finite sequences of repeated gambles. We are motivated by the fact that the strategy to maximize long-run wealth for sequences of repeated additive gambles (where gains and losses are independent of current wealth) is different to the strategy for repeated multiplicative gambles (where gains and losses are proportional to current wealth.) Accurate measurement of risk preferences would be needed to tell whether an agent is employing the optimal strategy or not. To generalize the types of gambles our agents face we use the Yeo-Johnson transformation, a tool borrowed from feature engineering for time series analysis, to construct a family of gambles that interpolates smoothly between the additive and multiplicative cases. We then analyze the optimal strategy for this family, both analytically and numerically. We find that it becomes increasingly difficult to distinguish the risk preferences of agents as their wealth increases. This is because agents with different risk preferences eventually make the same decisions for sufficiently high wealth. We believe that these findings are informative for the effective design of experiments to measure risk preferences in humans. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm

Authors: D. P. Aguillard, T. Albahri, D. Allspach, A. Anisenkov, K. Badgley, S. Baeßler, I. Bailey, L. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, E. Barzi, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, S. Braun, M. Bressler, G. Cantatore, R. M. Carey, B. C. K. Casey , et al. (166 additional authors not shown)

Abstract: We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable… ▽ More We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable beam, and improved knowledge of the magnetic field weighted by the muon distribution, $\tildeω'^{}_p$, and of the anomalous precession frequency corrected for beam dynamics effects, $ω_a$. From the ratio $ω_a / \tildeω'^{}_p$, together with precisely determined external parameters, we determine $a_μ= 116\,592\,057(25) \times 10^{-11}$ (0.21 ppm). Combining this result with our previous result from the 2018 data, we obtain $a_μ\text{(FNAL)} = 116\,592\,055(24) \times 10^{-11}$ (0.20 ppm). The new experimental world average is $a_μ(\text{Exp}) = 116\,592\,059(22)\times 10^{-11}$ (0.19 ppm), which represents a factor of 2 improvement in precision. △ Less

Submitted 4 October, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 8 pages, 3 figures

Report number: FERMILAB-PUB-23-385-AD-CSAID-PPD

Journal ref: Phys. Rev. Lett. 131, 161802 (2023)

Inhomogeneous Cosmology using General Relativistic Smoothed Particle Hydrodynamics coupled to Numerical Relativity

Authors: Spencer J. Magnall, Daniel J. Price, Paul D. Lasky, Hayley J. Macpherson

Abstract: We perform three-dimensional simulations of homogeneous and inhomogeneous cosmologies via the coupling of a numerical relativity code for spacetime evolution and smoothed particle hydrodynamics (SPH) code. Evolution of a flat dust and radiation dominated Friedmann-Lemaître-Roberston-Walker (FLRW) spacetime shows an agreement of exact solutions with residuals on the order $10^{-6}$ and $10^{-3}$ re… ▽ More We perform three-dimensional simulations of homogeneous and inhomogeneous cosmologies via the coupling of a numerical relativity code for spacetime evolution and smoothed particle hydrodynamics (SPH) code. Evolution of a flat dust and radiation dominated Friedmann-Lemaître-Roberston-Walker (FLRW) spacetime shows an agreement of exact solutions with residuals on the order $10^{-6}$ and $10^{-3}$ respectively, even at low grid resolutions. We demonstrate evolution of linear perturbations of density, velocity and metric quantities to the FLRW with residuals of only $10^{-2}$ compared to exact solutions. Finally, we demonstrate the evolution of non-linear perturbations of the metric past shell-crossing, such that dark matter halo formation is possible. We show that numerical relativistic smoothed particle hydrodynamics is a viable method for understanding non-linear effects in cosmology. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 15 pages, 12 figures, submitted to PRD

Second-Order Coherence Across the Brillouin Lasing Threshold

Authors: E. A. Cryer-Jenkins, G. Enzian, L. Freisem, N. Moroney, J. J. Price, A. Ø. Svela, K. D. Major, M. R. Vanner

Abstract: Brillouin-Mandelstam scattering is one of the most accessible nonlinear optical phenomena and has been widely studied since its theoretical discovery one hundred years ago. The scattering mechanism is a three-wave mixing process between two optical fields and one acoustic field and has found a broad range of applications spanning microscopy to ultra-narrow-linewidth lasers. Building on the success… ▽ More Brillouin-Mandelstam scattering is one of the most accessible nonlinear optical phenomena and has been widely studied since its theoretical discovery one hundred years ago. The scattering mechanism is a three-wave mixing process between two optical fields and one acoustic field and has found a broad range of applications spanning microscopy to ultra-narrow-linewidth lasers. Building on the success of utilizing this nonlinearity at a classical level, a rich avenue is now being opened to explore Brillouin scattering within the paradigm of quantum optics. Here, we take a key step in this direction by employing quantum optical techniques yet to be utilized for Brillouin scattering to characterize the second-order coherence of Stokes scattering across the Brillouin lasing threshold. We use a silica microsphere resonator and single-photon counters to observe the expected transition from bunched statistics of thermal light below the lasing threshold to Poissonian statistics of coherent light above the threshold. Notably, at powers approaching the lasing threshold, we also observe super-thermal statistics, which arise due to instability and a ``flickering'' in and out of lasing as the pump field is transiently depleted. The statistics observed across the transition, including the ``flickering'', are a result of the full nonlinear three-wave mixing process and cannot be captured by a linearized model. These measurements are in good agreement with numerical solutions of the three-wave Langevin equations and are well demarcated by analytical expressions for the instability and the lasing thresholds. These results demonstrate that applying second-order-coherence and photon-counting measurements to Brillouin scattering provides new methods to advance our understanding of Brillouin scattering itself and progress toward quantum-state preparation and characterization of acoustic modes. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: Main (8 pages, 2 figures) + Supplementary (6 pages, 1 figures), Submitted

Beam Spin Asymmetry Measurements of Deeply Virtual $π^0$ Production with CLAS12

Authors: A. Kim, S. Diehl, K. Joo, V. Kubarovsky, P. Achenbach, Z. Akbar, J. S. Alvarado, Whitney R. Armstrong, H. Atac, H. Avakian, C. Ayerbe Gayoso, L. Barion, M. Battaglieri, I. Bedlinskiy, B. Benkel, A. Bianconi, A. S. Biselli, M. Bondi, F. Bossù, S. Boiarinov, K. T. Brinkmann, W. J. Briscoe, W. K. Brooks, S. Bueltmann, V. D. Burkert , et al. (132 additional authors not shown)

Abstract: The new experimental measurements of beam spin asymmetry were performed for the deeply virtual exclusive $π^0$ production in a wide kinematic region with the photon virtualities $Q^2$ up to 8 GeV$^2$ and the Bjorken scaling variable $x_B$ in the valence regime. The data were collected by the CEBAF Large Acceptance Spectrometer (CLAS12) at Jefferson Lab with longitudinally polarized 10.6 GeV electr… ▽ More The new experimental measurements of beam spin asymmetry were performed for the deeply virtual exclusive $π^0$ production in a wide kinematic region with the photon virtualities $Q^2$ up to 8 GeV$^2$ and the Bjorken scaling variable $x_B$ in the valence regime. The data were collected by the CEBAF Large Acceptance Spectrometer (CLAS12) at Jefferson Lab with longitudinally polarized 10.6 GeV electrons scattered on an unpolarized liquid-hydrogen target. Sizable asymmetry values indicate a substantial contribution from transverse virtual photon amplitudes to the polarized structure functions.The interpretation of these measurements in terms of the Generalized Parton Distributions (GPDs) demonstrates their sensitivity to the chiral-odd GPD $\bar E_T$, which contains information on quark transverse spin densities in unpolarized and polarized nucleons and provides access to the proton's transverse anomalous magnetic moment. Additionally, the data were compared to a theoretical model based on a Regge formalism that was extended to the high photon virtualities. △ Less

Submitted 15 July, 2023; originally announced July 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2210.14557

Tidal truncation of circumplanetary disks fails above a critical disk aspect ratio

Authors: Rebecca G. Martin, Philip J. Armitage, Stephen H. Lubow, Daniel J. Price

Abstract: We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodyn… ▽ More We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodynamics code until a steady-state is reached. We find that the most significant dependence of the truncation boundary is on the disk aspect ratio $H/R$. Circumplanetary disks are efficiently truncated for $H/R \lesssim 0.2$. For $H/R \simeq 0.3$, up to about half of the injected mass, depending on the injection radius, flows outwards through the decretion disk and escapes. As expected from analytic arguments, the conditions ($H/R$ and Shakura-Sunyaev $α$) required for tidal truncation are independent of planet mass. A simulation with larger $α=0.1$ shows stronger outflow than one with $α=0.01$, but the dependence on transport efficiency is less important than variations of $H/R$. Our results suggest two distinct classes of circumplanetary disks: tidally truncated thin disks with dust-poor outer regions, and thicker actively decreting disks with enhanced dust-to-gas ratios. Applying our results to the PDS 70c system, we predict a largely truncated circumplanetary disk, but it is possible that enough mass escapes to support an outward flow of dust that could explain the observed disk size. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

Dust Formation in Common Envelope Binary Interaction -- I: 3D Simulations Using the Bowen Approximation

Authors: Miguel González-Bolívar, Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Mansi Kasliwal

Abstract: We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried… ▽ More We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried out simulations using both an ideal gas and a tabulated equations of state, with the latter considering the recombination energy of the envelope. We found that the dust-driven wind leads to a relatively small increase in the unbound gas, with the effect being smaller for the tabulated equation of state simulations and for the more massive primary. Dust acceleration does contribute to envelope expansion with only a slightly elongated morphology, if we believe the results from the tabulated equation of state as more reliable. The Bowen opacities in the outer envelopes of the two models, at late times, are large enough that the photosphere of the post-inspiral object is about ten times larger compared to the same without accounting for the dust opacities. As such, the prediction of the appearance of the transient would change substantially if dust is included. △ Less

Submitted 13 April, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 14 Pages, 13 Figures. Accepted for publication in MNRAS

An Integral-based Approach for the Vector Potential in Smoothed Particle Magnetohydrodynamics

Authors: Terrence S. Tricco, Daniel J. Price

Abstract: A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derive… ▽ More A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derived, but are not found to solve the MHD equations in the continuum limit. Tests are performed using a hybrid approach instead, whereby the equations of motion based on the magnetic field instead of vector potential are used. Test results experience the same numerical instability as with the Price (2010) formulation. We conclude that this new formulation is non-viable. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 7 pages, 5 figures, Proceedings of the "17th SPHERIC International Workshop", Rhodes, Greece, Jun 27-29, 2023

The Automatic Identification and Tracking of Coronal Flux Ropes -- Part I: Footpoints and Fluxes

Authors: Andreas Wagner, Emilia K. J. Kilpua, Ranadeep Sarkar, Daniel J. Price, Anshu Kumari, Farhad Daei, Jens Pomoell, Stefaan Poedts

Abstract: Investigating the early-stage evolution of an erupting flux rope from the Sun is important to understand the mechanisms of how it looses its stability and its space weather impacts. Our aim is to develop an efficient scheme for tracking the early dynamics of erupting solar flux ropes and use the algorithm to analyse its early-stage properties. The algorithm is tested on a data-driven simulation of… ▽ More Investigating the early-stage evolution of an erupting flux rope from the Sun is important to understand the mechanisms of how it looses its stability and its space weather impacts. Our aim is to develop an efficient scheme for tracking the early dynamics of erupting solar flux ropes and use the algorithm to analyse its early-stage properties. The algorithm is tested on a data-driven simulation of an eruption that took place in active region AR12473. We investigate the modelled flux rope's footpoint movement and magnetic flux evolution and compare with observational data from the Solar Dynamics Observatory's Atmospheric Imaging Assembly in the 211 $\unicode{x212B}$ and 1600 $\unicode{x212B}$ channels. To carry out our analysis, we use the time-dependent data-driven magnetofrictional model (TMFM). We also perform another modelling run, where we stop the driving of the TMFM midway through the flux rope's rise through the simulation domain and evolve it instead with a zero-beta magnetohydrodynamic (MHD) approach. The developed algorithm successfully extracts a flux rope and its ascend through the simulation domain. We find that the movement of the modelled flux rope footpoints showcases similar trends in both TMFM and relaxation MHD run: they recede from their respective central location as the eruption progresses and the positive polarity footpoint region exhibits a more dynamic behaviour. The ultraviolet brightenings and extreme ultraviolet dimmings agree well with the models in terms of their dynamics. According to our modelling results, the toroidal magnetic flux in the flux rope first rises and then decreases. In our observational analysis, we capture the descending phase of toroidal flux. In conclusion, the extraction algorithm enables us to effectively study the flux rope's early dynamics and derive some of its key properties such as footpoint movement and toroidal magnetic flux. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A81 (2023)

TauPETGen: Text-Conditional Tau PET Image Synthesis Based on Latent Diffusion Models

Authors: Se-In Jang, Cristina Lois, Emma Thibault, J. Alex Becker, Yafei Dong, Marc D. Normandin, Julie C. Price, Keith A. Johnson, Georges El Fakhri, Kuang Gong

Abstract: In this work, we developed a novel text-guided image synthesis technique which could generate realistic tau PET images from textual descriptions and the subject's MR image. The generated tau PET images have the potential to be used in examining relations between different measures and also increasing the public availability of tau PET datasets. The method was based on latent diffusion models. Both… ▽ More In this work, we developed a novel text-guided image synthesis technique which could generate realistic tau PET images from textual descriptions and the subject's MR image. The generated tau PET images have the potential to be used in examining relations between different measures and also increasing the public availability of tau PET datasets. The method was based on latent diffusion models. Both textual descriptions and the subject's MR prior image were utilized as conditions during image generation. The subject's MR image can provide anatomical details, while the text descriptions, such as gender, scan time, cognitive test scores, and amyloid status, can provide further guidance regarding where the tau neurofibrillary tangles might be deposited. Preliminary experimental results based on clinical [18F]MK-6240 datasets demonstrate the feasibility of the proposed method in generating realistic tau PET images at different clinical stages. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Strong Interaction Physics at the Luminosity Frontier with 22 GeV Electrons at Jefferson Lab

Authors: A. Accardi, P. Achenbach, D. Adhikari, A. Afanasev, C. S. Akondi, N. Akopov, M. Albaladejo, H. Albataineh, M. Albrecht, B. Almeida-Zamora, M. Amaryan, D. Androić, W. Armstrong, D. S. Armstrong, M. Arratia, J. Arrington, A. Asaturyan, A. Austregesilo, H. Avagyan, T. Averett, C. Ayerbe Gayoso, A. Bacchetta, A. B. Balantekin, N. Baltzell, L. Barion , et al. (419 additional authors not shown)

Abstract: This document presents the initial scientific case for upgrading the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) to 22 GeV. It is the result of a community effort, incorporating insights from a series of workshops conducted between March 2022 and April 2023. With a track record of over 25 years in delivering the world's most intense and precise multi-GeV electron… ▽ More This document presents the initial scientific case for upgrading the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) to 22 GeV. It is the result of a community effort, incorporating insights from a series of workshops conducted between March 2022 and April 2023. With a track record of over 25 years in delivering the world's most intense and precise multi-GeV electron beams, CEBAF's potential for a higher energy upgrade presents a unique opportunity for an innovative nuclear physics program, which seamlessly integrates a rich historical background with a promising future. The proposed physics program encompass a diverse range of investigations centered around the nonperturbative dynamics inherent in hadron structure and the exploration of strongly interacting systems. It builds upon the exceptional capabilities of CEBAF in high-luminosity operations, the availability of existing or planned Hall equipment, and recent advancements in accelerator technology. The proposed program cover various scientific topics, including Hadron Spectroscopy, Partonic Structure and Spin, Hadronization and Transverse Momentum, Spatial Structure, Mechanical Properties, Form Factors and Emergent Hadron Mass, Hadron-Quark Transition, and Nuclear Dynamics at Extreme Conditions, as well as QCD Confinement and Fundamental Symmetries. Each topic highlights the key measurements achievable at a 22 GeV CEBAF accelerator. Furthermore, this document outlines the significant physics outcomes and unique aspects of these programs that distinguish them from other existing or planned facilities. In summary, this document provides an exciting rationale for the energy upgrade of CEBAF to 22 GeV, outlining the transformative scientific potential that lies within reach, and the remarkable opportunities it offers for advancing our understanding of hadron physics and related fundamental phenomena. △ Less

Submitted 24 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: Updates to the list of authors; Preprint number changed from theory to experiment; Updates to sections 4 and 6, including additional figures

Report number: JLAB-PHY-23-3840

Observational Signatures of Circumbinary Discs I: Kinematics

Authors: Josh Calcino, Daniel J. Price, Christophe Pinte, Himanshi Garg, Brodie J. Norfolk, Valentin Christiaens, Hui Li, Richard Teague

Abstract: We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by i) a central cavity, ii) spiral arms both in and outside of their central cavities, iii) non-localised perturbations… ▽ More We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by i) a central cavity, ii) spiral arms both in and outside of their central cavities, iii) non-localised perturbations in their iso-velocity curves, iv) asymmetry between the lines of maximum speed of the blue and red-shifted wings and v) asymmetry between the area of the blue and red-shifted wings. We provide quantitative metrics for the last two criteria that can be used, in conjunction with the morphological criteria, to signal whether a protoplanetary disc is likely to be a circumbinary disc. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 26 pages, 26 figures, accepted to MNRAS

CENSUS-HWR: a large training dataset for offline handwriting recognition

Authors: Chetan Joshi, Lawry Sorenson, Ammon Wolfert, Dr. Mark Clement, Dr. Joseph Price, Dr. Kasey Buckles

Abstract: Progress in Automated Handwriting Recognition has been hampered by the lack of large training datasets. Nearly all research uses a set of small datasets that often cause models to overfit. We present CENSUS-HWR, a new dataset consisting of full English handwritten words in 1,812,014 gray scale images. A total of 1,865,134 handwritten texts from a vocabulary of 10,711 words in the English language… ▽ More Progress in Automated Handwriting Recognition has been hampered by the lack of large training datasets. Nearly all research uses a set of small datasets that often cause models to overfit. We present CENSUS-HWR, a new dataset consisting of full English handwritten words in 1,812,014 gray scale images. A total of 1,865,134 handwritten texts from a vocabulary of 10,711 words in the English language are present in this collection. This dataset is intended to serve handwriting models as a benchmark for deep learning algorithms. This huge English handwriting recognition dataset has been extracted from the US 1930 and 1940 censuses taken by approximately 70,000 enumerators each year. The dataset and the trained model with their weights are freely available to download at https://censustree.org/data.html. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Effects of optimisation parameters on data-driven magnetofrictional modelling of active regions

Authors: A. Kumari, D. J. Price, F. Daei, J. Pomoell, E. K. J. Kilpua

Abstract: Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing th… ▽ More Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing the eruption dynamics. We present a detailed study of the effects of optimisation of the non-inductive electric field on the TMFM of AR12473. We aim to study the effects of varying the non-inductive electric field on the data-driven coronal simulations, for two alternative parametrisations. By varying parameters controlling the strength of the non-inductive electric field, we wish to explore the changes in flux rope formation and their early evolution and other parameters, for instance, axial flux and magnetic field magnitude.The non-inductive electric field component in the photosphere is critical for energising and introducing twist to the coronal magnetic field, thereby allowing unstable configurations to be formed. We estimated this component using an approach based on optimising the injection of magnetic energy. However, the flux rope formation, evolution and eruption time varies depending on the values of the optimisation parameters. The flux rope is formed and has overall similar evolution and properties with a large range of non-inductive electric fields needed to determine the non-inductive electric field component that is critical for energising and introducing twist to the coronal magnetic field. This study shows that irrespective of non-inductive electric field values, flux ropes are formed and erupted, which indicates that data-driven TMFM can be used to estimate flux rope properties early in their evolution without employing a lengthy optimisation process. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures, 1 table. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 675, A80 (2023)

The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations

Authors: Cristiano Longarini, Philip J. Armitage, Giuseppe Lodato, Daniel J. Price, Simone Ceppi

Abstract: Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed p… ▽ More Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS, 20 pages

First measurement of hard exclusive $π^- Δ^{++}$ electroproduction beam-spin asymmetries off the proton

Authors: S. Diehl, N. Trotta, K. Joo, P. Achenbach, Z. Akbar, W. R. Armstrong, H. Atac, H. Avakian, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, F. Benmokhtar, A. Bianconi, A. S. Biselli, F. Bossu, K. -T. Brinkmann, W. J. Briscoe, D. Bulumulla, V. Burkert, R. Capobianco, D. S. Carman, J. C. Carvajal , et al. (120 additional authors not shown)

Abstract: The polarized cross section ratio $σ_{LT'}/σ_{0}$ from hard exclusive $π^{-} Δ^{++}$ electroproduction off an unpolarized hydrogen target has been extracted based on beam-spin asymmetry measurements using a 10.2 GeV / 10.6 GeV incident electron beam and the CLAS12 spectrometer at Jefferson Lab. The study, which provides the first observation of this channel in the deep-inelastic regime, focuses on… ▽ More The polarized cross section ratio $σ_{LT'}/σ_{0}$ from hard exclusive $π^{-} Δ^{++}$ electroproduction off an unpolarized hydrogen target has been extracted based on beam-spin asymmetry measurements using a 10.2 GeV / 10.6 GeV incident electron beam and the CLAS12 spectrometer at Jefferson Lab. The study, which provides the first observation of this channel in the deep-inelastic regime, focuses on very forward-pion kinematics in the valence regime, and photon virtualities ranging from 1.5 GeV$^{2}$ up to 7 GeV$^{2}$. The reaction provides a novel access to the $d$-quark content of the nucleon and to $p \rightarrow Δ^{++}$ transition generalized parton distributions. A comparison to existing results for hard exclusive $π^{+} n$ and $π^{0} p$ electroproduction is provided, which shows a clear impact of the excitation mechanism, encoded in transition generalized parton distributions, on the asymmetry. △ Less

Submitted 21 June, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Optimal interruption of P. vivax malaria transmission using mass drug administration

Authors: Md Nurul Anwar, Roslyn I. Hickson, Somya Mehra, David J. Price, James M. McCaw, Mark B. Flegg, Jennifer A. Flegg

Abstract: \textit{Plasmodium vivax} is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as hypnozoites, they reside in the liver following an initial exposure, before activating later to cause further infect… ▽ More \textit{Plasmodium vivax} is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as hypnozoites, they reside in the liver following an initial exposure, before activating later to cause further infections, referred to as relapses. As around 79-96$\%$ of infections are attributed to relapses, we expect it will be highly impactful to apply treatment to target the hypnozoite reservoir to eliminate \textit{P. vivax}. Treatment with a radical cure to target the hypnozoite reservoir is a potential tool to control or eliminate \textit{P. vivax}. We have developed a multiscale mathematical model as a system of integro-differential equations that captures the complex dynamics of \textit{P. vivax} hypnozoites and the effect of hypnozoite relapse on disease transmission. Here, we use our model to study the anticipated effect of radical cure treatment administered via a mass drug administration (MDA) program. We implement multiple rounds of MDA with a fixed interval between rounds, starting from different steady-state disease prevalences. We then construct an optimisation model to obtain the optimal MDA interval. We also incorporate mosquito seasonality in our model to study its effect on the optimal treatment regime. We find that the effect of MDA interventions is temporary and depends on the pre-intervention disease prevalence (and choice of model parameters) as well as the number of MDA rounds under consideration. We find radical cure alone may not be enough to lead to \textit{P. vivax} elimination under our mathematical model (and choice of model parameters) since the prevalence of infection eventually returns to pre-MDA levels. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 39 pages, 14 figures

MSC Class: 92D30

Confirmation and Keplerian motion of the gap-carving protoplanet HD 169142 b

Authors: Iain Hammond, Valentin Christiaens, Daniel J. Price, Claudia Toci, Christophe Pinte, Sandrine Juillard, Himanshi Garg

Abstract: We present the re-detection of a compact source in the face-on protoplanetary disc surrounding HD 169142, using VLT/SPHERE data in YJH bands. The source is found at a separation of 0.''319 ($\sim$37 au) from the star. Three lines of evidence argue in favour of the signal tracing a protoplanet: (i) it is found in the annular gap separating the two bright rings of the disc, as predicted by theory; (… ▽ More We present the re-detection of a compact source in the face-on protoplanetary disc surrounding HD 169142, using VLT/SPHERE data in YJH bands. The source is found at a separation of 0.''319 ($\sim$37 au) from the star. Three lines of evidence argue in favour of the signal tracing a protoplanet: (i) it is found in the annular gap separating the two bright rings of the disc, as predicted by theory; (ii) it is moving at the expected Keplerian velocity for an object at $\sim$37 au in the 2015, 2017 and 2019 datasets; (iii) we also detect a spiral-shaped signal whose morphology is consistent with the expected outer spiral wake triggered by a planet in the gap, based on dedicated hydrodynamical simulations of the system. The YJH colours we extracted for the object are consistent with tracing scattered starlight, suggesting that the protoplanet is enshrouded in a significant amount of dust, as expected for a circumplanetary disc or envelope surrounding a gap-clearing Jovian-mass protoplanet. △ Less

Submitted 23 February, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures

Kinematic and thermal signatures of the directly imaged protoplanet candidate around Elias 2-24

Authors: Christophe Pinte, Iain Hammond, Daniel J. Price, Valentin Christiaens, Sean M. Andrews, Gaël Chauvin, Laura M. Pérez, Sebastián Jorquera, Himanshi Garg, Brodie J. Norfolk, Josh Calcino, Mickaël Bonnefoy

Abstract: We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require loc… ▽ More We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models show this could either be a $\gtrsim 5$M$_\mathrm{Jup}$, or a lower mass ( $\gtrsim 2$M$_\mathrm{Jup}$) but accreting proto-planet. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Dynamics of dust grains in turbulent molecular clouds. Conditions for decoupling and limits of different numerical implementations

Authors: Benoît Commerçon, Ugo Lebreuilly, Daniel J. Price, Francesco Lovascio, Guillaume Laibe, Patrick Hennebelle

Abstract: Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of… ▽ More Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of two fundamentally different numerical implementations for the evolution of dust and gas mixtures in turbulent molecular clouds. We carried out a suite of numerical experiments using two different schemes. First, we used a monofluid formalism in the terminal velocity approximation (TVA) on a Eulerian grid. Second, we used a two-fluid scheme, in which the dust dynamics is handled with Lagrangian super-particles, and the gas dynamics on a Eulerian grid. The monofluid results are in good agreement with the theoretical critical size for decoupling. We report dust dynamics decoupling for Stokes number St>0.1, that is, dust grains of $s>4~μ$m in size. We find that the TVA is well suited for grain sizes of 10 $μ$m in molecular clouds, in particular in the densest regions. However, the maximum dust enrichment measured in the low-density material where St>1 is questionable. In the Lagrangian dust experiments, we show that the results are affected by the numerics for all dust grain sizes. At St<<1, the dust dynamics is largely affected by artificial trapping in the high-density regions, leading to spurious variations of the dust concentration. At St>1, the maximum dust enrichment is regulated by the grid resolution used for the gas dynamics. The results of previous similar numerical work should therefore be revisited with respect to the limitations we highlight in this study. Dust enrichment of submicron dust grains is unlikely to occur in the densest parts of molecular clouds. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 22 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 671, A128 (2023)

First CLAS12 measurement of DVCS beam-spin asymmetries in the extended valence region

Authors: CLAS Collaboration, G. Christiaens, M. Defurne, D. Sokhan, P. Achenbach, Z. Akbar, M. J. Amaryan, H. Atac, H. Avakian, C. Ayerbe Gayoso, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth, F. Bossù, S. Boiarinov, K. -Th. Brinkmann , et al. (146 additional authors not shown)

Abstract: Deeply virtual Compton scattering (DVCS) allows one to probe Generalized Parton Distributions (GPDs) describing the 3D structure of the nucleon. We report the first measurement of the DVCS beam-spin asymmetry using the CLAS12 spectrometer with a 10.2 and 10.6 GeV electron beam scattering from unpolarised protons. The results greatly extend the $Q^2$ and Bjorken-$x$ phase space beyond the existing… ▽ More Deeply virtual Compton scattering (DVCS) allows one to probe Generalized Parton Distributions (GPDs) describing the 3D structure of the nucleon. We report the first measurement of the DVCS beam-spin asymmetry using the CLAS12 spectrometer with a 10.2 and 10.6 GeV electron beam scattering from unpolarised protons. The results greatly extend the $Q^2$ and Bjorken-$x$ phase space beyond the existing data in the valence region and provide over 2000 new data points measured with unprecedented statistical uncertainty, setting new, tight constraints for future phenomenological studies. △ Less

Submitted 2 December, 2022; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: Revised Figure 4 and discussion around the number of effective ANNs after Bayesian reweighting

Report number: JLAB-PHY-22-3752

Beam-Based Diagnostics of Electric Guide Fields and Lattice Parameters for Run-1 of the Muon g-2 Storage Ring at Fermilab

Authors: David Tarazona, Martin Berz, Jason Crnkovic, Kim-Siang Khaw, Kyoko Makino, James Mott, Joseph Price, Michael Syphers, Vladimir Tishchenko

Abstract: A portion of the Muon g-2 Storage Ring electric system, which provides vertical beam focusing, exhibited an unexpected time dependence that produced a characteristic evolution of the stored beam during Run-1 of the Muon g-2 Experiment at Fermilab (E989). A method to reconstruct the Run-1 electric guide fields has been developed, which is based on a numerical model of the muon storage ring and opti… ▽ More A portion of the Muon g-2 Storage Ring electric system, which provides vertical beam focusing, exhibited an unexpected time dependence that produced a characteristic evolution of the stored beam during Run-1 of the Muon g-2 Experiment at Fermilab (E989). A method to reconstruct the Run-1 electric guide fields has been developed, which is based on a numerical model of the muon storage ring and optimization algorithms supported by COSY INFINITY. This method takes beam profile measurements from the Muon g-2 straw tracking detectors as input, and it produces a full reconstruction of the time-dependent fields. The fields can then be used for the reproduction of detailed beam tracking simulations and the calculation of ring lattice parameters for acceptance studies and systematic error evaluations. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: IPAC 2022

Report number: FERMILAB-CONF-22-683-AD

Hot Jupiter engulfment by a red giant in 3D hydrodynamics

Authors: Mike Y. M. Lau, Matteo Cantiello, Adam S. Jermyn, Morgan MacLeod, Ilya Mandel, Daniel J. Price

Abstract: Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been s… ▽ More Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been studied extensively as it may account for observed rapidly rotating and chemically enriched giant stars. We perform 3D hydrodynamical simulations of hot Jupiter engulfment by a 1 solar mass, 4 solar radii early red giant. Our "global" simulations simultaneously resolve the stellar envelope and planetary structure, modelling the hot Jupiter as a polytropic gas sphere. We find that approximately 90% of the hot Jupiter's mass is ablated in the convective part of the giant envelope, which would enhance the surface lithium abundance by 0.1 dex. The hot Jupiter is disrupted by a combination of ram pressure and tidal forces near the base of the convective envelope, with the deepest material penetrating to the radiative zone. The star experiences modest spin-up (~1 km/s), although engulfing a more massive companion could produce a rapidly rotating giant. Drag heating near the surface could exceed the unperturbed stellar luminosity and power an optical transient. For the amount of unbound ejecta recorded in the simulation, H-recombination could also power a transient that is around ten times the pre-engulfment luminosity, for several days. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 14 pages, 11 figures

A multidimensional study of the structure function ratio $σ_{LT'}/σ_{0}$ from hard exclusive $π^+$ electro-production off protons in the GPD regime

Authors: S. Diehl, A. Kim, K. Joo, P. Achenbach, Z. Akbar, M. J. Amaryan, H. Atac, H. Avagyan, C. Ayerbe Gayoso, L. Baashen, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth, F. Bossu, S. Boiarinov, K. -Th. Brinkmann, W. J. Briscoe, S. Bueltmann , et al. (129 additional authors not shown)

Abstract: A multidimensional extraction of the structure function ratio $σ_{LT'}/σ_{0}$ from the hard exclusive $\vec{e} p \to e^\prime n π^+$ reaction above the resonance region has been performed. The study was done based on beam-spin asymmetry measurements using a 10.6 GeV incident electron beam on a liquid-hydrogen target and the CLAS12 spectrometer at Jefferson Lab. The measurements focus on the very f… ▽ More A multidimensional extraction of the structure function ratio $σ_{LT'}/σ_{0}$ from the hard exclusive $\vec{e} p \to e^\prime n π^+$ reaction above the resonance region has been performed. The study was done based on beam-spin asymmetry measurements using a 10.6 GeV incident electron beam on a liquid-hydrogen target and the CLAS12 spectrometer at Jefferson Lab. The measurements focus on the very forward regime ($t/Q^{2}$ $\ll$ 1) with a wide kinematic range of $x_{B}$ in the valence regime (0.17 $<$ $x_{B}$ $<$ 0.55), and virtualities $Q^{2}$ ranging from 1.5 GeV$^{2}$ up to 6 GeV$^{2}$. The results and their comparison to theoretical models based on Generalized Parton Distributions demonstrate the sensitivity to chiral-odd GPDs and the directly related tensor charge of the nucleon. In addition, the data is compared to an extension of a Regge formalism at high photon virtualities. It was found that the Regge model provides a better description at low $Q^{2}$, while the GPD model is more appropriate at high $Q^{2}$. △ Less

Submitted 7 February, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: arXiv admin note: text overlap with arXiv:2007.15677

First Measurement of $Λ$ Electroproduction off Nuclei in the Current and Target Fragmentation Regions

Authors: T. Chetry, L. El Fassi, W. K. Brooks, R. Dupré, A. El Alaoui, K. Hafidi, P. Achenbach, K. P. Adhikari, Z. Akbar, W. R. Armstrong, M. Arratia, H. Atac, H. Avakian, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth , et al. (129 additional authors not shown)

Abstract: We report results of $Λ$ hyperon production in semi-inclusive deep-inelastic scattering off deuterium, carbon, iron, and lead targets obtained with the CLAS detector and the Continuous Electron Beam Accelerator Facility 5.014~GeV electron beam. These results represent the first measurements of the $Λ$ multiplicity ratio and transverse momentum broadening as a function of the energy fraction~($z$)… ▽ More We report results of $Λ$ hyperon production in semi-inclusive deep-inelastic scattering off deuterium, carbon, iron, and lead targets obtained with the CLAS detector and the Continuous Electron Beam Accelerator Facility 5.014~GeV electron beam. These results represent the first measurements of the $Λ$ multiplicity ratio and transverse momentum broadening as a function of the energy fraction~($z$) in the current and target fragmentation regions. The multiplicity ratio exhibits a strong suppression at high~$z$~and~an enhancement at~low~$z$. The measured transverse momentum broadening is an order of magnitude greater than that seen for light mesons. This indicates that the propagating entity interacts very strongly with the nuclear medium, which suggests that propagation of diquark configurations in the nuclear medium takes place at least part of the time, even at high~$z$. The trends of these results are qualitatively described by the Giessen Boltzmann-Uehling-Uhlenbeck transport model, particularly for the multiplicity ratios. These observations will potentially open a new era of studies of the structure of the nucleon as well as of strange baryons. △ Less

Submitted 1 April, 2023; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 14 pages, 6 figures

Solar Energetic Particle Time Series Analysis with Python

Authors: Christian Palmroos, Jan Gieseler, Nina Dresing, Diana E. Morosan, Eleanna Asvestari, Aleksi Yli-Laurila, Daniel J. Price, Saku Valkila, Rami Vainio

Abstract: Solar Energetic Particles (SEPs) are charged particles accelerated within the solar atmosphere or the interplanetary space by explosive phenomena such as solar flares or Coronal Mass Ejections (CMEs). Once injected into the interplanetary space, they can propagate towards Earth, causing space weather related phenomena. For their analysis, interplanetary in-situ measurements of charged particles ar… ▽ More Solar Energetic Particles (SEPs) are charged particles accelerated within the solar atmosphere or the interplanetary space by explosive phenomena such as solar flares or Coronal Mass Ejections (CMEs). Once injected into the interplanetary space, they can propagate towards Earth, causing space weather related phenomena. For their analysis, interplanetary in-situ measurements of charged particles are key. The recently expanded spacecraft fleet in the heliosphere not only provides much-needed additional vantage points, but also increases the variety of missions and instruments for which data loading and processing tools are needed. This manuscript introduces a series of Python functions that will enable the scientific community to download, load, and visualize charged particle measurements of the current space missions that are especially relevant to particle research as time series or dynamic spectra. In addition, further analytical functionality is provided that allows the determination of SEP onset times as well as their inferred injection times. The full workflow, which is intended to be run within Jupyter Notebooks and can also be approachable for Python laymen, will be presented with scientific examples. All functions are written in Python, with the source code publicly available at GitHub under a permissive license. Where appropriate, available Python libraries are used, and their application is described. △ Less

Submitted 14 December, 2022; v1 submitted 19 October, 2022; originally announced October 2022.

Journal ref: Front. Astron. Space Sci. 9:1073578 (2022)