Abstract visibility
Abstract text size

Papers for Monday, Aug 22 2022

Papers with local authors

Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold, Richard Anantua, Silke Britzen, Chiara Ceccobello, Koushik Chatterjee, Yongjun Chen, Nicholas S. Conroy, Geoffrey B. Crew, Alejandro Cruz-Osorio, Yuzhu Cui, Sheperd S. Doeleman, Razieh Emami, Joseph Farah, Christian M. Fromm, Peter Galison, Boris Georgiev, Luis C. Ho, David J. James, Britton Jeter, Alejandra Jimenez-Rosales, Jun Yi Koay, Carsten Kramer, Thomas P. Krichbaum, Sang-Sung Lee, Michael Lindqvist, Ivan Martí-Vidal, Karl M. Menten, Yosuke Mizuno, James M. Moran, Monika Moscibrodzka, Antonios Nathanail, Joey Neilsen, Chunchong Ni, Jongho Park, Vincent Pietu, Luciano Rezzolla, Angelo Ricarte, Bart Ripperda, Lijing Shao, Fumie Tazaki, Kenji Toma, Pablo Torne, Jonathan Weintroub, Maciek Wielgus, Feng Yuan, et al.
0 votes
Paper 4 — arXiv:2208.09004
0 votes
Paper 4 — arXiv:2208.09004

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~\mu{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

Tejas Prasanna, Matthew S. B. Coleman, Matthias J. Raives, Todd A. Thompson

16 pages, 10 figures

0 votes
Paper 11 — arXiv:2208.09042
0 votes
Paper 11 — arXiv:2208.09042

In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and super-luminous supernovae. Since the flow is non-relativistic at early times, and because the Alfv\'en radius is much larger than the proto-magnetar radius, spindown is millions of times more efficient than the typically-used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spindown. For polar magnetic field strengths $B_0\gtrsim5\times10^{14}$ G, the spindown timescale can be shorter than than the Kelvin-Helmholtz timescale. For $B_0\gtrsim10^{15}$ G, it is of order seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of $B_0$, $P_{\star0}$, and mass ($M$). Proto-magnetars born with $B_0$ greater than $\simeq1.3\times10^{15}\,{\rm\,G}\,(P_{\star0}/{400\,\rm\,ms})^{-1.4}(M/1.4\,{\rm M}_\odot)^{2.2}$ spin down to periods $> 1$ s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spindown. Spindown is more efficient for lower $M$ and for larger $P_{\star0}$. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348-5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars.

David T. Chuss, Brandon S. Hensley, Alan J. Kogut, Jordan A. Guerra, Hayley C. Nofi, Javad Siah

8 pages, 6 figures, submitted to ApJ

0 votes
Paper 12 — arXiv:2208.09049
0 votes
Paper 12 — arXiv:2208.09049

We use archival COBE/DIRBE data to construct a map of polycyclic aromatic hydrocarbon (PAH) emission in the $\lambda$-Orionis region. The presence of the 3.3 $\mu$m PAH feature within the DIRBE 3.5 $\mu$m band and corresponding lack of signficant PAH spectral features in the adjacent DIRBE bands (1.25 $\mu$m, 2.2 $\mu$m, and 4.9 $\mu$m) enable estimation of the PAH contribution to the 3.5 $\mu$m data. Having the shortest wavelength of known PAH features, the 3.3 $\mu$m feature probes the smallest PAHs, which are also the leading candidates for carriers of anomalous microwave emission (AME). We use this map to investigate the association between the AME and the emission from PAH molecules. We find that the spatial correlation in $\lambda$-Orionis is higher between AME and far-infrared dust emission (as represented by the DIRBE 240 $\mu$m map) than it is between our PAH map and AME. This finding, in agreement with previous studies using PAH features at longer wavelengths, is in tension with the hypothesis that AME is due to spinning PAHs. However, the expected correlation between mid-infrared and microwave emission could potentially be degraded by different sensitivities of each emission mechanism to local environmental conditions even if PAHs are the carriers of both.

All other papers

A.K. Hughes, D.J. Sand, A. Seth, J. Strader, C. Lidman, K. Voggel, A. Dumont, D. Crnojević, M. Mateo, N. Caldwell, D.A. Forbes, S. Pearson, P. Guhathakurta, E. Toloba

24 pages, 13 figures, submitted to AAS Journals

We present new radial velocity measurements from the Magellan and the Anglo-Australian Telescopes for 174 previously known and 122 newly confirmed globular clusters (GCs) around NGC 5128, the nearest accessible massive early-type galaxy at D=3.8 Mpc. Remarkably, 28 of these newly confirmed GCs are at projected radii >50' ($\gtrsim 54$ kpc), extending to $\sim 130$ kpc, in the outer halo where few GCs had been confirmed in previous work. We identify several subsets of GCs that spatially trace halo substructures that are visible in red giant branch star maps of the galaxy. In some cases, these subsets of GCs are kinematically cold, and may be directly associated with and originate from these specific stellar substructures. From a combined kinematic sample of 645 GCs, we see evidence for coherent rotation at all radii, with a higher rotation amplitude for the metal-rich GC subpopulation. Using the tracer mass estimator, we measure a total enclosed mass of $2.5\pm0.3 \times 10^{12} M_{\odot}$ within $\sim 120$ kpc, an estimate that will be sharpened with forthcoming dynamical modeling. The combined power of stellar mapping and GC kinematics makes NGC 5128 an ongoing keystone for understanding galaxy assembly at mass scales inaccessible in the Local Group.

I. Agudo, L. Amati, T. An, F. E. Bauer, S. Benetti, M. G. Bernardini, R. Beswick, K. Bhirombhakdi, T. de Boer, M. Branchesi, S. J. Brennan, M. D. Caballero-García, E. Cappellaro, N. Castro Rodríguez, A. J. Castro-Tirado, K. C. Chambers, E. Chassande-Mottin, S. Chaty, T.-W. Chen, A. Coleiro, S. Covino, F. D'Ammando, P. D'Avanzo, V. D'Elia, A. Fiore, A. Flörs, M. Fraser, S. Frey, C. Frohmaier, L. Galbany, C. Gall, H. Gao, J. García-Rojas, G. Ghirlanda, S. Giarratana, J. H. Gillanders, M. Giroletti, B. P. Gompertz, M. Gromadzki, K. E. Heintz, Y.-D. Hu, M. E. Huber, A. Inkenhaag, L. Izzo, Z. P. Jin, P. G. Jonker, D. A. Kann, E. C. Kool, R. Kotak, G. Leloudas, A. J. Levan, C.-C. Lin, J. D. Lyman, E. A. Magnier, K. Maguire, I. Mandel, B. Marcote, D. Mata Sánchez, et al. (41 additional authors not shown)

By the ENGRAVE collaboration (engrave-eso.org). 33 pages, 12 figures, submitted to A&A. Comments are welcome!

We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns.

Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold

Generating images from very long baseline interferometric observations poses a difficult, and generally not unique, inversion problem. This problem is simplified by the introduction of constraints, some generic (e.g., positivity of the intensity) and others motivated by physical considerations (e.g., smoothness, instrument resolution). It is further complicated by the need to simultaneously address instrumental systematic uncertainties and sparse coverage in the u-v plane. We report a new Bayesian image reconstruction technique in the parameter estimation framework Themis that has been developed for the Event Horizon Telescope. This has two key features: first, the full Bayesian treatment of the image reconstruction makes it possible to generate a full posterior for the images, permitting a rigorous and quantitative investigation into the statistical significance of image features. Second, it is possible to seamlessly incorporate directly modeled features simultaneously with image reconstruction. We demonstrate this second capability by incorporating a narrow, slashed ring in reconstructions of simulated M87 data in an attempt to detect and characterize the photon ring. We show that it is possible to obtain high-fidelity photon ring sizes, enabling mass measurements with accuracies of 2%-5% that are essentially insensitive to astrophysical uncertainties, and creating opportunities for precision tests of general relativity.

Brian Welch, Dan Coe, Erik Zackrisson, S.E. de Mink, Swara Ravindranath, Jay Anderson, Gabriel Brammer, Larry Bradley, Jinmi Yoon, Patrick Kelly, Jose M. Diego, Rogier Windhorst, Adi Zitrin, Paola Dimauro, Yolanda Jimenez-Teja, Abdurro'uf, Mario Nonino, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Matthew B. Bayliss, Alex Benitez, Tom Broadhurst, Rachana Bhatawdekar, Marusa Bradac, Gabriel Caminha, Wenlei Chen, Jan Eldridge, Ebraheem Farag, Michael Florian, Brenda Frye, Seiji Fujimoto, Sebastian Gomez, Alaina Henry, Tiger Y.-Y Hsiao, Taylor A. Hutchison, Bethan L. James, Meridith Joyce, Intae Jung, Gourav Khullar, Rebecca L. Larson, Guillaume Mahler, Nir Mandelker, Stephan McCandliss, Takahiro Morishita, Rosa Newshore, Colin Norman, Kyle O'Connor, Pascal A. Oesch, Masamune Oguri, et al. (13 additional authors not shown)

Submitted to ApJ. Data products, lens models, and analysis code will be available online at this https URL

The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $\zphot = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (\JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$\mu$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $\mu > 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system. The new imaging also supports the previous photometric redshift estimate, solidifying the interpretation of Earendel as a stellar source within the first billion years of the universe. Fitting stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Allowing for two stars with different temperatures can yield better fits to the photometry at the expense of increased free parameters, leading to non-unique well-fitting temperatures of $T_{{\rm eff}} \gtrsim 20000$ K and $T_{{\rm eff}} \lesssim 13000$ K for the hotter and cooler components, respectively. Follow-up observations, including \JWST\ NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe.

Michael J. Greener, Alfonso Aragón-Salamanca, Michael Merrifield, Thomas Peterken, Elizaveta Sazonova, Roan Haggar, Dmitry Bizyaev, Joel R. Brownstein, Richard R. Lane, Kaike Pan

15 pages, 8 figures, accepted for publication by MNRAS

We investigate archaeologically how the metallicity in both stellar and gaseous components of spiral galaxies of differing masses evolve with time, using data from the SDSS-IV MaNGA survey. For the stellar component, we can measure this evolution directly by decomposing the galaxy absorption-line spectra into populations of different ages and determining their metallicities. For the gaseous component, we can only measure the present-day metallicity directly from emission lines. However, there is a well-established relationship between gas metallicity, stellar mass and star formation rate which does not evolve significantly with redshift; since the latter two quantities can be determined directly for any epoch from the decomposition of the absorption-line spectra, we can use this relationship to infer the variation in gas metallicity over cosmic time. Comparison of present-day values derived in this way with those obtained directly from the emission lines confirms the validity of the method. Application of this approach to a sample of 1619 spiral galaxies reveals how the metallicity of these systems has changed over the last 10 billion years since cosmic noon. For lower-mass galaxies, both stellar and gaseous metallicity increase together, as one might expect in well-mixed fairly isolated systems. In higher-mass systems, the average stellar metallicity has not increased in step with the inferred gas metallicity, and actually decreases with time. Such disjoint behaviour is what one might expect if these more massive systems have accreted significant amounts of largely pristine gas over their lifetimes, and this material has not been well mixed into the galaxies.

Hinna Shivkumar, Amruta D. Jaodand, Arvind Balasubramanian, Christoffer Fremling, Alessandra Corsi, Anastasios Tzanidakis, Samaya Nissanke, Mansi Kasliwal, Murray Brightman, Geert Raaijmakers, Kristin Kruse Madsen, Fiona Harrison, Dario Carbone, Nayana A.J., Jean-Michel Désert, Igor Andreoni

20 pages, 14 figures, Submitted to ApJ

We present optical, radio and X-ray observations of a rapidly-evolving transient AT2019wxt (PS19hgw), discovered during the search for an electromagnetic (EM) counterpart to the gravitational-wave (GW) trigger S191213g (LIGO Scientific Collaboration & Virgo Collaboration 2019a). Although S191213g was not confirmed as a significant GW event in the off-line analysis of LIGO-Virgo data, AT2019wxt remained an interesting transient due its peculiar nature. The optical/NIR light curve of AT2019wxt displayed a double-peaked structure evolving rapidly in a manner analogous to currently know ultra-stripped supernovae (USSNe) candidates. This double-peaked structure suggests presence of an extended envelope around the progenitor, best modelled with two-components: i) early-time shock-cooling emission and ii) late-time radioactive $^{56}$Ni decay. We constrain the ejecta mass of AT2019wxt at $M_{ej} \approx{0.20 M_{\odot}}$ which indicates a significantly stripped progenitor that was possibly in a binary system. We also followed-up AT2019wxt with long-term Chandra and Jansky Very Large Array observations spanning $\sim$260 days. We detected no definitive counterparts at the location of AT2019wxt in these long-term X-ray and radio observational campaigns. We establish the X-ray upper limit at $9.93\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and detect an excess radio emission from the region of AT2019wxt. However, there is little evidence for SN1993J- or GW170817-like variability of the radio flux over the course of our observations. A substantial host galaxy contribution to the measured radio flux is likely. The discovery and early-time peak capture of AT2019wxt in optical/NIR observation during EMGW follow-up observations highlights the need of dedicated early, multi-band photometric observations to identify USSNe.

Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe

