Showing 1–50 of 86 results for author: Frost, J

Unveiling accretion in the massive YSO G033.3891. Spatial and kinematic constraints from the CO bandhead emission

Authors: E. Koumpia, D. Sun, M. Koutoulaki, J. D. Ilee, W. -J. de Wit, R. D. Oudmaijer, A. J. Frost

Abstract: The inner parts of the hot discs surrounding massive young stellar objects (MYSOs) are still barely explored due to observational limitations in terms of angular resolution, scarcity of diagnostic lines and the embedded and rare nature of these targets. We present the first K-band spectro-interferometric observations toward the MYSO G033.3891, which based on former kinematic evidence via the CO ba… ▽ More The inner parts of the hot discs surrounding massive young stellar objects (MYSOs) are still barely explored due to observational limitations in terms of angular resolution, scarcity of diagnostic lines and the embedded and rare nature of these targets. We present the first K-band spectro-interferometric observations toward the MYSO G033.3891, which based on former kinematic evidence via the CO bandhead emission is known to host an accreting disc. Using the high spectral resolution mode (R$\sim$4000) of the GRAVITY/VLTI, we spatially resolve the emission of the inner dusty disc and the crucial gaseous interface between the star and the dusty disc. Using detailed modelling on the K-band dust continuum and tracers known to be associated with the ionised and molecular gaseous interface (Br$γ$, CO), we report on the smallest scales of accretion/ejection. The new observations in combination with our geometric and kinematic models employed to fit former high spectral resolution observations on the source (R$\sim$30,000; CRIRES/VLTI) allow us to constrain the size of the inner gaseous disc both spatially and kinematically via the CO overtone emission at only 2 au. Our models reveal that both Br$γ$ and CO emissions are located well within the dust sublimation radius (5~au) as traced by the hot 2.2~$μ$m dust continuum. Our paper provides the first case study where the tiniest scales of gaseous accretion around the MYSO G033.3891 are probed both kinematically and spatially via the CO bandhead emission. This analysis of G033.3891 stands as only the second instance of such investigation within MYSOs, underscoring the gradual accumulation of knowledge regarding how massive young stars gain their mass, while further solidifying the disc nature of accretion at the smallest scales of MYSOs. △ Less

Submitted 27 August, 2024; originally announced August 2024.

Comments: Accepted for publication in A&A (9 pages, 5 figures, 2 tables)

Charge Transport and Defects in Sulfur-Deficient Chalcogenide Perovskite BaZrS$_3$

Authors: Garima Aggarwal, Adeem Saeed Mirza, Stefania Riva, Corrado Comparotto, Robert J. W. Frost, Soham Mukherjee, Monica Morales-Masis, Håkan Rensmo, Jonathan Staaf Scragg

Abstract: Exploring the conduction mechanism in the chalcogenide perovskite BaZrS$_3$ is of significant interest due to its potential suitability as a top absorber layer in silicon-based tandem solar cells and other optoelectronic applications. Theoretical and experimental studies anticipate native ambipolar doping in BaZrS$_3$, although experimental validation remains limited. This study reveals a transiti… ▽ More Exploring the conduction mechanism in the chalcogenide perovskite BaZrS$_3$ is of significant interest due to its potential suitability as a top absorber layer in silicon-based tandem solar cells and other optoelectronic applications. Theoretical and experimental studies anticipate native ambipolar doping in BaZrS$_3$, although experimental validation remains limited. This study reveals a transition from highly insulating behavior to n-type conductivity in BaZrS$_3$ through annealing in an S-poor environment. BaZrS$_3$ thin films are synthesized $\textit{via}$ a two step process: co-sputtering of Ba-Zr followed by sulfurization at 600 $^{\circ}$C, and subsequent annealing in high vacuum. UV-Vis measurement reveal a red-shift in the absorption edge concurrent with sample color darkening after annealing. The increase in defect density with vacuum annealing, coupled with low activation energy and n-type character of defects, strongly suggests that sulfur vacancies (V$_{\mathrm{S}}$) are responsible, in agreement with theoretical predictions. The shift of the Fermi level towards conduction band minimum, quantified by Hard X-ray Photoelectron Spectroscopy (Ga K$α$, 9.25 keV), further corroborates the induced n-type of conductivity in annealed samples. Our findings indicate that vacuum annealing induces V$_{\mathrm{S}}$ defects that dominate the charge transport, thereby making BaZrS$_3$ an n-type semiconductor under S-poor conditions. This study offers crucial insights into understanding the defect properties of BaZrS$_3$, facilitating further improvements for its use in solar cell applications. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 19 pages (single column), 6 figures in main manuscript and 8 figures in supplementary information

Anharmonic phonons with Gaussian processes

Authors: Keerati Keeratikarn, Jarvist Moore Frost

Abstract: We provide a method for calculating anharmonic lattice dynamics, by building a surrogate model based on Gaussian Processes (GPs). Due to the underlying Gaussian form of a GP, the model is infinitely differentiable. This allows us to train the model trained directly on forces (the derivative of PESs) reducing the evaluations required for a given accuracy. We can extend this differentiation to direc… ▽ More We provide a method for calculating anharmonic lattice dynamics, by building a surrogate model based on Gaussian Processes (GPs). Due to the underlying Gaussian form of a GP, the model is infinitely differentiable. This allows us to train the model trained directly on forces (the derivative of PESs) reducing the evaluations required for a given accuracy. We can extend this differentiation to directly calculate second and third order force-constants using automatic differentiation (AD). For the five model materials we study, we find that the force-constants are in close agreement with a standard finite-displacement approach. Our method appears to be linear scaling in the number of atoms at predicting both second and third-order (anharmonic) force-constants. △ Less

Submitted 6 June, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: 12 pages, 5 figures

Anisotropic electron-phonon interactions in 2D lead-halide perovskites

Authors: Jaco J. Geuchies, Johan Klarbring, Lucia Di Virgillio, Shuai Fu, Sheng Qu, Guangyu Liu, Hai Wang, Jarvist M. Frost, Aron Walsh, Mischa Bonn, Heejae Kim

Abstract: Two-dimensional hybrid organic-inorganic metal halide perovskites offer enhanced stability for perovskite-based applications. Their crystal structure's soft and ionic nature gives rise to strong interactions between charge carriers and ionic rearrangements. Here, we investigate the interaction of photo-generated electrons and ionic polarizations in single-crystal 2D perovskite butylammonium lead i… ▽ More Two-dimensional hybrid organic-inorganic metal halide perovskites offer enhanced stability for perovskite-based applications. Their crystal structure's soft and ionic nature gives rise to strong interactions between charge carriers and ionic rearrangements. Here, we investigate the interaction of photo-generated electrons and ionic polarizations in single-crystal 2D perovskite butylammonium lead iodide, varying the inorganic lammelae thickness in the 2D single crystals. We determined the directionality of the transition dipole moments of the relevant phonon modes (in the 0.3-3 THz range) by angle-and-polarization dependent THz transmission measurements. We find a clear anisotropy of the in-plane photoconductivity, with a 10% reduction along the axis parallel with the transition dipole moment of the most strongly coupled phonon. Detailed calculations, based on Feynman polaron theory, indicate that the anisotropy originates from directional electron-phonon interactions. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 50 pages, 5 figures main text, 13 figures in supplementary information

A magnetic massive star has experienced a stellar merger

Authors: A. J. Frost, H. Sana, L. Mahy, G. Wade, J. Barron, J. -B. Le Bouquin, A. Mérand, F. R. N. Schneider, T. Shenar, R. H. Barbá, D. M. Bowman, M. Fabry, A. Farhang, P. Marchant, N. I. Morrell, J. V. Smoker

Abstract: Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary sys… ▽ More Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Full paper with supplementary materials. 59 pages, 18 figures

Journal ref: Science, 11 Apr 2024, Vol 384, Issue 6692, pp. 214-217

Optimization of Ray-tracing Simulations to Confirm Performance of the GP-SANS Instrument at the High-Flux Isotope Reactor

Authors: James M. Rogers, Matthew J. Frost, Lisa DeBeer-Schmitt

Abstract: The CG-2 beamline at the High Flux Isotope Reactor (HFIR) exhibits a notable discrepancy between observed count rates and the count rates we would expect based on a Monte-Carlo neutron ray-trace simulation. These simulations consistently predict count rates approximately five times greater than those observed in four separate experimental runs involving different instrument configurations. This di… ▽ More The CG-2 beamline at the High Flux Isotope Reactor (HFIR) exhibits a notable discrepancy between observed count rates and the count rates we would expect based on a Monte-Carlo neutron ray-trace simulation. These simulations consistently predict count rates approximately five times greater than those observed in four separate experimental runs involving different instrument configurations. This discrepancy suggests that certain factors are causing losses in measurements that are not adequately accounted for in the simulation, in particular guide reflectivity or misalignment. To investigate these discrepancies, a high-dimensional simulation parameter approach is applied in order to understand the losses. Region of Interest (ROI) groups along the instrument are assigned to different surfaces of the guide components within the simulation. This allows the parameters of those guide components to be varied as a group to minimize the complexity of the search space. The result is an optimization of simulation parameters using an iterative scheme that aims to minimize the difference between experimentally measured count rates and simulated count rates across all tested collimator combinations. This proposed methodology holds the potential to reveal previously unrecognized sources of intensity loss in the CG-2 beamline at HFIR and improve the accuracy of simulations, leading to enhanced understanding and performance of the beamline for various scientific applications. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 14 pages, 19 figures

Analysis Facilities White Paper

Authors: D. Ciangottini, A. Forti, L. Heinrich, N. Skidmore, C. Alpigiani, M. Aly, D. Benjamin, B. Bockelman, L. Bryant, J. Catmore, M. D'Alfonso, A. Delgado Peris, C. Doglioni, G. Duckeck, P. Elmer, J. Eschle, M. Feickert, J. Frost, R. Gardner, V. Garonne, M. Giffels, J. Gooding, E. Gramstad, L. Gray, B. Hegner , et al. (41 additional authors not shown)

