Papers for Wednesday, Apr 21 2021

The wavelet analysis technique is a powerful tool and is widely used in broad disciplines of engineering, technology, and sciences. In this work, we present a novel scheme of constructing continuous wavelet functions, in which the wavelet functions are obtained by taking the first derivative of smoothing functions with respect to the scale parameter. Due to this wavelet constructing scheme, the inverse transforms are only one-dimensional integrations with respect to the scale parameter, and hence the continuous wavelet transforms constructed in this way are more ready to use than the usual scheme. We then apply the Gaussian-derived wavelet constructed by our scheme to computations of the density power spectrum for dark matter, the velocity power spectrum and the kinetic energy spectrum for baryonic fluid. These computations exhibit the convenience and strength of the continuous wavelet transforms. The transforms are very easy to perform, and we believe that the simplicity of our wavelet scheme will make continuous wavelet transforms very useful in practice.

Generating renewable energy on Mars is technologically challenging. Firstly, because compared to Earth, key energy resources such as solar and wind are weak as a result of very low atmospheric pressure and low solar irradiation. Secondly, because of the harsh environmental conditions, the required high degree of automation and the exceptional effort and costs to transport material to the planet. Like on Earth, it is crucial to combine complementary resources for an effective renewable energy solution. In this work, we present the result of a design synthesis exercise, a 10 kW microgrid solution, based on a pumping kite power system and photovoltaic solar modules to power the construction as well as the subsequent use of a Mars habitat. To buffer unavoidable energy fluctuations and balance seasonal and diurnal resource variations, the two energy systems are combined with a compressed gas storage system and lithium-sulfur batteries. The airborne wind energy solution was selected because of its low weight-to-wing-surface-area ratio, compact packing volume and high capacity factor which enables it to endure strong dust storms in an airborne parking mode. The surface area of the membrane wing is 50 m2 and the mass of the entire system, including the kite control unit and ground station, is 290 kg. The performance of the microgrid is assessed by computational simulation using available resource data for a chosen deployment location on Mars. The projected costs of the system are 8.95 million Euro, excluding transportation to Mars.

Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($\chi_\mathrm{eff}$), the precession spin parameter ($\chi_\mathrm{p}$), and the chirp mass ($m_\mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${\bar \chi}_\mathrm{typ}\equiv \overline{(\chi_\mathrm{eff}^2+\chi_\mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_\mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${\bar \chi}_\mathrm{typ}\sim 0.6$ at higher $m_\mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${\bar \chi}_\mathrm{typ}$ as a function of $m_\mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${\bar \chi}_\mathrm{typ}$ at low $m_\mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $\sim 2 \sigma$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $\sim 15$--$30\,M_\odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $\sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${\bar \chi}_\mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $\gtrsim O(100)$.

The Simons Observatory (SO) is a Cosmic Microwave Background (CMB) experiment to observe the microwave sky in six frequency bands from 30GHz to 290GHz. The Observatory -- at $\sim$5200m altitude -- comprises three Small Aperture Telescopes (SATs) and one Large Aperture Telescope (LAT) at the Atacama Desert, Chile. This research note describes the design and current status of the LAT along with its future timeline.

Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way. A local wake is predicted to trail the orbit of the LMC, while a large-scale over-density is predicted to exist across a large area of the northern Galactic hemisphere. Here we present the detection of both the local wake and Northern over-density (hereafter the "collective response") in an all-sky star map of the Galaxy based on 1301 stars at 60<R_gal<100 kpc. The location of the wake is in good agreement with an N-body simulation that includes the dynamical effect of the LMC on the Milky Way halo. The density contrast of the wake and collective response are both stronger in the data than in the simulation. The detection of a strong local wake is independent evidence that the Magellanic Clouds are on their first orbit around the Milky Way. The wake traces the path of the LMC, which will provide insight into the orbit of the LMC, which in turn is a sensitive probe of the mass of the LMC and the Milky Way. These data demonstrate that the outer halo is not in dynamical equilibrium, as is often assumed. The morphology and strength of the wake could be used to test the nature of dark matter and gravity.

In many astrophysical applications, the cost of solving a chemical network represented by a system of ordinary differential equations (ODEs) grows significantly with the size of the network, and can often represent a significant computational bottleneck, particularly in coupled chemo-dynamical models. Although standard numerical techniques and complex solutions tailored to thermochemistry can somewhat reduce the cost, more recently, machine learning algorithms have begun to attack this challenge via data-driven dimensional reduction techniques. In this work, we present a new class of methods that take advantage of machine learning techniques to reduce complex data sets (autoencoders), the optimization of multi-parameter systems (standard backpropagation), and the robustness of well-established ODE solvers to to explicitly incorporate time-dependence. This new method allows us to find a compressed and simplified version of a large chemical network in a semi-automated fashion that can be solved with a standard ODE solver, while also enabling interpretability of the compressed, latent network. As a proof of concept, we tested the method on an astrophysically-relevant chemical network with 29 species and 224 reactions, obtaining a reduced but representative network with only 5 species and 12 reactions, and a x65 speed-up.

We present the discovery of an extreme flaring event from Proxima Cen by ASKAP, ALMA, HST, TESS, and the du Pont Telescope that occurred on 2019 May 1. In the millimeter and FUV, this flare is the brightest ever detected, brightening by a factor of >1000 and >14000 as seen by ALMA and HST, respectively. The millimeter and FUV continuum emission trace each other closely during the flare, suggesting that millimeter emission could serve as a proxy for FUV emission from stellar flares and become a powerful new tool to constrain the high-energy radiation environment of exoplanets. Surprisingly, optical emission associated with the event peaks at a much lower level with a time delay. The initial burst has an extremely short duration, lasting for <10 sec. Taken together with the growing sample of millimeter M dwarf flares, this event suggests that millimeter emission is actually common during stellar flares and often originates from short burst-like events.

Gravitationally lensed curved arcs provide a wealth of information about the underlying lensing distortions. Extracting precise lensing information from extended sources is a key component in many studies aiming to answer fundamental questions about the Universe. To maintain accuracy with increased precision, it is of vital importance to characterize and understand the impact of degeneracies inherent in lensing observables. In this work, we present a formalism to describe the gravitational lensing distortion effects resulting in curved extended arcs based on the eigenvectors and eigenvalues of the local lensing Jacobian and their directional differentials. We identify a non-local and non-linear extended deflector basis that inherits these local properties. Our parameterization is tightly linked to observable features in extended sources and allows one to accurately extract the lensing information of extended images without imposing an explicit global deflector model. We quantify what degeneracies can be broken based on specific assumptions on the local lensing nature and assumed intrinsic source shape. Our formalism is applicable from the weak linear regime, the semi-linear regime all the way up to the highly non-linear regime of highly magnified arcs of multiple images. The methodology and implementation presented in this work provides a framework to assessing systematics, to guide inference efforts in the right choices in complexity based on the data at hand, and to quantify the lensing information extracted in a model-independent way.

Stellar streams produced from dwarf galaxies provide direct evidence of the hierarchical formation of the Milky Way. Here, we present the first comprehensive study of the "LMS-1" stellar stream, that we detect by searching for wide streams in the Gaia EDR3 dataset using the STREAMFINDER algorithm. This stream was recently discovered by Yuan et al. (2020). We detect LMS-1 as a $60\deg$ long stream to the north of the Galactic bulge, at a distance of $\sim 20$ kpc from the Sun, together with additional components that suggest that the overall stream is completely wrapped around the inner Galaxy. Using spectroscopic measurements from LAMOST, SDSS and APOGEE, we infer that the stream is very metal poor (${\rm \langle [Fe/H]\rangle =-2.1}$) with a significant metallicity dispersion ($\sigma_{\rm [Fe/H]}=0.4$), and it possesses a large radial velocity dispersion (${\rm \sigma_v=20 \pm 4\,km\,s^{-1}}$). These estimates together imply that LMS-1 is a dwarf galaxy stream. The orbit of LMS-1 is close to polar, with an inclination of $75\deg$ to the Galactic plane. Both the orbit and metallicity of LMS-1 are remarkably similar to the globular clusters NGC 5053, NGC 5024 and the stellar stream "Indus". These findings make LMS-1 an important contributor to the stellar population of the inner Milky Way halo.

Some of barred galaxies, including the Milky Way, host a boxy/peanut/X-shaped bulge (BPX-shaped bulge). Previous studiessuggested that the BPX-shaped bulge can either be developed by bar buckling or by vertical inner Lindblad resonance (vILR)heating without buckling. In this paper, we study the observable consequence of an BPX-shaped bulge built up quickly after barformation via vILR heating without buckling, using anN-body/hydrodynamics simulation of an isolated Milky Way-like galaxy.We found that the BPX-shaped bulge is dominated by stars born prior to bar formation. This is because the bar suppresses starformation, except for the nuclear stellar disc (NSD) region and its tips. The stars formed near the bar ends have higher Jacobienergy, and when these stars lose their angular momentum, their radial action increases to conserve Jacobi energy. This preventsthem from reaching the vILR to be heated to the BPX region. By contrast, the NSD forms after the bar formation. From thissimulation and general considerations, we expect that the age distributions of the NSD and BPX-shaped bulge formed withoutbar buckling do not overlap each other. Then, the transition age between these components betrays the formation time of the bar, and is testable in future observations of the Milky Way and extra-galactic barred galaxies

The kinematics of MgII absorbers is the key to understanding the origin of cool, metal-enriched gas clouds in the circumgalactic medium of massive quiescent galaxies. Exploiting the fact that the line-of-sight velocity distribution of those clouds is the only unknown for predicting the redshift-space distortion (RSD) of MgII absorbers from their 3D real-space distribution around galaxies, we develop a novel method to infer the cool cloud kinematics from the redshift-space galaxy-cloud cross-correlation function $\xi^{s}$. We measure $\xi^{s}$ for ~$10^4$ MgII absorbers around a sample of luminous red galaxies at $0.4{<}z{<}0.8$. By comparing with the RSD of galaxies, we discover that $\xi^{s}$ does not exhibit a strong Fingers-of-God effect, but is heavily truncated at relative velocity ~$300\,km/s$. We reconstruct both the redshift and real-space cloud number density distributions for absorbers that reside within the LRG halos, $\xi^{s}_{1h}$ and $\xi_{1h}$, respectively. Thus, for given model of cloud kinematics, we can predict $\xi^{s}_{1h}$ from the reconstructed $\xi_{1h}$, and compare to the observed $\xi^{s}_{1h}$ in a self-consistently manner. We consider four types of cloud kinematics, including an isothermal model with a single velocity dispersion, a satellite infall model in which cool clouds reside in the subhalos, a cloud accretion model in which cool clouds follow the cosmic dark matter accretion, and a tired wind model in which cool clouds originate from the galactic wind-driven bubbles. We find that all the four models provide statistically good fits to the RSD data, but only the tired wind model can reproduce the observed truncation by propagating an ancient wind bubble at ~$250\,km/s$ on scales ~$400\,kpc/h$. Our method provides an exciting path to decoding the dynamical origin of metal absorbers from the small-scale RSD modelling with upcoming spectroscopic surveys.

Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is particularly interesting if we consider that gravitational instabilities could be a fundamental first step for the formation of planetesimals and planets. In this Letter, we model the rotation curve obtained by CO data of Elias 2-27 with a theoretical rotation curve including both the disk self-gravity and the star contribution to the gravitational potential. We compare this model with a purely Keplerian one and with a simple power-law function. We find that (especially for the $^{13}$CO isotopologue) the rotation curve is better described by considering not only the star, but also the disk self-gravity. We are thus able to obtain for the first time a dynamical estimate of the disk mass of $0.08\pm0.04\,M_{\odot}$ and the star mass of $0.46\pm0.03\,M_{\odot}$ (in the more general case), the latter being comparable with previous estimates. From these values, we derive that the disk is 17$\%$ of the star mass, meaning that it could be prone to gravitational instabilities. This result would strongly support the hypothesis that the two spiral arms are generated by gravitational instabilities.

THESEUS, one of the two space mission concepts being studied by ESA as candidates for next M5 mission within its Comsic Vision programme, aims at fully exploiting Gamma-Ray Bursts (GRB) to solve key questions about the early Universe, as well as becoming a cornerstone of multi-messenger and time-domain astrophysics. By investigating the first billion years of the Universe through high-redshift GRBs, THESEUS will shed light on the main open issues in modern cosmology, such as the population of primordial low mass and luminosity galaxies, sources and evolution of cosmic re-ionization, SFR and metallicity evolution up to the "cosmic dawn" and across Pop-III stars. At the same time, the mission will provide a substantial advancement of multi-messenger and time-domain astrophysics by enabling the identification, accurate localisation and study of electromagnetic counterparts to sources of gravitational waves and neutrinos, which will be routinely detected in the late '20s and early '30s by the second and third generation Gravitational Wave (GW) interferometers and future neutrino detectors, as well as of all kinds of GRBs and most classes of other X/gamma-ray transient sources. In all these cases, THESEUS will provide great synergies with future large observing facilities in the multi-messenger domain. A Guest Observer programme, comprising Target of Opportunity (ToO) observations, will expand the science return of the mission, to include, e.g., solar system minor bodies, exoplanets, and AGN.

At peak, long-duration gamma-ray bursts are the most luminous sources of electromagnetic radiation known. Since their progenitors are massive stars, they provide a tracer of star formation and star-forming galaxies over the whole of cosmic history. Their bright power-law afterglows provide ideal backlights for absorption studies of the interstellar and intergalactic medium back to the reionization era. The proposed THESEUS mission is designed to detect large samples of GRBs at $z>6$ in the 2030s, at a time when supporting observations with major next generation facilities will be possible, thus enabling a range of transformative science. THESEUS will allow us to explore the faint end of the luminosity function of galaxies and the star formation rate density to high redshifts; constrain the progress of re-ionisation beyond $z\gtrsim6$; study in detail early chemical enrichment from stellar explosions, including signatures of Population III stars; and potentially characterize the dark energy equation of state at the highest redshifts.

THESEUS is a medium size space mission of the European Space Agency, currently under evaluation for a possible launch in 2032. Its main objectives are to investigate the early Universe through the observation of gamma-ray bursts and to study the gravitational waves electromagnetic counterparts and neutrino events. On the other hand, its instruments, which include a wide field of view X-ray (0.3-5 keV) telescope based on lobster-eye focusing optics and a gamma-ray spectrometer with imaging capabilities in the 2-150 keV range, are also ideal for carrying out unprecedented studies in time domain astrophysics. In addition, the presence onboard of a 70 cm near infrared telescope will allow simultaneous multi-wavelegth studies. Here we present the THESEUS capabilities for studying the time variability of different classes of sources in parallel to, and without affecting, the gamma-ray bursts hunt.

Multi-messenger astrophysics is becoming a major avenue to explore the Universe, with the potential to span a vast range of redshifts. The growing synergies between different probes is opening new frontiers, which promise profound insights into several aspects of fundamental physics and cosmology. In this context, THESEUS will play a central role during the 2030s in detecting and localizing the electromagnetic counterparts of gravitational wave and neutrino sources that the unprecedented sensitivity of next generation detectors will discover at much higher rates than the present. Here, we review the most important target signals from multi-messenger sources that THESEUS will be able to detect and characterize, discussing detection rate expectations and scientific impact.

The proposed THESEUS mission will vastly expand the capabilities to monitor the high-energy sky, and will exploit large samples of gamma-ray bursts to probe the early Universe back to the first generation of stars, and to advance multi-messenger astrophysics by detecting and localizing the counterparts of gravitational waves and cosmic neutrino sources. The combination and coordination of these activities with multi-wavelength, multi-messenger facilities expected to be operating in the thirties will open new avenues of exploration in many areas of astrophysics, cosmology and fundamental physics, thus adding considerable strength to the overall scientific impact of THESEUS and these facilities. We discuss here a number of these powerful synergies.

Aims: The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments within nearby star-forming galaxies from the PHANGS sample. Methods: Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of the resolved star formation main sequence (rSFMS), resolved Kennicutt-Schmidt relation (rKS), and resolved molecular gas main sequence (rMGMS), and we derived their slope and scatter at spatial resolutions from 100 pc to 1 kpc (under various assumptions). Results: All three relations were recovered at the highest spatial resolution (100 pc). Furthermore, significant variations in these scaling relations were observed across different galactic environments. The exclusion of non-detections has a systematic impact on the inferred slope as a function of the spatial scale. Finally, the scatter of the $\Sigma_\mathrm{mol. gas + stellar}$ versus $\Sigma_\mathrm{SFR}$ correlation is smaller than that of the rSFMS, but higher than that found for the rKS. Conclusions: The rMGMS has the tightest relation at a spatial scale of 100 pc (scatter of 0.34 dex), followed by the rKS (0.41 dex) and then the rSFMS (0.51 dex). This is consistent with expectations from the timescales involved in the evolutionary cycle of molecular clouds. Surprisingly, the rKS shows the least variation across galaxies and environments, suggesting a tight link between molecular gas and subsequent star formation. The scatter of the three relations decreases at lower spatial resolutions, with the rKS being the tightest (0.27 dex) at a spatial scale of 1 kpc. Variation in the slope of the rSFMS among galaxies is partially due to different detection fractions of $\Sigma_\mathrm{SFR}$ with respect to $\Sigma_\mathrm{stellar}$.

Accretion states, which are universally observed in stellar-mass black holes in X-ray binaries, are also anticipated in active galactic nuclei (AGN). This is the case at low luminosities, when the jet-corona coupling dominates the energy output in both populations. Previous attempts to extend this framework to a wider AGN population have been extremely challenging due to heavy hydrogen absorption of the accretion disc continuum and starlight contamination from the host galaxies. We present the luminosity-excitation diagram (LED), based on the [OIV]$_{25.9\mu m}$ and [NeII]$_{12.8\mu m}$ mid-infrared nebular line fluxes. This tool enables to probe the accretion disc contribution to the ionising continuum. When applied to a sample of 167 nearby AGN, the LED recovers the characteristic q-shaped morphology outlined by individual X-ray binaries during a typical accretion episode, allowing us to tentatively identify the main accretion states. The soft state would include broad-line Seyferts and about half of the Seyfert 2 population, showing highly excited gas and radio-quiet cores consistent with disc-dominated nuclei. The hard state mostly includes low-luminosity AGN ($\leq 10^{-3}\, \rm{L_{Edd}}$) characterised by low-excitation radio-loud nuclei and a negligible disc contribution. The remaining half of Seyfert 2 nuclei and the bright LINERs show low excitation at high accretion luminosities, and could be identified with the bright-hard and intermediate states. Their hosts show ongoing star formation in the central kiloparsecs. We discuss the above scenario, its potential links with the galaxy evolution picture and the possible presence of accretion state transitions in AGN, as suggested by the growing population of changing-look quasars.

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

The classification of the optical spectra of active galactic nuclei (AGN) into different types is well founded on AGN physics, but it involves some degree of human oversight and cannot be reliably scaled to large data sets. Machine learning (ML) tackles such a classification problem in a fast and reproducible way, but is often perceived as a black box. However, ML interpretability and explainability are active research areas in computer science, increasingly providing us with tools to alleviate this issue. We applied ML interpretability tools to a classifier trained to predict AGN type from spectra, to demonstrate the use of such tools in this context. We trained a support-vector machine on 3346 high-quality, low redshift AGN spectra from SDSS DR15 with an existing reliable classification as type 1, type 2, or intermediate type. On a selection of test-set spectra, we computed the gradient of the predicted class probability and we built saliency maps. We also visualized the high-dimensional space of AGN spectra using t-distributed stochastic neighbor embedding (t-SNE), showing where the spectra for which we computed a saliency map are located. Regions that affect the predicted AGN type often coincide with physically relevant features, such as spectral lines. t-SNE visualization shows good separability of type 1 and type 2 spectra, while intermediate-type spectra either lie in-between as expected or appear mixed with type 2 spectra. Saliency maps show why a given AGN type was predicted by our classifier, resulting in a physical interpretation in terms of regions of the spectrum that affected its decision, making it no longer a black box. These regions coincide with those used by human experts such as relevant spectral lines, and are even used in a similar way, with the classifier e.g. effectively measuring the width of a line by weighing its center and its tails oppositely.

The Central Molecular Zone (CMZ; the central ~500 pc of the Milky Way) hosts molecular clouds in an extreme environment of strong shear, high gas pressure and density, and complex chemistry. G0.253+0.016, also known as `the Brick', is the densest, most compact and quiescent of these clouds. High-resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed its complex, hierarchical structure. In this paper we compare the properties of recent hydrodynamical simulations of the Brick to those of the ALMA observations. To facilitate the comparison, we post-process the simulation and create synthetic ALMA maps of molecular line emission from eight molecules. We correlate the line emission maps to each other and to the mass column density, and find that HNCO is the best mass tracer of the eight emission lines. Additionally, we characterise the spatial structure of the observed and simulated cloud using the density probability distribution function (PDF), spatial power spectrum, fractal dimension, and moments of inertia. While we find good agreement between the observed and simulated data in terms of power spectra and fractal dimensions, there are key differences in terms of the density PDFs and moments of inertia, which we attribute to the omission of magnetic fields in the simulations. Models that include the external gravitational potential generated by the stars in the CMZ better reproduce the observed structure, highlighting that cloud structure in the CMZ results from the complex interplay between internal physics (turbulence, self-gravity, magnetic fields) and the impact of the extreme environment.

An interstellar communication network benefits from relay nodes placed in the gravitational lenses of stars. The signal gains are of order $10^{9}$ with optimal alignment, allowing for GBits connections at kW power levels with meter-sized probes over parsec distances. If such a network exists, there might be a node in our solar system: where is it? With some assumptions on the network topology, candidate sky positions can be calculated. Apparent positions are influenced by the parallax motion from the Earth's orbit around the Sun, and the (slow) drifts caused by proper motions of nearby stars. With Gaia astrometry, instantaneous positions can be determined with arcsec accuracy. These potential node locations can be observed in targeted

[Abridged] The study of disc kinematics has recently opened up as a promising method to detect unseen planets. However, a systematic, statistically meaningful analysis of such an approach remains missing. The aim of this work is to devise an automated, statistically robust technique to identify kinematical perturbations induced by the presence of planets in a gas disc, and to accurately infer their location. For this purpose, we produce hydro simulations of planet-disc interactions with different planet masses, 0.3, 1.0 and 3.0 $M_{Jup}$, at a radius of $R=100$ au in the disc, and perform radiative transfer calculations of CO to simulate observables for 13 planet azimuths. Using the DISCMINER package, we fit the synthetic data cubes with a Keplerian model of the channel-by-channel emission to study line profile differences, including deviations from Keplerian rotation. The detection technique, based on line centroid differences, captures localised planet-driven perturbations, and can distinguish them from axisymmetric velocity perturbations. The method can detect all three simulated planets, at all azimuths, with an average accuracy of $\pm3^\circ$ in azimuth and $\pm8$ au in radius. Owing to disc structure and line-of-sight projection effects, planets at azimuths close to $\pm45^\circ$ yield the highest velocity fluctuations, whereas those at limiting cases, $0^\circ$ and $\pm90^\circ$, drive the lowest. The observed peak velocities range within 40$-$70, 70$-$170 and 130$-$450 m s$^{-1}$ for 0.3, 1.0 and 3.0 $M_{Jup}$ planets. Our analysis indicates that the variance of peak velocities is boosted near planets due to organised gas motions prompted by their localised gravitational well. We propose an approach that exploits this velocity coherence to provide, for the first time, statistically significant detections of localised planet-driven perturbations in the gas disc kinematics.

IC 1459 is an early-type galaxy (ETG) with a rapidly counter-rotating stellar core, and is the central galaxy in a gas-rich group of spirals. In this work, we investigate the abundant ionized gas in IC 1459 and present new stellar orbital models to connect its complex array of observed properties and build a more complete picture of its evolution. Using the Multi-Unit Spectroscopic Explorer (MUSE), the optical integral field unit (IFU) on the Very Large Telescope (VLT), we examine the gas and stellar properties of IC 1459 to decipher the origin and powering mechanism of the galaxy's ionized gas. We detect ionized gas in a non-disk-like structure rotating in the opposite sense to the central stars. Using emission-line flux ratios and velocity dispersion from full-spectral fitting, we find two kinematically distinct regions of shocked emission-line gas in IC 1459, which we distinguished using narrow ($\sigma$ $\leq$ 155 km s$^{-1}$) and broad ($\sigma$ $>$ 155 km s$^{-1}$) profiles. Our results imply that the emission-line gas in IC 1459 has a different origin than that of its counter-rotating stellar component. We propose that the ionized gas is from late-stage accretion of gas from the group environment, which occurred long after the formation of the central stellar component. We find that shock heating and AGN activity are both ionizing mechanisms in IC 1459 but that the dominant excitation mechanism is by post-asymptotic giant branch stars from its old stellar population.

Superluminous supernovae (SLSNe) are luminous transients that can be detected to high redshifts with upcoming optical time-domain surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). An interesting open question is whether the properties of SLSNe evolve through cosmic time. To address this question, in this paper we model for the first time the multi-color light curves of all 21 Type I SLSNe from the Dark Energy Survey (DES) with a magnetar spin-down engine, implemented in the Modular Open Source Fitter for Transients (MOSFiT). With redshifts up to $z\approx 2$, this sample includes some of the highest-redshift SLSNe. We find that the DES SLSNe span a similar range of ejecta and magnetar engine parameters to previous samples of mostly lower-redshift SLSNe (spin period $P\approx 0.79-13.64$ ms, magnetic field $B\approx (0.03-7.12)\times10^{14}$ G, ejecta mass $M_{\rm ej}\approx 1.51-30.50$ M$_{\odot}$, and ejecta velocity $v_{\rm ej}\approx (0.56-1.45)\times 10^4$ km s$^{-1}$). The DES SLSN sample by itself exhibits the previously found negative correlation between $M_{\rm ej}$ and $P$, with a pronounced absence of SLSNe with low ejecta mass and rapid spin. Combining our results for the DES SLSNe with 60 previous SLSNe modeled in the same way, we find no evidence for redshift evolution in any of the key physical parameters.

The study of the quasi-periodic oscillations (QPOs) of X-ray flux observed in the stellar-mass black hole (BH) binaries or quasars can provide a powerful tool for testing the phenomena occurring in strong gravity regime. We thus fit the data of QPOs observed in the well known quasars as well as microquasars in the framework of the model of geodesic oscillations of Keplerian disks modified for the epicyclic oscillations of spinning test particles orbiting Kerr BHs. We show that the modified geodesic models of QPOs can explain the observational fixed data from the quasars and microquasars but not for all sources. We perform a successful fitting of the high frequency QPOs models of epicyclic resonance and its variants, relativistic precession and its variants, tidal disruption, as well as warped disc models, and discuss the corresponding constraints of parameters of the model, which are the spin $S$ of the test particle, mass $M$ and spin $a$ of the BH.

Here we report on X-ray observations of ten 17-60 keV sources discovered by the INTEGRAL satellite. The primary new information is sub-arcsecond positions obtained by the Chandra X-ray Observatory. In six cases (IGR J17040-4305, IGR J18017-3542, IGR J18112-2641, IGR J18434-0508, IGR J19504+3318, and IGR J20084+3221), a unique Chandra counterpart is identified with a high degree of certainty, and for five of these sources (all but J19504), Gaia distances or proper motions indicate that they are Galactic sources. For four of these, the most likely classifications are that the sources are magnetic Cataclysmic Variables (CVs). J20084 could be either a magnetic CV or a High Mass X-ray Binary. We classify the sixth source (J19504) as a likely Active Galactic Nucleus (AGN). In addition, we find likely Chandra counterparts to IGR J18010-3045 and IGR J19577+3339, and the latter is a bright radio source and probable AGN. The other two sources, IGR J12529-6351 and IGR J18013-3222 do not have likely Chandra counterparts, indicating that they are transient, highly variable, or highly absorbed.

Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from $\sim 5\times10^5$ up to $\sim 10^{9} M_\odot$, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [$\alpha$/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitor's lifetime is marked by several important timescales: its infall, star-formation quenching, and stream-formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than $\sim 2\times10^6 M_\odot$. Finally, we show that the main factor controlling the rate of phase-mixing, and therefore fading, of tidal streams from satellite galaxies in MW-mass hosts is non-adiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed.

Aims. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light received from off-axis sources, in particular from circumstellar disks. Methods. Built on the recent PACO algorithm dedicated to the detection of point-like sources, the proposed method is based on the local learning of patch covariances capturing the spatial fluctuations of the stellar leakages. From this statistical modeling, we develop a regularized image reconstruction algorithm (REXPACO) following an inverse problem approach based on a forward image formation model of the off-axis sources in the ADI sequences. Results. Injections of fake circumstellar disks in ADI sequences from the VLT/SPHERE-IRDIS instrument show that both the morphology and the photometry of the disks are better preserved by REXPACO compared to standard postprocessing methods like cADI. In particular, the modeling of the spatial covariances proves usefull in reducing typical ADI artifacts and in better disentangling the signal of these sources from the residual stellar contamination. The application to stars hosting circumstellar disks with various morphologies confirms the ability of REXPACO to produce images of the light distribution with reduced artifacts. Finally, we show how REXPACO can be combined with PACO to disentangle the signal of circumstellar disks from the signal of candidate point-like sources. Conclusions. REXPACO is a novel post-processing algorithm producing numerically deblurred images of the circumstellar environment. It exploits the spatial covariances of the stellar leakages and of the noise to efficiently eliminate this nuisance term.

Due to the recent growth of discoveries of strong gravitational lensing (SGL) systems, one can statistically study both lens properties and cosmological parameters from 161 galactic scale SGL systems. We analyze meVSL model with the velocity dispersion of lenses by adopting the redshift and surface mass density depending power-law mass model. Analysis shows that meVSL models with various dark energy models including $\Lambda$CDM, $\omega$CDM, and CPL provide the negative values of meVSL parameter, $b$ when we put the prior to the $\Omega_{m 0}$ value from Planck. These indicate the faster speed of light and the stronger gravitational force in the past. However, if we adopt the WMAP prior on $\Omega_{m0}$, then we obtain the null results on $b$ within 1-$\sigma$ CL for the different dark energy models.

Fast Radio Bursts (FRBs) are extremely energetic pulses of millisecond duration and unknown origin. In order to understand the phenomenon that emits these pulses, targeted and untargeted searches have been performed for multi-wavelength counterparts, including the optical. The objective of this work is to search for optical transients at the position of 8 well-localized FRBs, after the arrival of the burst on different time-scales (typically at one day, several months, and one year after FRB detection) in order to compare with known transient optical light curves. We used the Las Cumbres Observatory Global Telescope Network (LCOGT), which allows us to promptly take images owing to its network of twenty-three telescopes working around the world. We used a template subtraction technique on all the images we collected at different epochs. We have divided the subtractions into two groups, in one group we use the image of the last epoch as a template and in the other group we use the image of the first epoch as a template. We have searched for bright optical transients at the localizations of the FRBs (<1 arcsec) in the template subtracted images. We have found no optical transients, so we have set limiting magnitudes of optical counterparts. Typical limiting magnitudes in apparent (absolute) magnitudes for our LCOGT data are ~22 (-19) mag in the r-band. We have compared our limiting magnitudes with light curves of superluminous supernovae (SLSNe), type Ia supernovae (SNe), supernovae associated with gamma-ray bursts (GRB SNe), a kilonova, and tidal disruption events (TDEs). We rule out that FRBs are associated with SLSN at a confidence of ~99.9%. We can also rule out the brightest sub-types of type Ia SNe, GRB SNe and TDEs (under some conditions) at similar confidence, though we cannot exclude scenarios where FRBs are associated with the faintest sub-type of each of these transient classes.

We analyze the observed spatial, chemical, and dynamical distributions of local metal-poor stars, based on photometrically derived metallicity and distance estimates along with proper motions from the Gaia mission. Throughout the sample volume ($|Z| \leq 6$ kpc) along the Galactic prime meridian, we identify stellar populations with distinct properties in the metallicity versus rotational velocity space, and model their phase-space distributions. Our decomposition provides refined positions of individual populations in this phase space identified in the previous papers of this series. Stars associated with Gaia Sausage/Enceladus (GSE) exhibit a larger proportion of metal-poor stars at greater distances from the Galactic center ($R_{\rm GC}$), with a slope of $\Delta \langle {\rm [Fe/H]} \rangle/\Delta R_{\rm GC} = -0.05\pm0.02$ dex kpc$^{-1}$. This observed trend, along with a mild anti-correlation of a mean rotational velocity with metallicity ($\Delta \langle v_\phi \rangle / \Delta \langle {\rm [Fe/H]} \rangle \sim -10$ km s$^{-1}$ dex$^{-1}$), implies that more metal-rich stars in the inner region of the GSE progenitor were gradually stripped away, while the prograde orbit of the merger at infall became radialized by dynamical friction. The metal-rich GSE stars are causally disconnected in our decomposition from the Splash structure, whose stars are mostly found in prograde orbits and exhibit a more centrally concentrated distribution than GSE. The metal-weak thick disk shows a similar spatial distribution as the Splash, suggesting earlier dynamical heating of stars in the primordial disk of the Milky Way.

We consider the dark matter (DM) scenarios consisting of the mixture of WIMPs and PBHs and study how much fraction of the total DM can be PBHs. In such scenarios, PBHs can accrete the WIMPs and consequently enhance the heating and ionization in the intergalactic medium due to WIMP annihilations. We demonstrate that the CMB data can give the stringent bounds on the allowed PBH fraction which are comparable or even tighter than those from the gamma ray data depending on the DM masses. For instance, the MCMC likelihood analysis using the Planck CMB data leads to the bound on PBH DM fraction with respect to the total dark matter $f_{\rm PBH} \lesssim {\cal O}( 10^{-10}\sim 10^{-8})$ for the WIMP mass $m_{\chi}\sim {\cal O}(10\sim 10^3)$ GeV with the conventional DM annihilation cross section $\langle \sigma v \rangle=3 \times 10^{-26}~\rm cm^3/s $. We also investigate the feasibility of the global 21-cm signal measurement to provide the stringent constraints on the PBH fraction.

Void-defect is a possible origin of ferromagnetic feature on pure carbon materials. In our previous paper, void-defect on graphene-nanoribbon show highly polarized spin configuration. In this paper, we studied cases for graphene molecules by quantum theory, by astronomical observation and by laboratory experiment. Model molecules for the density functional theory are graphene molecules of C23 and C53 induced by a void-defect. They have carbon pentagon ring within a hexagon network. Single void has three radical carbons, holding six spins. Those spins make several spin-states, which affects to molecular structure and molecular vibration, finally to infrared spectrum. The stable spin state was triplet, not singlet. This suggests magnetic pure carbon molecule. It was a surprise that those molecules show close infrared spectrum with astronomically observed one, especially observed on carbon rich planetary nebulae. We could assign major band at 18.9 micrometer, and sub-bands at 6.6, 7.0, 7.6, 8.1, 8.5, 9.0 and 17.4 micrometer. Also, calculated spectrum roughly coincides with that of laboratory experiment by the laser-induced carbon plasma, which is an analogy of cosmic carbon creation in interstellar space.

Void-defect induced magnetism of graphene molecule was recently reported in our previous paper of this series study. This paper investigated the case of hydrogenated graphene molecule, in chemical term, polycyclic aromatic hydrocarbon (PAH). Molecular infrared spectrum obtained by density functional theory was compared with astronomical observation. Void-defect on PAH caused serious structure change. Typical example of C23H12 had two carbon pentagon rings among hexagon networks. Stable spin state was non-magnetic singlet state. This is contrary to pure carbon case of C23, which show magnetic triplet state. It was discussed that Hydrogen played an important role to diminish magnetism by creating an SP3-bond among SP2-networks. Such a structure change affected molecular vibration and finally to photoemission spectrum in infrared region. The dication-C23H12 showed featured bands at 3.2, 6.3, 7.7, 8.6, 11.2, and 12.7 micrometer. It was surprising that those calculated bands coincided well with astronomically observed bands in many planetary nebulae. To confirm our study, large size molecule of C53H18 was studied. Calculation reproduced again similar astronomical bands. Also, small size molecule of C12H8 showed good coincidence with the spectrum observed for young stars. This paper would be the first report to indicate the specific PAH in space.

Inferring the properties of exoplanets from their atmospheres, while confronting low resolution and low signal-to-noise in the context of the quantities we want to derive, poses rigorous demands upon the data collected from observation. Further compounding this challenge is that inferences of exoplanet properties are built from forward models, which can include errors due to incomplete or inaccurate assumptions in atmospheric physics and chemistry. The confluence of observational noise and model error makes developing techniques to identify predictive features that are robust to both low s/n and model error increasingly important for exoplanet science. We demonstrate how both issues can be addressed simultaneously by taking advantage of underutilized multivariate information already present in current atmospheric models, including thermodynamic statistics and reaction network structure. To do so, we provide a case study of the prediction of vertical mixing (parameterized as eddy diffusion) in hot Jupiter atmospheres and show how prediction efficacy depends on what model information is used - e.g. chemical species abundances, network statistics, and/or thermodynamic statistics. We also show how the variables with the most predictive power vary with planetary properties such as temperature and metallicity. Our results demonstrate how inferences built on single metrics do not have utility across all possible use cases. We also show how statistical measures derived from network analyses tend to be better predictors when accounting for the possibility of missing data or observational uncertainty. We discuss future directions applying multivariate and network model information as a framework for increasing confidence in inferences aimed at extracting features relevant to exoplanet atmospheres and future applications to the detection of life on terrestrial worlds.

It has been recently shown that the halo near the Sun contains several kinematic substructures associated to past accretion events. For the more distant halo, there is evidence of large-scale density variations - in the form of stellar clouds or overdensities. We study the link between the local halo kinematic groups and three of these stellar clouds: the Hercules-Aquila cloud, the Virgo Overdensity, and the Eridanus-Phoenix overdensity. We perform orbital integrations in a standard Milky Way potential of a local halo sample extracted from GaiaEDR3, with the goal of predicting the location of the merger debris elsewhere in the Galaxy. We specifically focus on the regions occupied by the three stellar clouds and compare their kinematic and distance distributions with those predicted from the orbits of the nearby debris. We find that the local halo substructures have families of orbits that tend to pile up in the regions where the stellar clouds have been found. The distances and velocities of the cloud's member stars are in good agreement with those predicted from the orbit integrations, particularly for Gaia-Enceladus stars. This is the dominant contributor of all three overdensities, with a minor part stemming from the Helmi streams and to an even smaller extent from Sequoia. The orbital integrations predict no asymmetries in the sky distribution of halo stars, and they pinpoint where additional debris associated with the local halo substructures may be located.

The pulsar light curves and energy spectra are explored in dissipative pulsar magnetospheres with the Aristotelian electrodynamics (AE), where particle acceleration is fullly balanced with radiation reaction. The AE magnetospheres with non-zero pair multiplicity are computed by a pseudo-spectral method in the co-moving frame. The dissipative region near the current sheet outside the LC is accurately captured by the high-resolution simulation. The pulsar light curves and spectra are computed by the test particle trajectory method including the influence of both the consistent accelerating electric field and radiation reaction. Our results can generally reproduce the double-peak light curves and the GeV cut-off energy spectra in agreement with the Fermi observations for the pair multiplicity $\kappa\gtrsim1$.

The results of a long-term multiwavelength study of the powerful flat spectrum radio quasar 3C 454.3 using {\it Fermi}-LAT and Swift XRT/UVOT data are reported. In the $\gamma$-ray band, {\it Fermi}-LAT observations show several major flares when the source flux was $>10^{-5}\:{\rm photon\:cm^{-2}\:s^{-1}}$; the peak $\gamma$-ray flux above $141.6$ MeV, $(9.22\pm1.96)\times10^{-5}\:{\rm photon\:cm^{-2}\:s^{-1}}$ observed on MJD 55519.33, corresponds to $2.15\times10^{50}\:{\rm erg\:s^{-1}}$ isotropic $\gamma$-ray luminosity. The analysis of Swift XRT and UVOT data revealed a flux increase, although with smaller amplitudes, also in the X-ray and optical/UV bands. The X-ray emission of 3C 454.3 is with a hard spectral index of $\Gamma_{\rm X}=1.16-1.75$, and the flux in the flaring states increased up to $(1.80\pm0.18)\times10^{-10}{\rm erg\:cm^{-2}\: s^{-1}}$. Through combining the analyzed data, it was possible to assemble 362 high-quality and quasi-simultaneous spectral energy distributions of 3C 454.3 in 2008-2018 which all were modeled within a one-zone leptonic scenario assuming the emission region is within the broad line region, involving synchrotron, synchrotron self-Compton and external Compton mechanisms. Such an extensive modeling is the key for constraining the underlying emission mechanisms in the 3C 454.3 jet and allows to derive the physical parameters of the jet and investigate their evolution in time. The modeling suggests that during the flares, along with the variation of emitting electron parameters, the Doppler boosting factor increased substantially implying that the emission in these periods has most likely originated in a faster moving region.

The circumgalactic medium (CGM) connects the gas between the interstellar medium (ISM) and the intergalactic medium, which plays an important role in galaxy evolution. We use the stellar mass-metallicity relationship to investigate whether sharing the CGM will affect the distribution of metals in galaxy pairs. The optical emission lines from the Sloan Digital Sky Survey Data Release (SDSS DR7) are used to measure the gas-phase metallicity. We find that there is no significant difference in the distribution of the metallicity difference between two members in star forming-star forming pairs ($\rm \Delta log(O/H)_{diff}$), metallicity offset from the best-fitted stellar mass-metallicity relationship of galaxies in pairs ($\rm \Delta log(O/H)_{MS}$), as compared to "fake" pairs. By looking at $\rm \Delta log(O/H)_{diff}$ and $\rm \Delta log(O/H)_{MS}$ as a function of the star formation rate (SFR), specific star formation rate (sSFR), and stellar mass ratio, no difference is seen between galaxies in pairs and control galaxies. From our results, the share of the CGM may not play an important role in shaping the evolution of metal contents of galaxies.

Sub-Neptunes (Rp~1.25-4 REarth) remain the most commonly detected exoplanets to date. However, it remains difficult for observations to tell whether these intermediate-sized exoplanets have surfaces and where their surfaces are located. Here we propose that the abundances of trace species in the visible atmospheres of these sub-Neptunes can be used as proxies for determining the existence of surfaces and approximate surface conditions. As an example, we used a state-of-the-art photochemical model to simulate the atmospheric evolution of K2-18b and investigate its final steady-state composition with surfaces located at different pressures levels (Psurf). We find the surface location has a significant impact on the atmospheric abundances of trace species, making them deviate significantly from their thermochemical equilibrium and "no-surface" conditions. This result arises primarily because the pressure-temperature conditions at the surface determine whether photochemically-produced species can be recycled back to their favored thermochemical-equilibrium forms and transported back to the upper atmosphere. For an assumed H2-rich atmosphere for K2-18b, we identify seven chemical species that are most sensitive to the existence of surfaces: ammonia (NH3), methane (CH4), hydrogen cyanide (HCN), acetylene (C2H2), ethane (C2H6), carbon monoxide (CO), and carbon dioxide (CO2). The ratio between the observed and the no-surface abundances of these species, can help distinguish the existence of a shallow surface (Psurf < 10 bar), an intermediate surface (10 bar < Psurf < 100 bar), and a deep surface (Psurf > 100 bar). This framework can be applied together with future observations to other sub-Neptunes of interest.

MQ Dra is a strongly magnetic Cataclysmic Variable whose white dwarf accretes material from its secondary star through a stellar wind at a low rate. TESS observations were made of MQ Dra in four sectors in Cycle 2 and show a short duration, high energy flare (~10^35 erg) which has a profile characteristic of a flare from the M5V secondary star. This is one of the few occasions where an energetic flare has been seen from a Polar. We find no evidence that the flare caused a change in the light curve following the event and consider whether a coronal mass ejection was associated with the flare. We compare the frequency of energetic flares from the secondary star in MQ Dra with M dwarf stars and discuss the overall flare rate of stars with rotation periods shorter than 0.2 d and how such fast rotators can generate magnetic fields with low differential rotation rates.

We report on the follow-up XMM-Newton observation of the persistent X-ray pulsar CXOU J225355.1+624336, discovered with the CATS@BAR project on archival Chandra data. The source was detected at $f_{\rm X}$(0.5-10 keV) = 3.4$\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, a flux level which is fully consistent with the previous observations performed with ROSAT, Swift, and Chandra. The measured pulse period $P$ = 46.753(3) s, compared with the previous measurements, implies a constant spin down at an average rate $\dot P = 5.3\times 10^{-10}$ s s$^{-1}$. The pulse profile is energy dependent, showing three peaks at low energy and a less structured profile above about 3.5 keV. The pulsed fraction slightly increases with energy. We described the time-averaged EPIC spectrum with four different emission models: a partially covered power law, a cut-off power law, and a power law with an additional thermal component (either a black body or a collisionally ionized gas). In all cases we obtained equally good fits, so it was not possible to prefer or reject any emission model on the statistical basis. However, we disfavour the presence of the thermal components, since their modeled X-ray flux, resulting from a region larger than the neutron star surface, would largely dominate the X-ray emission from the pulsar. The phase-resolved spectral analysis showed that a simple flux variation cannot explain the source variability and proved that it is characterized by a spectral variability along the pulse phase. The results of the XMM-Newton observation confirmed that CXOU J225355.1+624336 is a BeXB with a low-luminosity ($L_{\rm X} \sim 10^{34-35}$ erg s$^{-1}$), a limited variability, and a constant spin down. Therefore, they reinforce the source classification as a persistent BeXB.

The first appearance of radio type II burst emission at decameter-hectometer (DH) waves typically occurs in connection, and often simultaneously, with other types of radio emissions. As type II bursts are signatures of propagating shock waves that are associated with flares and coronal mass ejections (CMEs), a rich variety of radio emissions can be expected. However, sometimes DH type II bursts appear in the dynamic spectra without other or earlier radio signatures. One explanation for them could be that the flare-CME launch happens on the far side of the Sun, and the emission is observed only when the source gets high enough in the solar atmosphere. In this study we have analysed 26 radio type II bursts that started at DH waves and were well-separated ('isolated') from other radio emission features. These bursts were identified from all DH type II bursts observed in 1998-2016, and for 12 events we had observations from at least two different viewing angles with the instruments onboard Wind and STEREO satellites. We found that only 30% of the type II bursts had their source origin on the far side of the Sun, but also that no bursts originated from the central region of the Sun (longitudes E30 - W40). Almost all of the isolated DH type II bursts could be associated with a shock near the CME leading front, and only few were determined to be shocks near the CME flank regions. In this respect our result differs from earlier findings. Our analysis, which included inspection of various CME and radio emission characteristics, suggests that the isolated DH type II bursts could be a special subgroup within DH type II bursts, where the radio emission requires particular coronal conditions to form and to die out.

We studied the multi-wavelength timing and spectral properties of the high mass X-ray binary MAXI J1348$-$630 during two successive outbursts of April and June 2019 (from MJD 58596 to MJD 58710) using ALMA, Swift, Chandra and NuSTAR. We studied the first broadband spectrum for this black hole that includes fluxes in radio, optical, ultraviolet and X-ray energy bands. We have used the high-resolution ALMA data to study the millimetre radio emission from compact jets in band 3, 4, 6 and 7 (89.56$-$351.44 GHz). The ultraviolet and X-ray counterparts are found using UVOT ($\lambda=$2246 Angstrom) and XRT (0.9$-$9 keV) onboard the Swift observatory. The X-ray spectrum is intensively studied using NuSTAR data in the range of 3$-$70 keV. We have studied the interday variation of spectral parameters of the source using three NuSTAR observations (from MJD 58655 to MJD 58672) during the hard state of the outburst of June. We have detected low-frequency quasi-periodic oscillation using NICER data on MJD 58654 with a significance of 6.8$\sigma$ at frequency 0.82 Hz in the hard state with RMS variability of 7.6 per cent. Another QPO has been detected at 0.67 Hz on MJD 58655 using NuSTAR data with a significance of 5.58$\sigma$ and 2.1 per cent RMS. The hardness ratio shows significant variation during the outburst of April but remained almost constant during the outburst of June. The spectral evolution in the hardness intensity diagram is studied during the outbursts. The position of the source is measured by Chandra that gives RA=13h48m12.878s, Dec=$-$63$^{\circ}$16'28.85" with enhanced accuracy.

The widely used Novikov-Thorne relativistic thin disc equations are only valid down to the radius of the innermost stable circular orbit (ISCO). This leads to an undetermined boundary condition at the ISCO, known as the inner stress of the disc, which sets the luminosity of the disc at the ISCO and introduces considerable ambiguity in accurately determining the mass, spin and accretion rate of black holes from observed spectra. We resolve this ambiguity by self-consistently extending the relativistic disc solution through the ISCO to the black hole horizon by calculating the inspiral of an average disc particle subject to turbulent disc forces, using a new particle-in-disc technique. Traditionally it has been assumed that the stress at the ISCO is zero, with material plunging approximately radially into the black hole at close to the speed of light. We demonstrate that in fact the inspiral is less severe, with several (~4-17) orbits completed before the horizon. This leads to a small non-zero stress and luminosity at and inside the ISCO, with a local surface temperature at the ISCO between ~0.15 and 0.3 times the maximum surface temperature of the disc, in the case where no dynamically important net magnetic field is present. For a range of disc parameters we calculate the value of the inner stress/surface temperature, which is required when fitting relativistic thin disc models to observations. We resolve a problem in relativistic slim disc models in which turbulent heating becomes inaccurate and falls to zero inside the plunging region.

In this work we present a case that Microscopic Black Holes (MBH) of mass $10^{16}\ kg$ to $3 \cdot 10^{19}\ kg$ experience acceleration as they move within stellar material at low velocities. The accelerating forces are caused by the fact that a MBH moving through stellar material leaves a trail of hot rarified gas. The rarified gas behind a MBH exerts lower gravitational force on the MBH than the dense gas in front of it. The accelerating forces exceed the gravitational drag forces when MBH moves at Mach number $\mathcal{M}<\mathcal{M}_0$. The equilibrium Mach number $\mathcal{M}_0$ depends on MBH mass and stellar material characteristics. Our calculations open the possibility of MBH orbiting within stars including Sun at Mach number $\mathcal{M}_0$. At the end of this work we list some unresolved problems which result from our calculations.

CubeSats have become a meaningful option for deep-space exploration, but their autonomy must be increased to maximize the science return while limiting the complexity in operations. We present here a solution for an autonomous orbit determination in the context of a CubeSat cruising in deep space. The study case is a journey from Earth to Mars. An optical sensor at CubeSat standard is considered. The image processing is added to extract the direction of distant celestial bodies with 0.2 arcsec accuracy: it consists of a Multiple Cross-Correlation (MCC) algorithm that uses bright stars in the background of the images. Then, an Unscented Kalman Filter (UKF) is built to perform an asynchronous triangulation from the successive directions of the celestial bodies. The UKF meets the expected performance in contexts where linear approximations are not possible. The orbit reconstruction reaches a 3-sigma accuracy of 30 km in the middle of the Earth-Mars cruise. Additionally, the CPU cost of the filter is assessed at less than 1 second per iteration with a typical CubeSat hardware. It is ready for further improvements in terms of new observables associated with data fusion, quicker convergence and attitude control savings.

We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Parkes 64 m telescopes, spanning a frequency range from $\sim$150 MHz to 4 GHz. The pulsar is faint, with an estimated flux density ($S$) of $\sim$1 mJy at 400 MHz and a spectrum $S(\nu)\,\propto\,\nu^{-2.0 \pm 0.2}$, where $\nu$ is frequency. The DM-derived distance implies that it is also a low-luminosity source ($\sim$ 0.1 ${\rm mJy\,kpc^2}$ at 1400 MHz). The analysis of archival MWA observations reveals that the pulsar's mean flux density varies by up to a factor of $\sim$5-6 on timescales of several weeks to months. By combining MWA and uGMRT data, the pulsar position was determined to arcsecond precision. We also report on polarization properties detected in the MWA and Parkes bands. The pulsar's nondetection in previous pulsar and continuum imaging surveys, the observed high variability, and its detection in a small fraction of the survey data searched to date, all hint at a larger population of pulsars that await discovery in the southern hemisphere, with the MWA and the future low-frequency Square Kilometre Array.

Small scale jet-induced erosion experiments are useful for identifying the scaling of erosion with respect to the various physical parameters (gravity, grain size, gas velocity, gas density, grain density, etc.), and because they provide a data set for benchmarking numerical flow codes. We have performed experiments varying the physical parameters listed above (e.g., gravity was varied in reduced gravity aircraft flights). In all these experiments, a subsonic jet of gas impinges vertically on a bed of sand or lunar soil simulant forming a localized scour hole beneath the jet. Videography captures the erosion and scour hole formation processes, and analysis of these videos post-test identifies the scaling of these processes. This has produced important new insights into the physics of erosion. Based on these insights, we have developed an erosion rate model that can be applied to generalized situations, such as the erosion of soil beneath a horizontal gas flow on a planetary surface. This is important to lunar exploration because the rate of erosion beneath the rocket exhaust plume of a landing spacecraft will determine the amount of sand-blasting damage that can be inflicted upon surrounding hardware. Although the rocket exhaust plume at the exit of the nozzle is supersonic, the boundary layer on the lunar surface where erosion occurs is subsonic. The model has been benchmarked through comparison with the Apollo landing videos, which show the blowing lunar soil, and computational fluid dynamics simulations of those landings.

Stellar winds from fast-rotating Population III (Pop III) stars have long been suspected to make important contributions to early metal enrichment, as features in the nucleosynthesis of such `spinstars' are consistent with the chemical abundance patterns of some metal-poor stars in the local Universe. Particularly, stellar winds rich in light elements can provide another pathway towards explaining the carbon enhancement in carbon-enhanced metal-poor (CEMP) stars. In this work, we focus on the feedback of Pop III stellar winds combined with supernovae (SNe), and derive the resulting chemical signatures in the enriched medium. We explore a large parameter space of Pop III star formation and feedback with a semi-analytical model. The predicted pattern of carbon and iron abundances of second-generation stars agrees well with observations of CEMP-no stars ($[\rm Ba/Fe]<0$) at $[\rm Fe/H]\lesssim -3$ and $A(\mathrm{C})\lesssim 7$, under the optimistic assumption of significant mass loss by winds. In this scenario, carbon-rich but iron-free second-generation stars can form in systems dominated by enrichment from winds, gaining trace amounts of iron by accretion from the interstellar medium, to become the most iron-poor and carbon-enhanced stars seen in observations ($[\rm Fe/H]\lesssim -4$, $[\rm C/Fe]\gtrsim 2$). Wind feedback from Pop III spinstars deserves more detailed modelling in early cosmic structure formation.

The age of a supernova remnant (SNR) is, though undoubtedly one of the most important properties for study of its evolution, difficult to estimate reliably in most cases. In this study, we compare the dynamical and plasma ages of the SNRs and characteristic ages of their associated pulsars with the corresponding SNRs' ages that are generally thought to be reliable ($t_{\rm r}$): historical and light-echo ages of the SNRs, kinematic ages of the ejecta knots and kinematic ages of the associated neutron stars (NS). The kinematic age of ejecta knots or a NS is the time that they have taken to reach the current positions from the explosion center. We use all of the available 24 systems for which $t_{\rm r}$ is already available (historical, light-echo, and ejecta kinematic ages) or measurable (NS kinematic age). We estimate the NS kinematic ages for eight SNR-NS systems by determining quantitatively the geometric centers of the SNR shells. The obtained $t_{\rm r}$ ranges from 33 yr to $\approx 400$ kyr. We find that the two SNR ages, dynamical and plasma ages, are consistent with $t_{\rm r}$ within a factor of four, whereas the characteristic ages of the pulsars differ from $t_{\rm r}$ by more than a factor of four in some systems. Using the $t_{\rm r}$ summarized in this work, we present the initial spin periods of the associated pulsars, which are more strictly constrained than the previous works, as well.

We report the detection of a massive neutral gas outflow in the z=2.09 gravitationally lensed Dusty Star Forming Galaxy HATLASJ085358.9+015537 (G09v1.40), seen in absorption with the OH+(1_1-1_0) transition using spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations. The OH+ line, tracing diffuse atomic gas, is observed simultaneously with the CO(9-8) emission line and underlying dust continuum, providing kinematic and morphological information of the warm dense gas and dust in the host galaxy, respectively. Source plane reconstruction reveals a compact nuclear dust region from which the conical OH+ outflow emerges, away from the more extended and offset CO(9-8) and stellar distributions. The total atomic gas mass of the observed outflow is 6.7x10^9 M_sun>25 as massive as the molecular gas component within the galaxy. We derive outflow properties over a range of possible conical outflow inclinations (0.38 deg-64 deg), finding a neutral gas mass outflow rate between 83-25400 M_sun/yr, exceeding the star formation rate (788 M_sun/yr) if the inclination is >3.6 deg (mass-loading factor = 0.3-4.7). Kinetic energy and momentum fluxes are within 4.4-290x10^9 L_sun and 0.1-3.7x10^37 dyne, respectively (energy-loading factor = 0.013-16). If the inclination is small (<12 deg), star formation may be responsible for driving the neutral outflow via energy- or momentum-driven scenarios, depending on its orientation. If the inclination is >12 deg, unusually high coupling efficiencies to the ISM, or an additional energy source (e.g. an active galactic nucleus), is needed.

$\alpha$ Centauri A is the closest solar-type star to the Sun and offers the best opportunity to find and ultimately to characterize an Earth-sized planet located in its Habitable Zone (HZ). Here we describe initial results from an ALMA program to search for planets in the $\alpha$ Cen AB system using differential astrometry at millimeter wavelengths. Our initial results include new absolute astrometric measurements of the proper motion, orbital motion, and parallax of the $\alpha$ Cen system. These lead to an improved knowledge of the physical properties of both $\alpha$ Cen A and B. Our estimates of ALMA's relative astrometric precision suggest that we will ultimately be sensitive to planets of a few 10s of Earth mass in orbits from 1-3 AU, where stable orbits are thought to exist.

Using two sets of large $N$-body simulations, we study the origin of the correlations of halo assembly time ($z_{\rm f}$), concentration ($v_{\rm max}/v_{\rm 200}$) and spin ($\lambda$) with the large-scale density field around halos. We find that the correlations of the three halo properties with the large-scale density at $z=0$ are the secondary effects of their correlations with the initial linear density field. Using the linear density on different scales, we find two types of correlations. The $L$-correlation, which describes the correlation of halo properties with the mean linear over-density $\delta^{\rm i}_{\rm L}$ within the halo Lagrangian radius $R_{\rm L}$, is positive for both $z_{\rm f}$ and $v_{\rm max}/v_{\rm 200}$, and negative for $\lambda$. The $E$-correlation, which describes the correlations of halo properties with $\delta^{\rm i}(R/R_{\rm L}>1)$ for given $\delta^{\rm i}_{\rm L}$, shows trends opposite to the $L$-correlation. Both of the $E$- and $L$-correlations depend only weakly on halo mass, indicating a similar origin for halos of different masses. The dependence of the halo bias on the the three halo properties can be well explained by the competition of the $E$- and $L$-correlations and the correlationof the linear density field on different scales. These two types of correlation together can establish the complex halo-mass dependence of the clustering bias produced by the three halo properties.

The magnetic switchbacks observed recently by the Parker Solar Probe have raised the question about their nature and origin. One of the competing theories of their origin is the interchange reconnection in the solar corona. In this scenario, switchbacks are generated at the reconnection site between open and closed magnetic fields, and are either advected by an upflow or propagate as waves into the solar wind. In this paper we test the wave hypothesis, numerically modelling the propagation of a switchback, modelled as an embedded Alfv\'en wave packet of constant magnetic field magnitude, through the gravitationally stratified solar corona with different degrees of background magnetic field expansion. While switchbacks propagating in a uniform medium with no gravity are relatively stable, as reported previously, we find that gravitational stratification together with the expansion of the magnetic field act in multiple ways to deform the switchbacks. These include WKB effects, which depend on the degree of magnetic field expansion, and also finite-amplitude effects, such as the symmetry breaking between nonlinear advection and the Lorentz force. In a straight or radially expanding magnetic field the propagating switchbacks unfold into waves that cause minimal magnetic field deflections, while a super-radially expanding magnetic field aids in maintaining strong deflections. Other important effects are the mass uplift the propagating switchbacks induce and the reconnection and drainage of plasmoids contained within the switchbacks. In the Appendix, we examine a series of setups with different switchback configurations and parameters, which broaden the scope of our study.

Massive scalar fields provide excellent dark matter candidates, whose dynamics are often explored analytically and numerically using nonrelativistic Schr\"{o}dinger-Poisson (SP) equations in a cosmological context. In this paper, starting from the nonlinear and fully relativistic Klein-Gordon-Einstein (KGE) equations in an expanding universe, we provide a systematic framework for deriving the SP equations, as well as relativistic corrections to them, by integrating out `fast modes' and including nonlinear metric and matter contributions. We provide explicit equations for the leading-order relativistic corrections, which provide insight into deviations from the SP equations as the system approaches the relativistic regime. Upon including the leading-order corrections, our equations are applicable beyond the domain of validity of the SP system, and are simpler to use than the full KGE case in some contexts. As a concrete application, we calculate the mass-radius relationship of solitons in scalar dark matter and accurately capture the deviations of this relationship from the SP system towards the KGE one.

The progenitors of present-day galaxy clusters give important clues about the evolution of the large scale structure, cosmic mass assembly, and galaxy evolution. Simulations are a major tool for these studies since they are used to interpret observations. In this work, we introduce a set of "protocluster-lightcones", dubbed PCcones. They are mock galaxy catalogs generated from the Millennium Simulation with the L-GALAXIES semi-analytic model. These lightcones were constructed by placing a desired structure at the redshift of interest in the centre of the cone. This approach allows to adopt a set of observational constraints, such as magnitude limits and uncertainties in magnitudes and photometric redshifts (photo-zs), to produce realistic simulations of photometric surveys. We show that photo-zs obtained with PCcones are more accurate than those obtained directly with the Millennium Simulation, mostly due to the difference in how apparent magnitudes are computed. We apply PCcones in the determination of the expected accuracy of protocluster detection using photo-zs in the $z=1-3$ range in the wide-layer of HSC-SSP and the 10-year LSST forecast. With our technique, we expect to recover only $\sim38\%$ and $\sim 43\%$ of all massive galaxy cluster progenitors with more than 70\% of purity for HSC-SSP and LSST, respectively. Indeed, the combination of observational constraints and photo-z uncertainties affects the detection of structures critically for both emulations, indicating the need of spectroscopic redshifts to improve detection. We also compare our mocks of the Deep CFHTLS at $z<1.5$ with observed cluster catalogs, as an extra validation of the lightcones and methods.

Recent astronomical observations of both isomers E and Z of imines such as cyanomethanimine, ethanimine and 2-propyn-1-imine, have revealed that the abundances in the ISM of these isomers differ by factors of ~3-10. Several theories have been proposed to explain the observed behavior, but none of them successfully explains the [E]/[Z] ratios. In this work we present a detailed study of the kinetics of the one-step E-Z isomerization reactions of cyanomethanimine, ethanimine and 2-propyn-1-imine under interstellar conditions (in the 10-400 K temperature range). This reaction was previously thought to be non-viable in the ISM due to its associated high-energy barrier (about 13,000 K). In this Letter, we show that considering the multidimensional small curvature tunneling approximation, the tunneling effect enables the isomerization even at low temperatures. This is due to the fact that the representative tunneling energy lies in the vibrational ground state of the least stable isomer up to approximately 150 K, making the reaction constants of the isomerization from the least stable to the most stable isomer basically constant. The predicted [E]/[Z] ratios are almost the same as those reported from the astronomical observations for all imines observed. This study demonstrates that the [E]/[Z] ratio of imines in the ISM strongly depends on their relative stability.

We propose a new approach to constant-roll warm inflation as a generalization of constant-roll inflation. Based on the $\beta$-function formalism, it is shown that constant-roll warm inflation models with a natural end fall into universality classes defined by three different types of $\beta$-functions, under the assumption that radiation energy density is quasi-stable. Given that warm inflation is completely specified by the $\beta$-function and dissipation coefficient ratio, we investigate whether or not the inflation can physically be realized for enough number of e-foldings at the background level for some combinations of the $\beta$-functions and non-trivial dissipation coefficient ratios.

In gravitational-wave data analysis, we regularly work with a host of non-trivial prior probabilities on compact binary masses, redshifts, and spins. We must regularly manipulate these priors, computing the implied priors on a transformed basis of parameters or reweighting posterior samples from one prior to another. Here, I detail some common manipulations, presenting a table of Jacobians with which to transform priors between mass parametrizations, describing the conversion between source- and detector-frame priors, and deriving analytic expressions for priors on the "effective spin" parameters regularly invoked in gravitational-wave astronomy.

The axion-electron coupling $g_{ae}$ is a generic feature of non-hadronic axion models. This coupling may induce a variety of observable signatures, particularly in astrophysical environments. Here, we revisit the calculation of the axion-electron bremsstrahlung and provide a general formulation valid for a non-relativistic plasma with any level of degeneracy and for any axion mass. We apply our result to the Sun, red giant stars and white dwarfs. The most relevant impact of the corrections is in the core of red giant stars, where we find a reduction of $\sim25\%$ in the axion emissivity with respect to the calculation in the completely degenerate approximation, showing the need for a reliable prescription valid at intermediate degeneracy.

We investigate the null geodesic flow and in particular the light rings (LRs), fundamental photon orbits (FPOs) and shadows of a black hole (BH) immersed in a strong, uniform magnetic field, described by the Schwarzschilld-Melvin electrovacuum solution. The empty Melvin magnetic Universe contains a tube of planar LRs. Including a BH, for weak magnetic fields, the shadow becomes oblate, whereas the intrinsic horizon geometry becomes prolate. For strong magnetic fields (over-critical solutions), there are no LRs outside the BH horizon, a result explained using topological arguments. This feature, together with the light confining structure of the Melvin Universe yields \textit{panoramic shadows}, seen (almost) all around the equator of the observer's sky. Despite the lack of LRs, there are FPOs, including polar planar ones, which define the shadow edge. We also observe and discuss chaotic lensing, including in the empty Melvin Universe, and multiple disconnected shadows.

Void defect is a possible origin of ferromagnetic like feature of pure carbon material. Applying density functional theory to void defect induced graphene nano ribbon (GNR), a detailed relationship between multiple spin state and structure change was studied. An equitorial triangle of an initial initial void having six electrons is distorted to isosceles triangle by rebonding carbon atoms. Among possible spin states, the most stable state was Sz=2/2. The case of Sz=4/2 is remarkable that initial flat ribbon turned to three dimentional curled one having highly polarized spin configuration at ribbon edges. Total energy of Sz=4/2 was very close to that of Sz=2/2, which suggests coexistence of flat and curled ribbons. As a model of three dimensional graphite, bilayered AB stacked GNR was analyzed. Spin distribution was limited to the void created layer. Distributed void triangle show 60 degree clockwise rotation for differrent site void, which was consistent with experimental observation using the scanning tunneling microscope. (To be published on Journal of the Magnetic Society of Japan, 2021 )

Ultralight axion like particles (ALPs) of mass $m_a\in (10^{-21}\rm{eV}-10^{-22}\rm{eV})$ with axion decay constant $f_a\sim 10^{17}\rm{GeV}$ can be candidates for fuzzy dark matter (FDM). If celestial bodies like Earth and Sun are immersed in a low mass axionic FDM potential and if the ALPs have coupling with nucleons then the coherent oscillation of the axionic field results a long range axion hair outside of the celestial bodies. The range of the axion mediated Yukawa type fifth force is determined by the distance between the Earth and the Sun which fixes the upper bound of the mass of axion as $m_a\lesssim10^{-18}\rm{eV}$. The long range axionic Yukawa potential between the Earth and Sun changes the gravitational potential between them and contribute to the light bending and the Shapiro time delay. From the observational uncertainties of those experiments, we put an upper bound on the axion decay constant as $f_a\lesssim 9.85\times 10^{6}\rm{GeV}$, which is the stronger bound obtained from Shapiro time delay. This implies if ALPs are FDM, then they do not couple to nucleons.

Outer crusts of neutron stars and interiors of cool white dwarfs consist of bare atomic nuclei, arranged in a crystal lattice and immersed in a Fermi gas of degenerate electrons. We study electrostatic properties of such Coulomb crystals, taking into account the polarizability of the electron gas and considering different lattice structures, which can form the ground state. To take the electron background polarization into account, we use the linear response theory with the electron dielectric function given either by the Thomas-Fermi approximation or by the random-phase approximation (RPA). We compare the widely used nonrelativistic (Lindhard) version of the RPA approximation with the more general, relativistic (Jancovici) version. The results of the different approximations are compared to assess the importance of going beyond the Thomas-Fermi or Lindhard approximations. We also include contribution of zero-point vibrations of ions into the ground-state energy. We show that the bcc lattice forms the ground state for any charge number $Z$ of the atomic nuclei at the densities where the electrons are relativistic ($\rho\gtrsim10^6$ g cm$^{-3}$), while at lower densities the fcc and hcp lattices can form the ground state. The MgB$_2$-like lattice never forms the ground state at realistic densities in the crystallized regions of degenerate stars. The RPA corrections strongly affect the boundaries between the phases. As a result, transitions between different ground-state structures depend on $Z$ in a nontrivial way. The relativistic and quantum corrections produce less dramatic effects, moderately shifting the phase boundaries.

UCGretina, a GEANT4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton continuum. Both of these features are needed in order to accurately extract gamma-ray yields from spectra collected in in-beam gamma-ray spectroscopy measurements with beams traveling at $v/c \gtrsim 0.3$ at the National Superconducting Cyclotron Laboratory and the Facility for Rare Isotope Beams. In the process of developing the model, we determined that millimeter-scale layers of passive germanium surrounding the active volumes of the simulated crystals must be included in order to reproduce measured photopeak efficiencies. We adopted a simple model of effective passive layers and developed heuristic methods of determining passive-layer thicknesses by comparison of simulations and measurements for a single crystal and for the full array. Prospects for future development of the model are discussed.

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EoS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test we quantify the ability of current and future gravitational-wave observations to discriminate among families of realistic EoS which differ in particle content and ab-initio microscopic calculations. While the constraining power of GW170817 is limited, stringent constraints can be placed with approximately 10 coalescences detected by LIGO-Virgo at design sensitivity, but only for relatively stiff EoS which are already marginally in tension with GW170817. However, we show that even just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EoS with very strong statistical significance, and can discriminate among models which feature similar softness, hence constraining the properties of nuclear matter to unprecedented levels.

Guy Boistel questions the Sun / variable stars analogy at the heart of the work of Dutch astronomer Albert Brester. The problem of interdisciplinarity appears in the objectification of the Sun at the end of the 19th century: it is by mixing astronomy, physics and chemistry that Brester tried to apply knowledge of the solar atmosphere to the understanding of variable stars, a fairly confidential field of study still at that time.

The optimal lattice quantizer is the lattice which minimizes the (dimensionless) second moment $G$. In dimensions $1$ to $8$, it has been proven that the optimal lattice quantizer is one of the classical lattices, or there is good numerical evidence for this. In contrast, more than two decades ago, convincing numerical studies showed that in dimension $9$, a non-classical lattice is optimal. The structure and properties of this lattice depend upon a single positive real parameter $a$, whose value was only known approximately. Here, for $a^2 < 1/2$, we give an exact analytic description of this one-parameter family of lattices and their Voronoi cells, and calculate their second moment, which is a $19$th order polynomial in $a$. This allows us to determine the exact value of $a$ which minimizes $G$. It is an algebraic number, defined by the root of a $9$th order polynomial, with $a \approx 0.573223794$. We also show that for this value of $a$, the covariance matrix (second moment tensor) is proportional to the identity, consistent with a theorem of Zamir and Feder for optimal quantizers. The same method can be used for arbitrary one-parameter families of laminated lattices, so may provide a useful tool to identify optimal quantizers in other dimensions as well.

In general relativity (GR), gravitational waves (GWs) propagate the well-known plus and cross polarization modes which are the signature of a massless spin-2 field. However, diffraction of GWs caused by intervening objects along the line of sight can cause the apparent rise of additional polarizations due to GW-curvature interactions. In this paper, we continue the analysis by two of the authors of the present article, on lensing of gravitational waves beyond geometric optics. In particular, we calculate the lensing effect caused by a point-like lens, in the regime where its Schwarzschild radius $R_s$ is much smaller than the wavelength $\lambda$ of the signal, itself smaller than the impact parameter $b$. In this case, the curvature of spacetime induces distortions in the polarization of the wave such that effective scalar and vector polarizations may appear. We find that the amplitude of these apparent non-GR polarizations is suppressed by a factor $R_s\lambda/b^2$ with respect to the amplitude of the GR-like tensor modes. We estimate the probability to develop these extra polarization modes for a nearly monochromatic GW in the Pulsar Timing Arrays band traveling through a distribution of galaxies.