Submitted to MNRAS, 16 pages, 10 figures

Lyman-$\alpha$ (Ly$\alpha$) forest transmission towards high-$z$ quasars contains information on the state of the universe just after the epoch of reionization. Fluctuations in Ly$\alpha$ forest transmission are partially sourced from spatial fluctuations in the ultraviolet background (UVB), where the level of UVB fluctuations are set by the mean free path of ionizing photons ($\lambda_{\text{mfp}}$). The auto-correlation function of Ly$\alpha$ forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to $\lambda_{\text{mfp}}$. Recent measurements at $z \sim 6$ suggest an unexpected rapid evolution of $\lambda_{\text{mfp}}$ at $z>5.0$ which would leave a signature in the evolution of the auto-correlation function. For this forecast, we model mock Ly$\alpha$ forest data with properties similar to the XQR-30 extended data set at $5.4 \leq z \leq 6.0$. Due to cosmic variance, we investigate 100 mock data sets at each $z$. Additionally we look at ideal cases, where mock data matches model values of the auto-correlation function. For such ideal data with $\lambda_{\text{mfp}}=9.0$ cMpc at $z=6.0$, we recover $\lambda_{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of $\lambda_{\text{mfp}}$ from the stacking of quasar spectra beyond the Lyman limit. Hypothetical data with resolution over three times greater than XQR-30 data leads to a $\sim40\%$ reduction in the error bars over all $z$. The distribution of mock values of the auto-correlation function in this work is highly non-Gaussian for high-$z$, which should caution work with other statistics of the high-$z$ Ly$\alpha$ forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements.

Scott W. McIntosh, Phillip H. Scherrer, Leif Svalgaard, Robert J. Leamon

9 pages, 3 figures, submitted to Frontiers in Solar and Stellar Physics

Through meticulous daily observation of the Sun's large-scale magnetic field the Wilcox Solar Observatory (WSO) has catalogued two magnetic (Hale) cycles of solar activity. Those two (~22-year long) Hale cycles have yielded four ($\sim$11-year long) sunspot cycles (numbers 21 through 24). Recent research has highlighted the persistence of the ``Extended Solar Cycle'' (ESC) and its connection to the fundamental Hale Cycle - albeit through a host of proxies resulting from image analysis of the solar photosphere, chromosphere and corona. This short manuscript presents the correspondence of the ESC, the surface toroidal magnetic field evolution, and the evolution of the Hale Cycle. As Sunspot Cycle 25 begins, interest in observationally mapping the Hale and Extended cycles could not be higher given potential predictive capability that synoptic scale observations can provide.

Matthias Y. He, Eric B. Ford

13 pages, 6 figures, 1 table. Submitted to AAS Journals

A foundational idea in the theory of in situ planet formation is the "minimum mass extrasolar nebula" (MMEN), a surface density profile ($\Sigma$) of disk solids that is necessary to form the planets in their present locations. While most previous studies have fit a single power-law to all exoplanets in an observed ensemble, it is unclear whether most exoplanetary systems form from a universal disk template. We use an advanced statistical model for the underlying architectures of multi-planet systems to reconstruct the MMEN. The simulated physical and Kepler-observed catalogs allows us to directly assess the role of detection biases, and in particular the effect of non-transiting or otherwise undetected planets, on altering the inferred MMEN. We find that fitting a power-law of the form $\Sigma = \Sigma_0^* (a/a_0)^\beta$ to each multi-planet system results in a broad distribution of disk profiles; $\Sigma_0^* = 336_{-291}^{+727}$ g/cm$^2$ and $\beta = -1.98_{-1.52}^{+1.55}$ encompasses the 16th-84th percentiles of the marginal distributions in an underlying population, where $\Sigma_0^*$ is the normalization at $a_0 = 0.3$ AU. Around half of inner planet-forming disks have minimum solid masses of $\gtrsim 40 M_\oplus$ within 1 AU. While transit observations do not tend to bias the median $\beta$, they can lead to both significantly over- and under-estimated $\Sigma_0^*$ and thus broaden the inferred distribution of disk masses. Nevertheless, detection biases cannot account for the full variance in the observed disk profiles; there is no universal MMEN if all planets formed in situ. The great diversity of solid disk profiles suggests that a substantial fraction of planetary systems experienced a history of migration.

Maria V. Pruzhinskaya, Emille E. O. Ishida, Alexandra K. Novinskaya, Etienne Russeil, Alina A. Volnova, Konstantin L. Malanchev, Matwey V. Kornilov, Patrick D. Aleo, Vladimir S. Korolev, Vadim V. Krushinsky, Sreevarsha Sreejith, Emmanuel Gangler

23 pages with appendix, 12 figures, 2 tables, submitted to MNRAS

In order to explore the potential of adaptive learning techniques to big data sets, the SNAD team used Active Anomaly Discovery (AAD) as a tool to search for new supernova (SN) candidates in the photometric data from the first 9.4 months of the Zwicky Transient Facility survey - between 2018 March 17 and December 31 (58194 < MJD < 58483). We analysed 70 ZTF fields with high galactic latitude and visually inspected 2100 outliers. This resulted in 104 supernova-like objects found, 57 of them were reported to the Transient Name Server for the first time and 47 were previously mentioned in other catalogues either as supernovae with known types or as supernova candidates. We visually inspected the multi-colour light curves of the non-catalogued transients and performed their fit with different supernova models to assign it to a proper class: Ia, Ib/c, IIP, IIL, IIn. Moreover, we also identified unreported slow-evolving transients which are good superluminous SN candidates, and a few others non-catalogued objects, such as red dwarf flares and active galactic nuclei. Beyond confirming the effectiveness of human-machine integration underlying the AAD strategy, our results shed light on potential leaks in currently available pipelines and can help avoid similar losses in future large scale astronomical surveys. The algorithm enables directed search of any type of data and definition of anomaly chosen by the expert.

Joachim Saur, Stefan Duling, Alexandre Wennmacher, Clarissa Willmes, Lorenz Roth, Darrell F. Strobel, Frédéric Allegrini, Fran Bagenal, Scott J. Bolton, Bertrand Bonfond, George Clark, Randy Gladstone, T.K. Greathouse, Denis C. Grodent, Candice J. Hansen, W.S. Kurth, Glenn S. Orton, Kurt D. Retherford, Abigail M. Rymer, Ali H. Sulaiman

Accepted for publication with Geophys. Res. Let

We report results of Hubble Space Telescope observations from Ganymede's orbitally trailing side which were taken around the flyby of the Juno spacecraft on June 7, 2021. We find that Ganymede's northern and southern auroral ovals alternate in brightness such that the oval facing Jupiter's magnetospheric plasma sheet is brighter than the other one. This suggests that the generator that powers Ganymede's aurora is the momentum of the Jovian plasma sheet north and south of Ganymede's magnetosphere. Magnetic coupling of Ganymede to the plasma sheet above and below the moon causes asymmetric magnetic stresses and electromagnetic energy fluxes ultimately powering the auroral acceleration process. No clear statistically significant time variability of the auroral emission on short time scales of 100s could be resolved. We show that electron energy fluxes of several tens of mW m$^{-2}$ are required for its OI 1356 \AA$\;$ emission making Ganymede a very poor auroral emitter.

Norbert Magyar, Tom Van Doorsselaere

Accepted in ApJ

The origin and evolution of the 1/f power law observed in the energy spectrum of solar coronal and solar wind fluctuations at scales of around an hour is not entirely understood. Several existing theories aim at explaining it, involving both linear and nonlinear mechanisms. An often overlooked property of the solar corona and solar wind is their highly inhomogeneous nature. In this paper we investigate the linear evolution of pure Alfv\'en and surface Alfv\'en waves propagating through a plasma which is inhomogeneous across the magnetic field. The inhomogeneity is given by density, which we model to be two-dimensional colored noise, with power spectral slopes ranging from -2 to -1. Alfv\'en waves propagate independently on individual magnetic field lines, and eventually get completely out of phase through the process of phase mixing, leading to unrealistic spectra. When the coupling between the inhomogeneous background and the propagating waves is fully accounted for, transverse waves such as surface Alfv\'en waves (also referred to as kink or Alfv\'enic) appear, showing collective wave behavior of neighboring magnetic field lines with different Alfv\'en speeds. We show that the linear cascade of surface Alfv\'en wave energy, induced by phase mixing and resonant absorption, leads to a perpendicular wave energy spectrum which tends to the perpendicular power spectrum of the background density. Based on our model, we propose that a perpendicular density power spectrum of 1/f in the solar corona can induce, through linear processes, the 1/f spectrum of the fluctuations that is observed at the largest scales.

Franklin Wang, Jian Ge, Kevin Willis

14 pages, 22 Figures, 4 Tables; To be published in the Monthly Notices of the Royal Astronomical Society (MNRAS)

Although many near-Earth objects have been found by ground-based telescopes, some fast-moving ones, especially those near detection limits, have been missed by observatories. We developed a convolutional neural network for detecting faint fast-moving near-Earth objects. It was trained with artificial streaks generated from simulations and was able to find these asteroid streaks with an accuracy of 98.7% and a false positive rate of 0.02% on simulated data. This program was used to search image data from the Zwicky Transient Facility (ZTF) in four nights in 2019, and it identified six previously undiscovered asteroids. The visual magnitudes of our detections range from ~19.0 - 20.3 and motion rates range from ~6.8 - 24 deg/day, which is very faint compared to other ZTF detections moving at similar motion rates. Our asteroids are also ~1 - 51 m diameter in size and ~5 - 60 lunar distances away at close approach, assuming their albedo values follow the albedo distribution function of known asteroids. The use of a purely simulated dataset to train our model enables the program to gain sensitivity in detecting faint and fast-moving objects while still being able to recover nearly all discoveries made by previously designed neural networks which used real detections to train neural networks. Our approach can be adopted by any observatory for detecting fast-moving asteroid streaks.

The multi-messenger data of neutron star merger events are promising for constraining the Hubble constant. So far, GW170817 is still the unique gravitational wave event with multi-wavelength electromagnetic counterparts. In particular, its radio and X-ray emission have been measured in the past $3-5$ years. In this work, we fit the X-ray, optical, and radio afterglow emission of GW170817/GRB 170817A and find out that a relatively large viewing angle $\sim 0.5\, \rm rad$ is needed; otherwise, the late time afterglow data can not be well reproduced. Such a viewing angle has been taken as a prior in the gravitational wave data analysis, and the degeneracy between the viewing angle and the luminosity distance is broken. Finally, we have a Hubble constant $H_0=72.00^{+4.05}_{-4.13}\, \rm km\, s^{-1}\, Mpc^{-1}$, which is more consistent with that obtained by other local measurements. If rather similar values are inferred from multi-messenger data of future neutron star merger events, it will provide critical support to the existence of the Hubble tension.

In this paper, we introduce our open source implementation of the Coupled Discrete Unified Gas Kinetic Scheme (CDUGKS) of https://journals.aps.org/pre/abstract/10.1103/PhysRevE.98.053310, a phase space scheme capable of handling a wide range of flow regimes. We demonstrate its performance on several problems including a number of well known test problems from the astrophysical fluid dynamics literature such as the 1D Sod shock tube, 2D Kelvin-Helmholtz instability, 1D thermoacoustic wave, a triangular Gresho vortex, a sine wave velocity perturbation. For these problems, we show that the code can simulate flows ranging from the inviscid/Eulerian regime to the free-streaming regime, capturing shocks and emergent diffusive processes in the appropriate regimes. We also use a variety of Prandtl numbers to demonstrate the scheme's ability to simulate different thermal conductivities at fixed viscosity. The scheme is second-order accurate in space and time and, unlike many solvers, uses a time step that is independent of the mean free path of the gas. Our code (MP-CDUGKS) is public under a CC0 1.0 Universal license and is available on https://github.com/alvarozamora/CDUGKS

Hui-Fang Xue, Jia-Shu Niu, Jian-Ning Fu

26 pages, 5 figures, 3 tables

Stellar structure and evolution theory is one of the basis in modern astronomy. Stellar inner structures and their evolutionary states can be precisely tested by asteroseismology, since the inner information is brought to the stellar surface by the global oscillating waves and becomes observable. For stellar evolutionary speed (i.e. how long time scale does a star stay at a special evolution phase?), because of the insurmountable gap between the time scales of the evolutionary history of human civilization and a star, it can only be roughly tested by ensemble of stars in different evolutionary stages in most cases, and all the snapshots of these stars make up our global view of stellar evolution. The effect of stellar evolution on the structure and the corresponding global size of a pulsating star will lead to tiny period variations of its pulsation modes, which are the most valuable indicators of its evolutionary state and can be used to test the stellar evolution theory by a single star rather than ensemble of stars. Here, we report a High-Amplitude $\delta$ Scuti star AE Ursae Majoris, who locates in the post main-sequence (MS) evolutionary stage and its observed linear period variation rate can be practically ascribed to its evolutionary effect. The result tests the stellar evolution theory from the pre-MS to post-MS with an unprecedented precision by a single star, and the framework can be extended to other type of pulsating stars to perform precise evolutionary asteroseismology, which aims to test the current stellar evolution theory in different evolutionary stages, discover the discrepancies between the theory and observations, and ultimately build a complete and precise stellar evolution theory to backtrack the history of each of these stars.

A Marciniak (IRAP, LCAR), A Bonnamy (IRAP), S Zamith (LCAR), L Noguès (IRAP), O Coeur-Joly (IRAP), D Murat (IRAP), P Moretto-Capelle (LCAR), J.-M l'Hermite (LCAR), G Mulas (OAC, IRAP), C Joblin (IRAP)

Book Series: Astrophysics and Space Science Proceedings (ASSSP)

PIRENEA and PIRENEA 2 are experimental setups dedicated to the study of fundamental molecular processes involving species of astrochemical interest. The coupling of a VUV source to PIRENEA has allowed us to study the fragmentation pathways and stability of polycyclic aromatic hydrocarbons (PAHs) containing aliphatic bonds under conditions relevant for astrophysical photodissociation regions. PIRE-NEA 2 will open the possibility to extend these studies to larger systems such as PAH clusters, and more generally to study gas-nanograin-photon interactions.

P. Petit, T. Böhm, C.P. Folsom, F. Lignières, T. Cang

A&A, in press

Aims. The very weak magnetic field detected at the surface of Vega hints at a widespread population of weakly magnetic stars of A and B spectral types. We contribute here to gather more clues about the origin of this magnetism by investigating the long-term stability of the field geometry of this prototypical star. Methods. We use spectropolarimetric data collected as part of a long-term campaign, with more than 2,000 observations spread between 2008 and 2018. Using various sub-sets extracted from the whole time series, we reconstruct several maps of the large-scale surface magnetic field. Results. We confirm that the polarimetric signal is modulated according to a $\sim 0.68$ d period, which we interpret as the stellar rotation period. The surface magnetic field is organized in a complex geometry. We confirm the existence of a very localized, polar magnetic spot previously reported for Vega, with a radial field strength of about -5 G. We show that the surface of the star is also covered by a dipole, with a polar strength close to 9 G and a dipole obliquity close to $90^\circ$. Both magnetic structures are remarkably stable over one decade. The available data suggest that smaller-scale magnetic spots may not be limited to the polar region, although the poor reliability of their reconstruction does not allow us to firmly conclude about their temporal evolution.

Roldán A. Cala, José F. Gómez, Luis F. Miranda, Lucero Uscanga, Shari L. Breen, Joanne R. Dawson, Itziar de Gregorio-Monsalvo, Hiroshi Imai, Hai-Hua Qiao, Olga Suárez

17 pages, 7 figures, 5 tables. Accepted for publication in MNRAS

The evolution of asymptotic giant branch stars from the spherical symmetry into the diverse shapes of planetary nebulae (PNe) is a topic of intensive research. Young PNe provide a unique opportunity to characterize the onset of this transitional phase. In particular, OH maser-emitting PNe (OHPNe) are considered nascent PNe. In fact, only 6 OHPNe have been confirmed to date. In order to identify and characterize more OHPNe, we processed the unpublished continuum data of the interferometric follow-up of the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH). We then matched the interferometric positions of OH maser and radio continuum emission, considering the latter as a possible tracer of free-free emission from photoionized gas, characteristic of PNe. We report 8 objects with a positive coincidence, 4 of which are classified as candidate OHPNe here for the first time (IRAS 16372-4808, IRAS 17494-2645, IRAS 18019-2216 and OH 341.6811+00.2634). Available evidence strongly indicates that they are evolved stars, while the comparison with confirmed OHPNe indicates that they are likely to be PNe. Their final confirmation as bona fide PNe, however, requires optical/infrared spectroscopy. The obtained spectral indices of the radio continuum emission (between $\simeq$ 0.4 - 1.3) are consistent with partially optically thick free-free emission from photoionized gas. Also, they cluster in the same region of a WISE colour-colour diagram as that of the confirmed OHPNe ($9.5 \lesssim [3.4]-[22] \lesssim 13.5$, and $4.0 \lesssim [4.6]-[12] \lesssim 7.0$), thus this diagram could help to identify more OHPNe candidates in the future.

