Papers for Friday, Sep 09 2022

We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improved zeropoints for the ERS filters and show that alternate zeropoints derived by the community also improve the calibration significantly. We also find that the detector offsets appear to be time variable by up to at least 0.1 mag.

This paper evaluates the quality of CME analysis that has been undertaken with the rare Faraday rotation observation of an eruption. Exploring the capability of the FETCH instrument hosted on the MOST mission, a four-satellite Faraday rotation radio sounding instrument deployed between the Earth and the Sun, we discuss the opportunities and challenges to improving the current analysis approaches.

We measure the host galaxy properties of five quasars with $z\sim 1.6 - 3.5$ selected from SDSS and AEGIS, which fall within the JWST/HST CEERS survey area. A PSF library is constructed based on stars in the full field-of-view of the data and used with a 2-dimensional image modeling tool galight to decompose the quasar and its host with multi-band filters available for HST ACS+WFC3 and JWST NIRCAM (12 filters covering HST F606W to JWST F444W). As demonstrated, JWST provides the first capability to detect quasar hosts at $z>3$ and enables spatially-resolved studies of the underlying stellar populations at $z\sim2$ within morphological structures (spiral arms, bar) not possible with HST. Overall, we find quasar hosts to be disk-like, lack merger signatures, and have sizes generally more compact than typical star-forming galaxies at their respective stellar, thus in agreement with results at lower redshifts. The fortuitous face-on orientation of SDSSJ1420+5300A at $z = 1.646$ enables us to find higher star formation and younger ages in the central $2-4$ kpc region relative to the outskirts, which may help explain the relatively compact nature of quasar hosts and pose a challenge to AGN feedback models.

Gamma-ray bursts (GRBs) are flashes of high-energy radiation arising from energetic cosmic explosions. Bursts of long (>2 s) duration are produced by the core-collapse of massive stars, those of short (< 2 s) duration by the merger of two neutron stars (NSs). A third class of events with hybrid high-energy properties was identified, but never conclusively linked to a stellar progenitor. The lack of bright supernovae rules out typical core-collapse explosions, but their distance scales prevent sensitive searches for direct signatures of a progenitor system. Only tentative evidence for a kilonova has been presented. Here we report observations of the exceptionally bright GRB211211A that classify it as a hybrid event and constrain its distance scale to only 346 Mpc. Our measurements indicate that its lower-energy (from ultraviolet to near-infrared) counterpart is powered by a luminous (~1E42 erg/s) kilonova possibly formed in the ejecta of a compact binary merger.

Elongated bar-like features are ubiquitous in galaxies, occurring at the centers of approximately two-thirds of spiral disks. Due to gravitational interactions between the bar and the other components of galaxies, it is expected that angular momentum and matter will redistribute over long (Gyr) timescales in barred galaxies. Previous work ignoring the gas phase of galaxies has conclusively demonstrated that bars should slow their rotation over time due to their interaction with dark matter halos. We have performed a simulation of a Milky Way-like galactic disk hosting a strong bar which includes a state-of-the-art model of the interstellar medium and a live dark matter halo. In this simulation the bar pattern does not slow down over time, and instead remains at a stable, constant rate of rotation. This behavior has been observed in previous simulations using more simplified models for the interstellar gas, but the apparent lack of secular evolution has remained unexplained. We propose that the gas phase of the disk and the dark matter halo act in concert to stabilize the bar pattern speed and prevent the bar from slowing down or speeding up. We find that in a Milky Way-like disk, a gas fraction of only about 5% is necessary for this mechanism to operate. Our result naturally explains why nearly all observed bars rotate rapidly and is especially relevant for our understanding of how the Milky Way arrived at its present state.

Exoplanets display incredible diversity, from planetary system architectures around Sun-like stars that are very different to our Solar System, to planets orbiting post-main sequence stars or stellar remnants. Recently the B-star Exoplanet Abundance STudy (BEAST) reported the discovery of at least two super-Jovian planets orbiting massive stars in the Sco Cen OB association. Whilst such massive stars do have Keplerian discs, it is hard to envisage gas giant planets being able to form in such hostile environments. We use N-body simulations of star-forming regions to show that these systems can instead form from the capture of a free-floating planet, or the direct theft of a planet from one star to another, more massive star. We find that this occurs on average once in the first 10Myr of an association's evolution, and that the semimajor axes of the hitherto confirmed BEAST planets (290 and 556au) are more consistent with capture than theft. Our results lend further credence to the notion that planets on more distant (>100au) orbits may not be orbiting their parent star.

The gas surface density profile of protoplanetary disks is one of the most fundamental physical properties to understand planet formation. However, it is challenging to determine the surface density profile observationally, because the H$_2$ emission cannot be observed in low-temperature regions. We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the \co line toward the protoplanetary disk around TW Hya and discovered extremely broad line wings due to the pressure broadening. In conjunction with a previously reported optically thin CO isotopologue line, the pressure broadened line wings enabled us to directly determine the midplane gas density for the first time. The gas surface density at $\sim5$ au from the central star reaches $\sim 10^3\ {\rm g\ cm^{-2}}$, which suggests that the inner region of the disk has enough mass to form a Jupiter-mass planet. Additionally, the gas surface density drops at the inner cavity by $\sim2$ orders of magnitude compared to outside the cavity. We also found a low CO abundance of $\sim 10^{-6}$ with respect to H$_2$, even inside the CO snowline, which suggests conversion of CO to less volatile species. Combining our results with previous studies, the gas surface density jumps at $r\sim 20$ au, suggesting that the inner region ($3<r<20$ au) might be the magnetorotational instability dead zone. This study sheds light on direct gas-surface-density constraint without assuming the CO/H$_2$ ratio using ALMA.

We present a sample of four emission-line galaxies at $z=6.11-6.35$ that were serendipitously discovered using the commissioning data for the JWST/NIRCam wide-field slitless spectroscopy (WFSS) mode. One of them (at $z=6.11$) has been reported previously while the others are new discoveries. These sources are selected by the secure detections of both [O III] $\lambda$5007 and H$\alpha$ lines with other fainter lines tentatively detected in some cases (e.g., [O II] $\lambda$3727, [O III] $\lambda$4959 and [N II] $\lambda$6583). In the [O III]/H$\beta$ - [N II]/H$\alpha$ Baldwin-Phillips-Terlevich diagram, these galaxies occupy the same parameter space as that of $z\sim2$ star-forming galaxies, indicating that they have been enriched rapidly to sub-solar metallicities ($\sim$0.6 $Z_{\odot}$), similar to galaxies with comparable stellar masses at much lower redshifts. The detection of strong H$\alpha$ lines suggests a higher ionizing photon production efficiency within galaxies in the early Universe. We find brightening of the [O III] $\lambda$5007 line luminosity function (LF) from $z=3$ to 6, and no or weak redshift evolution of the H$\alpha$ line LF from $z=2$ to 6. Both LFs are under-predicted at $z\sim6$ by a factor of $\sim$10 in certain cosmological simulations. This further indicates a global Ly$\alpha$ photon escape fraction of 5-7% at $z\sim6$, much lower than previous estimates through the comparison of the UV-derived star-formation rate density and Ly$\alpha$ luminosity density. Our sample recovers $88^{+164}_{-57}$% of $z=6.0-6.6$ galaxies in the survey volume with stellar masses greater than $5\times10^8$ $M_{\odot}$, suggesting the ubiquity of strong H$\alpha$ and [O III] line emitters in the Epoch of Reionization, which will be further uncovered in the era of JWST.

Dipper stars are a classification of young stellar objects that exhibit dimming variability in their light curves, dropping in brightness by 10-50%, likely induced by occultations due to circumstellar disk material. This variability can be periodic, quasi-periodic, or aperiodic. Dipper stars have been discovered in young stellar associations via ground-based and space-based photometric surveys. We present the detection and characterization of the largest collection of dipper stars to date: 293 dipper stars, including 234 new dipper candidates. We have produced a catalog of these targets, which also includes young stellar variables that exhibit predominately bursting-like variability and symmetric variability (equal parts bursting and dipping). The total number of catalog sources is 414. These variable sources were found in a visual survey of TESS light curves, where dipping-like variability was observed. We found a typical age among our dipper sources of <5 Myr, with the age distribution peaking at ~2 Myr, and a tail of the distribution extending to ages older than 20 Myr. Regardless of the age, our dipper candidates tend to exhibit infrared excess, which is indicative of the presence of disks. TESS is now observing the ecliptic plane, which is rich in young stellar associations, so we anticipate many more discoveries in the TESS dataset. A larger sample of dipper stars would enhance the census statistics of light curve morphologies and dipper ages.

