Papers for Monday, Apr 10 2023

Compensated isocurvature perturbations (CIPs) are perturbations to the primordial baryon density that are accompanied by dark-matter-density perturbations so that the total matter density is unperturbed. Such CIPs, which may arise in some multi-field inflationary models, can be long-lived and only weakly constrained by current cosmological measurements. Here we show that the CIP-induced modulation of the electron number density interacts with the electron-temperature fluctuation associated with primordial adiabatic perturbations to produce, via the Biermann-battery mechanism, a magnetic field in the post-recombinaton Universe. This magnetic field may be larger than that (produced at second order in the adiabatic-perturbation amplitude) in the standard cosmological model and may provide seeds for galactic dynamos.

Hot giant planets like MASCARA-1 b are expected to have thermally inverted atmospheres, that makes them perfect laboratory for the atmospheric characterization through high-resolution spectroscopy. Nonetheless, previous attempts of detecting the atmosphere of MASCARA-1 b in transmission have led to negative results. In this paper we aim at the detection of the optical emission spectrum of MASCARA-1 b. We used the high-resolution spectrograph PEPSI to observe MASCARA-1 (spectral type A8) near the secondary eclipse of the planet. We cross-correlated the spectra with synthetic templates computed for several atomic and molecular species. We obtained the detection of FeI, CrI and TiI in the atmosphere of MASCARA-1 b with a S/N ~7, 4 and 5 respectively, and confirmed the expected systemic velocity of ~13 km/s and the radial velocity semi-amplitude of MASCARA-1 b of ~200 km/s. The detection of Ti is of particular importance in the context of the recently proposed Ti cold-trapping below a certain planetary equilibrium temperature. We confirm the presence of an the atmosphere around MASCARA-1 b through emission spectroscopy. We conclude that the atmospheric non detection in transmission spectroscopy is due to the high gravity of the planet and/or to the overlap between the planetary track and its Doppler shadow.

We investigate two-temperature accretion flows onto strongly magnetized compact stars. Matter is accreted in the form of an accretion disc upto the disc radius ($r_{\rm d}$), where, the magnetic pressure exceeds both the gas and ram pressure and thereafter the matter is channelled along the field lines onto the poles. We solve the equations of motion self-consistently along the field lines, incorporating radiative processes like bremsstrahlung, synchrotron and inverse-Comptonization. For a given set of constants of motion, the equations of motion do not produce unique transonic solution. Following the second law of thermodynamics the solution with the highest entropy is selected and thereby eliminating the degeneracy in solution. We study the properties of these solutions and obtain corresponding spectra as a function of the magnetic field ($B_*$), spin period ($P$) and accretion rate of the star ($\dot{M}$). A primary shock is always formed just near the surface. The enhanced radiative processes in this post-shock region slows down the matter and it finally settles on the surface of the star. This post-shock region contributes to $\gtrsim 99.99\%$ of the total luminosity obtained from the accretion flow. It is still important to study the full accretion flow because secondary shocks may be present for some combination of $B_*$, $P$ and ${\dot{M}}$ in addition to primary shocks. We find that secondary shocks, if present, produce an extended emission at higher energies in the spectra.

In 2021 and 2022 the hydrogen comae of three long period comets, C/2020 S3 (Erasmus), C/2021 A1 (Leonard) and C/2021 O3 (PanSTARRS) were observed with the Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliosphere Observer (SOHO) satellite. SWAN obtains nearly daily full-sky images of the hydrogen Lyman-alpha distribution of the interstellar hydrogen as it passes through the solar system yielding information about the solar wind and solar ultraviolet fluxes that eats away at it by ionization and charge exchange. The hydrogen comae of comets, when of sufficient brightness, are also observed. Water production rates have been calculated over time for each of these comets. Of particular interest are comet C/2021 O3 (PanSTARRS) which apparently disintegrated a few days before its perihelion at 0.28 au and C/2021 A1 (Leonard) which also disintegrated beginning about 20 days after its perihelion peak. The behavior of comet C/2020 S3 (Erasmus) was more typical without dramatic fading, but still was asymmetric about perihelion, with a more rapid turn on before perihelion and more extended activity well after perihelion.