C.-H. Dahlqvist, J. Milli, O. Absil, F. Cantalloube, L. Matra, E. Choquet, C. del Burgo, J. P. Marshall, M. Wyatt, S. Ertel

In the past decade, HCI surveys provided new insights about the frequency and properties of substellar companions at separation larger than 5 au. In this context, our study aims to detect and characterise potential exoplanets and brown dwarfs within debris disks, by considering the SHARDDS survey, which gathers 55 Main Sequence stars with known bright debris disk. We rely on the AutoRSM framework to perform an in-depth analysis of the targets, via the computation of detection maps and contrast curves. A clustering approach is used to divide the set of targets in multiple subsets, in order to reduce the computation time by estimating a single optimal parametrisation for each considered subset. The use of Auto-RSM allows to reach high contrast at short separations, with a median contrast of 10-5 at 300 mas, for a completeness level of 95%. Detection maps generated with different approaches are used along with contrast curves, to identify potential planetary companions. A new planetary characterisation algorithm, based on the RSM framework, is developed and tested successfully, showing a higher astrometric and photometric precision for faint sources compared to standard approaches. Apart from the already known companion of HD206893 and two point-like sources around HD114082 which are most likely background stars, we did not detect any new companion around other stars. A correlation study between achievable contrasts and parameters characterising HCI sequences highlights the importance of the strehl, wind speed and wind driven halo to define the quality of high contrast images. Finally, planet detection and occurrence frequency maps are generated and show, for the SHARDDS survey, a high detection rate between 10 and 100 au for substellar companions with mass >10MJ.