We study molecular outflows in a sample of 25 nearby (z< 0.17, d<750 Mpc) ULIRG systems (38 individual nuclei) as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) survey, using ~400 pc (0.1-1.0" beam FWHM) resolution ALMA CO(2-1) observations. We used a spectro-astrometry analysis to identify high-velocity (> 300 km/s) molecular gas disconnected from the galaxy rotation, which we attribute to outflows. In 77% of the 26 nuclei with $\log L_{IR}/L_{\odot}>11.8$, we identifid molecular outflows with an average $v_{out}= 490$ km/s, outflow masses $1-35 \times 10^7$ $M_{\odot}$, mass outflow rates $\dot{M}_{out}=6-300$ $M_{\odot}$ yr$^{-1}$, mass-loading factors $\eta = \dot{M}_{out}/SFR = 0.1-1$, and an average outflow mass escape fraction of 45%. The majority of these outflows (18/20) are spatially resolved with radii of 0.2-0.9 kpc and have short dynamical times ($t_{dyn}=R_{out}/v_{out}$) in the range 0.5-2.8 Myr. The outflow detection rate is higher in nuclei dominated by starbursts (SBs, 14/15=93%) than in active galactic nuclei (AGN, 6/11=55%). Outflows perpendicular to the kinematic major axis are mainly found in interacting SBs. We also find that our sample does not follow the $\dot{M}_{out}$ versus AGN luminosity relation reported in previous works. In our analysis, we include a sample of nearby main-sequence galaxies (SFR = 0.3-17 $M_{\odot}$ yr$^{-1}$) with detected molecular outflows from the PHANGS-ALMA survey to increase the $L_{IR}$ dynamic range. Using these two samples, we find a correlation between the outflow velocity and the SFR, as traced by $L_{IR}$ ($v_{out} \propto SFR^{0.25\pm0.01})$, which is consistent with what was found for the atomic ionised and neutral phases. Using this correlation, and the relation between $M_{out}/R_{out}$ and $v_{out}$, we conclude that these outflows are likely momentum-driven.

We present new broadband X-ray observations of the $z \sim 2.5$ lensed quasar 2MASS J1042+1641, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$ to provide coverage of the X-ray spectrum over the 0.3$-$40 keV bandpass in the observed frame, corresponding to the $\sim$1$-$140 keV band in the rest-frame of 2MASS J1042+1641. The X-ray data show clear evidence for strong (but still Compton-thin) X-ray absorption, $N_{\rm{H}} \sim 3-4 \times 10^{23}$ cm$^{-2}$, in addition to significant reprocessing by Compton-thick material that must lie away from our line-of-sight to the central X-ray source. We test two different interpretations for the latter: first that the reprocessing occurs in a classic AGN torus, as invoked in unification models, and second that the reprocessing occurs in the accretion disc. Both models can successfully reproduce the observed spectra, and both imply that the source is viewed at moderately low inclinations ($i < 50^{\circ}$) despite the heavy line-of-sight absorption. Combining the X-ray data with infrared data from $WISE$, the results seen from 2MASS J1042+1641 further support the recent suggestion that large X-ray and IR surveys may together be able to identify good lensed quasar candidates in advance of detailed imaging studies.

NGC 6441 is a bulge globular cluster with an unusual horizontal branch morphology and a rich population of RR Lyrae (RRL) and Type II Cepheid (T2C) variables that is unexpected for its relatively high metallicity. We present near-infrared (NIR, $JHK_s$) time-series observations of 42 RRL, 8 T2Cs, and 10 eclipsing binary candidate variables in NGC 6441. The multi-epoch observations were obtained using the FLAMINGOS-2 instrument on the 8-m Gemini South telescope. Multi-band data are used to investigate pulsation properties of RRL and T2C variables including their light curves, instability strip, period--amplitude diagrams, and period--luminosity and period--wesenheit relations (PLRs and PWRs) in $JHK_s$ filters. NIR pulsation properties of RRL variables are well fitted with theoretical models that have canonical helium content and mean-metallicity of NGC 6441 ([Fe/H]$=-0.44\pm0.07$~dex). The helium-enhanced RRL models predict brighter NIR magnitudes and bluer colors than the observations of RRL in the cluster. This suggests that these RRL variables in NGC 6441 are either not significantly helium enhanced as previously thought or the impact of such enhancement is smaller in NIR than the predictions of the pulsation models. We also use theoretical calibrations of RRL period--luminosity--metallicity (PLZ) relations to simultaneously estimate the mean reddening, $E(J-K_s)=0.26\pm0.06$~mag, and the distance, $d=12.67\pm0.09$~kpc, to NGC 6441. Our mean reddening value and the distance are consistent with the independent estimates from the bulge reddening map based on red clump stars and the latest Gaia data, respectively. Our distance and reddening values provide a very good agreement between the PLRs of T2Cs in NGC 6441 and those for RRL and T2Cs in Galactic globular clusters that span a broad range of metallicity.

New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history ($\mathcal{S}(Z,z)$) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for $\mathcal{S}(Z,z)$. Variations of this function can be easily interpreted, because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behaviour well. As an example application, we investigate the effect of systematic variations in the $\mathcal{S}(Z,z)$ parameters on the predicted mass distribution of locally merging binary black holes (BBH). Our main findings are: I) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and II) the low mass end is least affected by these variations. This is promising as it increases our chances to constrain the physics that governs the formation of these objects.

We present a strong lensing analysis of the cluster PSZ1 G311.65-18.48, based on Hubble Space Telescope imaging, archival VLT/MUSE spectroscopy, and Chandra X-ray data. This cool-core cluster (z=0.443) lenses the brightest lensed galaxy known, dubbed the "Sunburst Arc" (z=2.3703), a Lyman continuum (LyC) emitting galaxy multiply-imaged 12 times. We identify in this field 14 additional strongly-lensed galaxies to constrain a strong lens model, and report secure spectroscopic redshifts of four. We measure a projected cluster core mass of M(<250 kpc)=2.93+0.01/-0.02x10^14M_sun. The two least-magnified but complete images of the Sunburst Arc's source galaxy are magnified by ~13x, while the LyC clump is magnified by ~4-80x. We present time delay predictions and conclusive evidence that a discrepant clump in the Sunburst Arc, previously claimed to be a transient, is not variable, thus strengthening the hypothesis that it results from an exceptionally high magnification. A source plane reconstruction and analysis of the Sunburst Arc finds its physical size to be 1x2 kpc, and that it is resolved in three distinct directions in the source plane, 0, 40, and 75 degrees (east of North). We place an upper limit of r <~ 50 pc on the source plane size of unresolved clumps, and r<~ 32 pc for the LyC clump. Finally, we report that the Sunburst Arc is likely in a system of two or more galaxies separated by <~6 kpc in projection. Their interaction may drive star formation and could play a role in the mechanism responsible for the leaking LyC radiation.

We report direction detection constraints on the presence of hidden photon dark matter with masses between 20-30 ueV using a cryogenic emitter-receiver-amplifier spectroscopy setup designed as the first iteration of QUALIPHIDE (QUantum LImited PHotons In the Dark Experiment). A metallic dish sources conversion photons from hidden photon kinetic mixing onto a horn antenna which is coupled to a C-band kinetic inductance traveling wave parametric amplifier, providing for near quantum-limited noise performance. We demonstrate a first probing of the kinetic mixing parameter "chi" to just above 10^-12 for the majority of hidden photon masses in this region. These results not only represent stringent constraints on new dark matter parameter space but are also the first demonstrated use of wideband quantum-limited amplification for astroparticle applications

When the expansion of the universe is dominated by a perfect fluid with equation of state parameter $w$ and a sound speed $c_s$ satisfying $w = c_s^2 \le 1$, the Hubble parameter $H$ and time $t$ satisfy the bound $Ht \ge 1/3$. There has been recent interest in ``ultra-slow" expansion laws with $Ht < 1/3$ (sometimes described as ``fast expanding" models). We examine various models that can produce ultra-slow expansion: scalar fields with negative potentials, barotropic fluids, braneworld models, or a loitering phase in the early universe. Scalar field models and barotropic models for ultra-slow expansion are unstable to evolution toward $w = 1$ or $w \rightarrow \infty$ in the former case and $w \rightarrow \infty$ in the latter case. Braneworld models can yield ultra-slow expansion but require an expansion law beyond the standard Friedman equation. Loitering early universe models can produce a quasi-static expansion phase in the early universe but require an exotic negative-density component. These results suggest that appeals to an ultra-slow expansion phase in the early universe should be approached with some caution, although the loitering early universe may be worthy of further investigation. These results do not apply to ultra-slow contracting models.

The spatial distribution and physical sizes of star forming clumps at the smallest scales provide valuable information on hierarchical star formation (SF). In this context, we report the sites of ongoing SF at ~120 pc along the interacting galaxies in Stephan's Quintet (SQ) compact group using AstroSat-UVIT and JWST data. Since ultraviolet radiation is a direct tracer of recent SF, we identified star forming clumps in this compact group from the FUV imaging which we used to guide us to detect star forming regions on JWST IR images. The FUV imaging reveals star forming regions within which we detect smaller clumps from the higher spatial resolution images of JWST, likely produced by PAH molecules and dust ionised by FUV emission from young massive stars. This analysis reveals the importance of FUV imaging data in identifying star forming regions in the highest spatial resolution IR imaging available.

The VERITAS Imaging Atmospheric Cherenkov Telescope array (IACT) was augmented in 2019 with high-speed focal plane electronics to create a new Stellar Intensity Interferometry (SII) observational capability (VERITAS-SII, or VSII). VSII operates during bright moon periods, providing high angular resolution observations ( < 1 mas) in the B photometric band using idle telescope time. VSII has already demonstrated the ability to measure the diameters of two B stars at 416 nm (Bet CMa and Eps Ori) with < 5% accuracy using relatively short (5 hours) exposures. The VSII instrumentation was recently improved to increase instrumental sensitivity and observational efficiency. This paper describes the upgraded VSII instrumentation and documents the ongoing improvements in VSII sensitivity. The report describes VSII's progress in extending SII measurements to dimmer magnitude stars and improving the VSII angular diameter measurement resolution to better than 1%.