The formation of the first stars marks a watershed moment in the history of our universe. As the first luminous structures, these stars (also known as Population III, or Pop III stars) seed the first galaxies and begin the process of reionization. We construct an analytic model to self-consistently trace the formation of Pop III stars inside minihalos in the presence of the fluctuating ultraviolet background, relic dark matter-baryon relative velocities from the early universe, and an X-ray background, which largely work to suppress cooling of gas and delay the formation of this first generation of stars. We demonstrate the utility of this framework in a semi-analytic model for early star formation that also follows the transition between Pop III and Pop II star formation inside these halos. Using our new prescription for the criteria allowing Pop III star formation, we follow a population of dark matter halos from $z=50$ through $z=6$ and examine the global star formation history, finding that each process defines its own key epoch: (i) the stream velocity dominates at the highest redshifts ($z\gtrsim30$), (ii) the UV background sets the tone at intermediate times ($30\gtrsim z\gtrsim15$), and (iii) X-rays control the end of Pop III star formation at the latest times ($z\lesssim 15$). In all of our models, Pop III stars continue to form down to $z\sim 7-10$, when their supernovae will be potentially observable with forthcoming instruments. Finally, we identify the signatures of variations in the Pop III physics in the global 21-cm spin-flip signal of atomic hydrogen.

The interaction of supernova ejecta with a surrounding circumstellar medium (CSM) generates a strong shock which can convert the ejecta kinetic energy into observable radiation. Given the diversity of potential CSM structures (arising from diverse mass loss processes such as late-stage stellar outbursts, binary interaction, and winds), the resulting transients can display a wide range of light curve morphologies. We provide a framework for classifying the transients arising from interaction with a spherical CSM shell. The light curves are decomposed into five consecutive phases, starting from the onset of interaction and extending through shock breakout and subsequent shock cooling. The relative prominence of each phase in the light curve is determined by two dimensionless quantities representing the CSM-to-ejecta mass ratio $\eta$, and a breakout parameter $\xi$. These two parameters define four light curve morphology classes, where each class is characterized by the location of shock breakout and the degree of deceleration as the shock sweeps up the CSM. We compile analytic scaling relations connecting the luminosity and duration of each light curve phase to the physical parameters. We then run a grid of radiation hydrodynamics simulations for a wide range of ejecta and CSM parameters to numerically explore the landscape of interaction light curves, and to calibrate and confirm the analytic scalings. We connect our theoretical framework to several case studies of observed transients, highlighting the relevance in explaining slow-rising and superluminous supernovae, fast blue optical transients, and double-peaked light curves.

We illustrate the effect of boundary conditions on the evolution of structure in Fuzzy Dark Matter. Scenarios explored include the evolution of single, ground-state equilibrium solutions of the Schr\"odinger-Poisson system, the relaxation of a Gaussian density fluctuation, mergers of two equilibrium configurations, and the random merger of many solitons. For comparison, each scenario is evolved twice, with isolation boundary conditions and periodic boundary conditions, the two commonly used to simulate isolated systems and structure formation, respectively. Replacing isolation boundary conditions by periodic boundary conditions changes the domain topology and dynamics of each scenario, by affecting the outcome of gravitational cooling. With periodic boundary conditions, the ground-state equilibrium solution and Gaussian fluctuation each evolve toward the single equilibrium solitonic core of the isolated case, but surrounded by a tail, unlike the isolated versions. The case of head-on, binary mergers illustrates additional effects, caused by the pull suffered by the system due to the infinite network of periodic images along each direction of the domain. Binary merger with angular momentum is the first scenario we found in which the tail has a polynomial profile when using a periodic domain. Finally, the 3D merger of many, randomly-placed solitonic cores of different mass makes a solitonic core surrounded by a tail with power-law-like density profile, for periodic boundary conditions, while producing a core with a much sharper fall-off in the isolated case. This suggests that the conclusion of earlier work that the ground-state equilibrium solution is an attractor for the asymptotic state is true even in 3D and for general circumstances, but only if gravitational cooling is able to carry mass and energy off to infinity, which isolation boundary conditions allow, but periodic ones do not.

