Papers for Monday, Aug 29 2022

Our view of the last-scattering surface in the cosmic microwave background (CMB) is obscured by secondary anisotropies, sourced by scattering, extragalactic emission and gravitational processes between recombination and observation. Whilst it is established that non-Gaussianity from the correlation between the integrated-Sachs-Wolfe (ISW) effect and gravitational lensing can significantly bias primordial non-Gaussianity (PNG) searches, recent work by Hill (2018) has suggested that other combinations of secondary anisotropies can also produce significant biases. Building on that work, we use the WebSky and Sehgal et al.(2010) simulations to perform an extensive examination of possible biases to PNG measurements for the local, equilateral and orthogonal shapes. For a Planck-like CMB experiment, without foreground cleaning, we find significant biases from cosmic infrared background (CIB)-lensing and thermal Sunyaev-Zel'dovich (tSZ)-lensing bispectra for the local and orthogonal templates, and from CIB and tSZ bispectra for the equilateral template. For future experiments, such as the Simons Observatory, biases from correlations between the ISW effect and the tSZ and CIB will also become important. Finally we investigate whether foreground-cleaning techniques are able to suppress these biases sufficiently. We find that the majority of these biases are effectively suppressed by the internal-linear-combination method and the total bias to Planck-like and SO-like experiments is less than the $1\,\sigma$ statistical error. However, the small total bias arises from the cancellation of several $1\,\sigma$ biases for Planck-like experiments and $2\,\sigma$ biases for SO-like. As this cancellation is likely sensitive to the precise modelling, to ensure robustness against these biases explicit removal methods should be used, likely at the cost of decreased constraining power.

We present a search for sterile neutrino dark matter decay signals in the Milky Way's dark matter halo by considering the entirety of the Swift-XRT data archive. After filtering the raw archive, we analyze a $\sim$77 Ms data set containing the full field of view, as well as a $\sim$41 Ms data set with point-sources excised using the Swift-XRT Point Source catalog. We report non-detections of emission lines across the 3--6 keV continuum in both data sets, including at and around 3.5 keV. The point-sources excised data set is found to have higher sensitivity to faint dark matter decay signals due to its freedom from point-source contamination and is thus used to set constraints. Non-detections across the total data set's continuum are used to constrain the sterile neutrino dark matter parameter space, marginally strengthening existing X-ray constraints. Non-detections at $\sim$3.5 keV in data subsets grouped by angular distance from the galactic center are used to constrain the 3.5 keV line's galactic intensity profile, providing the strongest constraints to date across $\sim$1/4 of the galaxy.

The Pierre Auger Observatory (PAO) and Telescope Array (TA) collaborations report significant differences in the observed energy spectra of ultra-high-energy cosmic rays (UHECRs) above 30 EeV. In this work, we present a joint fit of TA and PAO data using the rigidity-dependent maximum energy model, and including full marginalization over all relevant parameters. We show that the presence of a local astrophysical source in the Northern Hemisphere, which is only visible by the TA experiment, can reconcile PAO and TA measurements up to the highest energies. We demonstrate that the presence of that local source is favored at the 5.6$\sigma$ level compared to the scenario where both experiments observe the same UHECR flux from a cosmological source distribution. We also quantify that the astrophysical explanation can describe the current data better than a scenario where the differences in the observations are explained by experimental systematics (i.e., energy-dependent shifts). Having tested different mass compositions emitted from the local source, we conclude that the data are best described by a source lying at a distance of about 14 Mpc that emits cosmic rays dominated by the silicon mass group; we also discuss possible source candidates.

A longstanding problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models based on the traditional Chandrasekhar DF formula have been proposed, yet they suffer from fundamental ambiguities in the definition of some terms in Chandrasekhar's formula when applied to real galaxies, as well as difficulty in evaluating continuous quantities from (spatially) discrete simulation data. In this work we present a new sub-grid dynamical friction estimator based on the discrete nature of $N$-body simulations, which also avoids the ambiguously-defined quantities in Chandrasekhar's formula. We test our estimator in the GIZMO code and find that it agrees well with high-resolution simulations where DF is fully captured, with negligible additional computational cost. We also compare it with a Chandrasekhar estimator and discuss its applications in real galactic simulations.

The Transiting Exoplanet Survey Satellite (TESS) is focusing on relatively bright stars and has found thousands of planet candidates. However, mainly because of the low spatial resolution of its cameras ($\approx$ 21 arcsec/pixel), TESS is expected to detect several false positives (FPs); hence, vetting needs to be done. Here, we present a follow-up program of TESS candidates orbiting solar-analogue stars that are in the all-sky PLATO input catalogue. Using Gaia photometry and astrometry we built an absolute colour-magnitude diagram and isolated solar-analogue candidates' hosts. We performed a probabilistic validation of each candidate using the VESPA software and produced a prioritized list of objects that have the highest probability of being genuine transiting planets. Following this procedure, we eliminated the majority of FPs and statistically vetted 23 candidates. For this remaining set, we performed a stellar neighbourhood analysis using Gaia Early Data Release 3 and centroid motion tests, greatly enhancing the on-target probability of 12 of them. We then used publicly available high-resolution imaging data to confirm their transit source and found five new, fully validated planets. For the remaining candidates, we propose on-off photometry to further refine the list of genuine candidates and prepare for the subsequent radial velocity follow-up.

Context: In recent years, stellar intensity interferometry has seen renewed interest from the astronomical community because it can be efficiently applied to Cherenkov telescope arrays. Aims: We have investigated the accuracy that can be achieved in reconstructing stellar sizes by fitting the visibility curve measured on the ground. The large number of expected available astronomical targets, the limited number of nights in a year, and the likely presence of multiple baselines will require careful planning of the observational strategy to maximise the scientific output. Methods: We studied the trend of the error on the estimated angular size, considering the uniform disk model, by varying several parameters related to the observations, such as the total number of measurements, the integration time, the signal-to-noise ratio, and different positions along the baseline. Results: We found that measuring the value of the zero-baseline correlation is essential to obtain the best possible results. Systems that can measure this value directly or for which it is known in advance will have better sensitivity. We also found that to minimise the integration time, it is sufficient to obtain a second measurement at a baseline half-way between 0 and that corresponding to the first zero of the visibility function. This function does not have to be measured at multiple positions. Finally, we obtained some analytical expressions that can be used under specific conditions to determine the accuracy that can be achieved in reconstructing the angular size of a star in advance. This is useful to optimise the observation schedule.