We develop a Bayesian model that jointly constrains receiver calibration, foregrounds and cosmic 21cm signal for the EDGES global 21\,cm experiment. This model simultaneously describes calibration data taken in the lab along with sky-data taken with the EDGES low-band antenna. We apply our model to the same data (both sky and calibration) used to report evidence for the first star formation in 2018. We find that receiver calibration does not contribute a significant uncertainty to the inferred cosmic signal (<1%), though our joint model is able to more robustly estimate the cosmic signal for foreground models that are otherwise too inflexible to describe the sky data. We identify the presence of a significant systematic in the calibration data, which is largely avoided in our analysis, but must be examined more closely in future work. Our likelihood provides a foundation for future analyses in which other instrumental systematics, such as beam corrections and reflection parameters, may be added in a modular manner.

Compact groups have high galaxy densities and low velocity dispersions, and their group members have experienced numerous and frequent interactions during their lifetimes. They provide a unique environment to study the evolution of galaxies. We examined the galaxies types and HI contents in groups to make a study on the galaxy evolution in compact groups. We used the group crossing time as an age indicator for galaxy groups. Our sample is derived from the Hickson Compact Group catalog. We obtained group morphology data from the Hyper-Leda database and the IR classification based on Wide-Field Infrared Survey Explorer (WISE) fluxes from Zucker et al. (2016). By cross-matching the latest released ALFALFA 100% HI source catalog and supplemented by data found in literature, we obtained 40 galaxy groups with HI data available. We confirmed that the weak correlation between HI mass fraction and group crossing time found by Ai & Zhu (2018) in SDSS groups also exists in compact groups. We also found that the group spiral galaxy fraction is correlated with the group crossing time, but the actively star-forming galaxy fraction is not correlated with the group crossing time. These results seem to fit with the hypothesis that the sequential acquisition of neighbors from surrounding larger-scale structures has affected the morphology transition and star formation efficiency in compact groups.

We summarize evidence that multiple supernovae exploded within 100 pc of Earth in the past few Myr. These events had dramatic effects on the heliosphere, compressing it to within ~20 au. We advocate for cross-disciplinary research of nearby supernovae, including on interstellar dust and cosmic rays. We urge for support of theory work, direct exploration, and study of extrasolar astrospheres.

Observations of helium in exoplanet atmospheres may reveal the presence of large gaseous envelopes, and indicate ongoing atmospheric escape. Orell-Miquel et al. (2022) used CARMENES to measure a tentative detection of helium for the sub-Neptune GJ 1214b, with a peak excess absorption reaching over 2% in transit depth at 10830 Angstroms. However, several non-detections of helium had previously been reported for GJ 1214b. One explanation for the discrepancy was contamination of the planetary signal by overlapping telluric absorption- and emission lines. We used Keck/NIRSPEC to observe another transit of GJ 1214b at 10830 Angstroms, at a time of minimal contamination by telluric lines, and did not observe planetary helium absorption. Accounting for correlated noise in our measurement, we place an upper limit on the excess absorption size of 1.22% (95% confidence). We find that the discrepancy between the CARMENES and NIRSPEC observations is unlikely to be caused by using different instruments or stellar activity. It is currently unclear whether the difference is due to correlated noise in the observations, or variability in the planetary atmosphere.

The acceleration of the solar coronal plasma to supersonic speeds is one of the most fundamental yet unresolved problem in heliophysics. Despite the success of Parker's pioneering theory on an isothermal solar corona, the realistic solar wind is observed to be non-isothermal, and the decay of its temperature with radial distance usually can be fitted to a polytropic model. In this work, we use Parker Solar Probe data from the first nine encounters to estimate the polytropic index of solar wind protons. We show that the polytropic index varies between 1.25 and $5/3$ and depends strongly on solar wind speed, faster solar wind on average displaying a smaller polytropic index. We comprehensively analyze the 1D spherically symmetric solar wind model with polytropic index $\gamma \in [1,5/3]$. We derive a closed algebraic equation set for transonic stellar flows, i.e. flows that pass the sound point smoothly. We show that an accelerating wind solution only exists in the parameter space bounded by $C_0/C_g < 1$ and $(C_0/C_g)^2 > 2(\gamma-1)$ where $C_0$ and $C_g$ are the surface sound speed and one half of the escape velocity of the star, and no stellar wind exists for $\gamma > 3/2$. With realist solar coronal temperatures, the observed solar wind with $\gamma \gtrsim 1.25$ cannot be explained by the simple polytropic model. We show that mechanisms such as strong heating in the lower corona that leads to a thick isothermal layer around the Sun and large-amplitude Alfv\'en wave pressure are necessary to remove the constraint in $\gamma$ and accelerate the solar wind to high speeds.

The vertical distribution of trace gases in planetary atmospheres can be obtained with observations of the atmosphere's thermal emission. Inverting radio observations to recover the atmospheric structure, however, is non-trivial, and the solutions are degenerate. We propose a modeling framework to prescribe a vertical distribution of trace gases that combines a thermo-chemical equilibrium model {based on a vertical temperature structure and compare these results to models where ammonia can vary between pre-defined pressure nodes}. To this means we retrieve nadir brightness temperatures and limb-darkening parameters, together with their uncertainties, from the Juno Microwave Radiometer (MWR). We then apply this framework to MWR observations during Juno's first year of operation (Perijove passes 1 - 12) and to longitudinally-averaged latitude scans taken with the upgraded Very Large Array (VLA) (de Pater 2016,2019a). We use the model to constrain the distribution of ammonia between -60$^{\circ}$ and 60$^{\circ}$ latitude and down to 100 bar. We constrain the ammonia abundance to be $340.5^{+34.8}_{-21.2}$ ppm ($2.30^{+0.24}_{-0.14} \times$ solar abundance), and find a depletion of ammonia down to a depth of $\sim$ 20 bar, which supports the existence of processes that deplete the atmosphere below the ammonia and water cloud layers. At the equator we find an increase of ammonia with altitude, while the zones and belts in the mid-latitudes can be traced down to levels where the atmosphere is well-mixed. The latitudinal variation in the ammonia abundance appears to be opposite to that shown at higher altitudes, which supports the existence of a stacked-cell circulation model.

It was recently proposed that exotic particles can trigger a new stellar instability which is analogous to the e-e+ pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ~400 keV--2 MeV and coupling with photons of g_{ag}~10^{-5} GeV^{-1}. It is found that the 56Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10--20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7--2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function.

Shock waves propagating in collisionless heliospheric and astrophysical plasmas have been studied extensively over the decades. One prime motivation is to understand the nonthermal particle acceleration at shocks. Although the theory of diffusive shock acceleration (DSA) has long been the standard for cosmic-ray acceleration at shocks, plasma physical understanding of particle acceleration remains elusive. In this review, we discuss nonthermal electron acceleration mechanisms at quasi-perpendicular shocks, for which substantial progress has been made in recent years. The discussion presented in this review is restricted to the following three specific topics. The first is stochastic shock drift acceleration (SSDA), which is a relatively new mechanism for electron injection into DSA. The basic mechanism, related in-situ observations and kinetic simulations results, and how it is connected with DSA will be discussed. Second, we discuss shock surfing acceleration (SSA) at very high Mach number shocks relevant to young supernova remnants (SNRs). While the original proposal under the one-dimensional assumption is unrealistic, SSA has now been proven efficient by a fully three-dimensional kinetic simulation. Finally, we discuss the current understanding of the magnetized Weibel-dominated shock. Spontaneous magnetic reconnection of self-generated current sheets within the shock structure is an interesting consequence of Weibel-generated strong magnetic turbulence. We argue that high Mach number shocks with both Alfven and sound Mach numbers exceeding 20-40 will likely behave as a Weibel-dominated shock. Despite a number of interesting recent findings, the relative roles of SSDA, SSA, and magnetic reconnection for electron acceleration at collisionless shocks and how the dominant particle acceleration mechanisms change depending on shock parameters remain to be answered.

In this paper, we discussed the multiple vector fields during the inflation era and the inflationary magnetogenesis with multiple vector fields. Instead of a single coupling function in single vector field models, the coupling matrix between vector fields and scalar field which drive the inflation is introduced. The dynamical equations for multiple vector fields are obtained and applied to the inflation era. We discussed three cases for the double-field model. In no mutual-coupling case, one can find that both electric and magnetic spectrum can be scale-invariant at the end of inflation, meanwhile, the strong coupling problem can be avoided. The effect of mutual-coupling between different vector fields is also discussed. We found that weak mutual-coupling can lead to the slightly blue spectrum of the magnetic field. On the other hand, in the strong mutual-coupling case, the scale-invariant magnetic spectrum can also be obtained but the energy density of electromagnetic fields either lead to the backreaction problem or is diluted by inflation.

We report 23 stars having Galactocentric velocities larger than $450~\mathrm{km\,s}^{-1}$ in the final data release of the APOGEE survey. This sample was generated using space velocities derived by complementing the high quality radial velocities from the APOGEE project in Sloan Digital Sky Survey's Data Release 17 (DR17) with distances and proper motions from Gaia early Data Release 3 (eDR3). We analyze the observed kinematics and derived dynamics of these stars, considering different potential models for the Galaxy. We find that three stars could be unbound depending on the adopted potential, but in general all of the stars show typical kinematics of halo stars. The APOGEE DR17 spectroscopic results and Gaia eDR3 photometry are used to assess the stellar parameters and chemical properties of the stars. All of the stars belong to the red giant branch, and, in general, they follow the abundance pattern of typical halo stars. There are a few exceptions that would deserve further analysis through high-resolution spectroscopy. In particular, we identify a high velocity Carbon-Enhanced Metal-Poor (CEMP) star, with Galactocentric velocity of 482 km\,s$^{-1}$. We do not confirm any hypervelocity star in the sample, but this result is very sensitive to the adopted distances, and less sensitive to the Galactic potential.