A broadband two-layer anti-reflection (AR) coating was developed for use on a sapphire half-wave plate (HWP) and an alumina infrared (IR) filter for cosmic microwave background (CMB) polarimetry. Measuring tiny CMB B-mode signals requires maximizing the number of photons reaching the detectors and minimizing spurious polarization due to reflection with an off-axis incident angle. However, a sapphire HWP and an alumina IR filter have high refractive indices of about 3.1, and an AR coating must be applied to them. Thermally sprayed mullite and Duroid 5880LZ were selected in terms of index and coefficient of thermal expansion for use at cryogenic temperatures. With these materials, the reflectivity was reduced to about 2% at 90/150 GHz and <1% at 220/280 GHz. The design, fabrication, and optical performance evaluation of the AR coatings are described. The coatings were used in a current ground-based CMB experiment called the Simons Array. They could also be applied to next-generation CMB experiments, such as the Simons Observatory.

Kunio Kaneta, Hye-Sung Lee, Jiheon Lee, Jaeok Yi

Despite its dominance in the present universe's energy budget, dark energy is the least understood component in the universe. Although there is a popular model for the dynamical dark energy, the quintessence scalar, the investigation is limited because of its highly elusive character. We present a model where the quintessence is gauged by an Abelian gauge symmetry. The quintessence is promoted to be a complex scalar whose real part is the dark energy field while the imaginary part is the longitudinal component of a new gauge boson. It brings interesting characters to dark energy physics. We study the general features of the model, including how the quintessence behavior is affected and how its gauge interaction is constrained by the solicited dark energy properties.

Yao Liu, Gesa H.-M. Bertrang, Mario Flock, Giovanni P. Rosotti, Ewine F. van Dishoeck, Yann Boehler, Stefano Facchini, Can Cui, Sebastian Wolf, Min Fang

13 pages, 7 figures, accepted for publication in Science China Physics, Mechanics & Astronomy