The early phases of the observed evolution of the supernovae (SNe) are expected to be dominated by the shock breakout and ``flash" ionization of the surrounding circumstellar medium. This material arises from the last stages of the evolution of the progenitor, such that photometry and spectroscopy of SNe at early times can place vital constraints on the latest and fastest evolutionary phases leading up to stellar death. These signatures are erased by the expansion of the ejecta within ~5 days after explosion. Here we present the earliest constraints, to date, on the polarization of ten transients discovered by the Zwicky Transient Facility (ZTF), between June 2018 and August 2019. Rapid polarimetric followup was conducted using the Liverpool Telescope RINGO3 instrument, including 3 SNe observed within <1 day of detection by the ZTF. The limits on the polarization within the first 5 days of explosion, for all SN types, is generally <2%, implying early asymmetries are limited to axial ratios >0.65 (assuming an oblate spheroidal configuration). We also present polarimetric observations of the Type I Superluminous SN 2018bsz and Type II SN 2018hna, observed around and after maximum light.

The remarkably tight relationship between optical color and stellar mass-to-light ratio ($M_*/L$) in galaxies is widely used for efficient stellar mass estimates. However, it has remained unclear whether this low scatter comes from a natural order in the galaxy population, or whether it is driven by simple relationships in the models used to describe them. In this work we investigate the origins of the relationship by contrasting the derived relationship from a simple 4-parameter physical model with a more sophisticated 14-dimensional Prospector-$\alpha$ model including nonparametric star formation histories. We apply these models to 63,430 galaxies at $0.5<z<3$ from the 3D-HST survey and fit the results with a hierarchical Bayesian model (HBM) for the population distribution in the $(g-r)$--$\log(M/L_g)$ plane. We find that Prospector-$\alpha$ infers systematically higher $M_*/L$ by 0.12 dex, a result of the nonparametric star formation history producing older ages. Prospector-$\alpha$ also infers systematically redder rest-frame $(g-r)$ by 0.06 dex owing to the nebular emission. Surprisingly, we observe similar average color--$M_*/L$ relationships for the two models due to the combined effects of the $M_*/L$ and $(g-r)$ offsets. Nevertheless, Prospector-$\alpha$ produces a much looser color-M/L relationship with a scatter of 0.28 dex compared to the simple model of 0.12 dex. Also, unlike the simple model, the Prospector-$\alpha$ model shows a substantial redshift evolution in the relationship due to stellar aging. Finally, we demonstrate that the HBM produces substantial shrinkage in the individual posteriors for faint galaxies, an important first step towards using the observed galaxy population directly to inform the priors in galaxy SED-fitting.

Anti-solar differential rotation profiles have been found for decades in numerical simulations of convective envelopes of solar-type stars. These profiles are characterized by a slow equator and fast poles (i.e., reversed with respect to the Sun) and have been found in simulations for high Rossby numbers (slow rotators). Rotation profiles like this have been reported observationally in evolved stars, but have never been unambiguously observed for cool solar-type stars on the main sequence. In this context, detecting this regime in main-sequence solar-type stars would improve our understanding of their magnetorotational evolution. The goal of this study is to identify the most promising cool main-sequence stellar candidates for anti-solar differential rotation in the \textit{Kepler} sample. First, we introduce a new theoretical formula to estimate fluid Rossby numbers, $Ro_{\rm f}$, of main-sequence solar-type stars, from observational quantities, and taking the influences of the internal structure and metallicity into account. We obtain a list of the most promising stars that are likely to show anti-solar differential rotation. We identify two samples: one at solar metallicity, including 14 targets, and another for other metallicities, including 8 targets. We find that the targets with the highest $Ro_{\rm f}$ are likely to be early-G or late-F stars at about log$_{10}g=4.37$~dex. We conclude that cool main-sequence stellar candidates for anti-solar differential rotation exist in the \textit{Kepler} sample. The most promising candidate is KIC~10907436, and two other particularly interesting candidates are the solar analog KIC~7189915 and the seismic target KIC~12117868. Future characterization of these 22 stars is expected to help us understand how dynamics can impact magnetic and rotational evolution of old solar-type stars at high Rossby number.

MWC656 has been reported as classical Be star with a black hole companion. Revisited spectral variability properties render this unlikely, with a hot subdwarf more probable.

We present 0.9 mm continuum and CO (3-2) line emission observations retrieved from the Atacama Large Millimeter/submillimeter Array (ALMA) archive toward the high-mass star formation region IRAS 16076-5134. We identify fourteen dense cores with masses between 0.3 to 22 M$_{\odot}$. We find an ensemble of filament-like CO (3-2) ejections from -62 to +83 km s$^{-1}$ that appear to arise radially from a common central position, close to the dense core MM8. The ensemble of filaments, has a quasi-isotropic distribution in the plane of the sky. The radial velocity of several filaments follow a linear velocity gradient, incresing from a common origin. Considering the whole ensemble of filaments, we estimate its total mass to be 138 and 216 M$_{\odot}$ from its CO emission, for 70 K and 140 K respectively. Also, assuming a constant velocity expansion of the filaments (of 83 km s$^{-1}$) we estimate the dynamical age of the outflowing material (3500 years), its momentum (~10$^{4}$ M$_{\odot}$ km s$^{-1}$) and its kinetic energy (~10$^{48-49}$ erg). The morphology and kinematics presented by the filaments suggest the presence of a dispersal outflow with explosive characteristics in IRAS 16076-5134. In addition, we make a raw estimate of the lower limit of the frequency rate of the explosive dispersal outflows in the Galaxy (one every 110 years) considering constant star formation rate and efficiency with respect to the galactocentric radius of the Galaxy. This may imply a comparable rate of dispersal outflows and supernovae (approximately one every 50 years), which may be important for the energy budget of the Interstellar Medium and the link between dispersal outflows and high-mass star formation.