Intermediate-mass black holes (with $\geq\!10^5\,M_\odot$) are promising candidates for the origin of supermassive black holes (with $\sim\!10^9\,M_\odot$) in the early universe (redshift $z\sim6$). Chon & Omukai (2020) firstly pointed out the direct collapse black hole (DCBH) formation in metal-enriched atomic-cooling halos (ACHs), which relaxes the DCBH formation criterion. On the other hand, Hirano et al. (2021) showed that the magnetic effects promote the DCBH formation in metal-free ACHs. We perform a set of magnetohydrodynamical simulations to investigate star formation in the magnetized ACHs with metallicities $Z/Z_\odot = 0$, $10^{-5}$, and $10^{-4}$. Our simulations show that the mass accretion rate onto the protostars becomes lower in metal-enriched ACHs than that of metal-free ACHs. However, many protostars form from gravitationally and thermally unstable metal-enriched gas clouds. Under such circumstances, the magnetic field rapidly increases as the magnetic field lines wind up due to the spin of protostars. The region with the amplified magnetic field expands outwards due to the orbital motion of protostars and the rotation of the accreting gas. The amplified magnetic field extracts the angular momentum from the accreting gas, promotes the coalescence of the low-mass protostars, and increases the mass growth rate of the primary protostar. We conclude that the magnetic field amplification is always realized in the metal-enriched ACHs regardless of the initial magnetic field strength, which relaxes the DCBH formation criterion. In addition, we find a qualitatively different trend from the previous unmagnetized simulations in that the mass growth rate is maximal for the extremely metal-poor ACHs with $Z/Z_\odot = 10^{-5}$.

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

The Telescope Array (TA) is the largest hybrid cosmic ray detector in the Northern Hemisphere, which observes primary particles in the energy range from 2 PeV to 100 EeV. The main TA detector consists of 507 plastic scintillation counters on a 1.2-km spacing square grid and fluorescence detectors at three stations overlooking the sky above the surface detector array. The TA Low energy Extension (TALE) hybrid detectors, which consists of ten fluorescence telescopes, and 80 infill surface detectors with 400-m and 600-m spacing, has continued to provide stable observations since its construction completion in 2018. The TAx4, a plan to quadruple the detection area of TA is also ongoing. About half of the planned detectors have been deployed, and the current TAx4 continues to operate stably as a hybrid detector. I review the present status of the TA experiment and the recent results on the cosmic-ray anisotropy, mass composition and energy spectrum.

Modifications of General Relativity generically contain additional degrees of freedom that can mediate forces between matter particles. One of the common manifestations of a fifth force in alternative gravity theories is a difference between the gravitational potentials felt by relativistic and non-relativistic particles, also known as ``the gravitational slip''. In contrast, a fifth force between dark matter particles, due to dark sector interaction, does not cause a gravitational slip, making the latter a possible smoking gun of modified gravity. In this article, we point out that a force acting on dark matter particles, as in models of coupled quintessence, would also manifest itself as a measurement of an effective gravitational slip by cosmological surveys of large-scale structure. This is linked to the fact that redshift-space distortions due to peculiar motion of galaxies do not provide a measurement of the true gravitational potential if dark matter is affected by a fifth force. Hence, it is extremely challenging to distinguish a dark sector interaction from a modification of gravity with cosmological data alone. Future observations of gravitational redshift from galaxy surveys can help to break the degeneracy between these possibilities, by providing a direct measurement of the distortion of time. We discuss this and other possible ways to resolve this important question.

After becoming ionized, low-density astrophysical plasmas will begin a process of slow recombination. Models for this still have significant uncertainties. The recombination cannot normally be observed in isolation, because the ionization follows the evolutionary time scale of the ionizing source. Laboratory experiments are unable to reach the appropriate conditions because of the required very long time scales. The extended nebula around the very late helium flash (VLTP) star V4334 Sgr provides a unique laboratory for this kind of study. The sudden loss of the ionizing UV radiation after the VLTP event has allowed the nebula to recombine free from other influences. More than 290 long slit spectra taken with FORS1/2 at the ESO VLT between 2007 and 2022 are used to follow the time evolution of lines of H, He, N, S, O, Ar. Hydrogen and helium lines, representing most of the ionized mass, do not show significant changes. A small increase is seen in [N II] (+2.8 %/yr; significance 2.7 sigma), while we see a decrease in [O III] (-1.96 %/yr; 2.0 sigma). The [S II] lines show a change of +3.0 %/yr; 1.6 sigma). The lines of [S III] and of Ar III] show no significant change. For [S III], the measurement differs from the predicted decrease by 4.5 sigma. A possible explanation is that the fraction of [S IV] and higher is larger than expected. Such an effect could provide a potential solution for the sulfur anomaly in planetary nebulae.

We have been developing the SOI pixel detector ``INTPIX'' for space use and general purpose applications such as the residual stress measurement of a rail and high energy physics experiments. INTPIX is a monolithic pixel detector composed of a high-resistivity Si sensor, a SiO2 insulator, and CMOS pixel circuits utilizing Silicon-On-Insulator (SOI) technology. We have considered the possibility of using INTPIX to observe X-ray polarization in space. When the semiconductor detector is used in space, it is subject to radiation damage resulting from high-energy protons. Therefore, it is necessary to investigate whether INTPIX has high radiation tolerance for use in space. The INTPIX8 was irradiated with 6 MeV protons up to a total dose of 2 krad at HIMAC, National Institute of Quantum Science in Japan, and evaluated the degradation of the performance, such as energy resolution and non-uniformity of gain and readout noise between pixels. After 500 rad irradiation, which is the typical lifetime of an X-ray astronomy satellite, the degradation of energy resolution at 14.4 keV is less than 10%, and the non-uniformity of readout noise and gain between pixels is constant within 0.1%.

Context: The Einasto model has become one of the most popular models for describing the density profile of dark matter haloes. There have been relatively few comprehensive studies on the dynamical structure of the Einasto model, mainly because only a limited number of properties can be calculated analytically. Aims: We want to systematically investigate the photometric and dynamical structure of the family of Einasto models over the entire model parameter space. Methods: We used the SpheCow code to explore the properties of the Einasto model. We systematically investigated how the most important properties change as a function of the Einasto index $n$. We considered both isotropic models and radially anisotropic models with an Osipkov-Merritt orbital structure. Results: We find that all Einasto models with $n<\tfrac12$ have a formal isotropic or Osipkov-Merritt distribution function that is negative in parts of phase space, and hence cannot be supported by such orbital structures. On the other hand, all models with larger values of $n$ can be supported by an isotropic orbital structure, or by an Osipkov-Merritt anisotropy, as long as the anisotropy radius is larger than a critical value. This critical anisotropy radius is a decreasing function of $n$, indicating that less centrally concentrated models allow for a larger degree of radial anisotropy. Conclusions: Studies of the structure and dynamics of models for galaxies and dark matter haloes should not be restricted to completely analytical models. Numerical codes such as SpheCow can help open up the range of models that are systematically investigated. This applies to the Einasto model discussed here, but also to other proposed models for dark matter haloes, including different extensions to the Einasto model.

The Fermi Gamma-ray Space Telescope was launched more than 13 years ago and since then it has dramatically changed our knowledge of the gamma-ray sky. With more than three billions photons from the whole sky, collected in the energy range between 20 MeV and more than 300 GeV, and beyond 6,000 detected sources, LAT observations have been crucial to improving our understanding of particle acceleration and gamma-ray production in astrophysical sources. In this proceeding, I will review recent science highlights from the LAT. I will focus on the recent source catalog release, as well as on the main transient phenomena seen with the LAT with multi-wavelength and multi-messenger connection.

The photometric and spectroscopic investigations of ten contact binaries were presented for the first time. It is discovered that the mass ratios of all the ten targets are smaller than 0.15, they are extremely low mass ratio contact binaries. Seven of them are deep contact binaries, two are medium contact binaries, while only one is a shallow contact system. Five of them show O'Connell effect, and a dark spot on one of the two components can lead to a good fit of the asymmetric light curves. The orbital period studies of the ten binaries reveal that they all exhibit long-term period changes, six of them are increasing, while the others are shrinking. The LAMOST spectra were analyzed by the spectral subtraction method, and all the ten targets exhibit excess emissions in the H$_\alpha$ line, indicating chromospheric activity. The evolutionary states of the two components of the ten binaries were studied, and it is found that their evolutionary states are identical to those of the other contact binaries. Based on the study of the relation between orbital angular momentum and total mass, we discovered the ten systems may be at the late evolutionary stage of a contact binary. The initial masses of the two components and the ages of them were obtained. By calculating the instability parameters, we found that the ten contact binaries are relatively stable at present.