Power-law distributions have been studied as a significant characteristic of non-linear dissipative systems. Since discovering the power-law distribution of solar flares that was later extended to nano-flares and stellar flares, it has been widely accepted that different scales of flares share the same physical process. Here, we present the newly developed Semi-Automated Jet Identification Algorithm (SAJIA) and its application for detecting more than 1200 off-limb solar jets during Solar Cycle 24. Power-law distributions have been revealed between the intensity/energy and frequency of these events, with indices found to be analogous to those for flares and coronal mass ejections (CMEs). These jets are also found to be spatially and temporally modulated by the solar cycle forming a butterfly diagram in their latitudinal-temporal evolution, experiencing quasi-annual oscillations in their analysed properties, and very likely gathering in certain active longitudinal belts. Our results show that coronal jets display the same nonlinear behaviour as that observed in flares and CMEs, in solar and stellar atmospheres, strongly suggesting that they result from the same nonlinear statistics of scale-free processes as their counterparts in different scales of eruptive events. Although these jets, like flares and other large-scale dynamic phenomena, are found to be significantly modulated by the solar cycle, their corresponding power-law indices still remain similar.

The gravitational wave observations have revealed four emerging peaks in the binary black hole mass distribution suggesting an overproduction of binaries clustered around specific mass values. Although the presence of the first and the third peaks has been attributed to binary black hole formation in star clusters or due to the evolution of stellar binaries in isolation, the second peak, because it lacks significance in the primary mass distribution, has received relatively less attention. In this article, we report that confidence in the second peak depends on the mass parameter we choose to model the population on. Unlike primary mass, when modelled on the chirp mass this peak is significant. We discuss the disparity as a consequence of mass asymmetry in the observations that cluster at the second peak. Finally, we report this asymmetry to be part of a potential trend in the mass ratio distribution which is manifested as a function of the chirp mass, but not as a function of primary mass, when we include the observation GW190814 in our modelling. Chirp mass is not a parameter of astrophysical relevance. Features present in the chirp mass, but not in the primary mass, are relatively difficult to explain and expected to garner significant interest.