We present a computational method to identify glitches in gravitational-wave data that occur nearby gravitational-wave signals. We flag any excess in the data surrounding a signal and compute the probability of such an excess occurring in Gaussian noise. We validate that the probabilities reported by this tool are self-consistent in colored Gaussian noise as well as data that contains a gravitational-wave event after subtracting the signal using the best-fit template. Furthermore, we compare our glitch identification results for events from LIGO-Virgo's third observing run against the list of events that required glitch mitigation. Finally, we discuss how the precise, automated information about the data quality surrounding gravitational-wave events this tool provides can be used to improve astrophysical analyses of these events.

The low-energy gamma-ray domain is an important window for the study of the high energy Universe. Here matter can be observed in extreme physical conditions and during powerful explosive events. However, observing gamma-rays from faint sources is extremely challenging with current instrumentation. With techniques used at present collecting more signal requires larger detectors, leading to an increase in instrumental background. For the leap in sensitivity that is required for future gamma-ray missions use must be made of flux concentrating telescopes. Fortunately, gamma-ray optics such as Laue or Fresnel lenses, based on diffraction, make this possible. Laue lenses work with moderate focal lengths (tens to a few hundreds of metres), but provide only rudimentary imaging capabilities. On the other hand, Fresnel lenses offer extremely good imaging, but with a very small field of view and a requirement for focal lengths $\sim$10$^8$ m. This chapter presents the basic concepts of these optics and describes their working principles, their main properties and some feasibility studies already conducted.

The 1950 discovery of activity emanating from asteroid (4015) Wilson-Harrington prompted astronomers to realize comet-like activity is not limited to comets. Since then < 30 active asteroids have been discovered, yet they hold clues about fundamental physical and chemical processes in the solar system. Around half of the activity is attributed to sublimation, highlighting asteroids as a "volatile reservoir" - a dynamical group of minor planets that harbor volatiles. Centaurs, found between the orbits of Jupiter and Neptune, were first recognized in 1977 and represent another reservoir. Active Centaurs are also rare, with < 20 known. Understanding the solar system volatile distribution has broad implications, from informing space exploration programs to illuminating how planetary systems form with volatiles prerequisite to life as we know it. We set out to increase the number of known active objects to enable their study as populations. In this dissertation I present (1) our pipeline that extracts images of known minor planets for presentation to Citizen Scientists, (2) our proof-of-concept demonstrating Dark Energy Camera images are ideal for activity detection (Chandler et al. 2018), (3) how we identified a potential new recurrent activity mechanism (Chandler et al. 2019), (4) a Centaur activity discovery, plus a novel technique to estimate which species are sublimating (Chandler et al. 2020), (5) how our project enabled us to classify an object as a member of the main-belt comets (Chandler et al. 2021), a rare (<10) active asteroid subset known for sublimation-driven activity, (6) the identification of a Quasi-Hilda comet and a dynamical pathway that may explain the presence of some active asteroids (Chandler et al. 2022), and (7) our NASA Partner Citizen Science project Active Asteroids (this http URL), including initial results.

It has long been uncertain how effective the dust-corrected H$\alpha$-to-UV luminosity ratio ($L(\rm H\alpha)/L(\rm UV)$) is in probing bursty SFHs for high-redshift galaxies. To address this issue, we present a statistical analysis of the resolved distribution of star-formation-rate surface density ($\Sigma_{\rm SFR}$) as well as stellar age and their correlations with the globally measured $L(\rm H\alpha)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37<z<1.70$ and $2.09<z<2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $\Sigma_{\rm SFR}$ and stellar age maps. We also analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck-LRIS, which includes 124 star-forming galaxies at redshifts $1.4<z<2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm H\alpha)/L(\rm UV)$. We find no evidence that galaxies with higher $L(\rm H\alpha)/L(\rm UV)$ are undergoing a burst of star formation based on the distribution of $\Sigma_{\rm SFR}$ and stellar age as well as the strengths of Siiv$\lambda\lambda1393,1402$ and Civ$\lambda\lambda1548,1550$ P-Cygni features from massive stars. Our results suggest that the $L(\rm H\alpha)/L(\rm UV)$ ratio is not a reliable tracer of bursty SFHs for typical star-forming galaxies at high redshift. We also study the variations observed in the strength of the nebular Heii$\lambda1640$ emission between the $L(\rm H\alpha)/L(\rm UV)$ subsamples. We find that variations in the upper-mass limit of the IMF cannot fully account for all the observed difference in the HeII emission between the $L(\rm H\alpha)/L(\rm UV)$ subsamples, and that another source of He$^{+}$ ionizing photons, such as X-ray binaries, may be needed to explain such a difference.

Coagulation of particles occurs when two particles collide and stick together. In the Martian atmosphere, dust coagulation would increase the effective particle size, as small particles accrete to larger particles. Murphy et al. (1990) showed that Brownian coagulation of dust in the Martian atmosphere was not significant, due to the low dust particle mixing ratios, while Montmessin et al. (2002) and Fedorova et al. (2014) showed that it mostly involves particle radii smaller than 0.1 um. However, the effects of coagulation have never been explored in 3D, during a global dust storm, i.e. in presence of larger numbers of small particles. Here we revisit this issue by using the NASA Ames Mars Global Climate Model (MGCM) to investigate the temporal and spatial changes in dust particle sizes during the 2018 global storm due to coagulation and the overall impact of these processes on Mars' climate. Our parameterization for coagulation includes the effect of Brownian motion, Brownian diffusion enhancement, and gravitational collection. We show that Brownian motion and Brownian diffusion enhancement dominate gravitational collection. The impact of coagulation is significant during the storm, with coagulation rates increased by a factor of 10 compared to non-storm conditions. The effective particle radius can be increased by a factor of 2 due to coagulation, leading to a 20K colder atmosphere above 30 km altitude. Overall, our parameterization improves the representation of the decay phase of the storm relative to observations. The process remains significant outside the storm period if large numbers of submicron-sized particles are involved. It may be possible, in GCMs, to lift larger amounts of submicron-sized particles from the surface without excess dust buildup in the atmosphere, thus improving the agreement with some of the observations without diverging from the observed column opacities.