Turbulent motions are believed to regulate angular momentum transport and influence dust evolution in protoplanetary disks. Measuring the strength of turbulence is challenging through gas line observations because of the requirement for high spatial and spectral resolution data, and an exquisite determination of the temperature. In this work, taking the well-known HD 163296 disk as an example, we investigated the contrast of gaps identified in high angular resolution continuum images as a probe for the level of turbulence. With self-consistent radiative transfer models, we simultaneously analyzed the radial brightness profiles along the disk major and minor axes, and the azimuthal brightness profiles of the B67 and B100 rings. By fitting all the gap contrasts measured from these profiles, we constrained the gas-to-dust scale height ratio $\Lambda$ to be $3.0_{-0.8}^{+0.3}$, $1.2_{-0.1}^{+0.1}$ and ${\ge}\,6.5$ for the D48, B67 and B100 regions, respectively. The varying gas-to-dust scale height ratios indicate that the degree of dust settling changes with radius. The inferred values for $\Lambda$ translate into a turbulence level of $\alpha_{\rm turb}\,{<}\,3\times10^{-3}$ in the D48 and B100 regions, which is consistent with previous upper limits set by gas line observations. However, turbulent motions in the B67 ring are strong with $\alpha_{\rm turb}\,{\sim}1.2\,{\times}\,10^{-2}$. Due to the degeneracy between $\Lambda$ and the depth of dust surface density drops, the turbulence strength in the D86 gap region is not constrained.

Afterglows of gamma-ray bursts (GRBs) are emitted from expanding forward shocks, which are expected to have magnetic field much stronger than the interstellar field, although the origin of the field is a long-standing problem. Two field amplification mechanisms, plasma kinetic instabilities and magnetohydrodynamic instabilities, have been discussed so far. The coherence length scales of the fields amplified by these two processes are different by 7-10 orders of magnitudes, and the polarimetric observations may distinguish them. We construct a semi-analytic model of the forward shock afterglow polarization under the assumption of hydrodynamic-scale turbulent magnetic field. We perform numerical calculations of synchrotron polarization for the isotropic turbulence and the zero viewing angle. We find that the polarization degrees are ~1-3% when the field coherence length scale in the fluid comoving frame is comparable to the thickness of the shocked regions. This range of polarization degree is comparable to that of the observed late-phase optical afterglows. Our model also shows that the radio polarization degrees are comparable to the optical ones on average but can be higher than the optical ones at some time intervals. The polarization angles are shown to vary randomly and continuously. These polarimetric properties are clearly different from the case of plasma kinetic instability. Simultaneous polarimetric observations of GRB afterglows at the radio and optical bands have recently started, which will help us constrain the magnetic field amplification mechanism.

V. S. Dhillon, M. R. Kennedy, R. P. Breton, C. J. Clark, D. Mata Sánchez, G. Voisin, E. Breedt, A. J. Brown, M. J. Dyer, M. J. Green, P. Kerry, S. P. Littlefair, T. R. Marsh, S. G. Parsons, I. Pelisoli, D. I. Sahman, J. F. Wild, M. H. van Kerkwijk, B. W. Stappers

10 pages, 4 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light curves, however, do not show this asymmetry, demonstrating that whatever mechanism is responsible for it varies on timescales of a decade or less. We fit the symmetrical HiPERCAM light curves with a direct-heating model to derive the system parameters, finding an orbital inclination of $55.9^{+4.8}_{-4.1}$ degrees, in good agreement with radio-eclipse constraints. We find that approximately half of the pulsar's spin-down energy is converted to optical luminosity, resulting in temperatures ranging from approximately $5150^{+190}_{-190}$ K on the day side to $2750^{+130}_{-150}$ K on the night side of the companion star. The companion star is close to filling its Roche lobe ($f_{\rm RL} =0.88^{+0.02}_{-0.02}$) and has a mass of $0.039^{+0.010}_{-0.011}$ M$_{\odot}$, giving a mean density of $20.24^{+0.59}_{-0.44}$ g cm$^{-3}$ and an apsidal motion constant in the range $0.0036 < k_2 < 0.0047$. The companion mass and mean density values are consistent with those of brown dwarfs, but the apsidal motion constant implies a significantly more centrally-condensed internal structure than is typical for such objects.

X. D. Jia, J. P. Hu, J.Yang, B. B. Zhang, F. Y. Wang (NJU)

14 pages, 7 figures, 3 tables, accepted for publication in MNRAS

Gamma-ray bursts (GRBs) are the most explosive phenomena and can be used to study the expansion of Universe. In this paper, we compile a long GRB sample for the $E_{\mathrm{iso}}$-$E_{\mathrm{p}}$ correlation from Swift and Fermi observations. The sample contains 221 long GRBs with redshifts from 0.03 to 8.20. From the analysis of data in different redshift intervals, we find no statistically significant evidence for the redshift evolution of this correlation. Then we calibrate the correlation in six sub-samples and use the calibrated one to constrain cosmological parameters. Employing a piece-wise approach, we study the redshift evolution of dark energy equation of state (EOS), and find that the EOS tends to be oscillating at low redshift, but consistent with $-1$ at high redshift. It hints a dynamical dark energy at $2\sigma$ confidence level at low redshift.

The Optical-to-Near-infrared variability time delay have already been reported for a small number (about 7) of AGNs and has been firmly established only for 5 of them. The time delay is probably increasing with the IR wavelengths. The most naturally this time delay can be interpreted by the model where IR emission is attributed to circumnuclear dust heated by the nuclear radiation. In given model a suggestion on narrowness of the near-infrared (NIR) emission region is quite natural, as far as the dust can be not saved on distances from the nucleus closer then some critical value, on which it is reached the sublimation temperature for graphite particles (Barvainis, 1987). For NGC 4151 case it has been shown that the NIR region has a form of thin ring or torus. The radius of this ring correlates with level of the nucleus activity (Oknyanskij et al. 1999). This dependency of radius of the NIR emission region from luminosity reveals itself as under object variability (as in the case of NGC4151), and also when objects with high and low luminosity are considered. We assume that the observed time delays allow us to derive a redshift independent luminosity distances to AGNs and estimate a Hubble constant. Some problems of using this strategy for the Hubble constant determination are discussed.

It is attempted to derive the general relativistic (GR) equation of motion for planet and its solution solely by the special relativity (SR) techniques. The motion of a planet relative to the sun and that of the sun to the planet are solved independently in special relativistic framework using the perturbation theory in the celestial mechanics. The solution reveals a nature of the structure of the spacetime under the gravitation of the sun, and then its effect on the planet's motion is examined. When the motion thus examined are compared with the one obtained by the general relativity theory in PN approximation, both are different concerning the mean motion and the radius of the orbit but exactly the same as for the perihelion precession.

R. Guadarrama, Eduard I. Vorobyov, Christian Rab, Manuel Güdel

19 pages, 18 figures

We study the influence of different metallicities on the physical, thermal, and chemical properties of protoplanetary disks, and in particular on the formation and destruction of carbon-based molecules. With the thermo-chemical code ProoDIMO we investigate the impact of lower metallicities on the radiation field, disk temperature, and the abundance of different molecules (H$_2$O, CH$_4$, CO, CO$_2$, HCN, CN, HCO$^+$ and N$_2$H$^+$). We use a fiducial disk model as a reference model and produce two models with lower metallicity. The resulting influence on different chemical species is studied by analyzing their abundance distribution throughout the disk and their vertical column density. Furthermore, the formation and destruction reactions of the chemical species are studied. The results show a relation between the metallicity of the disk and the strength of the stellar radiation field inside the disk. As the metallicity decreases the radiation field is able to penetrate deeper regions of the disk. As a result, there is a stronger radiation field overall in the disk with lower metallicity which also heats up the disk. This triggers a series of changes in the chemical formation and destruction efficiencies for different chemical species. In most cases, the available species abundances change and have greater values compared to scaled-down abundances by constant factors. Metallicity has a clear impact on the snowline of the molecules studied here as well. As metallicity decreases the snowlines are pushed further out and existing snow rings shrink in size.

Erik Høg (Niels Bohr Institute - Copenhagen)

62 pages (=1+22+19+20, including abstracts and references)