What makes the study of exoplanetary atmospheres so hard is the extraction of its tiny signal from observations, usually dominated by telluric absorption, stellar spectrum and instrumental noise. The High Resolution Spectroscopy has emerged as one of the leading techniques for detecting atomic and molecular species (Birkby 2018), but although it is particularly robust against contaminant absorption in the Earth's atmosphere, the non-stationary stellar spectrum -- in the form of either Doppler shift or distortion of the line profile during planetary transits -- creates a non-negligible source of noise that can alter or even prevent the detection. Recently, significant improvements have been achieved by using 3D, radiative hydrodynamical (RHD) simulations for the star and Global Circulation Models (GCM) for the planet (e.g., Chiavassa & Brogi 2019, Flowers et al. 2019). However, these numerical simulations have been computed independently so far, while acquired spectra are the result of the natural coupling at each phase along the planet orbit. With our work, we aim at generating emission spectra of G,F, and K-type stars and Hot Jupiters and coupling them at any phase of the orbit. This approach is expected to be particularly advantageous for those molecules that are present in both the atmospheres (e.g., CO) and form in the same region of the spectrum, resulting in mixed and overlapped spectral lines. We also present the analysis of transmission spectra of the Hot Saturn WASP-20b, observed in the K-band of the recently upgraded spectrograph CRIRES+ at a resolution ~92, 000 during the first night of the Science Verification of the instrument and that led to a tentative detection of H2O.

Models of exoplanet atmospheres based on Parker-wind density and velocity profiles are a common choice in fitting spectroscopic observations tracing planetary atmospheric escape. Inferring atmospheric properties using these models often results in a degeneracy between the temperature and the mass-loss rate, and thus provides weak constraints on either parameter. We present a framework that can partially resolve this degeneracy by placing more stringent constraints on the expected thermospheric temperature. We use the photoionization code Cloudy within an iterative scheme to compute the temperature structure of a grid of 1D Parker wind models, including the effects of radiative heating/cooling, as well as the hydrodynamic effects (expansion cooling and heat advection). We constrain the parameter space by identifying models that are not self-consistent through a comparison of the simulated temperature in the He 10830 $\r{A}$ line-forming region to the temperature assumed in creating the models. We demonstrate this procedure on models based on HD 209458 b. By investigating the Parker wind models with an assumed temperature between 4000 and 12000 K, and a mass-loss rate between $10^{8}$ and $10^{11}$ g s$^{-1}$, we are able to rule out a large portion of this parameter space. Furthermore, we fit the models to previous observational data and combine both constraints to find a preferred thermospheric temperature of $T=8200^{+1200}_{-1100}$ K and a mass-loss rate of $\dot{M}=10^{9.84^{+0.24}_{-0.27}}$ g s$^{-1}$ assuming a fixed atmospheric composition and no gas pressure confinement by the stellar wind. Using the same procedure, we constrain the temperatures and mass-loss rates of WASP-69 b, WASP-52 b, HAT-P-11 b, HAT-P-18 b and WASP-107 b.

Star-forming molecular clouds are characterised by the ubiquity of intertwined filaments. The filaments have been observed in both high- and low-mass star-forming regions, and are thought to split into collections of sonic fibres. The locations where filaments converge are termed hubs, and these are associated with the young stellar clusters. However, the observations of filamentary structures within hubs at distances require a high angular resolution that limits the number of such studies conducted so far. The integral shaped filament of the Orion A molecular cloud is noted for harbouring several hubs within which no filamentary structures have been observed so far. The goal of our study is to investigate the nature of the filamentary structures within one of these hubs, which is the chemically rich hub OMC-2 FIR 4, and to analyse their emission with high density and shock tracers. We observed the OMC-2 FIR 4 proto-cluster using Band 6 of the ALMA in Cycle 4 with an angular resolution of ~0.26"(100 au). We analysed the spatial distribution of dust, the shock tracer SiO, and dense gas tracers (i.e., CH$_{3}$OH, CS, and H$^{13}$CN). We also studied gas kinematics using SiO and CH3OH maps. Our observations for the first time reveal interwoven filamentary structures within OMC-2 FIR 4 that are probed by several tracers. Each filamentary structure is characterised by a distinct velocity as seen from the emission peak of CH$_{3}$OH lines. They also show transonic and supersonic motions. SiO is associated with filaments and also with multiple bow-shock features. In addition, for the first time, we reveal a highly collimated SiO jet (~1$^{\circ}$) with a projected length of ~5200 au from the embedded protostar VLA15. Our study shows that multi-scale observations of these regions are crucial for understanding the accretion processes and flow of material that shapes star formation.

Puffy disc is a numerical model, expected to capture the properties of the accretion flow in X-ray black hole binaries in the luminous, mildly sub-Eddington state. We fit the kerrbb and kynbb spectral models in XSPEC to synthetic spectra of puffy accretion discs, obtained in general relativistic radiative magnetohydrodynamic simulations, to see if they correctly recover the black hole spin and mass accretion rate assumed in the numerical simulation. We conclude that neither of the two models is capable of correctly interpreting the puffy disc parameters, which highlights a necessity to develop new, more accurate, spectral models for the luminous regime of accretion in X-ray black hole binaries. We propose that such spectral models should be based on the results of numerical simulations of accretion.

Classical Cepheids (DCEPs) are the first fundamental step in the calibration of the cosmological distance ladder. Furthermore, they represent powerful tracers in the context of Galactic studies. We have collected high-resolution spectroscopy with UVES@VLT for a sample of 65 DCEPs. The majority of them are the faintest DCEPs ever observed in the Milky Way. For each target, we derived accurate atmospheric parameters, radial velocities, and abundances for 24 different species. The resulting iron abundances range between +0.3 and $-$1.1 dex with the bulk of stars at [Fe/H]$\sim-0.5$ dex. Our sample includes the most metal-poor DCEPs observed so far with high-resolution spectroscopy. We complement our sample with literature data obtaining a complete sample of 637 DCEPs and use Gaia Early Data Release 3 (EDR3) photometry to determine the distance of the DCEPs from the Period-Wesenheit-Metallicity relation. Our more external data trace the Outer arm (at Galactocentric radius ($R_{GC})\sim$16--18 kpc) which appears significantly warped. We investigate the metallicity gradient of the Galactic disc using this large sample, finding a slope of $-0.060 \pm 0.002$ dex kpc$^{-1}$, in very good agreement with previous results based both on DCEPs and open clusters. We also report a possible break in the gradient at $R_{GC}$=9.25 kpc with slopes of $-0.063 \pm 0.007$ and $-0.079 \pm 0.003$ dex kpc$^{-1}$ for the inner and outer sample, respectively. The two slopes differ by more than 1 $\sigma$. A more homogeneous and extended DCEPs sample is needed to further test the plausibility of such a break.

The stable secondary-to-primary flux ratios of cosmic rays (CRs), represented by the boron-to-carbon ratio (B/C), are the main probes of the Galactic CR propagation. However, the fluorine-to-silicon ratio (F/Si) predicted by the CR diffusion coefficient inferred from B/C is significantly higher than the latest measurement of AMS-02. This anomaly is commonly attributed to the uncertainties of the F production cross sections. In this work, we give a careful test to this interpretation. We consider four different cross-section parametric models. Each model is constrained by the latest cross-section data. We perform combined fits to the B/C, F/Si, and cross-section data with the same propagation framework. Two of the cross-section models have good overall goodness of fit with $\chi^2/n_{d.o.f.}\sim1$. However, the goodness of fit of the cross-section part is poor with $\chi^2_{\rm{cs}}/n_{\rm{cs}}\gtrsim2$ for these models. The best-fitted F production cross sections are systematically larger than the measurements, while the fitted cross sections for B production are systematically lower than the measurements. This indicates that the F anomaly can hardly be interpreted by neither the random errors of the cross-section measurements nor the differences between the existing cross-section models. We then propose that the spatially dependent diffusion model could help to explain B/C and F/Si consistently. In this model, the average diffusion coefficient of the Ne-Si group is expected to be larger than that of the C-O group.

Pulsar wind nebula (PWN) Boomerang and the associated supernova remnant (SNR) G106.3+2.7 are among candidates for the ultra-high-energy (UHE) gamma-ray counterparts published by LHAASO. Although the centroid of the extended source, LHAASO J2226+6057, deviates from the pulsar's position by about $0.3^\circ$, the source partially covers the PWN. Therefore, we cannot totally exclude the possibility that a part of the UHE emission comes from the PWN. Previous studies mainly focus on whether the SNR is a PeVatron, while neglecting the energetic PWN. Here, we explore the possibility of the Boomerang Nebula being a PeVatron candidate by studying its X-ray radiation. By modelling the diffusion of relativistic electrons injected in the PWN, we fit the radial profiles of the X-ray surface brightness and the photon index. The solution with a magnetic field $B=140\mu$G can well reproduce the observed profiles and implies a severe suppression of IC scattering of electrons. Therefore, a proton component need be introduced to account for the UHE emission, in light of recent LHAASO's measurement on Crab Nebula, if future observations reveal part of the UHE emission originating from the PWN. In this sense, Boomerang Nebula would be a hadronic PeVatron.

The Andromeda Galaxy (M31) is an object of ongoing study with the Ultraviolet Imaging Telescope (UVIT) on AstroSat. UVIT FUV and NUV photometry is carried out here for a set of 239 clusters in the NE disk and bulge of M31 which overlap with the HST/PHAT survey. Padova stellar models were applied to derive ages, masses, metallicities and extinctions for 170 clusters. The ages show a narrow peak at $\sim$4 Myr and a broad peak around 100 Myr. log(Z/Z$_{\odot}$) values are mostly between $-0.3$ and $+0.3$. The 7 clusters in the bulge are low metallicity and high mass. Most clusters are in the spiral arms and have metallicities in the range noted above. The youngest clusters are mostly high metallicity and are concentrated along the brightest parts of the spiral arms. The UVIT FUV and NUV data are sensitive to young stars and detect a new metal rich peak in star formation in the disk at age $\sim4$ Myr.