In order to understand the formation and evolution of galaxies fully, it is important to study their three-dimensional gravitational potential for a large sample of galaxies. Since polar-ring galaxies (PRGs) provide useful laboratories for this investigation, we have started our detailed study of a sample of known PRGs by using the data set obtained by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). During the course of this study, we have discovered a new PRG, identified as SDSS J095351.58+012036.1. Its photometric redshift is estimated as z ~ 0.2. The polar ring structure in this PRG appears to be almost perpendicular to the disk of its host galaxy without any disturbed features. Therefore, this PRG will provide us with useful information on the formation of such an undisturbed polar structure. We discuss its photometric properties in detail.

The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of $35\times 35$ arcminutes and will measure the Stokes parameters $I$, $q$, and $u$ in a single exposure in the SDSS-r broadband filter. The design goal for the instrument is to achieve an overall polarimetric measurement accuracy of 0.1 % over the entire field of view. We present here the complete polarimetric modeling of the instrument, characterizing the amount and sources of instrumental polarization. To accurately retrieve the real Stokes parameters of a source from the measured values, we have developed a calibration method for the instrument. Using this calibration method and simulated data, we demonstrate how to correct instrumental polarization and obtain 0.1 % accuracy in the degree of polarization, $p$. Additionally, we tested and validated the calibration method by implementing it on a table-top WALOP-like test-bed polarimeter in the laboratory.

For a flat $\Lambda$CDM universe, the dipole of the luminosity distance can be utilized to measure the Hubble parameter. It is here shown that this is not the case in more general settings where curvature and cosmic backreaction is permitted. This implies that a discordance between $H(z)$ measurements obtained using such dipole luminosity distance data and "true"/actual $H(z)$ data obtained from e.g. cosmic chronometers is a signal of curvature and/or cosmic backreaction. \newline\indent By considering mock future gravitational wave measurements of the Hubble parameter obtained through the dipole luminosity distance, it is shown that already a $1\%$ curvature could in principle just barely show up in the determination. However, for realistic mock data generation using models with as much as 5 \% curvature, parameter estimates do not yield reliable measures of inconsistency between the false $H(z)$ measurements and true measurements of $H(z)$. At the same time, cosmic backreaction is hard to detect even if it makes up $10\%$ of the "energy budget" in the current universe, even when considering a highly idealized situation with low errors. The results concerning backreaction are based on specific "scaling solutions" to the backreaction problem and the study shows that the possibility of detecting a signal of backreaction through the dipole of the luminosity distance depends strongly on the particular backreaction model.

One of the outstanding mysteries surrounding the rich diversity found in supernova remnants (SNRs) is the recent discovery of over-ionized or recombining plasma from a number of dynamically evolved objects. To help decipher its formation mechanism, we have developed a new simulation framework capable of modeling the time evolution of the ionization state of the plasma in an SNR. The platform is based on a one-dimensional hydrodynamics code coupled to a fully time-dependent non-equilibrium ionization calculation, accompanied by a spectral synthesis code to generate space-resolved broadband X-ray spectra for SNRs at arbitrary ages. We perform a comprehensive parametric survey to investigate the effects of different circumstellar environments on the ionization state evolution in SNRs up to a few 10,000 years. A two-dimensional parameter space, spanned by arrays of interstellar medium (ISM) densities and mass-loss rates of the progenitor, is used to create a grid of models for the surrounding environment, in which a core-collapse explosion is triggered. Our results show that a recombining plasma can be successfully reproduced in the case of a young SNR (a few 100 to 1,000 years old) expanding fast in a spatially extended low-density wind, an old SNR (> a few 1,000 years) expanding in a dense ISM, or an old SNR broken out from a confined dense wind region into a tenuous ISM. Finally, our models are confronted with observations of evolved SNRs, and an overall good agreement is found except for a couple of outliers.

(Abridged) In this paper, we present the results of an exceptional repeating X-ray nuclear transient, eRASSt J045650.3-203750 (hereafter J0456-20), uncovered by SRG/eROSITA in a quiescent galaxy at redshift of z~0.077. The main results are: 1) J0456-20 cycles through four distinctive phases: an X-ray rising phase leading into an X-ray plateau phase which lasts for ~2 months. This is terminated by a rapid X-ray flux drop phase during which the X-ray flux can drastically drop by more than a factor of 100 within 1 week followed by an X-ray faint state for about two months before it starts the X-ray rising phase again; 2) the X-ray spectra are generally soft in the rising phase with a photon index >3.0, and become harder as the X-ray flux increases. There is evidence of a multi-colour disk with inner region temperature of $T_\text{in}=70$ eV at the beginning of the X-ray rising phase. The high quality XMM-Newton data suggest that a warm and hot corona could be responsible for the X-ray emission, through inverse Comptonisation of soft disk seed photons, during the plateau phase and at the bright end of the rising phase; 3) J0456-20 shows only moderate UV variability and no significant optical variability; 4) radio emission is only detected (as yet) in the X-ray plateau phase, and shows a rapid decline on a time-scale of 2 weeks. We conclude that J0456-20 is likely a repeating nuclear transient with a tentative recurrence time of ~223 days. We discuss several possibilities to explain J0456-20's observational properties, and currently favour a repeating partial tidal disruption event (TDE) as the most likely scenario. The long-term X-ray evolution is explained as a transition between a thermal disk-dominated soft state and a steep power-law state, implying that the corona can be formed within a few months and destroyed within a few weeks.