Michael Perryman has interviewed some of the scientists and project leaders in the Hipparcos and Gaia missions, the interviews with photos of the persons are given at his site: https://www.michaelperryman.co.uk . Michael has also written essays -- 84 to date ! -- about results from the Gaia mission and they are placed at his site. Three of the interviews are with me and transcriptions, co-authored with Michael, are provided below with the titles: #1. An interview about astronomy and astrometry up to 1980. #2. An interview about the revival of astrometry after 1980. #3. The billion-star astrometry after 1990. The third interview begins in 1990 when I had the first ideas for a Hipparcos successor. In 1992 I made a detailed design with direct imaging on CCD detectors in a satellite proposal called Roemer. In 1993 a supposedly better option was proposed with the acronym GAIA where the capital "I" stood for Interferometer. In 1998, however, interferometry was shown to be unsuited for the purpose and we returned to the original idea from 1992 for the further development. The name was later changed to Gaia -- for the sake of continuity.

Tristan Cantat-Gaudin, Morgan Fouesneau, Hans-Walter Rix, Anthony G. A. Brown, Alfred Castro-Ginard, Ronald Drimmel, David W. Hogg, Andrew R. Casey, Shourya Khanna, Semyeong Oh, Adrian M. Price Whelan, Vasily Belokurov, Andrew K. Saydjari, Gregory M. Green

to be submitted to A&A

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky. On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project.

A. Evans (Keele University, UK), T. R. Geballe (Gemini Observatory/NSF's NOIRLab), C. E. Woodward (Minnesota Institute for Astrophysics), D. P. K. Banerjee (Physical Research Laboratory, Ahmedabad, India), R. D. Gehrz (Minnesota Institute for Astrophysics), S. Starrfield (Arizona State University), M. Shahbandeh (Florida State University)

Accepted for publication in MNRAS on 2022 August 1. 19 pages, 10 figures

We present infrared spectroscopy of the 2019 eruption of the recurrent nova V3890 Sgr, obtained over the period 5.1-46.3 days after the eruption. The spectrum of the red giant became more prominent as the flux declined, and by day 46.3 dominated the spectrum. Hydrogen and helium emission lines consisted of a narrow component superposed on a broad pedestal. The full width at half maximum of the narrow components declined with time $t$ as the eruption progressed, as ${t}^{-0.74}$, whereas those of the broad components remained essentially constant. Conversely, the line fluxes of the narrow components of Pa,$\beta$ remained roughly constant, while those of the broad components declined by a factor $\sim30$ over a period of $\lesssim25$~days. The behaviour of the broad components is consistent with them arising in unencumbered fast-flowing ejecta perpendicular to the binary plane, in material that was ejected in a short $\sim3.3$-day burst. The narrow components arise in material that encounters the accumulated circumstellar material. The outburst spectra were rich in coronal lines. There were two coronal line phases, one that originated in gas ionised by supersoft X-ray source, the other in shocked gas. From the relative fluxes of silicon and sulphur coronal lines on day 23.4 - when the emitting gas was shocked - we deduce that the temperature of the coronal gas was $9.3\times10^5$~K, and that the abundances are approximately solar.

Omar Contigiani, Henk Hoekstra, Margot M. Brouwer, Andrej Dvornik, Maria Cristina Fortuna, Cristóbal Sifón, Ziang Yan, Mohammadjavad Vakili

The masses of galaxy clusters can be measured using data obtained exclusively from wide photometric surveys in one of two ways: directly from the amplitude of the weak lensing signal or, indirectly, through the use of scaling relations calibrated using binned lensing measurements. In this paper, we build on a recently proposed idea and implement an alternative method based on the radial profile of the satellite distribution. This technique relies on splashback, a feature associated with the apocenter of recently accreted galaxies that offers a clear window into the phase-space structure of clusters without the use of velocity information. We carry out this dynamical measurement using the stacked satellite distribution around a sample of luminous red galaxies in the fourth data release of the Kilo-Degree Survey and validate our results using abundance-matching and lensing masses. To illustrate the power of this measurement, we combine dynamical and lensing mass estimates to robustly constrain scalar-tensor theories of gravity at cluster scales. Our results exclude departures from General Relativity of order unity. We conclude the paper by discussing the implications for future data sets. Because splashback mass measurements scale only with the survey volume, stage-IV photometric surveys are well-positioned to use splashback to provide high-redshift cluster masses.

P. Mas-Buitrago, E. Solano, A. González-Marcos, C. Rodrigo, E. L. Martín, J. A. Caballero, F. Jiménez-Esteban, P. Cruz, A. Ederoclite, J. Ordieres-Meré, A. Bello-García, R. A. Dupke, A. J. Cenarro, D. Cristóbal-Hornillos, C. Hernández-Monteagudo, C. López-Sanjuan, A. Marín-Franch, M. Moles, J. Varela, H. Vázquez Ramió, J. Alcaniz, L. Sodré Jr., R. E. Angulo

Accepted in A&A

Ultracool dwarfs (UCDs) comprise the lowest mass members of the stellar population and brown dwarfs, from M7 V to cooler objects with L, T, and Y spectral types. Most of them have been discovered using wide-field imaging surveys, for which the Virtual Observatory (VO) has proven to be of great utility. We aim to perform a search for UCDs in the entire Javalambre Photometric Local Universe Survey (J-PLUS) second data release (2176 deg$^2$) following a VO methodology. We also explore the ability to reproduce this search with a purely machine learning (ML)-based methodology that relies solely on J-PLUS photometry. We followed three different approaches based on parallaxes, proper motions, and colours, respectively, using the VOSA tool to estimate the effective temperatures. For the ML methodology, we built a two-step method based on principal component analysis and support vector machine algorithms. We identified a total of 7827 new candidate UCDs, which represents an increase of about 135% in the number of UCDs reported in the sky coverage of the J-PLUS second data release. Among the candidate UCDs, we found 122 possible unresolved binary systems, 78 wide multiple systems, and 48 objects with a high Bayesian probability of belonging to a young association. We also identified four objects with strong excess in the filter corresponding to the Ca II H and K emission lines and four other objects with excess emission in the H$\alpha$ filter. With the ML approach, we obtained a recall score of 92% and 91% in the test and blind test, respectively. We consolidated the proposed search methodology for UCDs, which will be used in deeper and larger upcoming surveys such as J-PAS and Euclid. We concluded that the ML methodology is more efficient in the sense that it allows for a larger number of true negatives to be discarded prior to analysis with VOSA, although it is more photometrically restrictive.

Josh Hoffmann, David Sloan

24 pages, 25 figures

There are many single field inflationary models that are consistent with the recent Planck 2018 measurements of the spectral index $n_s$ and tensor-to-scalar ratio $r$. Despite good agreement with observational data some of these models suffer from having unregularized potentials which would produce a collapsing universe shortly after the end of inflation. In this paper we show that how one chooses to correct the behaviour potential towards the end of inflation can have a significant effect on the inflationary predictions of the model, specifically in the case of quartic hilltop and radiatively corrected Higgs inflation.

D. Alvarez-Castillo

7 pages, 1 figure. Presented at the International Conference on Precision Physics and Fundamental Physical Constants 2021

The study of the variation of fundamental constants through time or in localized regions of space is one of the goals of the Cosmic Ray Extremely Distributed Observatory which consists of multiple detectors over the Earth. In this paper, the various effects which can be potentially identified through cosmic rays detections by CREDO are presented.

Filip Huško (1), Cedric G. Lacey (1), ((1) - ICC, Durham)

Submitted to MNRAS

We use SWIFT, a smoothed particle hydrodynamics code, to simulate the evolution of bubbles inflated by active galactic nuclei (AGN) jets, as well as their interactions with the ambient intracluster medium (ICM). These jets inflate lobes that turn into bubbles after the jets are turned off (at $t=50$ Myr). Almost all of the energy injected into the jets is transferred to the ICM very quickly after they are turned off, with roughly $70$ per cent of it in thermal form and the rest in kinetic. At late times ($t>500$ Myr) we find the following: 1) the bubbles draw out trailing filaments of low-entropy gas, similar to those recently observed, 2) the action of buoyancy and the uplift of the filaments dominates the energetics of both the bubbles and the ICM and 3) almost all of the originally injected energy is in the form of gravitational potential energy, with the bubbles containing $15$ per cent of it, and the rest contained in the ICM. These findings indicate that feedback proceeds mainly through the displacement of gas to larger radii. We find that the uplift of these filaments permanently changes the thermodynamic properties of the ICM by reducing the central density and increasing the central temperature (within $30$ kpc). We propose that jet feedback proceeds not only through the heating of the ICM (which can delay cooling), but also through the uplift-related reduction of the central gas density. The latter also delays cooling, on top of reducing the amount of gas available to cool.