Abstract: This white paper presents the current status of the R&D for Analysis Facilities (AFs) and attempts to summarize the views on the future direction of these facilities. These views have been collected through the High Energy Physics (HEP) Software Foundation's (HSF) Analysis Facilities forum, established in March 2022, the Analysis Ecosystems II workshop, that took place in May 2022, and the WLCG/HS… ▽ More This white paper presents the current status of the R&D for Analysis Facilities (AFs) and attempts to summarize the views on the future direction of these facilities. These views have been collected through the High Energy Physics (HEP) Software Foundation's (HSF) Analysis Facilities forum, established in March 2022, the Analysis Ecosystems II workshop, that took place in May 2022, and the WLCG/HSF pre-CHEP workshop, that took place in May 2023. The paper attempts to cover all the aspects of an analysis facility. △ Less

Submitted 15 April, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Improved Limits on $n \rightarrow n'$ Transformation from the Spallation Neutron Source

Authors: Francisco M. Gonzalez, Cary Rock, Leah J. Broussard, Lisa DeBeer-Schmitt, Matthew J. Frost, Lawrence Heilbronn, Erik B. Iverson, Yuri Kamyshkov, Michael Kline, David Milstead, Devyn Powers, James Rogers, Valentina Santoro, Shaun Vavra

Abstract: Conversions between neutrons $n$ and Dark Matter candidate sterile neutrons $n'$ have been proposed as a mechanism for Baryon Number $\mathcal{B}$ violation. In the case that there is a small mass difference $Δ{m}$ between the $n$ and the $n'$ states, oscillations can be induced by compensating for $Δ{m}$ with a magnetic field. A search for such neutron oscillations was performed at the Spallation… ▽ More Conversions between neutrons $n$ and Dark Matter candidate sterile neutrons $n'$ have been proposed as a mechanism for Baryon Number $\mathcal{B}$ violation. In the case that there is a small mass difference $Δ{m}$ between the $n$ and the $n'$ states, oscillations can be induced by compensating for $Δ{m}$ with a magnetic field. A search for such neutron oscillations was performed at the Spallation Neutron Source by looking for anomalous neutron transmission through a strongly absorbing cadmium wafer inside of a $6.6$~T magnet. The approach described here saw no regenerated neutrons above background, which provides an improved limit for neutron - sterile neutrons transformations for a range of $Δ{m}$ between $0.1$~neV and $1000$~neV. △ Less

Submitted 24 February, 2024; originally announced February 2024.

Comments: 18 pages, 15 figures

The onset of stellar multiplicity in massive star formation: A search for low-mass companions of massive young stellar objects with $L'$-band adaptive optics imaging

Authors: E. Bordier, W. -J. de Wit, A. J. Frost, H. Sana, T. Pauwels, E. Koumpia