In the past three decades, many stars orbiting about the supermassive black hole (SMBH) at the Galactic Centre (Sgr A*) were identified. Their orbital nature can give stringent constraints for the mass of the SMBH. In particular, the star S2 has completed at least one period since our first detection of its position, which can provide rich information to examine the properties of the SMBH, and the astrophysical environment surrounding the SMBH. Here, we report an interesting phenomenon that if a significant amount of dark matter or stellar mass is distributed around the SMBH, the precession speed of the S2 stellar orbit could be `slow down' by at most 27\% compared with that without dark matter surrounding the SMBH, assuming the optimal dark matter scenario. We anticipate that future high quality observational data of the S2 stellar orbit or other stellar orbits can help reveal the actual mass distribution near the SMBH and the nature of dark matter.

Our objective is to critically assess the X-ray flux variability as a tool for measuring the black hole (BH) mass in active galactic nuclei (AGN). We aim to establish a prescription for estimating BH masses based on measurements of the normalised excess variance from X-ray data. We discuss the minimum requirements in terms of the light-curve duration and X-ray signal-to-noise ratio (S/N) to enable a reliable determination that is comparable to what can be derived from the continuum and emission line reverberation studies. We used the light curves of local Seyfert from the Nuclear Spectroscopic Telescope Array hard X-ray mission (NuSTAR), to compute the normalised excess variance (NXV) in the 3-10 and 10-20 keV bands, thus extending the analysis to an energy band higher than 10 keV. The excess variance measurements were then combined with independent BH mass estimates from the literature to establish the MBH versus NXV relation for different samples and weigh its accuracy in terms of the light-curve duration and X-ray S/N. We find that it is possible to accurately measure the BH mass in AGN using excess variance measurements in the 3-10 and the 10-20 keV bands, however, strong quality requirements should be applied. The minimum necessary S/N and duration of the light curves used to compute the excess variance ought to be 3 and approximately 100 ks, respectively. We provide a linear relationship between the normalised excess variance and the black hole mass that can be used to estimate the latter, with an average uncertainty of the order of 0.4 to 0.25 dex (depending on the adopted light-curve segment duration).

Ultra-diffuse galaxies (UDGs) are as faint as dwarf galaxies but whose sizes are similar to those of spiral galaxies. A variety of formation mechanisms have been proposed, some of which could result in different disk thicknesses. In this study, we measure the radial profile of the HI scale height (h_g) and flaring angle (h_g/R) of AGC 242019 through the joint Poisson-Boltzmann equation based on its well spatially-resolved HI gas maps. The mean HI scale height of AGC 242019 is <h_g> \approx 537.15 \pm 89.4 pc, and the mean flaring angle is <h_g/R> \approx 0.19 \pm 0.03. As a comparison, we also derive the disk thickness for a sample of 14 dwarf irregulars. It is found that the HI disk of AGC 242019 has comparable thickness to dwarfs. This suggests that AGC 242019 is unlikely to experience much stronger stellar feedback than dwarf galaxies, which otherwise leads to a thicker disk for this galaxy.

Magnetic reconnection preferentially takes place at the intersection of two separatrices or two quasi-separatrix layers, which can be quantified by the squashing factor Q, whose calculation is computationally expensive due to the need to trace as many field lines as possible. We developed a method (FastQSL) optimized for obtaining Q and the twist number in a 3D data cube. FastQSL utilizes the hardware acceleration of the graphic process unit (GPU) and adopts a step-size adaptive scheme for the most computationally intensive part: tracing magnetic field lines. As a result, it achieves a computational efficiency of 4.53 million Q values per second. FastQSL is open source, and user-friendly for data import, export, and visualization.

We present the results of analyses of the ground level enhancements (GLEs) of cosmic ray (CR) events of 29 September 1989; 15 April 2001 and 20 January 2005. This involve examination of hourly raw CR counts of an array of neutron monitors (NMs) spread across different geographical latitudes and longitudes. Using awk script and computer codes implemented in R-software, the pressure corrected raw data plots of the NMs were grouped into low-, mid-, and, high-latitudes. The results show both similarities and differences in the structural patterns of the GLE signals. In an attempt to explain why the CR count during the decay phase of GLEs is always higher than the count before peak, we interpreted all counts prior to the peak as coming from direct solar neutrons and those in the decay phase including the peak as coming from secondary CR neutrons generated by the interactions of primary CRs with the atoms and molecules in the atmosphere. We identified NMs that detected these primary neutrons and found that they are close in longitudes. Previous authors seemingly identified these two species as impulsive and gradual events. Although there are a number of unexplained manifestations of GLE signals, some of the results suggest that geomagnetic rigidity effectively determines the intensity of CRs at low- and mid-latitudes. Its impact is apparently insignificant in high-latitude regions. Nevertheless, the results presented should be validated before making any firm statements. Principally, the contributions of the ever-present and intractable CR diurnal anisotropies to GLE signals should be accounted for in future work.

Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K \sigma'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $\sigma'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a S\'ersic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $\sigma'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $\sigma'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a S\'ersic profile-based $R$ estimate, then a S\'{e}rsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex.

The first experimental study of the low-temperature kinetics of the gas-phase reaction of NH2 with formaldehyde (CH2O) has been performed. This reaction has previously been suggested as a source of formamide (NH2CHO) in interstellar environments. A pulsed Laval nozzle equipped with laser-flash photolysis and laser-induced fluorescence spectroscopy was used to create and monitor the temporal decay of NH2 in the presence of CH2O. No loss of NH2 could be observed via reaction with CH2O and we place an upper-limit on the rate coefficient of <6x10-12 cm3 molecule-1 s-1 at 34K. Ab initio calculations of the potential energy surface were combined with RRKM calculations to predict a rate coefficient of 6.2x10-14 cm3 molecule-1 s-1 at 35K, consistent with the experimental results. The presence of a significant barrier, 18 kJ mol-1, for the formation of formamide as a product, means that only the H-abstraction channel producing NH3 + CHO, in which the transfer of an H-atom can occur by quantum mechanical tunnelling through a 23 kJ mol-1 barrier, is open at low temperatures. These results are in contrast with a recent theoretical study which suggested that the reaction could proceed without a barrier and was therefore a viable route to gas-phase formamide formation. The calculated rate coefficients were used in an astrochemical model which demonstrated that this reaction produces only negligible amounts of gas-phase formamide under interstellar and circumstellar conditions. The reaction of NH2 with CH2O is therefore not an important source of formamide at low temperatures in interstellar environments.

We report a unique conjugate observation of fast flows and associated current sheet disturbances in the near-Earth magnetotail by MMS (Magnetospheric Multiscale) and Cluster preceding a positive bay onset of a small substorm at ~14:10 UT, Sep. 8, 2018. MMS and Cluster were located both at X ~-14 RE. A dipolarization front (DF) of a localized fast flow was detected by Cluster and MMS, separated in the dawn-dusk direction by ~4 RE, almost simultaneously. Adiabatic electron acceleration signatures revealed from comparison of the energy spectra confirm that both spacecraft encounter the same DF. We analyzed the change in the current sheet structure based on multi-scale multi-point data analysis. The current sheet thickened during the passage of DF, yet, temporally thinned subsequently associated with another flow enhancement centered more on the dawnward side of the initial flow. MMS and Cluster observed intense perpendicular and parallel current in the off-equatorial region mainly during this interval of the current sheet thinning. Maximum field-aligned currents both at MMS and Cluster are directed tailward. Detailed analysis of MMS data showed that the intense field-aligned currents consisted of multiple small-scale intense current layers accompanied by enhanced Hall-currents in the dawn-dusk flow-shear region. We suggest that the current sheet thinning is related to the flow bouncing process and/or to the expansion/activation of reconnection. Based on these mesoscale and small-scale multipoint observations, 3D evolution of the flow and current-sheet disturbances was inferred preceding the development of a substorm current wedge.

We present high-resolution (~2-3"; ~0.1 pc) radio observations of the Galactic center cloud M0.10-0.08 using the Very Large Array at K and Ka band (~25 and 36 GHz). The M0.10-0.08 cloud is located in a complex environment near the Galactic center Radio Arc and the adjacent M0.11-0.11 molecular cloud. From our data, M0.10-0.08 appears to be a compact molecular cloud (~3 pc) that contains multiple compact molecular cores (5+; <0.4 pc). In this study we detect a total of 15 molecular transitions in M0.10-0.08 from the following molecules: NH3, HC3N, CH3OH, HC5N, CH3CN, and OCS. We have identified more than sixty 36 GHz CH3OH masers in M0.10-0.08 with brightness temperatures above 400 K and 31 maser candidates with temperatures between 100-400 K. We conduct a kinematic analysis of the gas using NH3 and detect multiple velocity components towards this region of the Galactic center. The bulk of the gas in this region has a velocity of 51.5 km/s (M0.10-0.08) with a lower velocity wing at 37.6 km/s. We also detect a relatively faint velocity component at 10.6 km/s that we attribute to being an extension of the M0.11-0.11 cloud. Analysis of the gas kinematics, combined with past X-ray fluorescence observations, suggests M0.10-0.08 and M0.11-0.11 are located in the same vicinity of the Galactic center and could be physically interacting.

The Near-Infrared Spectrograph (NIRSpec) is one the four focal plane instruments on the James Webb Space Telescope (JWST) which was launched on December 25, 2021. We present the in-flight status and performance of NIRSpec's detector system as derived from the instrument commissioning data. The instrument features two 2048 x 2048 HAWAII-2RG sensor chip assemblies (SCAs) that are operated at a temperature of about 42.8 K and are read out via a pair of SIDECAR ASICs. NIRSpec supports "Improved Reference Sampling and Subtraction" (IRS2) readout mode that was designed to meet NIRSpec's stringent noise requirements and to reduce 1/f and correlated noise. In addition, NIRSpec features subarrays optimized for bright object time series observations, e.g. for the observation of exoplanet transit around bright host stars. We focus on the dark signal as well as the read and total noise performance of the detectors.

The large galactic scales are connected to the many orders of magnitude smaller supermassive black hole (SMBH) scales by an episodic cycle of feeding and feedback. Active galactic nuclei (AGN) are powered by accretion onto SMBH and the majority of AGN energy, in near-Eddington regime, is produced in thin sub-pc accretion discs. Currently, it is very difficult to model processes that occur on vastly different scales, ranging from the circumnuclear gas reservoirs at tens to hundreds of parsecs, down to the accretion disc scales at <0.01 pc. While sub-grid prescriptions used in large-scale or cosmological simulations are able to reproduce large-scale feedback, we propose using a more realistic model in parsec-scale simulations, where it is important to get accurate timescales to understand how feedback affects gas dynamics and star formation in the vicinity of the AGN. To test our approach we use a sub-resolution thin accretion disc model, coupled to the SMBH, in a set of hydrodynamical simulations of a retrograde collision between a gas ring and a molecular cloud in an environment similar to the Galactic centre using the SPH code Gadget-3. The disc-mediated feeding of the SMBH is relatively smooth and delayed compared to an instantaneous feeding prescription. While the reduction of accretion due to feedback is present in both accretion disc and instantaneous feeding simulations, a clear central cavity appears only in accretion disc runs - hinting that a less volatile accretion phase could have a greater impact on the surrounding gas.

We analyze 5108 AFGKM stars with at least five high precision radial velocity points as well as Gaia and Hipparcos astrometric data utilizing a novel pipeline developed in previous work. We find 914 radial velocity signals with periods longer than 1000\,d. Around these signals, 167 cold giants and 68 other types of companions are identified by combined analyses of radial velocity, astrometry, and imaging data. Without correcting for detection bias, we estimate the minimum occurrence rate of the wide-orbit brown dwarfs to be 1.3\%, and find a significant brown dwarf valley around 40 $M_{\rm Jup}$. We also find a power-law distribution in the host binary fraction beyond 3 au similar to that found for single stars, indicating no preference of multiplicity for brown dwarfs. Our work also reveals nine sub-stellar systems (GJ 234 B, GJ 494 B, HD 13724 b, HD 182488 b, HD 39060 b and c, HD 4113 C, HD 42581 d, HD 7449 B, and HD 984 b) that have previously been directly imaged, and many others that are observable at existing facilities. Depending on their ages we estimate that an additional 10-57 sub-stellar objects within our sample can be detected with current imaging facilities, extending the imaged cold (or old) giants by an order of magnitude.

We report the discovery of a new example of the rare class of highly magnetized, rapidly rotating, helium enhanced, early B stars that produce anomalously wide hydrogen emission due to a centrifugal magnetosphere (CM). The star is Trumpler 16-26, a B1.5 V member of the Trumpler 16 open cluster. A CM was initially suspected based on hydrogen Brackett series emission observed in SDSS/APOGEE $H$-band spectra. Similar to the other stars of this type, the emission was highly variable and at all times remarkable due to the extreme velocity separations of the double peaks (up to 1300 km s$^{-1}$.) Another clue lay in the TESS lightcurve, which shows two irregular eclipses per cycle when phased with the likely 0.9718115 day rotation period, similar to the behavior of the well known CM host star $\sigma$ Ori E. To confirm a strong magnetic field and rotation-phase-locked variability, we initiated a follow-up campaign consisting of optical spectropolarimetry and spectroscopy. The associated data revealed a longitudinal magnetic field varying between $-3.1$ and $+1.6$ kG with the period found from photometry. The optical spectra confirmed rapid rotation ($v \sin i=195$ km s$^{-1}$), surface helium enhancement, and wide, variable hydrogen emission. Tr16-26 is thus confirmed as the 20$^{\rm th}$ known, the fourth most rapidly rotating, and the faintest CM host star yet discovered. With a projected dipole magnetic field strength of $B_{\rm d}>11$ kG, Tr16-26 is also among the most magnetic CM stars.

In cold and shielded environments, molecules freeze out on dust grain surfaces to form ices such as H2O, CO, CO2, CH4, CH3OH, and NH3. In protoplanetary disks, the exact radial and vertical ice extension depend on disk mass, geometry, and stellar UV irradiation. The goal of this work is to present a computationally efficient method to compute ice and bare-grain opacities in protoplanetary disk models consistently with the chemistry and to investigate the effect of ice opacities on the physico-chemical state and optical appearance of the disk. A matrix of Mie efficiencies is pre-calculated for different ice species and thicknesses, from which the position dependent opacities of icy grains are then interpolated. This is implemented in the PRODIMO code by a self-consistent solution of ice opacities and the local composition of ices, which are obtained from our chemical network. Locally, the opacity can change significantly, for example, an increase by a factor of more than 200 in the midplane, especially at UV and optical wavelengths, due to ice formation. This is mainly due to changes in the size distribution of dust grains resulting from ice formation. However, since the opacity only changes in the optically thick regions of the disk, the thermal disk structure does not change significantly. For the same reason, the spectral energy distributions computed with our disk models with ice opacities generally show only faint ice emission features at far-IR wavelengths. The ice absorption features are only seen in the edgeon orientation. The assumption made on how the ice is distributed across the grain size distribution influences the far-IR and millimeter slope of the SED. The ice features and their strengths are influenced by the ice power law and the type of chemistry. Our models predict stronger ice features for observations that can spatially resolve the disk, particularly in absorption.

The SVOM mission under development will carry four instruments, and in particular the coded-mask telescope named ECLAIRs, with a large field of view of about 2 sr, operating in the 4-150 keV energy band. The trigger software on board ECLAIRs will search for high-energy transients such as gamma-ray bursts and peculiar behaviour (e.g. strong outbursts) from known X-ray sources, in order to repoint the satellite to perform follow-up observations with the onboard narrow field of view instruments. The image trigger, one of the two algorithms implemented in the software on board ECLAIRs, produces images over periods of exposure ranging from 20 seconds to 20 minutes during which the Earth can cross the field of view. The CXB and contributions from known X-ray sources are expected to dominate the ECLAIRs astrophysical and instrumental background and must be taken into account and corrected prior to coded-mask image deconvolution in order to optimise the sensitivity to faint transients. To correct these background components, we implemented and studied a traditional fitting method and a new method based on wavelet decomposition of the detector image. In order to study and to assess the performance of these methods, we performed a one-year simulation of the image trigger on board ECLAIRs. From the images produced during this realistic observation scenario of the SVOM mission, we also defined a way to analyse the sky images to search for new sources. We present the algorithms behind the image trigger on board SVOM/ECLAIRs. We show that the wavelet method we implemented provides similar results in terms of cleaning performance compared to the traditional fitting method, and has the benefit of not requiring any assumption on the shape of the background on the detector. We also calibrate the detection threshold to be adaptive and based on the quality of the reconstructed sky image.

Forecasting solar energetic particles (SEPs), and identifying flare/CMEs from active regions (ARs) that will produce SEP events in advance is extremely challenging. We investigate the magnetic field environment of AR 12673, including the AR's magnetic configuration, the surrounding field configuration in the vicinity of the AR, the decay index profile, and the footpoints of Earth-connected magnetic field, around the time of four eruptive events. Two of the eruptive events are SEP-productive (2017 September 4 at 20:00~UT and September 6 at 11:56~UT), while two are not (September 4 at 18:05~UT and September 7 at 14:33~UT). We analysed a range of EUV and white-light coronagraph observations along with potential field extrapolations and find that the CMEs associated with the SEP-productive events either trigger null point reconnection that redirects flare-accelerated particles from the flare site to the Earth-connected field and/or have a significant expansion (and shock formation) into the open Earth-connected field. The rate of change of the decay index with height indicates that the region could produce a fast CME ($v >$ 1500~km~s$^{-1}$), which it did during events two and three. The AR's magnetic field environment, including sites of open magnetic field and null points along with the magnetic field connectivity and propagation direction of the CMEs play an important role in the escape and arrival of SEPs at Earth. Other SEP-productive ARs should be investigated to determine whether their magnetic field environment and CME propagation direction are significant in the escape and arrival of SEPs at Earth.

We present numerical simulation results for the propagation of Alfv\'{e}n waves in the charge starvation regime. This is the regime where the plasma density is below the critical value required to supply the current for the wave. We analyze a conservative scenario where Alfv\'{e}n waves pick up charges from the region where the charge density exceeds the critical value and advect them along at a high Lorentz factor. The system consisting of the Alfv\'{e}n wave and charges being carried with it, which we call charge-carrying Alfv\'{e}n wave (CC-AW), moves through a medium with small, but non-zero, plasma density. We find that the interaction between CC-AW and the stationary medium has a 2-stream like instability which leads to the emergence of a strong electric field along the direction of the unperturbed magnetic field. The growth rate of this instability is of order the plasma frequency of the medium encountered by the CC-AW. Our numerical code follows the system for hundreds of wave periods. The numerical calculations suggest that the final strength of the electric field is of order a few percent of the Alfv\'{e}n wave amplitude. Little radiation is produced by the sinusoidally oscillating currents associated with the instability during the linear growth phase. However, in the nonlinear phase, the fluctuating current density produces strong EM radiation near the plasma frequency and limits the growth of the instability.

In this article, we study Coleman bounce in weakly nonlocal theories which are motivated from string field theory. The kinetic term is extended via an infinite series of high-order derivatives, which comes into play at an energy scale $M$, without introducing any new states or ghosts in the mass spectrum. We calculate the bubble nucleation in thin-wall approximation, treating the system in semi-classical manner. We find that the effect of nonlocal scale $M$ in the theory is to suppress the vacuum tunneling rate from false to true vacuum compared to the standard local bouncing scenario. Likewise, we show that as we move further away from the bubble wall, the effects of nonlocality gets reduced and this suppression is significant only around the wall of the nucleated bubble. From our investigations, we conclude that the main effect is due to the fact that the nonlocality smears the solution of the local bubble profile. However, the energy of the bubble wall remains unaffected by the microscopic nonlocal behavior of the theory in the thin-wall approximation. We also discuss the cases for Lee-Wick theories and applications of our result to cosmology.

This work demonstrates that nontopological solitons with large global charges and masses, even above the Planck scale, can form in the early universe and dominate the dark matter abundance. In solitosynthesis, solitons prefer to grow as large as possible under equilibrium dynamics when an initial global charge asymmetry is present. Their abundance is set by when soliton formation via particle fusion freezes out, and their charges are set by the time it takes to accumulate free particles. This work improves the estimation of both quantities, and in particular shows that much larger-charged solitons form than previously thought. The results are estimated analytically and validated numerically by solving the coupled Boltzmann equations. Without solitosynthesis, phase transitions can still form solitons from particles left inside false-vacuum pockets and determine their present-day abundance and properties. Even with zero charge asymmetry, solitons formed in this way can have very large charges on account of statistical fluctuations in the numbers of (anti)particles inside each pocket.

This paper presents a novel encoding method for fine time data of a tapped delay line (TDL) time-to-digital Converter (TDC). It is based on divide-and-conquer strategy, and has the advantage of significantly reducing logic resource utilization while retaining low dead-time performance. Furthermore, the problem of high bubble depth in advanced devices can be resolved with this method. Four examples are demonstrated, which were implemented in a Xilinx Artix-7 Field Programmable Gate Array (FPGA) device, and encoding method presented in this paper was employed to encode fine time data for normal TDL TDC, a half-length delay line TDC, and double-edge and four-edge wave union TDCs. Compared with TDCs from the latest published papers that adopt traditional encoders, the logic utilization of TDCs in this paper were reduced by a factor of 45% to 70% in different situations, while the encoding dead time can be restricted in one clock cycle. Acceptable resolutions of the demonstrated TDCs were also obtained, proving the functionality of the encoding method.

Pulsar timing offers an independent avenue to test general relativity and alternative gravity theories. This requires an understanding of how metric polarizations beyond the familiar transverse tensor ones imprint as a stochastic gravitational wave background and correlate the arrival time of radio pulses from a pair of millisecond pulsars. In this work, we focus on an isotropic stochastic gravitational wave background and present a straightforward, self-contained formalism for obtaining the power spectrum and the overlap reduction function, the relevant physical observable in a pulsar timing array, for generic gravitational degrees of freedom featuring both transverse and longitudinal modes off the light cone. We additionally highlight our consideration of finite pulsar distances, which we find significant in two ways: first, making all the modes well defined, and second, keeping the small scale power that is contained by pulsars of subdegree separations in the sky. We discuss this for tensor, vector, and scalar polarizations, for each one focusing on the angular power spectrum and the overlap reduction function for an isotropic stochastic gravitational wave background. Our results pave the road for an efficient numerical method for examining the gravitational wave induced spatial correlations across millisecond pulsars in a pulsar timing array.

One of the most common assumptions has been made that the pressure inside the star is isotropic in nature. However, the pressure is locally anisotropic in nature which is a more realistic case. In this study, we investigate certain properties of anisotropic neutron stars with the scalar pressure anisotropy model. Different perfect fluid conditions are tested within the star with the relativistic mean-field model equation of states (EOSs). The anisotropic neutron star properties such as mass ($M$), radius ($R$), compactness ($C$), Love number ($k_2$), dimensionless tidal deformability ($\Lambda$), and the moment of inertia ($I$) are calculated. The magnitude of the quantities as mentioned above increases (decreases) with the positive (negative) value of anisotropy except $k_2$ and $\Lambda$. The Universal relation $I-$Love$-C$ is calculated with almost 58 EOSs spans from relativistic to non-relativistic cases. We observed that the relations between them get weaker when we include anisotropicity. With the help of the GW170817 tidal deformability limit and radii constraints from different approaches, we find that the anisotropic parameter is less than 1.0 if one uses the BL model. Using the universal relation and the tidal deformability bound given by the GW170817, we put a theoretical limit for the canonical radius, $R_{1.4}=10.74_{-1.36}^{+1.84}$ km, and the moment of inertia, $I_{1.4} = 1.77_{-0.09}^{+0.17}\times10^{45}$ g cm$^2$ with 90% confidence limit for isotropic stars. Similarly, for anisotropic stars with $\lambda_{\rm BL}=1.0$, the values are $R_{1.4}=11.74_{-1.54}^{+2.11}$ km, $I_{1.4} = 2.40_{-0.08}^{+0.17} \times10^{45}$ g cm$^2$ respectively.

We study the dependence of primordial nuclear abundances on fundamental nuclear observables such as binding energies, scattering lengths, neutron lifetime, \textit{etc.} by varying these quantities. The numerical computations were performed with four publicly available codes, thus facilitating an investigation of the model-dependent (systematic) uncertainties on these dependences. Indeed deviations of the order of a few percent are found. Moreover, accounting for the temperature dependence of the sensitivity of the rates to some relevant parameters leads to a reduction of the sensitivity of the final primordial abundances, which in some cases is appreciable. These effects are considered to be relevant for studies of the dependence of the nuclear abundances on fundamental parameters such as quark masses or couplings underlying the nuclear parameters studied here.

Assuming unitarity, locality, causality, and Lorentz invariance of the, otherwise unknown, UV completion, we derive a new set of constraints on the effective field theory coefficients for the most general, ghost-free Generalized Proca and Proca Nuevo massive vector models. For the Generalized Proca model, we include new interactions that had not been previously considered in the context of positivity bounds and find these additional terms lead to a widened parameter space for the previously considered interactions. Although, the Generalized Proca and Proca Nuevo models are inequivalent, we find interesting analogues between the coefficients parameterizing the two models and the roles they play in the positivity bounds.