We have conducted a search for open clusters in the vicinity of Galactic classical Cepheids based on high-quality astrometry from the ESA mission Gaia's third data release (DR3) to improve the calibration of the Leavitt law. Our approach requires no prior knowledge of existing clusters, allowing us to both detect new host clusters and cross-check previously reported associations. Our Gold sample consists of 34 Cepheids residing in 28 open clusters, including 27 fundamental mode and 7 overtone Cepheids. Three new bona fide cluster Cepheids are reported (V0378 Cen, ST Tau and GH Lup) and the host cluster identifications for VW Cru, IQ Nor and SX Vel is corrected. The fraction of Cepheids occurring in open clusters within 2 kpc of the Sun is $f_{\mathrm{CC,2kpc}}=0.092 \pm 0.009$. Combining cluster and field Cepheids, we calibrate the Galactic Cepheid Leavitt law for a variety of individual optical and near-infrared passbands, as well as using reddening-free Wesenheit magnitudes based on Gaia and HST photometry, while solving for the residual offset applicable to Cepheid parallaxes, $\Delta \varpi_{\text{Cep}}$. The most direct comparison of our results with the SH0ES distance ladder yields excellent ($0.3\sigma$) agreement for both the absolute magnitude of a 10 d Solar metallicity Cepheid in the NIR HST Wesenheit magnitudes, $M_H^W=-5.915\pm0.017$ mag, and the residual parallax offset, $\varpi_{\text{Cl}}=-13\pm5\,\mu$as. Using 26 cluster Cepheids and 225 MW Cepheids with Gaia DR3 astrometry and photometry, we determine the absolute magnitude of a 10 d fundamental mode Cepheid in the Gaia Wesenheit magnitude at Solar metallicity $M_G^W =-5.996\pm 0.019$ mag and $\Delta \varpi_{\text{Cep}}=-19 \pm 3\,\mu$as. These results mark the currently most accurate absolute calibrations of the Cepheid luminosity scale based purely on observations of MW Cepheids.

Thea Kozakis, João M. Mendonça, Lars A. Buchhave

19 pages, 7 figures, 4 tables; in press with A & A

Molecular oxygen (O2) paired with a reducing gas is regarded as a promising biosignature pair for the atmospheric characterization of terrestrial exoplanets. In circumstances when O2 may not be detectable in a planetary atmosphere (e.g., at mid-IR wavelengths) it has been suggested that ozone (O3), the photochemical product of O2, could be used as a proxy to infer the presence of O2. However, O3 production has a nonlinear dependence on O2 and is strongly influenced by the UV spectrum of the host star. To evaluate the reliability of O3 as a proxy for O2, we used Atmos, a 1D coupled climate/photochemistry code, to study the O2-O3 relationship for "Earth-like" habitable zone planets around a variety of stellar hosts (G0V-M5V) and O2 abundances. Overall, we found that the O2-O3 relationship differed significantly with stellar hosts and resulted in different trends for hotter stars (G0V-K2V) vs cooler stars (K5V-M5V). Planets orbiting hotter host stars counter-intuitively experience an increase in O3 when O2 levels are initially decreased from 100% Earth's present atmospheric level (PAL), with a maximum O3 abundance occurring at 25-55% PAL O2. As O2 abundance initially decreases, larger amounts of UV photons capable of O2 photolysis reach the lower (denser) regions of the atmosphere where O3 production is more efficient, resulting in these increased O3 levels. This effect does not occur for cooler host stars (K5V-M5V), as the weaker incident UV flux does not allow O3 formation to occur at dense enough regions of the atmosphere where the faster O3 production can outweigh a smaller source of O2 from which to create O3. Overall it will be extremely difficult (or impossible) to infer precise O2 levels from an O3 measurement, however, with information about the UV spectrum of the host star and context clues, O3 will provide valuable information about potential surface habitability of an exoplanet.

Juan Carlos Hidalgo, Luis E. Padilla, Gabriel German

11 pages, 5 figures. Comments are welcome

At times prior to Big Bang Nucleosynthesis, the universe could show a primordial structure formation period if dominated by a fast oscillating inflaton field during reheating. In this context, we have postulated a new mechanism of primordial black hole formation [L. E. Padilla, J. C. Hidalgo, and K. A. Malik, Phys. Rev. D, vol. 106, p. 023519, Jul 2022], that draws the analogy between an extended reheating era and the scalar field dark matter model, stipulating the gravitational collapse of inflaton halos and inflaton stars. In this paper we look at the requirements for the realization of this new mechanism. We show that a generic primordial power spectrum with a peak at small scales is most suitable for the production of a considerable number of PBHs. When such requirement is met, and if reheating lasts long enough, large populations of PBHs with $M_{\rm PBH}\sim 1~\mathrm{gram}$ may be produced. We find in particular, that the mass fraction of PBHs is orders of magnitude larger than that obtained when PBHs form via direct collapse in a universe dominated by radiation or pressure-less dust. Looking at observable implications of our findings, we explore the possibility that the PBHs component may dominate the energy density of the universe at some point after the end of reheating.

Aleks Scholz (St Andrews), Koraljka Muzic (Lisbon), Ray Jayawardhana (Cornell), Lyra Quinlan, James Wurster (St Andrews)

16 pages, 3 figures, 2 tables, preprint, submitted to PASP, comments welcome

Free-floating (or rogue) planets are planets that are liberated (or ejected) from their host systems. Although simulations predict their existence in substantial numbers, direct observational evidence for free-floating planets with masses below 5 Mjup is still lacking. Several cycle-1 observing programs with JWST aim to hunt for them in four different star-forming clusters. These surveys are designed to be sensitive to masses of 1-15 Mjup (assuming a hot-start formation), which corresponds to spectral types of early L to late T for the ages of these clusters. If the existing simulations are not widely off the mark, we show here that the planned programs are likely to find up to 10-20 giant rogue planets in moderate density clusters like NGC1333 or IC348, and several dozen to about 100 in high-density regions like NGC2024 and the Orion Nebula Cluster. This number is expected to be modulated with the environment, and could be substantially higher if stars form multiple giant planets at birth. In contrast, the number of free-floating brown dwarfs, formed from core collapse (i.e., like stars) is expected to be significantly lower, only about 0.25% of the number of stars, or 1-7 for the clusters considered here. Below 10 Mjup that number drops further by an order of magnitude. We also show that the planned surveys are not at risk of being significantly contaminated by field brown dwarfs in the foreground or background, after spectroscopic confirmation. Taken together, our results imply that if a population of L and T dwarfs were to be found in these JWST surveys, it is expected to be predominantly made up of rogue planets.

Giacomo Cordoni, Antonino P. Milone, Anna F. Marino, Michele Cignoni, Edoardo P. Lagioia, Marco Tailo, Marília Carlos, Emanuele Dondoglio, Sohee Jang, Anjana Mohandasan, Maria V. Legnardi

9 pages, 5 figures, published on Nature Communications DOI this https URL

The origin of young star clusters represents a major challenge for modern stellar astrophysics. While stellar rotation partially explains the colour spread observed along main-sequence turn-offs, i.e. where stars leave the main-sequence after the exhaustion of hydrogen in their core, and the multiple main sequences in the colour-magnitude diagrams of stellar systems younger than approximately 2 Gyr, it appears that an age difference may still be required to fulfill the observational constraints. Here we introduce an alternative approach that exploits the main-sequence turn-on, i.e. the point alongside the colour-magnitude diagram where pre-main-sequence stars join the main-sequence, to disentangle between the effects of stellar rotation and age to assess the presence, or lack thereof, of prolonged star formation in the approximately 40-Myr-old cluster NGC1818. Our results provide evidence for a fast star formation, confined within 8 Myr, thus excluding age differences as responsible for the extended main-sequence turn-offs, and leading the way to alternative observational perspectives in the exploration of stellar populations in young clusters.

Prakash Sarkar, Biswajit Pandey, Suman Sarkar

10 pages, 4 figures, comments are welcome

The filaments and sheets are the most striking visual patterns in the cosmic web. The maximum extent of these large-scale structures are difficult to determine due to their structural variety and complexity. We construct a volume-limited sample of galaxies in a cubic region from the SDSS, divide it into smaller subcubes and shuffle them around. We quantify the average filamentarity and planarity in the three-dimensional galaxy distribution as a function of the density threshold and compare them with those from the shuffled realizations of the original data. The analysis is repeated for different shuffling lengths by varying the size of the subcubes. The average filamentarity and planarity in the shuffled data show a significant reduction when the shuffling scales are smaller than the maximum size of the genuine filaments and sheets. We observe a statistically significant reduction in these statistical measures even at a shuffling scale of $\sim 130 \, \rm Mpc$, indicating that the filaments and sheets in three dimensions can extend up to this length scale. They may extend to somewhat larger length scales that are missed by our analysis due to the limited size of the SDSS data cube. However, our analysis safely infer that the size of the most extended filament and most giant sheet in the galaxy distribution must be smaller than $200 \, \rm Mpc$.