Exoplanets smaller than Neptune are common around red dwarf stars (M dwarfs), with those that transit their host star constituting the bulk of known temperate worlds amenable for atmospheric characterization. We analyze the masses and radii of all known small transiting planets around M dwarfs, identifying three populations: rocky, water-rich, and gas-rich. Our results are inconsistent with the previously known bimodal radius distribution arising from atmospheric loss of a hydrogen/helium envelope. Instead, we propose that a density gap separates rocky from water-rich exoplanets. Formation models that include orbital migration can explain the observations: Rocky planets form within the snow line, whereas water-rich worlds form outside it and later migrate inward.

The soft X-ray pulsar RX J1856.5-3754 is the brightest member of a small class of thermally-emitting, radio-silent, isolated neutron stars. Its X-ray spectrum is almost indistinguishable from a blackbody with $kT^\infty\approx 60$ eV, but evidence of harder emission above $\sim 1$ keV has been recently found. We report on a spectral and timing analysis of RX J1856.5-3754 based on the large amount of data collected by XMM-Newton in 2002--2022, complemented by a dense monitoring campaign carried out by NICER in 2019. Through a phase-coherent timing analysis we obtained an improved value of the spin-down rate $\dot{\nu}=-6.042(4)\times10^{-16}$ Hz s$^{-1}$, reducing by more than one order magnitude the uncertainty of the previous measurement, and yielding a characteristic spin-down field of $1.47\times10^{13}$ G. We also detect two spectral components above $\sim1$ keV: a blackbody-like one with $kT^\infty=138\pm13$ eV and emitting radius $31_{-16}^{+8}$ m, and a power law with photon index $\Gamma=1.4_{-0.4}^{+0.5}$. The power-law 2--8\,keV flux, $(2.5_{-0.6}^{+0.7})\times10{-15}$ erg cm$^{-2}$ s$^{-1}$, corresponds to an efficiency of $10^{-3}$, in line with that seen in other pulsars. We also reveal a small difference between the $0.1$--$0.3$ keV and $0.3$--$1.2$ keV pulse profiles, as well as some evidence for a modulation above $1.2$ keV. These results show that, notwithstanding its simple spectrum, \eighteen still has a non-trivial thermal surface distribution and features non-thermal emission as seen in other pulsars with higher spin-down power.

We present the analysis of three more planets from the KMTNet 2021 microlensing season. KMT-2021-BLG-0119Lb is a $\sim 6\, M_{\rm Jup}$ planet orbiting an early M-dwarf or a K-dwarf, KMT-2021-BLG-0192Lb is a $\sim 2\, M_{\rm Nep}$ planet orbiting an M-dwarf, and KMT-2021-BLG-0192Lb is a $\sim 1.25\, M_{\rm Nep}$ planet orbiting a very--low-mass M dwarf or a brown dwarf. These by-eye planet detections provide an important comparison sample to the sample selected with the AnomalyFinder algorithm, and in particular, KMT-2021-BLG-2294, is a case of a planet detected by-eye but not by-algorithm. KMT-2021-BLG-2294Lb is part of a population of microlensing planets around very-low-mass host stars that spans the full range of planet masses, in contrast to the planet population at $\lesssim 0.1\, $ au, which shows a strong preference for small planets.

Most stars and thus most planetary systems do not form in isolation. The larger star-forming environment affects protoplanetary disks in multiple ways: gravitational interactions with other stars truncate disks and alter the architectures of exoplanet systems; external irradiation from nearby high-mass stars truncates disks and shortens their lifetimes; and remaining gas and dust in the environment affects dynamical evolution (if removed by feedback processes) and provides some shielding for disks from external irradiation. The dynamical evolution of the region regulates when and how long various feedback mechanisms impact protoplanetary disks. Density is a key parameter that regulates the intensity and duration of UV irradiation and the frequency of dynamical encounters. The evolution of larger star-forming complexes may also play an important role by mixing populations. Observations suggest that clusters are not a single-age population but multiple populations with small age differences which may be key to resolving several timescale issues (i.e., proplyd lifetimes, enrichment). In this review, we consider stellar clusters as the ecosystems in which most stars and therefore most planets form. We review recent observational and theoretical results and highlight upcoming contributions from facilities expected to begin observations in the next five years. Looking further ahead, we argue that the next frontier is large-scale surveys of low-mass stars in more distant high-mass star-forming regions. The future of ecosystem studies is bright as faint low-mass stars in more distant high-mass star-forming regions will be routinely observable in the era of Extremely Large Telescopes (ELTs).

Recent TESS-based studies have suggested that the dayside of KELT-1b, a strongly-irradiated brown dwarf, is significantly brighter in visible light than what would be expected based on Spitzer observations in infrared. We observe eight eclipses of KELT-1b with CHEOPS (CHaracterising ExOPlanet Satellite) to measure its dayside brightness temperature in the bluest passband observed so far, and model the CHEOPS photometry jointly with the existing optical and NIR photometry from TESS, LBT, CFHT, and Spitzer. Our modelling leads to a self-consistent dayside spectrum for KELT-1b covering the CHEOPS, TESS, H , Ks, and Spitzer IRAC 3.6 and 4.5 ${\mu}$m bands, where our TESS, H , Ks, and Spitzer band estimates largely agree with the previous studies, but we discover a strong discrepancy between the CHEOPS and TESS bands. The CHEOPS observations yield a higher photometric precision than the TESS observations, but do not show a significant eclipse signal, while a deep eclipse is detected in the TESS band. The derived TESS geometric albedo of $0.36^{+0.12}_{-0.13}$ is difficult to reconcile with a CHEOPS geometric albedo that is consistent with zero because the two passbands have considerable overlap. Variability in cloud cover caused by the transport of transient nightside clouds to the dayside could provide an explanation for reconciling the TESS and CHEOPS geometric albedos, but this hypothesis needs to be tested by future observations.

Almost all meteorite impacts occur at oblique incidence angles, but the effect of impact angle on crater size is not well understood, especially for large craters. To improve oblique impact crater scaling, we present a suite of simulations of complex crater formation on Earth and the Moon over a range of impact angles, velocities and impactor sizes. We show that crater diameter is larger than predicted by existing scaling relationships for oblique impacts and for impacts steeper than 45$^{\circ}$ shows little dependence on obliquity. Crater depth, volume and diameter depend on impact angle in different ways such that relatively shallower craters are formed by more oblique impacts. Our simulation results have implications for how crater populations are determined from impactor populations and vice-versa. Our results suggest that existing approaches to account for impact obliquity may underestimate the number of craters larger than a given size by as much as 40%.

Estimating the distances to asymptotic giant branch (AGB) stars using optical measurements of their parallaxes is not straightforward because of the large uncertainties introduced by their dusty envelopes, their large angular sizes, and their surface brightness variability. We compared the parallaxes from Gaia DR3 with parallaxes measured with maser observations with very long baseline interferometry (VLBI) to determine a statistical correction factor for the DR3 parallaxes using a sub-sample of 33 maser-emitting oxygen-rich nearby AGB stars. We found that the nominal errors on the Gaia DR3 parallaxes are underestimated by a factor of 5.44 for the brightest sources ($G<8$ mag). Fainter sources ($8\leq G<12$) require a lower parallax error inflation factor of 2.74. We obtain a Gaia DR3 parallax zero-point offset of -0.077 mas for bright AGB stars. The offset becomes more negative for fainter AGB stars. We calculated the distances of about 200 AGB stars in the DEATHSTAR project using a Bayesian statistical approach on the corrected DR3 parallaxes and a prior based on the previously determined Galactic distribution of AGB stars. We found that the derived Gaia distances are associated with significant, asymmetrical errors for more than 40% of the sources in our sample. We performed radiative transfer modelling of the stellar and dust emission to determine the luminosity of the sources in the VLBI sub-sample based on the distances derived from maser parallaxes, and derived a new bolometric period-luminosity relation for Galactic oxygen-rich Mira variables. A new distance catalogue based on these results is provided for the sources in the DEATHSTAR sample. We found that a RUWE (re-normalised unit weight error) below 1.4 does not guarantee reliable distance estimates and we advise against the use of only the RUWE to measure the quality of Gaia DR3 astrometric data for individual AGB stars.

Among supernovae (SNe) of different luminosities, many double-peaked light curves (LCs) have been observed, representing a broad morphological variety. In this work, we investigate which of these can be modelled by assuming a double-peaked distribution of their radioactive material, as originally proposed for SN2005bf. The inner zone corresponds to the regular explosive nucleosynthesis and extends outwards, according to the usual scenario of mixing. The outer 56Ni-rich shell may be related to the effect of jet-like outflows that have interacted with more distant portions of the star before the arrival of the SN shock. As the outer layer is covered by matter that is optically less thick, its energy emerges earlier and generates a first peak of radiation. To investigate this scenario in more detail, we have applied our hydrodynamic code that follows the shock propagation through the progenitor star and takes into account the effect of the gamma-ray photons produced by the decay of the radioactive isotopes. We present a simple parametric model for the 56Ni abundance profile and explore the consequences on the LC of individually varying the quantities that define this distribution, setting our focus onto the stripped-envelope progenitors. In this first study, we are interested in the applicability of this model to SNe that have not been classified as superluminous, thus, we have selected our parameter space accordingly. Then, within the same mathematical prescription for the 56Ni-profile, we revisited the modelling process for a series of objects: SN2005bf, PTF2011mnb, SN2019cad, and SN2008D. In some cases, a decrease in the gamma ray opacity is required to fit the late time observations. We also discuss the other cases in which this scenario might be likely to explain the LC morphology.