Abstract: Given the high incidence of binaries among mature field massive stars, it is clear that multiplicity is an inevitable outcome of high-mass star formation. Understanding how massive multiples form requires the study of the birth environments of massive stars, covering the innermost to outermost regions. We aim to detect and characterise low-mass companions around massive young stellar objects (MYSO… ▽ More Given the high incidence of binaries among mature field massive stars, it is clear that multiplicity is an inevitable outcome of high-mass star formation. Understanding how massive multiples form requires the study of the birth environments of massive stars, covering the innermost to outermost regions. We aim to detect and characterise low-mass companions around massive young stellar objects (MYSOs) during and shortly after their formation phase. To investigate large spatial scales, we carried out an $L'$-band high-contrast direct imaging survey seeking low-mass companions (down to $L_{\text{bol}}\approx 10 L_{\odot}$, or late A-type) around thirteen previously identified MYSOs using the VLT/NACO instrument. From those images, we looked for the presence of companions on a wide orbit, covering scales from 300 to 56,000 au. Detection limits were determined for all targets and we tested the gravitational binding to the central object based on chance projection probabilities. We have discovered a total of thirty-nine potential companions around eight MYSOs, the large majority of which have never been reported to date. We derived a multiplicity frequency (MF) of $62\pm13$% and a companion fraction (CF) of $3.0\pm0.5$. The derived MF and CF are compared to other studies for similar separation ranges. The comparisons are effective for a fixed evolutionary stage spanning a wide range of masses and vice versa. We find an increased MF and CF compared to the previous studies targeting MYSOs, showing that the statement in which multiplicity scales with primary mass also extends to younger evolutionary stages. The separations at which the companions are found and their location with relation to the primary star allow us to discuss the implications for the massive star formation theories. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: 24 pages, 2 Tables, 5 Figures and 4 appendices (containing figures)

Journal ref: A&A 681, A85 (2024)

Light-yield response of liquid scintillators using 2--6 MeV tagged neutrons

Authors: N. Mauritzson, K. G. Fissum, J. R. M. Annand, H. Perrey, R. Al Jebali, A. Backis, R. Hall-Wilton, K. Kanaki, V. Maulerova-Subert, F. Messi, R. J. W. Frost, E. Rofors, J. Scherzinger

Abstract: Knowledge of the neutron light-yield response is crucial to the understanding of scintillator-based neutron detectors. In this work, neutrons from 2--6 MeV have been used to study the scintillation light-yield response of the liquid scintillators NE 213A, EJ 305, EJ 331 and EJ 321P using event-by-event waveform digitization. Energy calibration was performed using a GEANT model to locate the edge p… ▽ More Knowledge of the neutron light-yield response is crucial to the understanding of scintillator-based neutron detectors. In this work, neutrons from 2--6 MeV have been used to study the scintillation light-yield response of the liquid scintillators NE 213A, EJ 305, EJ 331 and EJ 321P using event-by-event waveform digitization. Energy calibration was performed using a GEANT model to locate the edge positions of the Compton distributions produced by gamma-ray sources. The simulated light yield for neutrons from a PuBe source was compared to measured recoil proton distributions, where neutron energy was selected by time-of-flight. This resulted in an energy-dependent Birks parametrization to characterize the non-linear response to the lower energy neutrons. The NE 213A and EJ 305 results agree very well with existing data and are reproduced nicely by the simulation. New results for EJ 331 and EJ 321P, where the simulation also reproduces the data well, are presented. △ Less

Submitted 24 November, 2023; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: 31 pages, 15 figures, 4 tables

Roadmap on Photovoltaic Absorber Materials for Sustainable Energy Conversion

Authors: James C. Blakesley, Ruy S. Bonilla, Marina Freitag, Alex M. Ganose, Nicola Gasparini, Pascal Kaienburg, George Koutsourakis, Jonathan D. Major, Jenny Nelson, Nakita K. Noel, Bart Roose, Jae Sung Yun, Simon Aliwell, Pietro P. Altermatt, Tayebeh Ameri, Virgil Andrei, Ardalan Armin, Diego Bagnis, Jenny Baker, Hamish Beath, Mathieu Bellanger, Philippe Berrouard, Jochen Blumberger, Stuart A. Boden, Hugo Bronstein , et al. (61 additional authors not shown)

Abstract: Photovoltaics (PVs) are a critical technology for curbing growing levels of anthropogenic greenhouse gas emissions, and meeting increases in future demand for low-carbon electricity. In order to fulfil ambitions for net-zero carbon dioxide equivalent (CO₂eq) emissions worldwide, the global cumulative capacity of solar PVs must increase by an order of magnitude from 0.9 TWp in 2021 to 8.… ▽ More Photovoltaics (PVs) are a critical technology for curbing growing levels of anthropogenic greenhouse gas emissions, and meeting increases in future demand for low-carbon electricity. In order to fulfil ambitions for net-zero carbon dioxide equivalent (CO₂eq) emissions worldwide, the global cumulative capacity of solar PVs must increase by an order of magnitude from 0.9 TWp in 2021 to 8.5 TWp by 2050 according to the International Renewable Energy Agency, which is considered to be a highly conservative estimate. In 2020, the Henry Royce Institute brought together the UK PV community to discuss the critical technological and infrastructure challenges that need to be overcome to address the vast challenges in accelerating PV deployment. Herein, we examine the key developments in the global community, especially the progress made in the field since this earlier roadmap, bringing together experts primarily from the UK across the breadth of the photovoltaics community. The focus is both on the challenges in improving the efficiency, stability and levelized cost of electricity of current technologies for utility-scale PVs, as well as the fundamental questions in novel technologies that can have a significant impact on emerging markets, such as indoor PVs, space PVs, and agrivoltaics. We discuss challenges in advanced metrology and computational tools, as well as the growing synergies between PVs and solar fuels, and offer a perspective on the environmental sustainability of the PV industry. Through this roadmap, we emphasize promising pathways forward in both the short- and long-term, and for communities working on technologies across a range of maturity levels to learn from each other. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 160 pages, 21 figures

Gaia uncovers difference in B and Be star binarity at small scales: evidence for mass transfer causing the Be phenomenon

Authors: Jonathan M. Dodd, René D. Oudmaijer, Isaac C. Radley, Miguel Vioque, Abigail J. Frost

Abstract: Be stars make up almost 20% of the B star population, and are rapidly rotating stars surrounded by a disc; however the origin of this rotation remains unclear. Mass transfer within close binaries provides the leading hypothesis, with previous detections of stripped companions to Be stars supporting this. Here, we exploit the exquisite astrometric precision of Gaia to carry out the largest to date… ▽ More Be stars make up almost 20% of the B star population, and are rapidly rotating stars surrounded by a disc; however the origin of this rotation remains unclear. Mass transfer within close binaries provides the leading hypothesis, with previous detections of stripped companions to Be stars supporting this. Here, we exploit the exquisite astrometric precision of Gaia to carry out the largest to date comparative study into the binarity of matched samples of nearby B and Be stars from the Bright Star Catalogue. By utilising new "proper motion anomaly" values, derived from Gaia DR2 and DR3 astrometric data alongside previous values calculated using Hipparcos and Gaia data, and the Gaia provided RUWE, we demonstrate that we can identify unresolved binaries down to separations of 0.02". Using these measures, we find that the binary fractions of B and Be stars are similar between 0.04 - 10", but the Be binary fraction is significantly lower than that of the B stars for separations below 0.04". As the separation range of these "missing" binaries is too large for mass transfer, and stripped companions are not retrieved by these measures, we suggest the companions migrate inwards via binary hardening within a triple system. This confirms statistically for the first time the hypothesis that binary interaction causes the Be phenomenon, with migration causing the dearth of Be binaries between 0.02 - 0.04". Furthermore, we suggest that triplicity plays a vital role in this migration, and thus in the formation of Be stars as a whole. △ Less

Submitted 10 October, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 11 pages, 4 figures; typos corrected

MWC 656 is unlikely to contain a black hole

Authors: S. Janssens, T. Shenar, N. Degenaar, J. Bodensteiner, H. Sana, J. Audenaert, A. J. Frost

Abstract: Context. MWC 656 was reported as the first known Be star with a black-hole (BH) companion in a 60 d period. The mass of the proposed BH companion is estimated to be between 4 - 7 MSun. This estimate is based on radial velocity (RV) measurements derived from the Fe ii 4583 emission line of the Be star disc and from the He ii 4686 emission line, assumed to be formed in a disc around the putative BH.… ▽ More Context. MWC 656 was reported as the first known Be star with a black-hole (BH) companion in a 60 d period. The mass of the proposed BH companion is estimated to be between 4 - 7 MSun. This estimate is based on radial velocity (RV) measurements derived from the Fe ii 4583 emission line of the Be star disc and from the He ii 4686 emission line, assumed to be formed in a disc around the putative BH. Aims. Using new high-resolution spectroscopic data, we investigate whether MWC 656 truly contains a BH. Methods. We used the cross-correlation method to calculate the RVs of both the Be star and the He ii 4686 emission line and we derive a new orbital solution. We also performed disentangling to look for the spectral signature of a companion. Results. We derive an orbital period of 59.028 +- 0.011 d and a mass ratio q = M_Heii/M_Be = 0.12 +- 0.03, much lower than the previously reported q = 0.41 +- 0.07. Adopting a mass of the Be star of M_Be = 7.8 +- 2.0MSun, the companion has a mass of 0.94 +- 0.34MSun. For the upper limit of M_Be = 16MSun and q = 0.15, the companion has a mass 2.4MSun. Performing disentangling on mock spectra shows that the spectral signature of a non-degenerate stellar companion with such a low mass cannot be retrieved using our data. Conclusions. Our measurements do not support the presence of a BH companion in MWC 656. The derived upper limit on the mass of the companion rather indicates that it is a neutron star, a white dwarf, or a hot helium star. Far-UV data will help to reject or confirm a hot helium-star companion. △ Less

Submitted 21 August, 2023; v1 submitted 16 August, 2023; originally announced August 2023.

Comments: 6 pages, 2 + 6 figures, 3 + 1 tables, accepted for publication in A&A

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

Surrogate model solver for impurity-induced superconducting subgap states

Authors: Virgil V. Baran, Emil J. P. Frost, Jens Paaske

Abstract: A simple impurity solver is shown to capture the impurity-induced superconducting subgap states in quantitative agreement with the numerical renormalization group and quantum Monte-Carlo simulations. The solver is based on the exact diagonalization of a single-impurity Anderson model with discretized superconducting reservoirs including only a small number of effective levels. Their energies and c… ▽ More A simple impurity solver is shown to capture the impurity-induced superconducting subgap states in quantitative agreement with the numerical renormalization group and quantum Monte-Carlo simulations. The solver is based on the exact diagonalization of a single-impurity Anderson model with discretized superconducting reservoirs including only a small number of effective levels. Their energies and couplings to the impurity $d$-level are chosen so as to best reproduce the Matsubara frequency dependence of the hybridization function. We provide a number of critical benchmarks and demonstrate the solvers efficiency in combination with the reduced basis method [Phys. Rev. B 107, 144503 (2023)] by calculating the phase diagram for an interacting three-terminal junction. △ Less

Submitted 13 December, 2023; v1 submitted 21 July, 2023; originally announced July 2023.

Comments: 19 pages, 20 figures; v2: new results on the phase diagram of a finite-gap interacting topological three-terminal Josephson junction

Report number: NBI QDEV CMT 2023

Journal ref: Phys. Rev. B 108, L220506 (2023)

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Multiplicity of northern bright O-type stars with optical long baseline interferometry

Authors: Cyprien Lanthermann, Jean-Baptiste Le Bouquin, Hugues Sana, Antoine Mérand, John D. Monnier, Karine Perraut, Abigail J. Frost, Laurent Mahy, Eric Gosset, Michael De Becker, Stefan Kraus, Narsireddy Anugu, Claire L. Davies, Jacob Ennis, Tyler Gardner, Aaron Labdon, Benjamin Setterholm, Theo ten Brummelaar, Gail H. Schaefer

Abstract: The study of the multiplicity of massive stars gives hints on their formation processes and their evolutionary paths, which are still not fully understood. Large separation binaries (>50 milliseconds of arc, mas) can be probed by adaptive-optics-assisted direct imaging and sparse aperture masking, while close binaries can be resolved by photometry and spectroscopy. However, optical long baseline i… ▽ More The study of the multiplicity of massive stars gives hints on their formation processes and their evolutionary paths, which are still not fully understood. Large separation binaries (>50 milliseconds of arc, mas) can be probed by adaptive-optics-assisted direct imaging and sparse aperture masking, while close binaries can be resolved by photometry and spectroscopy. However, optical long baseline interferometry is mandatory to establish the multiplicity of Galactic massive stars at the separation gap between 1 and 50 mas. In this paper, we aim to demonstrate the capability of the new interferometric instrument MIRC-X, located at the CHARA Array, to study the multiplicity of O-type stars and therefore probe the full range of separation for more than 120 massive stars (H<7.5 mag). We initiated a pilot survey of bright O-type stars (H<6.5mag) observable with MIRC-X. We observed 29 O-type stars, including two systems in average atmospheric conditions around a magnitude of H=7.5 mag. We systematically reduced the obtained data with the public reduction pipeline of the instrument. We analyzed the reduced data using the dedicated python software CANDID to detect companions. Out of these 29 systems, we resolved 19 companions in 17 different systems with angular separations between ~0.5 and 50 mas. This results in a multiplicity fraction fm=17/29=0.59+/-0.09, and an average number of companions fc=19/29=0.66+/-0.13. Those results are in agreement with the results of the SMASH+ survey in the Southern Hemisphere. Thirteen of these companions have been resolved for the first time, including the companion responsible for the nonthermal emission in Cyg OB2-5 A and the confirmation of the candidate companion of HD 47129 suggested by SMASH+. A large survey on more than 120 northern O-type stars (H<7.5) is possible with MIRC-X and will be fruitful. △ Less

Submitted 8 February, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 15 pages, 9 figures, 5 tables, accepted in A&A

Journal ref: A&A 672, A6 (2023)

First spatially resolved Na I and He I transitions towards an MYSO. Finding new tracers for the gaseous star/disc interface

Authors: Evgenia Koumpia, M. Koutoulaki, W. -J. de Wit, R. D. Oudmaijer, A. J. Frost, S. L. Lumsden, J. M. Pittard

Abstract: With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We pre… ▽ More With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We present the analysis of K-band spectro-interferometric observations toward the Massive Young Stellar Object IRAS 13481-6124, which is known to host an accreting dusty disc. Using GRAVITY on the VLTI, we trace the crucial au-scales of the warm inner interface between the star and the accretion dusty disc. We detect and spatially resolve the Na I doublet and He I transitions towards an object of this class for the first time. The new observations in combination with our geometric models allowed us to probe the smallest au-scales of accretion/ejection around an MYSO. We find that Na I originates in the disc at smaller radii than the dust disc and is more compact than any of the other spatially resolved diagnostics (Br$γ$, He I, and CO). Our findings suggest that Na I can be a new powerful diagnostic line in tracing the warm star/disc accreting interface of forming (massive) stars, while the similarities between He I and Br$γ$ point towards an accretion/ejection origin of He I △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS Letters; 5 pages, 3 Figures

Metal-silicate partitioning of W and Mo and the role of carbon in controlling their abundances in the Bulk Silicate Earth

Authors: E. S. Jennings, S. A. Jacobson, D. C. Rubie, Y. Nakajima, A. K. Vogel, L. A. Rose-Weston, D. J. Frost

Abstract: The liquid metal-liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf-W isotopic system is used to date early planetary evolution. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous exper… ▽ More The liquid metal-liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf-W isotopic system is used to date early planetary evolution. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous experiments to re-examine the systematics of Mo and W partitioning. W partitioning is particularly sensitive to silicate and metallic melt compositions and becomes more siderophile with increasing temperature. We show that W has a 6+ oxidation state in silicate melts over the full experimental fO2 range of $Δ$IW -1.5 to -3.5. Mo has a 4+ oxidation state and its partitioning is less sensitive to silicate melt composition, but also depends on metallic melt composition. DMo stays approximately constant with increasing depth in Earth. Both W and Mo become more siderophile with increasing C content of the metal, so we fit epsilon interaction parameters. W and Mo along with C are incorporated into a combined N-body accretion and core-mantle differentiation model. We show that W and Mo require the early accreting Earth to be sulfur-depleted and carbon-enriched so that W and Mo are efficiently partitioned into Earth's core and do not accumulate in the mantle. If this is the case, the produced Earth-like planets possess mantle compositions matching the BSE for all simulated elements. However, there are two distinct groups of estimates of the bulk mantle's C abundance in the literature: low (100 ppm), and high (800 ppm), and all models are consistent with the higher estimated carbon abundance. The low BSE C abundance would be achievable when the effects of the segregation of dispersed metal droplets produced in deep magma oceans by the disproportionation of Fe2+ to Fe3+ plus metallic Fe is considered. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: 61 pages including supplemental material

Journal ref: Geochimica et Cosmochimica Acta. 293, 40-69 (2021)

The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries

Authors: Hugues Sana, Michael Abdul-Masih, Gareth Banyard, Julia Bodensteiner, Dominic M. Bowman, Karan Dsilva, C. Eldridge, Matthias Fabry, Abigail J. Frost, Calum Hawcroft, Soetkin Janssens, Laurent Mahy, Pablo Marchant, Norbert Langer, Timothy Van Reeth, Koushik Sen, Tomer Shenar

Abstract: Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open quest… ▽ More Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: To appear in the proceedings of IAUS361: Massive stars near and far; 6 pages, 1 figure

Predicting polaron mobility in organic semiconductors with the Feynman variational approach

Authors: Bradley A. A. Martin, Jarvist Moore Frost

Abstract: We extend the Feynman variational method applied to the parabolic-band Fröhlich (continuum) large polaron~\cite{Feynman1955} to a Holstein (lattice) small polaron, with a parabolic-band. This new theory shows a discrete localisation as a function of coupling strength. Having build the theory with the same quasi-particle Lagrangian as the 1955 work, we can directly use the FHIP~\cite{Feynman1962} r… ▽ More We extend the Feynman variational method applied to the parabolic-band Fröhlich (continuum) large polaron~\cite{Feynman1955} to a Holstein (lattice) small polaron, with a parabolic-band. This new theory shows a discrete localisation as a function of coupling strength. Having build the theory with the same quasi-particle Lagrangian as the 1955 work, we can directly use the FHIP~\cite{Feynman1962} response theory to calculate DC mobility and complex conductivity. We show that we can take matrix elements from electronic structure calculations on real materials, by modelling charge-carrier mobility in crystalline Rubrene. Good agreement is found to measurement, with a predicted mobility of $μ= 47.72$~\si{cm^2 V^{-1} s^{-1}} at $300$~\si{K}. △ Less

Submitted 12 June, 2024; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, 9 figures, 72 equations

Technique for the measurement of intrinsic pulse-shape discrimination for organic scintillators using tagged neutrons

Authors: N. Mauritzson, K. G. Fissum, J. R. M. Annand, H. Perrey, R. J. W. Frost, R. Al Jebali, A. Backis, R. Hall-Wilton, K. Kanaki, V. Maulerova-Subert, C. Maurer, F. Messi, E. Rofors

Abstract: Fast-neutron/gamma-ray pulse-shape discrimination has been performed for the organic liquid scintillators NE 213A and EJ 305 using a time-of-flight based neutron-tagging technique and waveform digitization on an event-by-event basis. Gamma-ray sources and a Geant4-based simulation were used to calibrate the scintillation-light yield. The difference in pulse shape for the neutron and gamma-ray even… ▽ More Fast-neutron/gamma-ray pulse-shape discrimination has been performed for the organic liquid scintillators NE 213A and EJ 305 using a time-of-flight based neutron-tagging technique and waveform digitization on an event-by-event basis. Gamma-ray sources and a Geant4-based simulation were used to calibrate the scintillation-light yield. The difference in pulse shape for the neutron and gamma-ray events was analysed by integrating selected portions of the digitized waveform to produce a figure-of-merit for neutron/gamma-ray separation. This figure-of-merit has been mapped as a function of detector threshold and also of neutron energy determined from time-of-flight. It shows clearly that the well-established pulse-shape discrimination capabilities of NE 213A are superior to those of EJ 305. The extra information provided by the neutron-tagging technique has resulted in a far more detailed assessment of the pulse-shape discrimination capabilities of these organic scintillators. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: 22 pages, 12 figures, 2 tables

Round Trees and Conformal Dimension in Random Groups: low density to high density

Authors: Jordan Frost

Abstract: We investigate conformal dimension for the class of infinite hyperbolic groups in the Gromov density model $\mathcal{G}^d_{m,l}$ of random groups with $m \geq 2$ fixed generators, density $0 < d < 1/2$ and relator length $l \to \infty$. Our main result is a lower bound linear in $l$ at all densities $0 < d < 1/2$ achieved by building undistorted round trees coming directly from lower density Gromo… ▽ More We investigate conformal dimension for the class of infinite hyperbolic groups in the Gromov density model $\mathcal{G}^d_{m,l}$ of random groups with $m \geq 2$ fixed generators, density $0 < d < 1/2$ and relator length $l \to \infty$. Our main result is a lower bound linear in $l$ at all densities $0 < d < 1/2$ achieved by building undistorted round trees coming directly from lower density Gromov random groups. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 31 pages, 9 figures, comments welcome!

MSC Class: 20F67 (Primary) 20F05; 20F06; 20P05; 30L10 (Secondary)

Multiple phonon modes in Feynman path-integral variational polaron mobility

Authors: Bradley A. A. Martin, Jarvist Moore Frost

Abstract: The Feynman path-integral variational approach to the polaron problem\cite{Feynman1955}, along with the associated FHIP linear-response mobility theory\cite{Feynman1962}, provides a computationally amenable method to predict the frequency-resolved temperature-dependent charge-carrier mobility, and other experimental observables in polar semiconductors. We show that the FHIP mobility theory predict… ▽ More The Feynman path-integral variational approach to the polaron problem\cite{Feynman1955}, along with the associated FHIP linear-response mobility theory\cite{Feynman1962}, provides a computationally amenable method to predict the frequency-resolved temperature-dependent charge-carrier mobility, and other experimental observables in polar semiconductors. We show that the FHIP mobility theory predicts non-Drude transport behaviour, and shows remarkably good agreement with the recent diagrammatic Monte-Carlo mobility simulations of Mishchenko et al.\cite{Mishchenko2019} for the abstract Fröhlich Hamiltonian. We extend this method to multiple phonon modes in the Fröhlich model action. This enables a slightly better variational solution, as inferred from the resulting energy. We carry forward this extra complexity into the mobility theory, where it shows richer structure in the frequency and temperature dependent mobility, due to the different phonon modes activating at different energies. The method provides a computationally efficient and fully quantitative method of predicting polaron mobility and response in real materials. △ Less

Submitted 22 December, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 28 pages, 17 figures

On the origin of close massive binaries in the M17 star-forming region

Authors: E. Bordier, A. J. Frost, H. Sana, M. Reggiani, A. Mérand, A. Rainot, M. C. Ramírez-Tannus, W. J. de Wit

Abstract: Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be eluci… ▽ More Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 20 pages, 14 figures, 4 tables, 3 appendices

Journal ref: A&A 663, A26 (2022)

HR 6819 is a binary system with no black hole -- revisiting the source with infrared interferometry and optical integral field spectroscopy

Authors: A. J. Frost, J. Bodensteiner, Th. Rivinius, D. Baade, A. Merand, F. Selman, M. Abdul-Masih, G. Banyard, E. Bordier, K. Dsilva, C. Hawcroft, L. Mahy, M. Reggiani, T. Shenar, M. Cabezas, P. Hadrava, M. Heida, R. Klement, H. Sana

Abstract: Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stella… ▽ More Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 16 pages, 11 figures, subject of ESO press release eso2204

Journal ref: A&A, Volume 659, March 2022, L3, 12, Letters to the Editor

Planet Hunters TESS IV: A massive, compact hierarchical triple star system TIC 470710327

Authors: Nora L. Eisner, Cole Johnston, Silvia Toonen, Abigail J. Frost, Soetkin Janssens, Chris J. Lintott, Suzanne Aigrain, Hugues Sana, Michael Abdul-Masih, Karla Z. Arellano-Córdova, Paul G. Beck, Emma Bordier, Emily Canon, Ana Escorza, Mattias Fabry, Lars Hermansson, Steve Howell, Grant Miller, Shreeya Sheyte, Safaa Alhassan, Elisabeth M. L. Baeten, Frank Barnet, Stewart. J. Bean, Mikael Bernau, David M. Bundy , et al. (15 additional authors not shown)

Abstract: We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, w… ▽ More We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS (11 pages, 8 figures)

Reconfigurable Compute-In-Memory on Field-Programmable Ferroelectric Diodes

Authors: Xiwen Liu, John Ting, Yunfei He, Merrilyn Mercy Adzo Fiagbenu, Jeffrey Zheng, Dixiong Wang, Jonathan Frost, Pariasadat Musavigharavi, Giovanni Esteves, Kim Kisslinger, Surendra B. Anantharaman, Eric A. Stach, Roy H. Olsson III, Deep Jariwala

Abstract: The deluge of sensors and data generating devices has driven a paradigm shift in modern computing from arithmetic-logic centric to data centric processing. At a hardware level, this presents an urgent need to integrate dense, high-performance and low-power memory units with Si logic-processor units. However, data-heavy problems such as search and pattern matching also require paradigm changing inn… ▽ More The deluge of sensors and data generating devices has driven a paradigm shift in modern computing from arithmetic-logic centric to data centric processing. At a hardware level, this presents an urgent need to integrate dense, high-performance and low-power memory units with Si logic-processor units. However, data-heavy problems such as search and pattern matching also require paradigm changing innovations at the circuit and architecture level to enable compute in memory (CIM) operations. CIM architectures that combine data storage yet concurrently offer low-delay and small footprint are highly sought after but have not been realized. Here, we present Aluminum Scandium Nitride (AlScN) ferroelectric diode (FeD) memristor devices that allow for storage, search and neural network-based pattern recognition in a transistor-free architecture. Our devices can be directly integrated on top of Si processors in a scalable, back-end-of-line process. We leverage the field-programmability, non-volatility and non-linearity of FeDs to demonstrated circuit blocks that can support search operations in-situ memory with search delay times < 0.1 ns and a cell footprint < 0.12 um^2. In addition, we demonstrate matrix multiplication operations with 4-bit operation of the FeDs. Our results highlight FeDs as promising candidates for fast, efficient, and multifunctional CIM platforms. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Impact of Side Chain Hydrophilicity on Packing, Swelling and Ion Interactions in Oxy-bithiophene Semiconductors

Authors: Nicholas Siemons, Drew Pearce, Camila Cendra, Hang Yu, Sachetan M. Tuladhar, Rawad K. Hallani, Rajendar Sheelamanthula, Garrett S. LeCroy, Lucas Siemons, Andrew J. P. White, Iain Mcculloch, Alberto Salleo, Jarvist M. Frost, Alexander Giovannitti, Jenny Nelson

Abstract: Exchanging hydrophobic alkyl-based side chains to hydrophilic glycol-based side chains is a widely adopted method for improving mixed-transport device performance, despite the impact on solid state packing and polymer-electrolyte interactions being poorly understood. Presented here is a Molecular Dynamics (MD) force field for modelling alkoxylated and glycolated polythiophenes. The force field is… ▽ More Exchanging hydrophobic alkyl-based side chains to hydrophilic glycol-based side chains is a widely adopted method for improving mixed-transport device performance, despite the impact on solid state packing and polymer-electrolyte interactions being poorly understood. Presented here is a Molecular Dynamics (MD) force field for modelling alkoxylated and glycolated polythiophenes. The force field is validated against known packing motifs for their monomer crystals. MD simulations, coupled with X-ray Diffraction (XRD), show that alkoxylated polythiophenes will pack with a `tilted stack' and straight interdigitating side chains, whilst their glycolated counterpart will pack with a `deflected stack' and an s-bend side chain configuration. MD simulations reveal water penetration pathways into the alkoxylated and glycolated crystals - through the π-stack and through the lamellar stack respectively. Finally, the two distinct ways tri-ethylene glycol polymers can bind to cations are revealed, showing the formation of a meta-stable single bound state, or an energetically deep double bound state, both with a strong side chain length dependance. The minimum energy pathways for the formation of the chelates are identified, showing the physical process through which cations can bind to one or two side chains of a glycolated polythiophene, with consequences for ion transport in bithiophene semiconductors. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 10 pages, 4 figures

Journal ref: Advanced Materials 34 (2022)

Explaining Reinforcement Learning Policies through Counterfactual Trajectories

Authors: Julius Frost, Olivia Watkins, Eric Weiner, Pieter Abbeel, Trevor Darrell, Bryan Plummer, Kate Saenko

Abstract: In order for humans to confidently decide where to employ RL agents for real-world tasks, a human developer must validate that the agent will perform well at test-time. Some policy interpretability methods facilitate this by capturing the policy's decision making in a set of agent rollouts. However, even the most informative trajectories of training time behavior may give little insight into the a… ▽ More In order for humans to confidently decide where to employ RL agents for real-world tasks, a human developer must validate that the agent will perform well at test-time. Some policy interpretability methods facilitate this by capturing the policy's decision making in a set of agent rollouts. However, even the most informative trajectories of training time behavior may give little insight into the agent's behavior out of distribution. In contrast, our method conveys how the agent performs under distribution shifts by showing the agent's behavior across a wider trajectory distribution. We generate these trajectories by guiding the agent to more diverse unseen states and showing the agent's behavior there. In a user study, we demonstrate that our method enables users to score better than baseline methods on one of two agent validation tasks. △ Less

Submitted 18 March, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: Accepted at ICML HILL 2021 Workshop

ACM Class: I.2.6

Experimental Search for Neutron to Mirror Neutron Oscillations as an Explanation of the Neutron Lifetime Anomaly

Authors: L. J. Broussard, J. L. Barrow, L. DeBeer-Schmitt, T. Dennis, M. R. Fitzsimmons, M. J. Frost, C. E. Gilbert, F. M. Gonzalez, L. Heilbronn, E. B. Iverson, A. Johnston, Y. Kamyshkov, M. Kline, P. Lewiz, C. Matteson, J. Ternullo, L. Varriano, S. Vavra

Abstract: An unexplained $>4\,σ$ discrepancy persists between "beam" and "bottle" measurements of the neutron lifetime. A new model proposed that conversions of neutrons $n$ into mirror neutrons $n'$, part of a dark mirror sector, can increase the apparent neutron lifetime by $1\%$ via a small mass splitting $Δ{m}$ between $n$ and $n'$ inside the 4.6 T magnetic field of the National Institute of Standards a… ▽ More An unexplained $>4\,σ$ discrepancy persists between "beam" and "bottle" measurements of the neutron lifetime. A new model proposed that conversions of neutrons $n$ into mirror neutrons $n'$, part of a dark mirror sector, can increase the apparent neutron lifetime by $1\%$ via a small mass splitting $Δ{m}$ between $n$ and $n'$ inside the 4.6 T magnetic field of the National Institute of Standards and Technology Beam Lifetime experiment. A search for neutron conversions in a 6.6 T magnetic field was performed at the Spallation Neutron Source which excludes this explanation for the neutron lifetime discrepancy. △ Less

Submitted 26 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Fröhlich polaron effective mass and localization length in cubic materials: degenerate and anisotropic electronic bands

Authors: Bogdan Guster, Pedro Melo, Bradley A. A. Martin, Véronique Brousseau-Couture, Anna Miglio, Matteo Giantomassi, Michel Côté, Jarvist M. Frost, Matthieu J. Verstraete, Xavier Gonze

Abstract: Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fröhlich model: a simple (non-degenerate) parabolic isotropic electronic band… ▽ More Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fröhlich model: a simple (non-degenerate) parabolic isotropic electronic band coupled to one dispersionless longitudinal optical phonon branch. The Fröhlich model allows one to understand characteristics such as polaron formation energy, radius, effective mass and mobility. Real cubic materials, instead, have electronic band extrema that are often degenerate or anisotropic and present several phonon modes. In the present work, we address such issues. We keep the continuum hypothesis inherent to the large polaron Fröhlich model, but waive the isotropic and non-degeneracy hypotheses, and also include multiple phonon branches. For polaron effective masses, working at the lowest order of perturbation theory, we provide analytical results for the case of anisotropic electronic energy dispersion, with two distinct effective masses (uniaxial) and numerical simulations for the degenerate 3-band case, typical of III-V and II-VI semiconductor valence bands. We also deal with the strong-coupling limit, using a variational treatment: we propose trial wavefunctions for the above-mentioned cases, providing polaron radii and energies. Then, we evaluate the polaron formation energies, effective masses and localisation lengths using parameters representative of a dozen II-VI, III-V and oxide semiconductors, for both electron and hole polarons...In the non-degenerate case, we compare the perturbative approach with the Feynman path integral approach in characterisizing polarons in the weak coupling limit... △ Less

Submitted 26 September, 2021; originally announced September 2021.

Comments: 18 pages, 10 figures in main text; 11 pages, 1 figure in supplementary material

Geant4-based calibration of an organic liquid scintillator

Authors: N. Mauritzson, K. G. Fissum, H. Perrey, J. R. M. Annand, R. J. W. Frost, R. Hall-Wilton, R. Al Jebali, K. Kanaki, V. Maulerova, F. Messi, E. Rofors

Abstract: A light-yield calibration of an NE 213A organic liquid scintillator detector has been performed using both monoenergetic and polyenergetic gamma-ray sources. Scintillation light was detected in a photomultiplier tube, and the corresponding pulses were subjected to waveform digitization on an event-by-event basis. The resulting Compton edges have been analyzed using a Geant4 simulation of the detec… ▽ More A light-yield calibration of an NE 213A organic liquid scintillator detector has been performed using both monoenergetic and polyenergetic gamma-ray sources. Scintillation light was detected in a photomultiplier tube, and the corresponding pulses were subjected to waveform digitization on an event-by-event basis. The resulting Compton edges have been analyzed using a Geant4 simulation of the detector which models both the interactions of the ionizing radiation as well as the transport of scintillation photons. The simulation is calibrated and also compared to well-established prescriptions used to determine the Compton edges, resulting ultimately in light-yield calibration functions. In the process, the simulation-based method produced information on the gain and intrinsic pulse-height resolution of the detector. It also facilitated a previously inaccessible understanding of the systematic uncertainties associated with the calibration of the scintillation-light yield. The simulation-based method was also compared to well-established numerical prescriptions for locating the Compton edges. Ultimately, the simulation predicted as much as 17% lower light-yield calibrations than the prescriptions. These calibrations indicate that approximately 35% of the scintillation light associated with a given gamma-ray reaches the photocathode. It is remarkable how well two 50 year old prescriptions for calibrating scintillation-light yield in organic scintillators have stood the test of time. △ Less

Submitted 22 October, 2021; v1 submitted 24 September, 2021; originally announced September 2021.

Comments: 20 pages, 8 figures, 4 tables

Tying the geometrical traits of massive young stellar objects and their discs to a potential evolutionary sequence using infrared observations

Authors: A. J. Frost, R. D. Oudmaijer, S. L. Lumsden, W-J de Wit

Abstract: Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is li… ▽ More Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is limited. Here, we combine the results of two studies with the aim to attribute geometrical features to an evolutionary sequence for a sample of seven MYSOs. The time evolution is based on a near-IR spectral features, while the geometry is determined from a multi size-scale study of MYSOs. We find that MYSO discs with determined geometrical substructure turn out to be also spectroscopically more evolved. This implies that disc evolution and dispersal are occurring within MYSOs, similar to low-mass YSO disc evolution, despite their faster formation timescales. △ Less

Submitted 12 August, 2021; v1 submitted 11 August, 2021; originally announced August 2021.

The first interferometric survey in the K-band of massive YSOs. On the hot dust, ionised gas, and binarity at au scales

Authors: E. Koumpia, W. -J. de Wit, R. D. Oudmaijer, A. J. Frost, S. Lumsden, A. Caratti o Garatti, S. P. Goodwin, B. Stecklum, I. Mendigutıa, J. D. Ilee, M. Vioque

Abstract: Circumstellar discs are essential for high mass star formation, while multiplicity, in particular binarity, appears to be an inevitable outcome since the vast majority of massive stars (> 8 Msun) are found in binaries (up to 100%). We constrain the sizes of the dust and ionised gas (Brgamma) emission of the innermost regions towards a sample of six MYSOs, and provide high-mass binary statistics of… ▽ More Circumstellar discs are essential for high mass star formation, while multiplicity, in particular binarity, appears to be an inevitable outcome since the vast majority of massive stars (> 8 Msun) are found in binaries (up to 100%). We constrain the sizes of the dust and ionised gas (Brgamma) emission of the innermost regions towards a sample of six MYSOs, and provide high-mass binary statistics of young stars at 2-300 au scales using VLTI (GRAVITY, AMBER) observations. We determine the inner radius of the dust emission and place MYSOs with K-band measurements in a size-luminosity diagram for the first time, and compare our findings to T Tauris and Herbig AeBes. We also compare the observed K-band sizes to the sublimation radius predicted by three different disc scenarios. Lastly, we apply binary geometries to trace close binarity among MYSOs. The inner sizes of MYSOs, Herbig AeBe and T Tauri stars appear to follow a universal trend at which the sizes scale with the square-root of the stellar luminosity. The Brgamma emission originates from somewhat smaller and co-planar area compared to the 2.2 μm continuum emission. We discuss this new finding with respect to disc-wind or jet origin. Finally, we report an MYSO binary fraction of 17-25% at milli-arcsecond separations (2-300 au). The size-luminosity diagram indicates that the inner regions of discs around young stars scale with luminosity independently of the stellar mass. At the targeted scales (2-300 au), the MYSO binary fraction is lower than what was previously reported for the more evolved main sequence massive stars, which, if further confirmed, could implicate the predictions from massive binary formation theories. Lastly, we spatially resolve the crucial star/disc interface in a sample of MYSOs, showing that au-scale discs are prominent in high-mass star formation and similar to their low-mass equivalents. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: Accepted for publication in A&A. 26 pages, 16 Figures, 7 Tables, full abstract to appear in A&A

Journal ref: A&A 654, A109 (2021)

Multi-Pulse Terahertz Spectroscopy Unveils Hot Polaron Photoconductivity Dynamics in Metal-Halide Perovskites

Authors: Xijia Zheng, Thomas R. Hopper, Andrei Gorodetsky, Marios Maimaris, Weidong Xu, Bradley A. A. Martin, Jarvist M. Frost, Artem A. Bakulin

Abstract: The behavior of hot carriers in metal-halide perovskites (MHPs) present a valuable foundation for understanding the details of carrier-phonon coupling in the materials as well as the prospective development of highly efficient hot carrier and carrier multiplication solar cells. Whilst the carrier population dynamics during cooling have been intensely studied, the evolution of the hot carrier prope… ▽ More The behavior of hot carriers in metal-halide perovskites (MHPs) present a valuable foundation for understanding the details of carrier-phonon coupling in the materials as well as the prospective development of highly efficient hot carrier and carrier multiplication solar cells. Whilst the carrier population dynamics during cooling have been intensely studied, the evolution of the hot carrier properties, namely the hot carrier mobility, remain largely unexplored. To address this, we introduce a novel ultrafast visible pump - infrared push - terahertz probe spectroscopy (PPP-THz) to monitor the real-time conductivity dynamics of cooling carriers in methylammonium lead iodide. We find a decrease in mobility upon optically depositing energy into the carriers, which is typical of band-transport. Surprisingly, the conductivity recovery dynamics are incommensurate with the intraband relaxation measured by an analogous experiment with an infrared probe (PPP- IR), and exhibit a negligible dependence on the density of hot carriers. These results and the kinetic modelling reveal the importance of highly-localized lattice heating on the mobility of the hot electronic states. This collective polaron-lattice phenomenon may contribute to the unusual photophysics observed in MHPs and should be accounted for in devices that utilize hot carriers. △ Less

Submitted 19 August, 2021; v1 submitted 30 June, 2021; originally announced June 2021.

Comments: 28 pages, 4 figures, 77 references

Resolving the dynamical mass tension of the massive binary 9 Sagittarii

Authors: M. Fabry, C. Hawcroft, A. J. Frost, L. Mahy, P. Marchant, J-B. Le Bouquin, H. Sana

Abstract: Direct dynamical mass measurements of stars with masses above 30 M${}_\odot$ are rare. This is the result of the low yield of the upper initial mass function and the limited number of such systems in eclipsing binaries. Long-period, double-lined spectroscopic binaries that are also resolved astrometrically offer an alternative for obtaining absolute masses of stellar objects. 9 Sgr is one such lon… ▽ More Direct dynamical mass measurements of stars with masses above 30 M${}_\odot$ are rare. This is the result of the low yield of the upper initial mass function and the limited number of such systems in eclipsing binaries. Long-period, double-lined spectroscopic binaries that are also resolved astrometrically offer an alternative for obtaining absolute masses of stellar objects. 9 Sgr is one such long-period, high-mass binary. Unfortunately, a large amount of tension exists between its total dynamical mass inferred from radial velocity measurements and that from astrometric data. We obtained the astrometric orbit from VLTI/PIONIER and VLTI/GRAVITY interferometric measurements. Using archival and new spectroscopy, we performed a grid-based spectral disentangling search to constrain the semi-amplitudes of the radial velocity curves. We computed atmospheric parameters and surface abundances by adjusting \textsc{fastwind} atmosphere models and we compared our results with evolutionary tracks computed with the Bonn Evolutionary Code (BEC). Grid spectral disentangling of 9 Sgr supports the presence of a 53 M${}_\odot$ primary and a 39 M${}_\odot$ secondary. Comparison with BEC evolutionary tracks shows the components of 9 Sgr are most likely coeval with an age of roughly 1 Myr. Our analysis clears up the contradiction between mass and orbital inclination estimates reported in previous studies. We detect the presence of significant CNO-processed material at the surface of the primary, suggesting enhanced internal mixing compared to currently implemented in the BEC models. The present measurements provide a high-quality high-mass anchor to validate stellar evolution models and to test the efficiency of internal mixing processes. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 14 pages, 5 figures, 4 tables (+6 figures and 3 tables in appendices). Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 651, A119 (2021)

Giant Huang-Rhys Factor for Electron Capture by the Iodine Interstitial in Perovskite Solar Cells

Authors: Lucy D. Whalley, Puck van Gerwen, Jarvist M. Frost, Sunghyun Kim, Samantha N. Hood, Aron Walsh

Abstract: Improvement in the optoelectronic performance of halide perovskite semiconductors requires the identification and suppression of non-radiative carrier trapping processes. The iodine interstitial has been established as a deep level defect, and implicated as an active recombination centre. We analyse the quantum mechanics of carrier trapping. Fast and irreversible electron capture by the neutral io… ▽ More Improvement in the optoelectronic performance of halide perovskite semiconductors requires the identification and suppression of non-radiative carrier trapping processes. The iodine interstitial has been established as a deep level defect, and implicated as an active recombination centre. We analyse the quantum mechanics of carrier trapping. Fast and irreversible electron capture by the neutral iodine interstitial is found. The effective Huang-Rhys factor exceeds 300, indicative of the strong electron-phonon coupling that is possible in soft semiconductors. The accepting phonon mode has a frequency of 53 cm$^{-1}$ and has an associated electron capture coefficient of 10$^{-10}$cm$^3$s$^{-1}$. The inverse participation ratio is used to quantify the localisation of phonon modes associated with the transition. We infer that suppression of octahedral rotations is an important factor to enhance defect tolerance. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Accelerated hot-carrier cooling in MAPbI3 perovskite by pressure-induced lattice compression

Authors: Loreta A. Muscarella, Eline M. Hutter, Jarvist M. Frost, Gianluca G. Grimaldi, Jan Versluis, Huib J. Bakker, Bruno Ehrler

Abstract: Hot-carrier cooling (HCC) in metal halide perovskites in the high-density regime is significantly slower compared to conventional semiconductors. This effect is commonly attributed to a hot-phonon bottleneck but the influence of the lattice properties on the HCC behaviour is poorly understood. Using pressure-dependent transient absorption spectroscopy (fs-TAS) we find that at an excitation density… ▽ More Hot-carrier cooling (HCC) in metal halide perovskites in the high-density regime is significantly slower compared to conventional semiconductors. This effect is commonly attributed to a hot-phonon bottleneck but the influence of the lattice properties on the HCC behaviour is poorly understood. Using pressure-dependent transient absorption spectroscopy (fs-TAS) we find that at an excitation density below Mott transition, pressure does not affect the HCC. On the contrary, above Mott transition, HCC in methylammonium lead iodide (MAPbI3) is around two times as fast at 0.3 GPa compared to ambient pressure. Our electron-phonon coupling calculations reveal about two times stronger electron-phonon coupling for the inorganic cage mode at 0.3 GPa. However, our experiments reveal that pressure promotes faster HCC only above Mott transition. Altogether, these findings suggest a change in the nature of excited carriers in the high-density regime, providing insights on the electronic behavior of devices operating at such high charge-carrier density. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Unveiling the traits of massive young stellar objects through a multi-scale survey

Authors: A. J. Frost, R. D. Oudmaijer, W. J. de Wit, S. L. Lumsden

Abstract: The rarity and deeply embedded nature of stars with masses larger than 8 solar masses has limited our understanding of their formation. Previous work has shown that complementing spectral energy distributions with interferometric and imaging data can probe the circumstellar environments of massive young stellar objects (MYSOs) well. However, complex studies of single objects often use different ap… ▽ More The rarity and deeply embedded nature of stars with masses larger than 8 solar masses has limited our understanding of their formation. Previous work has shown that complementing spectral energy distributions with interferometric and imaging data can probe the circumstellar environments of massive young stellar objects (MYSOs) well. However, complex studies of single objects often use different approaches in their analysis. Therefore the results of these studies cannot be directly compared. This work aims to obtain the physical characteristics of a sample of MYSOs at ~0.01" scales, at 0.1" scales, and as a whole, which enables us to compare the characteristics of the sources. We apply the same multi-scale method and analysis to a sample of MYSOs. High-resolution interferometric data, near-diffraction-limited imaging data, and a multi-wavelength spectral energy distribution are combined. By fitting simulated observables derived from 2.5D radiative transfer models of disk-outflow-envelope systems to our observations, the properties of the MYSOs are constrained. We find that the observables of all the MYSOs can be reproduced by models with disk-outflow-envelope geometries, analogous to the Class I geometry associated with low-mass protostars. The characteristics of the envelopes and the cavities within them are very similar across our sample. On the other hand, the disks seem to differ between the objects, in particular with regards to what we interpret as evidence of complex structures and inner holes. This is comparable to the morphologies observed for low-mass young stellar objects. A strong correlation is found between the luminosity of the central MYSO and the size of the transition disk-like inner hole for the MYSOs, implying that photoevaporation or the presence of binary companions may be the cause. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Journal ref: A&A 648, A62 (2021)

Is HR 6819 a triple system containing a black hole? -- An alternative explanation

Authors: J. Bodensteiner, T. Shenar, L. Mahy, M. Fabry, P. Marchant, M. Abdul-Masih, G. Banyard, D. M. Bowman, K. Dsilva, A. J. Frost, C. Hawcroft, M. Reggiani, H. Sana

Abstract: HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate t… ▽ More HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Submitted to A&A, 13 pages (16 figures and 2 tables); 4 pages supplementary material (4 figures and 4 tables). Comments are welcome

Journal ref: A&A 641, A43 (2020)

Neural Parameter Allocation Search

Authors: Bryan A. Plummer, Nikoli Dryden, Julius Frost, Torsten Hoefler, Kate Saenko

Abstract: Training neural networks requires increasing amounts of memory. Parameter sharing can reduce memory and communication costs, but existing methods assume networks have many identical layers and utilize hand-crafted sharing strategies that fail to generalize. We introduce Neural Parameter Allocation Search (NPAS), a novel task where the goal is to train a neural network given an arbitrary, fixed par… ▽ More Training neural networks requires increasing amounts of memory. Parameter sharing can reduce memory and communication costs, but existing methods assume networks have many identical layers and utilize hand-crafted sharing strategies that fail to generalize. We introduce Neural Parameter Allocation Search (NPAS), a novel task where the goal is to train a neural network given an arbitrary, fixed parameter budget. NPAS covers both low-budget regimes, which produce compact networks, as well as a novel high-budget regime, where additional capacity can be added to boost performance without increasing inference FLOPs. To address NPAS, we introduce Shapeshifter Networks (SSNs), which automatically learn where and how to share parameters in a network to support any parameter budget without requiring any changes to the architecture or loss function. NPAS and SSNs provide a complete framework for addressing generalized parameter sharing, and can also be combined with prior work for additional performance gains. We demonstrate the effectiveness of our approach using nine network architectures across four diverse tasks, including ImageNet classification and transformers. △ Less

Submitted 15 March, 2022; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Accepted at ICLR 2022

Higgs self-coupling measurements using deep learning in the $b\bar{b}b\bar{b}$ final state

Authors: Jacob Amacker, William Balunas, Lydia Beresford, Daniela Bortoletto, James Frost, Cigdem Issever, Jesse Liu, James McKee, Alessandro Micheli, Santiago Paredes Saenz, Michael Spannowsky, Beojan Stanislaus

Abstract: Measuring the Higgs trilinear self-coupling $λ_{hhh}$ is experimentally demanding but fundamental for understanding the shape of the Higgs potential. We present a comprehensive analysis strategy for the HL-LHC using di-Higgs events in the four $b$-quark channel ($hh \to 4b$), extending current methods in several directions. We perform deep learning to suppress the formidable multijet background wi… ▽ More Measuring the Higgs trilinear self-coupling $λ_{hhh}$ is experimentally demanding but fundamental for understanding the shape of the Higgs potential. We present a comprehensive analysis strategy for the HL-LHC using di-Higgs events in the four $b$-quark channel ($hh \to 4b$), extending current methods in several directions. We perform deep learning to suppress the formidable multijet background with dedicated optimisation for BSM $λ_{hhh}$ scenarios. We compare the $λ_{hhh}$ constraining power of events using different multiplicities of large radius jets with a two-prong structure that reconstruct boosted $h \to bb$ decays. We show that current uncertainties in the SM top Yukawa coupling $y_t$ can modify $λ_{hhh}$ constraints by $\sim 20\%$. For SM $y_t$, we find prospects of $-0.8 < λ_{hhh} / λ_{hhh}^\text{SM} < 6.6$ at 68% CL under simplified assumptions for 3000~fb$^{-1}$ of HL-LHC data. Our results provide a careful assessment of di-Higgs identification and machine learning techniques for all-hadronic measurements of the Higgs self-coupling and sharpens the requirements for future improvement. △ Less

Submitted 12 October, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 36 pages, 15 figures + bibliography and appendices

Report number: IPPP/20/11

Journal ref: JHEP 12 (2020) 115

Sparse hierarchical representation learning on molecular graphs

Authors: Matthias Bal, Hagen Triendl, Mariana Assmann, Michael Craig, Lawrence Phillips, Jarvist Moore Frost, Usman Bashir, Noor Shaker, Vid Stojevic

Abstract: Architectures for sparse hierarchical representation learning have recently been proposed for graph-structured data, but so far assume the absence of edge features in the graph. We close this gap and propose a method to pool graphs with edge features, inspired by the hierarchical nature of chemistry. In particular, we introduce two types of pooling layers compatible with an edge-feature graph-conv… ▽ More Architectures for sparse hierarchical representation learning have recently been proposed for graph-structured data, but so far assume the absence of edge features in the graph. We close this gap and propose a method to pool graphs with edge features, inspired by the hierarchical nature of chemistry. In particular, we introduce two types of pooling layers compatible with an edge-feature graph-convolutional architecture and investigate their performance for molecules relevant to drug discovery on a set of two classification and two regression benchmark datasets of MoleculeNet. We find that our models significantly outperform previous benchmarks on three of the datasets and reach state-of-the-art results on the fourth benchmark, with pooling improving performance for three out of four tasks, keeping performance stable on the fourth task, and generally speeding up the training process. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 4 pages, 2 figures, accepted as a DLG 2019 workshop paper at KDD 2019

The Neutral Gas Properties of Extremely Isolated Early-Type Galaxies III

Authors: Trisha Ashley, Pamela M. Marcum, Mehmet Alpaslan, Michael N. Fanelli, James D. Frost

Abstract: We report on the neutral hydrogen gas content (21-cm emission) of eight extremely isolated early-type galaxies (IEGs) using the Green Bank Telescope. Emission is detected in seven of the eight objects. This paper is the third in a series that collectively present new HI observations for 20 IEGs. Among the 14 HI detections in our observations, eight exhibit a Gaussian-like HI line profile shape, fo… ▽ More We report on the neutral hydrogen gas content (21-cm emission) of eight extremely isolated early-type galaxies (IEGs) using the Green Bank Telescope. Emission is detected in seven of the eight objects. This paper is the third in a series that collectively present new HI observations for 20 IEGs. Among the 14 HI detections in our observations, eight exhibit a Gaussian-like HI line profile shape, four are double-peaked, one is triple-peaked, and another has a plateaued rectangular shape. Five additional IEGs observed in previous surveys were added to our analysis, bringing the total number of IEGs with HI observations to 25. Of these objects, emission is detected in 19 (76%). The 25 IEGs in our combined study have gas masses that are systematically larger than their luminosity-matched comparison galaxy counterparts. The isolated early-type galaxies presented here follow a trend of increasing gas-richness with bluer B-V colors. This correlation is also observed in a comparison sample drawn from the literature composed of loose group and field early-type galaxies. Two IEGs, KIG164 and KIG870, exhibit properties highly anomalous for spheroidal systems: luminous (M_B = -20.5, -20.1) and blue (B-V = 0.47, 0.48) respectively, with substantial neutral gas, M_HI = 4.1 & 5.5x10^9 M_(sun). Other IEG systems may represent early-type galaxies continuing to assemble via quiescent HI accretion from the cosmic web or relaxed merged systems. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: 11 pages; published in the Astronomical Journal

A multi-scale exploration of a massive young stellar object - a transition disk around G305.20+0.21?

Authors: A. J. Frost, R. D. Oudmaijer, W. J. de Wit, S. L. Lumsden

Abstract: The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplem… ▽ More The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplementary to such a dataset. By consolidating these observations with modelling, the features of a massive protostellar environment can be constrained. This work simultaneously fits the aforementioned observations and a spectral energy distribution (SED) with a 2.5D radiative transfer model, providing an extensive view of the physical characteristics of the accreting regions of the MYSO G305.20+0.21. The high-resolution observations were obtained using the Very Large Telescope's MIDI and VISIR instruments, producing visibilities in the N-band and near-diffraction-limited imaging in the Q-band respectively. A model including a central protostar with a luminosity of $\sim$5$\times$10$^{4}$L$_{\odot}$ surrounded by a low-density bipolar cavity, a flared 1M$_{\odot}$ disk and an envelope provides a sufficient fit all three types of observation. The need to include a disk in the model implies that this MYSO follows a scaled-up version of the low-mass star formation process. The weak silicate absorption feature within the SED requires low-density envelope cavities to be successfully fit and is an atypical characteristic in comparison to previously studied MYSOs. Additionally, the inner radius of the disk must be three times the dust sublimation radius to satisfy the MIDI visibilities. The low density, low extinction environment implies the object is a more evolved MYSO and this combined with large inner radius of the disk suggests that it could be an example of a transitional disk around an MYSO. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 16 pages, 8 figures

Journal ref: A&A 625, A44 (2019)

Review: Solid-state physics of halide perovskites

Authors: Jarvist Moore Frost

Abstract: Halide perovskite solar cells presented a unique opportunity to apply modern computational materials science techniques to an (initially) poorly understood new material. In this review, we recount the key understanding developed during the last five years, through a narrative review of research progress. The central enigma of the material is how it can be so defective, and yet work so well as a ph… ▽ More Halide perovskite solar cells presented a unique opportunity to apply modern computational materials science techniques to an (initially) poorly understood new material. In this review, we recount the key understanding developed during the last five years, through a narrative review of research progress. The central enigma of the material is how it can be so defective, and yet work so well as a photovoltaic. The physical properties of the material were understood through molecular and lattice dynamic calculations, revealing the material to show large dynamic responses on a wide range of time scales. Longer length scales in the material was simulated with effective classical potentials, showing that complex domains can be generated by the interacting molecular dipoles, generating structured features in the electrostatic potential of the lattice. Relativistic electronic structure reveals unique features in the bands, which may explain observed slow recombination, and could be used in high efficiency photovoltaics. The large dielectric response of the lattice leads to a strong drive for the formation of polarons, some device physics of which are discussed. These polarons offer a possible explanation for the observed slow cooling of photoexcitations in the material. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: 10 pages

Impact of non-parabolic electronic band structure on the optical and transport properties of photovoltaic materials

Authors: Lucy D. Whalley, Jarvist M. Frost, Benjamin J. Morgan, Aron Walsh

Abstract: For semiconductors used in photovoltaic devices, the effective mass approximation allows calculation of important material properties from first-principles calculations, including optical properties (e.g. exciton binding energies), defect properties (e.g. donor and acceptor levels) and transport properties (e.g. carrier mobilities). The conduction and valence bands of semiconductors are commonly a… ▽ More For semiconductors used in photovoltaic devices, the effective mass approximation allows calculation of important material properties from first-principles calculations, including optical properties (e.g. exciton binding energies), defect properties (e.g. donor and acceptor levels) and transport properties (e.g. carrier mobilities). The conduction and valence bands of semiconductors are commonly approximated as parabolic around their extrema, which gives a simple theoretical description, but ignores the complexity of real materials. In this work, we use density functional theory to assess the impact of band non-parabolicity on four common thin-film photovoltaic materials - GaAs, CdTe, Cu$_2$ZnSnS$_4$ and CH$_3$NH$_3$PbI$_3$ - at temperatures and carrier densities relevant for real-world applications. First, we calculate the effective mass at the band edges. We compare finite-difference, unweighted least-squares and thermally weighted least-squares approaches. We find that the thermally weighted least-squares method reduces sensitivity to the choice of sampling density. Second, we employ a Kane quasi-linear dispersion to quantify the extent of non-parabolicity, and compare results from different electronic structure theories to consider the effect of spin-orbit coupling and electron exchange. Finally, we focus on the halide perovskite CH$_3$NH$_3$PbI$_3$ as a model system to assess the impact of non-parabolicity on calculated electron transport and optical properties at high carrier concentrations. We find that at a concentration of 10$^{20}$ cm$^-3$ the optical effective mass increases by a factor of two relative to the low carrier-concentration value, and the polaron mobility decreases by a factor of three. Our work suggests that similar adjustments should be made to the predicted optical and transport properties of other semiconductors with significant band non-parabolicity. △ Less

Submitted 31 January, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: Updated after review process

Journal ref: Phys. Rev. B 99, 085207 (2019)

Perspective: Dielectric and Ferroic Properties of Halide Perovskite Solar Cells

Authors: Jacob N. Wilson, Jarvist M. Frost, Suzanne K. Wallace, Aron Walsh

Abstract: Halide perovskite semiconductors respond to electric fields in a way that varies across time and length scales. We discuss the microscopic processes that give rise to the macroscopic polarization of these materials, ranging from the optical and vibrational response to the transport of ions and electrons. The strong frequency dependence of the dielectric permittivity can be understood by separating… ▽ More Halide perovskite semiconductors respond to electric fields in a way that varies across time and length scales. We discuss the microscopic processes that give rise to the macroscopic polarization of these materials, ranging from the optical and vibrational response to the transport of ions and electrons. The strong frequency dependence of the dielectric permittivity can be understood by separating the static dielectric constant into its constituents, including the orientional polarization due to rotating dipoles, which connects theory with experimental observations. The controversial issue of ferroelectricity is addressed, where we highlight recent progress in materials and domain characterization, but emphasize the challenge associated with isolating spontaneous lattice polarization from other processes such as charged defect formation and transport. We conclude that CH$_3$NH$_3$PbI$_3$ exhibits many features characteristic of a ferroelastic electret, where a spontaneous lattice strain is coupled to long-lived metastable polarization states. △ Less

Submitted 3 December, 2018; v1 submitted 5 November, 2018; originally announced November 2018.

Acoustic Phonon Lifetimes Limit Thermal Transport in Methylammonium Lead Iodide

Authors: Aryeh Gold-Parker, Peter M. Gehring, Jonathan M. Skelton, Ian C. Smith, Dan Parshall, Jarvist M. Frost, Hemamala I. Karunadasa, Aron Walsh, Michael F. Toney

Abstract: Hybrid organic-inorganic perovskites (HOIPs) have become an important class of semiconductors for solar cells and other optoelectronic applications. Electron-phonon coupling plays a critical role in all optoelectronic devices, and although the lattice dynamics and phonon frequencies of HOIPs have been well studied, little attention has been given to phonon lifetimes. We report the first high-preci… ▽ More Hybrid organic-inorganic perovskites (HOIPs) have become an important class of semiconductors for solar cells and other optoelectronic applications. Electron-phonon coupling plays a critical role in all optoelectronic devices, and although the lattice dynamics and phonon frequencies of HOIPs have been well studied, little attention has been given to phonon lifetimes. We report the first high-precision measurements of acoustic phonon lifetimes in the hybrid perovskite methylammonium lead iodide (MAPI), using inelastic neutron spectroscopy to provide high energy resolution and fully deuterated single crystals to reduce incoherent scattering from hydrogen. Our measurements reveal extremely short lifetimes on the order of picoseconds, corresponding to nanometer mean free paths and demonstrating that acoustic phonons are unable to dissipate heat efficiently. Lattice-dynamics calculations using ab-initio third-order perturbation theory indicate that the short lifetimes stem from strong three-phonon interactions and a high density of low-energy optical phonon modes related to the degrees of freedom of the organic cation. Such short lifetimes have significant implications for electron-phonon coupling in MAPI and other HOIPs, with direct impacts on optoelectronic devices both in the cooling of hot carriers and in the transport and recombination of band edge carriers. These findings illustrate a fundamental difference between HOIPs and conventional photovoltaic semiconductors and demonstrate the importance of understanding lattice dynamics in the effort to develop metal halide perovskite optoelectronic devices. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: 31 pages, 14 figures

Journal ref: PNAS 115, 11905 (2018)

Bioadhesive Graft-Antenna for Stimulation and Repair of Peripheral Nerves

Authors: Ashour Sliow, Zhi Ma, Gaetano Gargiulo, David Mahns, Damia Mawad, Paul Breen, Marcus Stoodley, Jessica Houang, Rhiannon Kuchel, Giuseppe Tettamanzi, Richard D Tilley, Samuel J Frost, John Morley, Leonardo Longo, Antonio Lauto

Abstract: Peripheral nerve injuries are difficult to treat due to limited axon regeneration; brief electrical stimulation of injured nerves is an emerging therapy that can relieve pain and enhance regeneration. We report an original wireless stimulator based on a metal loop (diameter ~1 mm) that is powered by a transcranial magnetic stimulator (TMS). The loop can be integrated in a chitosan scaffold that fu… ▽ More Peripheral nerve injuries are difficult to treat due to limited axon regeneration; brief electrical stimulation of injured nerves is an emerging therapy that can relieve pain and enhance regeneration. We report an original wireless stimulator based on a metal loop (diameter ~1 mm) that is powered by a transcranial magnetic stimulator (TMS). The loop can be integrated in a chitosan scaffold that functions as a graft when applied onto transected nerves (graft-antenna). The graft-antenna was bonded to rat sciatic nerves by a laser without sutures; it did not migrate after implantation and was able to trigger steady compound muscle action potentials for 12 weeks (CMAP ~1.3 mV). Eight weeks post-operatively, axon regeneration was facilitated in transected nerves that were repaired with the graft-antenna and stimulated by the TMS for 1 hour/week. The graft-antenna is an innovative and minimally-invasive device that functions concurrently as a wireless stimulator and adhesive scaffold for nerve repair. △ Less

Submitted 8 July, 2018; originally announced July 2018.

Comments: 33 pages, 9 figures (including 1 supplementary figure), 5 tables (including 3 supplementary tables)