Ian J. M. Crossfield, Matej Malik, Michelle L. Hill, Stephen R. Kane, Bradford Foley, Alex S. Polanski, David Coria, Jonathan Brande, Yanzhe Zhang, Katherine Wienke, Laura Kreidberg, Nicolas B. Cowan, Diana Dragomir, Varoujan Gorjian, Thomas Mikal-Evans, Bjoern Benneke, Jessie L. Christiansen, Drake Deming, Farisa Y. Morales

ApJL in press. 16 pages, 12 figures, 10 eclipses, 1 bandpass. Models will be available at journal website

The increasing numbers of rocky, terrestrial exoplanets known to orbit nearby stars (especially M dwarfs) has drawn increased attention to the possibility of studying these planets' surface properties, and atmospheric compositions & escape histories. Here we report the detection of the secondary eclipse of the terrestrial exoplanet GJ1252b using the Spitzer Space Telescope's IRAC2 4.5 micron channel. We measure an eclipse depth of 149(+25/-32) ppm, corresponding to a day-side brightness temperature of 1410(+91/-125) K and consistent with the prediction for no atmosphere. Comparing our measurement to atmospheric models indicates that GJ1252b has a surface pressure of <10 bar, substantially less than Venus. Assuming energy-limited escape, even a 100 bar atmosphere would be lost in <1 Myr, far shorter than estimated age of 3.9+/-0.4 Gyr. The expected mass loss could be overcome by mantle outgassing, but only if the mantle's carbon content were >7% by mass - over two orders of magnitude greater than that found in Earth. We therefore conclude that GJ1252b has no significant atmosphere. Model spectra with granitoid or feldspathic surface composition, but with no atmosphere, are disfavored at >2 sigma. The eclipse occurs just +1.4(+2.8/-1.0) min after orbital phase 0.5, indicating e cos omega=+0.0025(+0.0049/-0.0018), consistent with a circular orbit. Tidal heating is therefore likely to be negligible to GJ1252b's global energy budget. Finally, we also analyze additional, unpublished TESS transit photometry of GJ1252b which improves the precision of the transit ephemeris by a factor of ten, provides a more precise planetary radius of 1.180+/-0.078 R_E, and rules out any transit timing variations with amplitudes <1 min.

Łukasz Bratek, Joanna Jałocha, Marek Kutschera

This is not the final published version

Radial steady-state accretion of polytropic matter is investigated under cylindrical symmetry in the Levi-Civita background metric. The model can be considered as a cylindrical analog of Bondi accretion in strong gravitational field. As a byproduct of this study, the issue of defining the line mass density is addressed and the role of the metric free parameters is discussed on the example of physical observables. The form of radial accretion equations is insensitive to the structure of the interior solution. Accordingly, the accretion solution analysis can be limited to a special Wilson form of Levi-Civita metric describing a structureless homogeneous string.

Łukasz Bratek, Joanna Jałocha

the final version of this manuscript belongs to the Special Issue "Ultra High Energy Photons"

(also inside: this manuscript introduces the reader to the argument against the existence of magnetic monopoles, which forms an essential part of Staruszkiewicz's Quantum Mechanics of the Electric Charge) Ultra-high energy (UHE) photons with energies exceeding $10^{18}\eV$ can potentially be observed. They are produced in various processes involving electrically charged particles. However, more exotic scenarios are also possible. UHE photons could be emitted in encounters of massive magnetically charged monopole--antimonopole pairs or in processes associated with monopoles accelerated to high energies, typically $10^{21}\eV$ or beyond. Observing UHE photons can pose constraints on the properties of magnetic monopoles. There are compelling theoretical reasons in favor of the presence of magnetic monopoles in nature. The predicted observational signatures of these particles are therefore searched for in dedicated experiments currently in operation. Despite these attempts, magnetic monopoles have yet to be empirically proved. There are also theoretical reasons why magnetic monopoles allowed by Dirac's theory might not be realized in nature in the form of isolated particles. Detection or non-detection of UHE photon signatures of magnetic monopoles would bring us closer to solving this fascinating puzzle.

A. D. V. Di Virgilio, G. Terreni, A. Basti, N. Beverini, G. Carelli, D. Ciampini, F. Fuso, E. Maccioni, P. Marsili, J. Kodet, K. U. Schreiber

10 pages, 6 figures

The absolute measurement of the Earth angular rotation rate with ground-based instruments becomes challenging if the 1 part in $10^9$ of precision has to be obtained. This threshold is important for fundamental physics and for geodesy, to investigate effects of General Relativity and Lorentz violation in the gravity sector and to provide the fast variation of the Earth rotation rate. High sensitivity Ring Laser Gyroscopes (RLG) are currently the only promising technique to achieve this task in the near future, but their precision has been so far limited by systematics related to the laser operation. In this paper we analyze two different sets of observations, each of them three days long. They were obtained from the G ring laser at the Geodetic Observatory Wettzell. The applied method has been developed for the GINGERINO ring laser in order to identify and extract the laser systematics. For the available data sets the residuals show mostly white noise behavior and the Allan deviation drops below 1 part in $10^9$ after an integration time of about $10^4$~s.

Commenting the 11-year sunspot cycle, Wolf (1859, MNRAS 19, 85-86) conjectured that "the variations of spot frequency depend on the influences of Venus, Earth, Jupiter, and Saturn". The high synchronization of our planetary system is already nicely revealed by the fact that the ratios of the planetary orbital radii are closely related to each other through a scaling-mirror symmetry equation (Bank and Scafetta, Front. Astron. Space Sci. 8, 758184, 2022). Reviewing the many planetary harmonics and the orbital invariant inequalities that characterize the planetary motions of the solar system from the monthly to the millennial time scales, we show that they are not randomly distributed but clearly tend to cluster around some specific values that also match those of the main solar activity cycles. In some cases, planetary models have even been able to predict the time-phase of the solar oscillations including the Schwabe 11-year sunspot cycle. We also stress that solar models based on the hypothesis that solar activity is regulated by its internal dynamics alone have never been able to reproduce the variety of the observed cycles. Although planetary tidal forces are weak, we review a number of mechanisms that could explain how the solar structure and the solar dynamo could get tuned to the planetary motions. In particular, we discuss how the effects of the weak tidal forces could be significantly amplified in the solar core by an induced increase in the H-burning. Mechanisms modulating the electromagnetic and gravitational large-scale structure of the planetary system are also discussed.

Jesper Leong (1), Theo F. Motta ((1) and (2)), Anthony W. Thomas (1), P. A. M. Guichon (3) ((1) CSSM and ARC Centre of Excellence for Dark Matter Particle Physics, Department of Physics, University of Adelaide, (2) Institut für Theoretische Physik, Justus-Liebig-Universität Giessen, (3) Irfu, CEA, Université Paris-Saclay)

16 pages, 10 figures

The possibility of new short-distance physics applicable inside the cores of NS is incorporated into the equation of state generated by the quark-meson coupling model. The contribution of this new physics to the energy density is taken to be proportional to the amount of overlap between the quark cores of the baryons involved. With no change to the properties of symmetric nuclear matter at saturation density, including an incompressibility compatible with data on giant monopole resonances, one can sustain neutron stars with a maximum mass $M_{max}>2.1$ M$_\odot$, even when hyperons are included.

One of the fundamental challenges in string theory is to derive realistic four-dimensional cosmological backgrounds from it and it has been recently shown that there are strict consistency conditions which must be satisfied in string compactifications, thus constraining its possible low-energy backgrounds. In this work, we focus on energy conditions as \textit{covariant and background independent} consistency requirements in order to classify possible backgrounds coming from low-energy string theory in two steps. Firstly, we show how supergravity actions typically obey many relevant energy conditions, under some reasonable assumptions. Remarkably, we find that the energy conditions are satisfied even in the presence of objects which individually violate them due to the tadpole cancellation condition. Thereafter, we list a set of required conditions for a higher-dimensional energy condition to imply the corresponding lower-dimensional one, thereby categorizing the allowed low-energy solutions. As for any no-go theorem, our aim is to highlight the assumptions which must be circumvented for deriving four-dimensional spacetimes that necessarily violate these energy conditions, with emphasis on cosmological backgrounds.