We study the topology of the network of ionized and neutral regions that characterized the intergalactic medium during the Epoch of Reionization. Our analysis uses the formalism of persistent homology, which offers a highly intuitive and comprehensive description of the ionization topology in terms of the births and deaths of topological features. Features are identified as $k$-dimensional holes in the ionization bubble network, whose abundance is given by the $k$th Betti number: $\beta_0$ for ionized bubbles, $\beta_1$ for tunnels, and $\beta_2$ for neutral islands. Using semi-numerical models of reionization, we investigate the dependence on the properties of sources and sinks of ionizing radiation. Of all topological features, we find that the tunnels dominate during reionization and that their number is easiest to observe and most sensitive to the astrophysical parameters of interest, such as the gas fraction and halo mass necessary for star formation. Seen as a phase transition, the importance of the tunnels can be explained by the entanglement of two percolating clusters and the fact that higher-dimensional features arise when lower-dimensional features link together. We also study the relation between the morphological components of the bubble network (bubbles, tunnels, islands) and those of the cosmic web (clusters, filaments, voids), describing a correspondence between the $k$-dimensional features of both. Finally, we apply the formalism to mock observations of the 21-cm signal. Assuming 1000 observation hours with HERA Phase II, we show that astrophysical models can be differentiated and confirm that persistent homology provides additional information beyond the power spectrum.

We measure the tidal alignment of the major axes of Luminous Red Galaxies (LRGs) from the Legacy Imaging Survey and use it to infer the artificial redshift-space distortion signature that will arise from an orientation-dependent, surface-brightness selection in the Dark Energy Spectroscopic Instrument (DESI) survey. Using photometric redshifts to down-weight the shape-density correlations due to weak lensing, we measure the intrinsic tidal alignment of LRGs. Separately, we estimate the net polarization of LRG orientations from DESI's fiber-magnitude target selection to be of order 10^-2 along the line of sight. Using these measurements and a linear tidal model, we forecast a 0.5% fractional decrease on the quadrupole of the 2-point correlation function for projected separations of 40-80 Mpc/h. We also use a halo catalog from the Abacus Summit cosmological simulation suite to reproduce this false quadrupole.

Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop new N-body codes, which allow for experimentation with various tidal models and numerical implementations. We present the open-source N-body code TIDYMESS, which stands for ``TIdal DYnamics of Multi-body ExtraSolar Systems''. This code implements a creep deformation law for the bodies, parametrized by their fluid Love numbers and fluid relaxation times. Due to tidal and centrifugal deformations, we approximate the general shape of a body to be an ellipsoid. We calculate the associated gravitational field to quadruple order, from which we derive the gravitational accelerations and torques. The equations of motion for the orbits, spins and deformations are integrated directly using a fourth-order integration method based on a symplectic composition. We implement a novel integration method for the deformations, which allows for a time step solely dependent on the orbits, and not on the spin periods or fluid relaxation times. This feature greatly speeds up the calculations, while also improving the consistency when comparing different tidal regimes. We demonstrate the capabilities and performance of TIDYMESS, particularly in the niche regime of parameter space where orbits are chaotic and tides become non-linear.

Recent sky surveys have discovered a large number of stellar substructures. It is highly likely that there are dark matter (DM) counterparts to these stellar substructures. We examine the implications of DM substructures for electron recoil (ER) direct detection (DD) rates in dual phase xenon experiments. We have utilized the results of the LAMOST survey and considered a few benchmark substructures in our analysis. Assuming that these substructures constitute $\sim 10\%$ of the local DM density, we study the discovery limits of DM-electron scattering cross sections considering one kg-year exposure and 1, 2, and 3 electron thresholds. With this exposure and threshold, it is possible to observe the effect of the considered DM substructure for the currently allowed parameter space. We also explore the sensitivity of these experiments in resolving the DM substructure fraction. For all the considered cases, we observe that DM having mass $\mathcal{O}(10)\,$MeV has a better prospect in resolving substructure fraction as compared to $\mathcal{O}(100)\,$MeV scale DM. We also find that within the currently allowed DM-electron scattering cross-section; these experiments can resolve the substructure fraction (provided it has a non-negligible contribution to the local DM density) with good accuracy for $\mathcal{O}(10)\,$MeV DM mass with one electron threshold.

The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction postulate then introduces fundamental uncertainties in the Cartesian phase space coordinates that propagate through classical three-body dynamics to macroscopic scales within the triple's lifetime. The spontaneous collapse of the wave function then drives the global chaotic behavior of the Universe through the superposition of triple systems (and probably multi-body systems). The paradox of infinitesimal granularity then arises from the superposition principle, which states that any multi-body system is composed of an ensemble of three-body problems.

Neutron stars are versatile in their application to studying various important aspects of fundamental physics, in particular strong-field gravity tests and the equation of state for super-dense nuclear matter at low temperatures. However, in many cases these two objectives are degenerate to each other. We discuss how pulsar timing and gravitational waves provide accurate measurements of neutron star systems and how to effectively break the degeneracy using tools like universal relations. We also present perspectives on future opportunities and challenges in the field of neutron star physics.

The irreducible upscattering of cold dark matter by cosmic rays opens up the intriguing possibility of detecting even light dark matter in conventional direct detection experiments or underground neutrino detectors. The mechanism also significantly enhances sensitivity to models with very large nuclear scattering rates, where the atmosphere and rock overburden efficiently stop standard non-relativistic dark matter particles before they could reach the detector. In this article, we demonstrate that cosmic-ray upscattering essentially closes the window for strongly interacting dark matter in the (sub-)GeV mass range. Arriving at this conclusion crucially requires a detailed treatment of both nuclear form factors and inelastic dark matter-nucleus scattering, as well as including the full momentum-transfer dependence of scattering amplitudes. We illustrate the latter point by considering three generic situations where such a momentum-dependence is particularly relevant, namely for interactions dominated by the exchange of light vector or scalar mediators, respectively, and for dark matter particles of finite size. As a final concrete example, we apply our analysis to a putative hexaquark state, which has been suggested as a viable baryonic dark matter candidate. Once again, we find that the updated constraints derived in this work close a significant part of otherwise unconstrained parameter space.

The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a function of mass ratio for eccentric non-spinning binary black holes. We produce $52$ NR simulations with mass ratios between $1:10$ and $1:1$ and initial eccentricities up to $0.7$. From these we extract quantities like gravitational wave energy and angular momentum fluxes and periastron advance, and assess their accuracy. To facilitate comparison, we develop tools to map between NR and SMR inspiral evolutions of eccentric binary black holes. We derive post-Newtonian accurate relations between different definitions of eccentricity. Based on these analyses, we introduce a new definition of eccentricity based on the (2,2)-mode of the gravitational radiation, which reduces to the Newtonian definition of eccentricity in the Newtonian limit. From the comparison between NR simulations and SMR results, we quantify the unknown next-to-leading order SMR contributions to the gravitational energy and angular momentum fluxes, and periastron advance. We show that in the comparable mass regime these contributions are subdominant and higher order SMR contributions are negligible.

Mentorship is critical to student academic success and persistence, especially for students from historically underrepresented (HU) groups. In a program designed to support the academic success of HU undergraduates in STEM who wish to pursue a PhD in those fields, students experience comprehensive support including financial aid, highly-engaged mentoring, dual faculty mentorship, professional development workshops, and summer research experiences. Scholars in this program, the Cal-Bridge program, consistently report that faculty mentorship is the most impactful feature. While mentorship was rated highly, preliminary evaluation indicated an early deficit in a sense of community among scholars. In response, faculty professional development and support for peer networking were implemented to expand and enhance the relationships that support scholar success. Here we present a promising multifaceted model of mentorship that can support the academic success of HU undergraduates.

We revisit the quark-mass density-dependent model -- a phenomenological equation of state for deconfined quark matter in the high-density low-temperature regime -- and show that thermodynamic inconsistencies that have plagued the model for decades, can be solved if the model is formulated in the canonical ensemble instead of the grand canonical one. Within the new formulation, the minimum of the energy per baryon occurs at zero pressure, and the Euler's relation is verified. Adopting a typical mass-formula, we first analyze in detail a simple model with one particle species. We show that a ``bag'' term that produces quark confinement naturally appears in the pressure (and not in the energy density) due to density dependence of the quark masses. Additionally, the chemical potential gains a new term as in other models with quark repulsive interactions. Then, we extend the formalism to the astrophysically realistic case of charge-neutral three-flavor quark matter in equilibrium under weak interactions, focusing on two different mass formulae: a flavor dependent and a flavor blind one. For these two models, we derive the equation of state and analyze its behavior for several parameter choices. We systematically analyze the parameter space and identify the regions corresponding to self-bound 2-flavor and 3-flavor quark matter, hybrid matter and causal behavior.