We investigate the roles of major and minor mergers during Brightest Cluster Galaxy (BCG) assembly using surface brightness profiles, line indices, and fundamental plane relations. Based on our own sample (Kluge et al.) and consistently reanalyzed SDSS data, we find that BCGs and luminous normal Ellipticals (LNEs) have similar central velocity dispersions, central absorption line strengths, and central surface brightnesses. However, BCGs are more luminous due to their much larger radial extent. These properties result in a flattening of the Faber-Jackson and Mg$_{\rm b}$-luminosity relations above 10$^{10.6}$ L$_{\odot,g'}$. We use this effect to estimate an amount of 60-80% of accreted and merged light in BCGs relative to LNEs, which agrees with results from cosmological simulations. We determine the contribution of this excess light (EL) at each radius from the difference between the surface flux profiles of BCGs and LNEs. It is small in the center but increases steeply to 50% at already $\sim$3 kpc radius. The shape of these profiles suggests that BCGs could be formed from LNEs in 3 major merger processes. This is also consistent with the mild increase of the S\'ersic indices from $n\approx4$ to $n\approx6$, as confirmed in merger simulations. We note that minor mergers cannot be the dominant origin of the BCG's EL because they deposit too few stars at intermediate radii $r\lesssim20$ kpc. The shape of the EL profile also explains a detected offset of 0.14 dex of the fundamental planes for BCGs and LNEs relative to each other.

Neutron stars (NSs) play essential roles in modern astrophysics. Magnetic fields and spin periods of newborn (zero age) NSs have large impact on the further evolution of NSs, which are however poorly explored in observation due to the difficulty of finding newborn NSs. In this work, we aim to infer the magnetic fields and spin periods (Bi and Pi) of zero-age NSs from the observed properties of NS population. We select non-accretion NSs (NANSs) whose evolution is solely determined by magnetic dipole radiation. We find that both Bi and Pi can be described by log-normal distribution and the fitting sensitively depends on our parameters.

Context. Recent observational data show that the Milky Way (MW) galaxy contains about 170 globular clusters (GCs). A fraction of them is likely formed in dwarf galaxies accreted onto the MW in the past, while the remaining of clusters are formed in-situ. Therefore, different parameters, including orbits, of the globular clusters is a valuable tool for studying the Milky Way evolution. However, since the evolution of the 3D mass distribution of the MW is poorly constrained, the orbits of the clusters are usually calculated in static potentials. Aims. In this work, we study the evolution of the GCs in several external potentials, where we aim to quantify the effects of the evolving galaxy potential on the orbits of the GCs. Methods. For the orbits calculation we used five MW-like potentials from IllustrisTNG-100 simulation. The orbits of 159 GCs were integrated using a high-order N-body parallel dynamic code phi-GPU, with initial conditions obtained from recent Gaia DR3 catalogues. Results. We provide a classification of the GCs orbits according to their 3D shapes and association with different components of the MW (disk, halo, bulge). We also found that the globular clusters in the external potentials have roughly similar energy-angular momentum distributions at the present time. However, both total energy and total angular momentum of the GCs are not conserved due to time-varying nature of the potentials. In some extreme cases, the total energy can change up to 40% (18 objects) over the last 5 Gyr of evolution. We found that the in-situ formed GCs are less affected by the evolution of the TNG potentials as compared to the clusters which are likely formed ex-situ. Therefore, our results suggest that time-varying potentials significantly affect the orbits of the GC, thus making it vital for understanding the formation of the MW.

The nonlinear interaction between electromagnetic waves and plasmas attracts significant attention in astrophysics because it can affect the propagation of Fast Radio Bursts (FRBs) -- luminous millisecond-duration pulses detected at radio frequency. The filamentation instability (FI) -- a type of nonlinear wave-plasma interaction -- is considered to be dominant near FRB sources, and its nonlinear development may also affect the inferred dispersion measure of FRBs. In this paper, we carry out fully kinetic particle-in-cell simulations of the FI in unmagnetized pair plasmas. Our simulations show that the FI generates transverse density filaments, and that the electromagnetic wave propagates in near vacuum between them, as in a waveguide. The density filaments keep merging until force balance between the wave ponderomotive force and the plasma pressure gradient is established. We estimate the merging timescale and discuss the implications of filament merging for FRB observations.

Reactions between cyano radical and aromatic hydrocarbons are believed to be important pathways for the formation of aromatic nitriles in the interstellar medium (ISM) including those identified in the Taurus molecular cloud (TMC-1). Aromatic nitriles might participate in the formation of polycyclic aromatic nitrogen containing hydrocarbons (PANHs) in Titan's atmosphere. Here, ab initio kinetics simulations reveal a high efficiency of $\rm \sim10^{-10}~cm^{3}~s^{-1}$ and the competition of the different products of 30-1800 K and $10^{-7}$-100 atm of the CN + toluene reaction. In the star-forming region of TMC-1 environment, the product yields of benzonitrile and tolunitriles for CN reacting with toluene may be approximately 17$\%$ and 83$\%$, respectively. The detection of main products, tolunitriles, can serve as proxies for the undetected toluene in the ISM due to their much larger dipole moments. The competition between bimolecular and unimolecular products is extremely intense under the warmer and denser PANH forming region of Titan's stratosphere. The computational results show that the fractions of tolunitriles, adducts, and benzonitrile are 19$\%$-68$\%$, 15$\%$-64$\%$ and 17$\%$, respectively, at 150-200 K and 0.0001-0.001 atm (Titan's stratosphere). Then, benzonitrile and tolunitriles may contribute to the formation of PANHs by consecutive $\rm C_{2}H$ additions. Kinetic information of aromatic nitriles for the CN + toluene reaction calculated here helps to explain the formation mechanism of polycyclic aromatic hydrocarbons (PAHs) or PANHs under different interstellar environments and constrains corresponding astrochemical models.

Atmospheric seeing is one of the most important parameters for evaluating and monitoring an astronomical site. Moreover, being able to predict the seeing in advance can guide observing decisions and significantly improve the efficiency of telescopes. However, it is not always easy to obtain long-term and continuous seeing measurements from a standard instrument such as differential image motion monitor (DIMM), especially for those unattended observatories with challenging environments such as Dome A, Antarctica. In this paper, we present a novel machine learning-based framework for estimating and predicting seeing at a height of 8 m at Dome A, Antarctica, using only the data from a multi-layer automated weather station (AWS). In comparison with DIMM data, our estimate has a root mean square error (RMSE) of 0.18 arcsec, and the RMSE of predictions 20 minutes in the future is 0.12 arcsec for the seeing range from 0 to 2.2 arcsec. Compared with the persistence, where the forecast is the same as the last data point, our framework reduces the RMSE by 37 percent. Our method predicts the seeing within a second of computing time, making it suitable for real-time telescope scheduling.

We present a complete analysis of the individual components of the ABC visual triple system HIP 32475. AB pair was discovered during the Hipparcos mission, with a separation of 412 mas. Later, in 2015, a third component was added to the system by discovering it at a small angular distance from B. In our analysis, we follow Al-Wardat's method for analyzing binary and multiple stellar systems, which is a computational spectrophotometric method. Using estimated parameters, the components' positions on the H-R diagram, evolutionary tracks, and isochrones are defined. Depending on the analysis, we estimate the age of the system as 1.259 Gyr with a metallicity of $Z=0.019$. The results show that component A started to evolve from the main sequence to the sub-giants stage, while components B and C are still in the main sequence stage.

[Abridged] The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around the young stars MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the discminer framework, in three bright CO lines, 12CO, 13CO and C18O $J=2-1$, to provide a comprehensive summary of the kinematic and column density substructures that planets might be actively sculpting in these discs. We find convincing evidence for the presence of four giant planets located at wide orbits in three of the discs in the sample: two around HD 163296, one in MWC 480, and one in AS 209. One of the planet candidates in HD 163296, P94, originally proposed by Izquierdo et al. (2022) using lower velocity resolution 12CO data, is confirmed and linked to localised velocity and line width perturbations in 13CO and C18O too. We highlight that line widths are also powerful tracers of planet-forming sites as they are sensitive to turbulent motions triggered by planet-disc interactions. In MWC 480, we identified non-axisymmetric line width enhancements around the radial separation of candidate planet-driven buoyancy spirals proposed by Teague et al. (2021), which we used to narrow the location of the planet candidate to an orbital radius of $R=245$ au and $\rm{PA}=193^\circ$. In the disc of AS 209, we found excess 12CO line widths centred at $R=210$ au, $\rm{PA}=151^\circ$, spanning around the immediate vicinity of the circumplanetary disc candidate proposed by Bae et al. (2022), which further supports its presence. Our simultaneous analysis of multiple tracers and observables aims to lay the groundwork for robust studies of molecular line properties focused on the search for young planets in discs.

We propose to use low-rank matrix approximation using the component-wise L1-norm for direct imaging of exoplanets. Exoplanet detection is a challenging task for three main reasons: (1) the host star is several orders of magnitude brighter than exoplanets, (2) the angular distance between exoplanets and star is usually very small, and (3) the speckles are very similar to the exoplanet signal both in shape and intensity. First, we empirically examine the statistical noise assumptions of the models, second, we evaluate the performance of the proposed L1 low-rank approximation (L1-LRA) algorithm based on visual comparisons and receiver operating characteristic (ROC) curves. We compare the results of the L1-LRA with the widely used truncated singular value decomposition (SVD) based on the L2 norm in two different annuli, one close to the star and one far away.

Magnetoacoustic oscillations are nowadays routinely observed in various regions of the solar corona. This allows them to be used as means of diagnosing plasma parameters and processes occurring in it. Plasma diagnostics, in turn, requires a sufficiently reliable MHD model to describe the wave evolution. In our paper, we focus on obtaining the exact analytical solution to the problem of the linear evolution of standing slow magnetoacoustic (MA) waves in coronal loops. Our consideration of the properties of slow waves is conducted using the infinite magnetic field assumption. The main contribution to the wave dynamics in this assumption comes from such processes as thermal conduction, unspecified coronal heating, and optically thin radiation cooling. In our consideration, the wave periods are assumed to be short enough so that the thermal misbalance has a weak effect on them. Thus, the main non-adiabatic process affecting the wave dynamics remains thermal conduction. The exact solution of the evolutionary equation is obtained using the Fourier method. This means that it is possible to trace the evolution of any harmonic of the initial perturbation, regardless of whether it belongs to entropy or slow mode. We show that the fraction of energy between entropy and slow mode is defined by the thermal conduction and coronal loop parameters. It is shown for which parameters of coronal loops it is reasonable to associate the full solution with a slow wave, and when it is necessary to take into account the entropy wave. Furthermore, we obtain the relationships for the phase shifts of various plasma parameters applicable to any values of harmonic number and thermal condition coefficient. In particular, it is shown that the phase shifts between density and temperature perturbations for the second harmonic of the slow wave vary between $\pi/2$ to 0, but are larger than for the fundamental harmonic.

Triumvirate is a Python/C++ package for measuring the three-point clustering statistics in large-scale structure (LSS) cosmological analyses. Given a catalogue of discrete particles (such as galaxies) with their spatial coordinates, it computes estimators of the multipoles of the three-point correlation function, also known as the bispectrum in Fourier space, in the tri-polar spherical harmonic (TripoSH) decomposition proposed by Sugiyama et al. (2019). The objective of Triumvirate is to provide efficient end-to-end measurement of clustering statistics which can be fed into downstream galaxy survey analyses to constrain and test cosmological models. To this end, it builds upon the original algorithms in the hitomi code developed by Sugiyama et al. (2018, 2019), and supplies a user-friendly interface with flexible input/output (I/O) of catalogue data and measurement results, with the built program configurable through external parameter files and tracked through enhanced logging and warning/exception handling. For completeness and complementarity, methods for measuring two-point clustering statistics are also included in the package.

We propose a new wide-field imaging method that exploits the Localized Surface Plasmon Resonance phenomenon to produce super-resolution images with an optical microscope equipped with a custom design polarization analyzer module. In this paper we describe the method and apply it to the analysis of low-energy carbon ion tracks implanted in a nuclear emulsion film. The result is then compared with the measurements of the same tracks carried out at an electronic microscope. The images set side by side show their close similarity. The resolution achieved with the current microscope setup is estimated to be about 60 nm.

When coupled to electromagnetism via a Chern-Simons interaction, axion-like particles (ALP) produce a rotation of the plane of linear polarization of photons known as cosmic birefringence. Recent measurements of cosmic birefringence obtained from the polarization of the cosmic microwave background (CMB) hint at the existence of an isotropic birefringence angle of $\beta\approx 0.3^\circ$, currently excluding $\beta=0$ with a statistical significance of $3.6\sigma$. Were such measurement to be confirmed as a cosmological signal, CMB information alone could constrain the ALP parameter space for masses $m_\phi\lesssim 10^{-27}$eV and axion-photon coupling constants $g_{\phi\gamma}\gtrsim 10^{-20}$GeV$^{-1}$.

Here we present the first optical photometric monitoring results of a sample of twelve newly discovered blazars from the ICRF - Gaia CRF astrometric link. The observations were performed from April 2013 until August 2019 using eight telescopes located in Europe. For a robust test for the brightness and colour variability, we use Abbe criterion and F-test. Moreover, linear fittings are performed to investigate the relation in the colour-magnitude variations of the blazars. Variability was confirmed in the case of 10 sources; two sources, 1429+249 and 1556+335 seem to be possibly variable. Three sources (1034+574, 1722+119, and 1741+597) have displayed large amplitude brightness change of more than one magnitude. We found that the seven sources displayed bluer-when-brighter variations, and one source showed redder-when-brighter variations. We briefly explain the various AGN emission models which can explain our results.

We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets.

The origin of the recently discovered new class of transients, X-ray quasi-periodic eruptions (QPEs), remains a puzzle. Due to their periodicity and association with active galactic nuclei (AGN), it is natural to relate these eruptions to stars or compact objects in tight orbits around supermassive black holes (SMBHs). In this paper, we predict the properties of emission from bow shocks produced by stars crossing AGN disks, and compare them to the observed properties of QPEs. We find that when a star's orbit is retrograde and has a low inclination ($\lesssim 20^\circ$) with respect to the AGN disk, the breakout emission from the bow shock can explain the observed duration ($\sim$ hours) and X-ray luminosity ($\sim$few$\times10^{42}~{\rm erg~s^{-1}}$) of QPEs. This model can further explain various observed features of QPEs, such as their complex luminosity evolution, the gradual decline of luminosity of the flares over several years, the evolution of the hardness ratio, the modulation of the luminosity during quiescent phases, and the preference of the central SMBHs to have low masses.

We applied principal component analysis (PCA) to the study of five ground level enhancement (GLE) of cosmic ray (CR) events. The nature of the multivariate data involved makes PCA a useful tool for this study. A subroutine program written and implemented in R software environment generated interesting principal components. Analysis of the results shows that the method can distinguish between neutron monitors (NMs) that observed Forbush decrease (FD) from those that observed GLE at the same time. The PCA equally assigned NMs with identical signal counts with the same correlation factor (r) and those with close r values equally have a close resemblance in their CR counts. The results further indicate that while NMs that have the same time of peak may not have the same r, most NMs that had the same r also had the same time of peak. Analyzing the second principal components yielded information on the differences between NMs having opposite but the same or close values of r. NMs that had the same r equally had the tendency of being in close latitude.

The recent discovery of a radio-emitting neutron star with an ultralong spin period of 76 s, PSR J0901-4046, raises a fundamental question on how such a slowly rotating magnetized object can be active in the radio band. A canonical magnetic field of $1.3\times10^{14}$ G estimated from the pulsar period and its time derivative is wholly insufficient for PSR J0901-4046 to operate. Consideration of a magnetic inclination angle of $10^\circ$ estimated from the pulse width gives a higher magnetic field of $1.5\times10^{15}$ G, which is still an order of magnitude lower than the necessary minimum of $2.5\times10^{16}$ G following from the death line for radio pulsars with magnetic fields exceeding the critical value $4.4\times10^{13}$ G. We show that if the observed microstructure of single pulses reflects relativistic beaming, the inferred surface magnetic field appears to be $3.2\times10^{16}$ G, and without this assumption it is no less than $2.7\times10^{16}$ G, which explains the existence of radio emission from PSR J0901-4046. This estimation makes PSR J0901-4046 a radio pulsar with the strongest magnetic field known and is a sign that PSR J0901-4046 slows down not by magnetic-dipole radiation, but rather by an electric current of 56 MA, when rotational energy is expended in accelerating charged particles over the polar cap.

A key assumption of the standard cosmological model is that the temperature of the cosmic microwave background (CMB) radiation scales with cosmological redshift $z$ as $T_{\rm CMB}(z) \propto (1+z)$ at all times after recombination at $z_\star \simeq 1090$. However, this assumption has only been precisely tested at $z \lesssim 3$. Here, we consider cosmological models with post-recombination reheating (PRR), in which the CMB monopole temperature abruptly increases due to energy injection after last scattering. Such a scenario can potentially resolve tensions between inferences of the current cosmic expansion rate (the Hubble constant, $H_0$). We consider an explicit model in which a metastable sub-component of dark matter (DM) decays to Standard Model photons, whose spectral energy distribution is assumed to be close to that of the CMB blackbody. A fit to Planck CMB anisotropy, COBE/FIRAS CMB monopole, and SH0ES distance-ladder measurements yields $H_0 = 71.2 \pm 1.1$ km/s/Mpc, matter fluctuation amplitude $S_8 = 0.774 \pm 0.018$, and CMB temperature increase $\delta T_{\rm CMB} = 0.109^{+0.033}_{-0.044}$ K, which is sourced by DM decay at $z \gtrsim 10$. However, matter density constraints from baryon acoustic oscillation and supernovae data highly constrain this scenario, with a joint fit to all datasets yielding $H_0 = 68.69 \pm 0.35$ km/s/Mpc, $S_8 = 0.8035 \pm 0.0081$, and $\delta T_{\rm CMB} < 0.0342$ K (95% CL upper limit). These bounds can be weakened if additional dark relativistic species are present in the early universe, yielding higher $H_0$. We conclude that current data disfavor models with significant PRR solely through its impact on background and linear-theory observables, completely independent of CMB spectral distortion constraints. However, a small amount of such energy injection could play a role in restoring cosmological concordance.

We investigate whether an accelerating universe can be realized as an asymptotic late-time solution of FLRW-cosmology with multi-field multi-exponential potentials. Late-time cosmological solutions exhibit a universal behavior which enables us to bound the rate of time variation of the Hubble parameter. In string-theoretic realizations, if the dilaton remains a rolling field, our bound singles out a tension in achieving asymptotic late-time cosmic acceleration. Our findings go beyond previous no-go theorems in that they apply to arbitrary multi-exponential potentials and make no specific reference to vacuum or slow-roll solutions. We also show that if the late-time solution approaches a critical point of the dynamical system governing the cosmological evolution, the criterion for cosmic acceleration can be generally stated in terms of a directional derivative of the potential.

Using a quantum-mechanical close-coupling method, we calculate cross sections for fine structure excitation and relaxation of Si and S atoms in collisions with atomic hydrogen. Rate coefficients are calculated over a range of temperatures for astrophysical applications. We determine the temperature-dependent critical densities for the relaxation of Si and S in collisions with H and compare these to the critical densities for collisions with electrons. The present calculation should be useful in modeling environments exhibiting the [S i] 25 {\mu}m and [S i] 57 {\mu}m far-infrared emission lines or where cooling of S and Si by collisions with H are of interest.

We calculate the $S$-factor for proton-proton fusion using chiral effective field theory interactions and currents. By performing order-by-order calculations with a variety of chiral interactions that are regularized and calibrated in different ways, we assess the uncertainty in the $S$-factor from the truncation of the effective field theory expansion and from the sensitivity of the $S$-factor to the short-distance axial current determined from three- and four-nucleon observables. We find that $S(0)=(4.100\pm0.019\mathrm{(syst)}\pm0.013\mathrm{(stat)}\pm0.008(g_A))\times10^{-23}~\mathrm{MeV\,fm}^2\,,$ where the three uncertainties arise, respectively, from the truncation of the effective field theory expansion, use of the two-nucleon axial current fit to few-nucleon observables and variation of the axial coupling constant within the recommended range.

We study spin and flavor oscillations of astrophysical neutrinos under the influence of external fields in curved spacetime. First, we consider spin oscillations in case of neutrinos gravitationally scattered off a rotating supermassive black hole surrounded by a thin magnetized accretion disk. We find that the gravitational interaction only does not result in the spin-flip of scattered ultrarelativistic neutrinos. Realistic magnetic fields lead to the significant reduction of the observed flux of neutrinos possessing reasonable magnetic moments. Second, we study neutrino flavor oscillations in stochastic gravitational waves (GWs). We derive the effective Hamiltonian for neutrinos interacting with a plane GW having an arbitrary polarization. Then, we consider stochastic GWs with arbitrary correlators of amplitudes. The equation for the density matrix for neutrino oscillations is solved analytically and the probabilities to detect certain neutrino flavors are derived. We find that the interaction of neutrinos, emitted by a core-collapsing supernova, with the stochastic GW background results in the several percent change of the neutrino fluxes. The observability of the predicted effects is discussed.

Nominally anhydrous minerals (NAMs) are the primary carriers of water in rocky planet mantles. Therefore, studying water solubilities of major NAMs in the mantle can help us estimate the water storage capacities of rocky planet mantles and indirectly constrain the actual water contents of their interiors. By using data science methods such as statistics and statistical learning algorithms, in this paper, current modeling studies on the mantle water storage capacities of Earth, Mars, and exoplanets have been introduced and summarized. Firstly, the thermodynamic model for mantle water storage capacity has been reviewed. Then, based on the two case studies on Earth and Mars, how to translate atomic-scale experimental data of water solubility and their measurement errors into planetary-scale models of mantle water storage capacity has been explored by using robust regression, Monte Carlo methods, and bootstrap aggregation algorithms. Thirdly, how the large sample data from the exoplanet observational campaigns can help us understand the statistical properties of the mantle water storage capacities of rocky exoplanets has been introduced. Finally, the application limitations of data science methods in mineral physics research have been discussed, and how to better combine statistics and statistical algorithms with mineral physics data research has been prospected.

Collisionless and weakly collisional plasmas often exhibit non-thermal quasi-equilibria. Among these quasi-equilibria, distributions with power-law tails are ubiquitous. It is shown that the statistical-mechanical approach originally suggested by Lynden-Bell (1967) can easily recover such power-law tails. Moreover, we show that, despite the apparent diversity of Lynden-Bell equilibria, a generic form of the equilibrium distribution at high energies is a `hard' power-law tail $\propto \varepsilon^{-2}$, where $\varepsilon$ is the particle energy. The shape of the `core' of the distribution, located at low energies, retains some dependence on the initial condition but it is the tail (or `halo') that contains most of the energy. Thus, a degree of universality exists in collisionless plasmas.

Using simple, intuitive arguments, we discuss the expected accuracy with which astrophysical parameters can be extracted from an observed gravitational wave signal. The observation of a chirp like signal in the data allows for measurement of the component masses and aligned spins, while measurement in three or more detectors enables good localization. The ability to measure additional features in the observed signal -- the existence or absence of power in i) the second gravitational wave polarization, ii) higher gravitational wave multipoles or iii) spin-induced orbital precession -- provide new information which can be used to significantly improve the accuracy of parameter measurement. We introduce the simple-pe algorithm which uses these methods to generate rapid parameter estimation results for binary mergers. We present results from a set of simulations, to illustrate the method, and compare results from simple-pe with measurements from full parameter estimation routines. The simple-pe routine is able to provide initial parameter estimates in a matter of CPU minutes, which could be used in real-time alerts and also as input to significantly accelerate detailed parameter estimation routines.