Ultralight bosons are a class of hypothetical particles that could potentially solve critical problems in fields ranging from cosmology to astrophysics and fundamental physics. If ultralight bosons exist, they form clouds around spinning black holes with sizes comparable to their particle Compton wavelength through superradiance, a well-understood classical wave amplification process that has been studied for decades. After these clouds form, they dissipate and emit continuous gravitational waves through the annihilation of ultralight bosons into gravitons. These gravitons could be detected with ground-based gravitational-wave detectors using continuous-wave searches. However, it is conceivable for other continuous-wave sources to mimic the emission from the clouds, which could lead to false detections. Here we investigate how one can use continuous waves from clouds formed around known merger remnants to alleviate this problem. In particular, we simulate a catalogue of merger remnants that form clouds around them and demonstrate with select "golden" merger remnants how one can perform a Bayesian cross-verification of the ultralight boson hypothesis that has the potential to rule out alternative explanations. Our proof-of-concept study suggest that, in the future, there is a possibility that a merger remnant exists close enough for us to perform the analysis and test the boson hypothesis if the bosons exist in the relevant mass range. Future research will focus on building more sophisticated continuous-wave tools to perform this analysis in practice.

We present exact solutions of test particle orbits spiralling inward from the innermost stable circular orbit (ISCO) of a Kerr black hole. Our results are valid for any allowed value of the angular momentum $a$-parameter of the Kerr metric. These solutions are of considerable physical interest. In particular, the radial 4-velocity of these orbits is both remarkably simple and, with the radial coordinate scaled by its ISCO value, universal in form, otherwise completely independent of the black hole spin.

Hawking's area theorem is one of the fundamental laws of black holes (BHs), which has been tested at a confidence level of $\sim 95\%$ with gravitational wave (GW) observations by analyzing the inspiral and ringdown portions of GW signals independently. In this work, we propose to carry out the test in a new way with the hierarchical triple merger (i.e., two successive BH mergers occurred sequentially within the observation window of GW detectors), for which the properties of the progenitor BHs and the remnant BH of the first coalescence can be reliably inferred from the inspiral portions of the two mergers. As revealed in our simulation, a test of the BH area law can be achieved at the significance level of $\gtrsim 3\sigma$ for the hierarchical triple merger events detected in LIGO/Virgo/KAGRA's O4/O5 runs. If the hierarchical triple mergers contribute a $\gtrsim 0.1\%$ fraction to the detected BBHs, a precision test of the BH area law with such systems is achievable in the near future. Our method also provides an additional criterion to establish the hierarchical triple merger origin of some candidate events.

Assuming a standard Maxwellian velocity distribution for the WIMPs in the halo of our Galaxy we use the null results of an exhaustive set of 6 direct detection experiments to calculate the maximal variation of the exclusion plot for each Wilson coefficient of the most general Galilean-invariant effective Hamiltonian for a WIMP of spin one half due to interferences. We consider 56 Wilson coefficients $c_i^{p,n}$ and $\alpha_i^{n,p}$ for WIMP-proton and WIMP-neutron contact interactions ${\cal O}_i^{p,n}$ and the corresponding long range interaction ${\cal O}_i^{p,n}/q^2$, parameterized by a massless propagator $1/q^2$. For each coupling we provide a different exclusion plot when the following set of operators is allowed to interfere: proton-neutron, i.e. $c_i^{p}$-$c_i^{n}$ or $\alpha_i^{p}$-$\alpha_i^{n}$; contact-contact or long range-long range, i.e. $c_i^{p,n}$-$c_j^{p,n}$ or $\alpha_i^{p,n}$-$\alpha_j^{p,n}$; contact-long range, i.e. $c_i^{p,n}$-$\alpha_j^{p,n}$. For each of the 56 Wilson coefficients $c_i^{p,n}$ and $\alpha_j^{p,n}$ and for the largest number of interfering operators the exclusion plot variation can reach 3 orders of magnitude and reduces to a factor as small as a few for the Wilson coefficients of the effective interactions where the WIMP couples to the nuclear spin, thanks to the combination of experiments using proton-odd and neutron-odd targets. Some of the conservative bounds require an extremely high level of cancellation, putting into question the reliability of the result. We analyze this issue in a systematic way, showing that it affects some of the couplings driven by the operators ${\cal O}_{1}$, ${\cal O}_{3}$, ${\cal O}_{11}$, ${\cal O}_{12}$ and ${\cal O}_{15}$, especially when interferences among contact and long range interactions are considered.

Using self-force methods, we consider the hyperbolic-type scattering of a pointlike particle carrying a scalar charge $Q$ off a Schwarzschild black hole. For given initial velocity and impact parameter, back-reaction from the scalar field modifies the scattering angle by an amount $\propto\! Q^2$, which we calculate numerically for a large sample of orbits (neglecting the gravitational self-force). Our results probe both strong-field and field-weak scenarios, and in the latter case we find a good agreement with post-Minkowskian expressions. The scalar-field self-force has a component tangent to the four-velocity that exchanges particle's mass with scalar-field energy, and we also compute this mass exchange as a function along the orbit. The expressions we derive for the scattering angle (in terms of certain integrals of the self-force along the orbit) can be used to obtain the gravitational self-force correction to the angle in the physical problem of a binary black hole with a large mass ratio. We discuss the remaining steps necessary to achieve this goal.

The tidal disruption of stars in the vicinity of massive black holes is discussed in the context of $\Lambda$-gravity. The latter provides an explanation to the Hubble tension as a possible consequence of two Hubble flows, the local and global ones. The bunch of notions which play role for the considered tidal effect are obtained, along with the rate of the disrupted stars. The role of pulsars is emphasized due to their ability to penetrate up to the horizon of the massive black hole as for them the tidal radius can reach the horizon. Tidal disruption mechanism also can lead to segregation of stars by their mean density vs the distance from the black hole, the denser stars surviving at shorter distances. The interplay of the central gravity field and the repulsive $\Lambda$-term increasing with radius and its certain observational consequences are discussed.

In this work we investigate the characteristics of the group velocity of obliquely propagating Alfv\'en waves in a dusty plasma typical of a stellar wind. The dispersion relation is derived with the aid of the kinetic theory for a magnetized dusty plasma consisting of electrons and ions, with distribution of momenta described by a Maxwellian function. The dust particles are considered to be immobile and have all the same size; they are electrically charged by absorption of plasma particles via inelastic collisions and by photoionization. We numerically solve the dispersion relation and calculate the components of group velocity (along and transverse to the magnetic field) for the normal modes, namely the compressional and shear Alfv\'en waves (CAW and SAW). The results show that the direction of the group velocity of CAWs is greatly modified with the wave-vector direction. On the other hand, SAWs will present group velocity propagating practically along the magnetic field. The changes in dust parameters, such as number density and equilibrium electrical charge, may significantly change the waves' characteristics. It is seen that for sufficiently high dust to ion number density ratio, the SAWs may present perpendicular group velocity propagating in opposite direction to the perpendicular phase velocity, in a small interval of wavenumber values; we also notice that this interval may change, or even vanish, when the flux of radiation incident on the dust is altered, changing the equilibrium electrical charge of the grains.

Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform ongoing efforts to construct an accurate description of the general-relativistic binary dynamics via an effective-one-body description or other approaches. Existing self-force methods are to a large extent specialised to bound, inspiral orbits. Here we derive the first-order conservative self-force correction to the scattering angle, show its agreement with recent post-Minkowsian results, and develop a technique for (numerical) self-force calculations that can efficiently tackle scatter orbits. In the method, the metric perturbation is reconstructed from a Hertz potential that satisfies (mode-by-mode) a certain inhomogeneous version of the Teukolsky equation. The crucial ingredient in this formulation are certain jump conditions that the (multipole modes of the) Hertz potential must satisfy along the worldline of the small body's orbit. We present a closed-form expression for these jumps, for an arbitrary geodesic orbit in Schwarzschild spacetime. To begin developing the numerical infrastructure, a scalar-field evolution code on a Schwarzschild background (in 1+1D) is developed. Following this, results for the conservative scalar self-force corrections to the scatter angle are calculated. We continue by constructing a Teukolsky evolution code on a Schwarzschild background. This produces numerically unstable solutions due to unphysical homogeneous solutions of the Teukolsky equation at the horizon and null infinity being seeded by numerical error. This can be resolved by a change of variables to a Regge-Wheeler-like field. We then present a full numerical implementation of this method for circular and scatter orbits in Schwarzschild.

Although the cosmic microwave background (CMB) is largely understood to be homogeneous and isotropic, the hemispherical asymmetry anomaly seems to breakdown the isotropy, since the difference between the power spectrum in the two hemispheres of the CMB is of the order of $10^{-2}$ at large angular scales. We argue that the existence of an anisotropic power spectrum can simply be explained by considering the existence of two distinct power spectra in the two hemispheres of the CMB. We achieve this by proposing a double vacuum structure for (single field) inflationary quantum fluctuations based on discrete spacetime transformations ($\mathcal{P}\mathcal{T}$) in a gravitational context, first in de Sitter and finally in quasi de Sitter. As a result we obtain inflationary quantum fluctuations that are produced in pairs with which we are able to reproduce the amplitude of the observed dipolar asymmetry at different scales of $ 10^{-4} {\rm Mpc^{-1}}\lesssim k\lesssim 1 {\rm Mpc^{-1}}$ fixing the pivot scale $k=0.05 \, {\rm Mpc}^{-1}$ for $N=55$ e-foldings of inflation. We also predict that a similar hemispherical asymmetry should arise for the primordial gravitational waves (PGWs) as well and we compute the power asymmetry of PGW spectra at various wave numbers. In our framework we do not introduce any new parameters.