Papers for Tuesday, Aug 23 2022

A red giant star is an evolved low- or intermediate-mass star that has exhausted its central hydrogen content, leaving a helium core and a hydrogen-burning shell. Oscillations of stars can be observed as periodic dimmings and brightenings in the optical light curves. In red giant stars, non-radial acoustic waves couple to gravity waves and give rise to mixed modes, which behave as pressure (p) modes in the envelope and gravity (g) modes in the core. These modes were previously used to measure the internal rotation of red giants, leading to the conclusion that purely hydrodynamical processes of angular momentum transport from the core are too inefficient. Magnetic fields could produce the additional required transport. However, due to the lack of direct measurements of magnetic fields in stellar interiors, very little is currently known about their properties. Asteroseismology can provide direct detection of magnetic fields because, like rotation, the fields induce shifts in the oscillation mode frequencies. Here we report the measurement of magnetic fields in the cores of three red giant stars observed with the Kepler satellite. The fields induce shifts that break the symmetry of dipole mode multiplets. We thus measure field strengths ranging from ~30 to ~100 kG in the vicinity of the hydrogen-burning shell and place constraints on the field tolopolgy.

In recent years, high-angular-resolution observations of the dust and gas in circumstellar discs have revealed a variety of morphologies, naturally triggering the question of whether these substructures are driven by forming planets. While it remains difficult to directly image embedded planets, a promising method to distinguish disc-shaping mechanisms is to study the gas kinematics as characterising deviations from Keplerian rotation can be used to probe underlying perturbations such as planets. Creating spiral structures, the latter can also be traced in the brightness temperature. Here we analyse the brightness temperatures and kinematics of a sample of 36 transition discs observed with ALMA to search for substructures possibly tracing companions. We use archival Band 6 and 7 observations of different CO isotopologues and fit Keplerian disc models to the velocity fields. After subtraction of an azimuthally averaged brightness temperature and Keplerian rotation model from the data, we find significant substructures in both residuals of eight sources. Other sources show tentative features, while half of our sample does not show any substructures that may be indicative of planet-disc interactions. For the first time, we compare the substructures from our analysis with various other indicators of planets. About 20% of discs show strong features such as arcs or spirals, possibly associated with the presence of planets, while the majority do not present as clear planet-driven signatures. Almost all discs that exhibit spirals in near-infrared scattered light show at least tentative features in the CO data. The present data are able to reveal only very massive bodies and a lack of features may suggest that, if there are planets, they are of lower mass (<1-3Mj) or located closer to the star within deep cavities. Dedicated observations and modelling efforts are needed to confirm such scenarios.

Young M-type binaries are particularly useful for precise isochronal dating by taking advantage of their extended pre-main sequence evolution. Orbital monitoring of these low-mass objects becomes essential in constraining their fundamental properties, as dynamical masses can be extracted from their Keplerian motion. Here, we present the combined efforts of the AstraLux Large Multiplicity Survey, together with a filler sub-programme from the SpHere INfrared Exoplanet (SHINE) project and previously unpublished data from the FastCam lucky imaging camera at the Nordical Optical Telescope (NOT) and the NaCo instrument at the Very Large Telescope (VLT). Building on previous work, we use archival and new astrometric data to constrain orbital parameters for 20 M-type binaries. We identify that eight of the binaries have strong Bayesian probabilities and belong to known young moving groups (YMGs). We provide a first attempt at constraining orbital parameters for 14 of the binaries in our sample, with the remaining six having previously fitted orbits for which we provide additional astrometric data and updated Gaia parallaxes. The substantial orbital information built up here for four of the binaries allows for direct comparison between individual dynamical masses and theoretical masses from stellar evolutionary model isochrones, with an additional three binary systems with tentative individual dynamical mass estimates likely to be improved in the near future. We attained an overall agreement between the dynamical masses and the theoretical masses from the isochrones based on the assumed YMG age of the respective binary pair. The two systems with the best orbital constrains for which we obtained individual dynamical masses, J0728 and J2317, display higher dynamical masses than predicted by evolutionary models.

Double white dwarf (DWD) mergers are possibly the leading formation channel of massive, rapidly rotating, high-field magnetic white dwarfs (HFMWDs). However, the direct link connecting a DWD merger to any observed HFMWD is still missing. We here show that the HFMWDs SDSS J221141.80+113604.4 (hereafter J2211+1136) and ZTF J190132.9+145808.7 (hereafter J1901+1458), might be DWD merger products. J2211+1136 is a $1.27\, M_\odot$ WD with a rotation period of $70.32$ s and a surface magnetic field of $15$ MG. J1901+1458 is a $1.35\, M_\odot$ WD with a rotation period of $416.20$ s, and a surface magnetic field in the range $600$-$900$ MG. With the assumption of single-star evolution, the currently measured WD masses and surface temperatures, the cooling ages of J2211+1136 and J1901+1458 are, respectively, $2.61$-$2.85$ Gyr and $10$-$100$ Myr. We hypothesize that these WDs are DWD merger products and compute the evolution of the post-merged configuration formed by a central WD surrounded by a disk. We show that the post-merger system evolves through three phases depending on whether accretion, mass ejection (propeller), or magnetic braking dominates the torque onto the central WD. We calculate the duration of these phases and find that the WD rotational age, i.e., the total time elapsed since the merger to the instant where the WD central remnant reaches the current measured rotation period, agrees with the cooling age. We infer the value of the accretion rate, the disk mass, and the value of the primary and secondary WD components of the DWD merger that lead to a post-merger evolution consistent with the observations.

The Solar Dynamics Observatory (SDO), a NASA multi-spectral decade-long mission that has been daily producing terabytes of observational data from the Sun, has been recently used as a use-case to demonstrate the potential of machine learning methodologies and to pave the way for future deep-space mission planning. In particular, the idea of using image-to-image translation to virtually produce extreme ultra-violet channels has been proposed in several recent studies, as a way to both enhance missions with less available channels and to alleviate the challenges due to the low downlink rate in deep space. This paper investigates the potential and the limitations of such a deep learning approach by focusing on the permutation of four channels and an encoder--decoder based architecture, with particular attention to how morphological traits and brightness of the solar surface affect the neural network predictions. In this work we want to answer the question: can synthetic images of the solar corona produced via image-to-image translation be used for scientific studies of the Sun? The analysis highlights that the neural network produces high-quality images over three orders of magnitude in count rate (pixel intensity) and can generally reproduce the covariance across channels within a 1% error. However the model performance drastically diminishes in correspondence of extremely high energetic events like flares, and we argue that the reason is related to the rareness of such events posing a challenge to model training.

The Epoch of Reionization Spectrometer (EoR-Spec) will be an instrument module for the Prime-Cam receiver on the CCAT-prime Collaboration's Fred Young Submillimeter Telescope (FYST), a 6-m primary mirror Crossed Dragone telescope. With its Fabry-Perot interferometer (FPI), EoR-Spec will step through frequencies between 210 and 420 GHz to perform line intensity mapping of the 158 $\mu$m [CII] line in aggregates of star-forming galaxies between redshifts of 3.5 and 8 to trace the evolution of structure in the universe during the epoch of reionization. Here we present the optical design of the module including studies of the optical quality and other key parameters at the image surface. In order to achieve the required resolving power (R$\sim$100) with the FPI, it is important to have a highly collimated beam at the Lyot stop of the system; the optimization process to achieve this goal with four lenses instead of three as used in other Prime-Cam modules is outlined. As part of the optimization, we test the effect of replacing some of the aspheric lenses with biconic lenses in this Crossed Dragone design and find that the biconic lenses tend to improve the image quality across the focal plane of the module.

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors' third observing run and an archival search on the 80 confident events reported in GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly 200 stars and planets. The observatory consists of a 60 cm, f/10 telescope paired to a long-slit spectrograph, yielding simultaneous, almost continuous coverage between 1203 {\AA} and 4000 {\AA}, with resolutions ranging from 6000 to 240. The efficient instrument provides throughputs > 4% (far-UV; FUV) and > 15% (near-UV; NUV), comparable to HST/COS and much better than HST/STIS, over the same spectral range. A key design feature is the LiF prism, which serves as a dispersive element and provides high throughput even after accounting for radiation degradation. The use of two delta-doped Electron-Multiplying CCD detectors with UV-optimized, single-layer anti-reflection coatings provides high quantum efficiency and low detector noise. From the Earth-Sun second Lagrangian point, UV-SCOPE will continuously observe planetary transits and stellar variability in the full FUV-to-NUV range, with negligible astrophysical background. All these features make UV-SCOPE the ideal instrument to study exoplanetary atmospheres and the impact of host stars on their planets. UV-SCOPE was proposed to NASA as a Medium Explorer (MidEx) mission for the 2021 Announcement of Opportunity. If approved, the observatory will be developed over a 5-year period. Its primary science mission takes 34 months to complete. The spacecraft carries enough fuel for 6 years of operations.

Non-thermal emission from clusters of galaxies at the high-energy X-ray regime has been searched with various instruments, but the detection significance of this emission has yet been found to be either marginal or controversial. Taking advantage of NuSTAR's unique capability to focus X-rays in the hard energy band, we present a detailed analysis of 238 ks NuSTAR observations of the merging galaxy cluster SPT-CL J2031-4037, searching for non-thermal inverse Compton emission. Our spectral analysis of SPT-CL J2031-4037 shows a possibility that the hard X-ray emission of the cluster can be described by a non-thermal component, though we cannot completely rule out a purely thermal origin for this hard emission. Including the statistical and systematic uncertainties, our best model fit yields a 20-80 keV non-thermal flux of $3.93_{-1.10}^{+1.24} \times 10^{-12}$ erg s$^{-1}$ cm$^{-2}$. The estimated non-thermal flux is comparable to those found in other galaxy clusters using NuSTAR and other X-ray instruments. Using this non-thermal flux with the existing radio data of the cluster, we estimate a volume-averaged magnetic field strength in the range of around 0.1-0.2 $\mu$G.

Prime-Cam is a first-generation instrument for the Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) Facility. The 850$~$GHz module for Prime-Cam will probe the highest frequency of all the instrument modules. We describe the parameter space of the 850$~$GHz optical system between the F$\lambda$ spacing, beam size, pixel sensitivity, and detector count. We present the optimization of an optical design for the 850$~$GHz instrument module for CCAT-prime. We further describe the development of the cryogenic RF chain design to accommodate $>$30 readout lines to read 41,400 kinetic inductance detectors (KIDs) within the cryogenic testbed.

We report on the effects of strong magnetic fields on neutrino emission in the modified Urca process. We show that the effect of Landau levels on the various Urca pairs affects the neutrino emission spectrum and leads to an angular asymmetry in the neutrino emission. For low magnetic fields the Landau levels have almost no effect on the cooling. However, as the field strength increases, the electron chemical potential increases resulting in a lower density at which Urca pairs can exist. For intermediate field strength there is an interesting interference between the Landau level distribution and the Fermi distribution. For high enough field strength, the entire electron energy spectrum is eventually confined to single Landau level producing dramatic spikes in the emission spectrum.

Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software "VapoRock" combines the MELTS liquid model (Ghiorso et al. 1995) with gas-species properties from multiple thermochemistry tables (e.g. Chase et al. 1998). VapoRock calculates abundances of 34 gaseous species in equilibrium with magmatic liquid in the system Si-Mg-Fe-Al-Ca-Na-K-Ti-Cr-O at desired temperatures and oxygen fugacities (fO2, or partial pressure of O2). Comparison with experiments shows that pressures and melt-oxide activities (which vary over many orders of magnitude) are reproduced within a factor of ~3, consistent with measurement uncertainties. We also benchmark against a wide selection of igneous rock compositions including bulk silicate Earth, predicting elemental vapor abundances that are comparable (Na, Ca, & Al) or more realistic (K, Si, Mg, Fe, & Ti) than those of the closed-source MAGMA code (with maximum deviations by factors of 30-300 for K). Calculated vapor abundances depend critically on liquid activities, and the MELTS model underpinning VapoRock was specifically calibrated on natural igneous liquids and has been extensively tested & refined over the last 3 decades. In contrast, MAGMA's underlying liquid model assumes ideal mixtures of liquid pseudo-species, which are incapable of capturing the non-ideal compositional interactions that typify the behavior of natural silicate melts. Using VapoRock, we finally explore how relative abundances of SiO and SiO2 provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constraining fO2 in outgassed exoplanetary mantles.

As a dedicated solar radioheliograph, the MingantU SpEctral RadioHeliograph (MUSER) has a maximum baseline of more than 3000 meters and a frequency range of 400 MHz -- 15 GHz. According to the classical radio interferometry theory, the non-coplanar baseline effect (i.e., w-term effect) would be considered and calibrated for such a radio instrument. However, little previous literature made the qualitative or quantitative analyses on w-term effects of solar radioheliograph in-depth. This study proposes a complete quantitative analysis of w-term effects for the MUSER. After a brief introduction of the MUSER, we systematically investigate the baseline variations over a year and analyze the corresponding variations of w-term. We further studied the effects of the w-term in the imaging for the specified extended source, i.e., the Sun. We discussed the possible effects of the w-term, such as image distortion and so on. The simulated results show that the w-term is an essential and unavoidable issue for solar radio imaging with high spatial resolution.

Gravitational wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernova (PISN). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections from merging BHs. We use \MESA\ to evolve single, non-rotating, massive helium cores with a metallicity of $Z = 10^{-5}$ until they either collapse to form a BH or explode as a PISN without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the $\pm 3\sigma$ uncertainty in our high resolution tabulated $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction rate probability distribution function. We extensively test the temporal and mass resolution to resolve the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower resolution models reach convergence. We establish a new lower edge of the upper mass gap as M\textsubscript{lower} $\simeq$\,60$^{+32}_{-14}$\,\Msun\ from the $\pm 3\sigma$ uncertainty in the $^{12}\text{C}(\alpha, \gamma) ^{16}\text{O}$ rate. We explore the effect of a larger 3-$\alpha$ rate on the lower edge of the upper mass gap, finding M\textsubscript{lower} $\simeq$\,69$^{+34}_{-18}$\,\Msun. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.

We report on the fast Nova Sgr 2016 N.4 being surprisingly trapped in a long-lasting and bright plateau (Delta I >= 10 mag above quiescence) six years past the nova eruption. Very few other novae experience a similar occurrence. We carried out an intensive observing campaign collecting daily BVRI photometry and monthly high-resolution optical spectroscopy, and observed the nova in ultraviolet and X-rays with Swift satellite at five distinct epochs. The bolometric luminosity radiated during the plateau is ~4200 Lsun (scaled to the distance of the Galactic Bulge), corresponding to stable nuclear burning on a 0.6 Msun white dwarf. A stable wind is blown off at FWZI~1600 km/s, with episodic reinforcement of a faster FWZI~3400 km/s mass loss, probably oriented along the polar directions. The collision of these winds could power the emission detected in X-rays. The burning shell has an outer radius of ~25 Rsun at which the effective temperature is ~7600 K, values similar to those of a F0 II/Ib bright giant. The Delta m < 1 mag variability displayed during the plateau is best described as chaotic, with the irregular appearance of quasi-periodic oscillations with a periodicity of 15-17 days. A limited amount of dust (~3x10^(-11) Msun) continuously condenses at T(dust)~1200 K in the outflowing wind, radiating L(dust)~52 Lsun.

The goal of this paper is to investigate the flux and spectral index distribution of FR II radio galaxy 4C 14.11. We focused on the distribution of flux and spectral indices over the lobes, as well as in their hot spots. For that purpose, we used publicly available observations of this radio galaxy given at 1450 and 8440 MHz. Particularly, we used Leahy's Atlas of radio-emitting double radio sources, Jodrell Bank Centre for Astrophysics in Manchester, as well as NASA/IPAC Extragalactic Database. We found that the non-thermal (synchrotron) radiation dominates over the areas of the lobes. Distinction between hot spots and rest of the lobes are much smaller in the spectral index than in the flux. We also found that over the inner parts of both lobes, spectral index $\alpha$ is flat in average and significantly higher than 1.2, indicating that 4C 14.11 is old AGN.

Orbital characteristics based on Gaia Early Data Release 3 astrometric parameters are analyzed for ${\sim} 1700$ $r$-process-enhanced (RPE; [Eu/Fe] $> +0.3$) metal-poor stars ([Fe/H] $\leq -0.8$) compiled from the $R$-Process Alliance, the GALactic Archaeology with HERMES (GALAH) DR3 survey, and additional literature sources. We find dynamical clusters of these stars based on their orbital energies and cylindrical actions using the \HDBSCAN ~unsupervised learning algorithm. We identify $36$ Chemo-Dynamically Tagged Groups (CDTGs) containing between $5$ and $22$ members; $17$ CDTGs have at least $10$ member stars. Previously known Milky Way (MW) substructures such as Gaia-Sausage-Enceladus, the Splashed Disk, the Metal-Weak Thick Disk, the Helmi Stream, LMS-1 (Wukong), and Thamnos are re-identified. Associations with MW globular clusters are determined for $7$ CDTGs; no recognized MW dwarf galaxy satellites were associated with any of our CDTGs. Previously identified dynamical groups are also associated with our CDTGs, adding structural determination information and possible new identifications. Carbon-Enhanced Metal-Poor RPE (CEMP-$r$) stars are identified among the targets; we assign these to morphological groups in the Yoon-Beers $A$(C)$_{c}$ vs. [Fe/H] Diagram. Our results confirm previous dynamical analyses that showed RPE stars in CDTGs share common chemical histories, influenced by their birth environments.

It is proposed that GLEAM-X~J162759.5-523504.3, the newly discovered radio transient with an unusually long spin-period (P$_s$ = 1091.1690s), can be identified to be a Radio Magnetar which has a dipolar surface magnetic field of $2.5 \times 10^{16}$~G. It is shown that - a) it is possible to anchor such a strong field at the core-crust boundary of a neutron star, and b) the energy of field dissipation can explain the observed luminosity (radio & X-ray) of this source.

The formyl cation (HCO+) is one of the most abundant ions in molecular clouds and plays a major role in the interstellar chemistry. For this reason, accurate collisional rate coefficients for the rotational excitation of HCO+ and its isotopes due to the most abundant perturbing species in interstellar environments are crucial for non-local thermal equilibrium models and deserve special attention. In this work, we determined the first hyperfine resolved rate coefficients of HC17O+ in collision with H2 (j=0). Indeed, despite no scattering calculations on its collisional parameters have been performed so far, the HC17O+ isotope assumes a prominent role for astrophysical modelling applications. Computations are based on a new four dimensional (4D) potential energy surface, obtained at the CCSD(T)-F12a/aug-cc-pVQZ level of theory. A test on the corresponding cross section values pointed out that, to a good approximation, the influence of the coupling between rotational levels of H2 can be ignored. For this reason, the H2 collider has been treated as a spherical body and an average of the potential based on five orientations of H2 has been employed for scattering calculations. State-to-state rate coefficients resolved for the HC17O+ hyperfine structure for temperature ranging from 5 to 100 K have been computed using recoupling techniques. This study provides the first determination of HC17O+/H2 inelastic rate coefficients directly computed from full quantum close-coupling equations, thus supporting the reliability of future radiative transfer modellings of HC17O+ in interstellar environments.

We present new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey for the two protoclusters, G12.42+0.50 and G19.88-0.53. The 3 mm continuum maps reveal seven cores in each of the two globally contracting protoclusters. These cores satisfy the radius-mass relation and the surface mass density criteria for high-mass star formation. Similar to their natal clumps, the virial analysis of the cores suggests that they are undergoing gravitational collapse ($\rm \alpha_{vir} << 2$). The clump to core scale fragmentation is investigated and the derived core masses and separations are found to be consistent with thermal Jeans fragmentation. We detect large-scale filamentary structures with velocity gradients and multiple outflows in both regions. Dendrogram analysis of the H$^{13}$CO$^{+}$ map identifies several branch and leaf structures with sizes $\sim$ 0.1 and 0.03 pc, respectively. The supersonic gas motion displayed by the branch structures is in agreement with the Larson power-law indicating that the gas kinematics at this spatial scale is driven by turbulence. The transition to transonic/subsonic gas motion is seen to occur at spatial scales of $\sim$0.1 pc indicating the dissipation of turbulence. In agreement with this, the leaf structures reveal gas motions that deviate from the slope of Larson's law. From the large-scale converging filaments to the collapsing cores, the gas dynamics in G12.42+0.50 and G19.88-0.53 show scale-dependent dominance of turbulence and gravity and the combination of these two driving mechanisms needs to be invoked to explain massive star formation in the protoclusters.

We present a convenient tool (ChaSES) which allows to search for extended structures in Chandra X-ray Observatory Advanced CCD Imaging Spectrometer (ACIS) images. The tool relies on DBSCAN clustering algorithm to detect regions with overdensity of photons compared to the background. Here we describe the design and functionality of the tool which we make publicly available on GitHub. We also provide online extensive examples of its applications to the real data.

Gamma-ray bursts (GRBs) are classified into long and short populations (i.e., LGRBs and SGRBs) based on the observed bimodal distribution of duration $T_{90}$. Multimessenger observations indicated that most SGRBs and LGRBs should be powered by ultrarelativistic jets launched from black hole (BH) hyperaccretion in compact object mergers and massive collapsars, respectively. However, the duration criterion sometimes cannot correctly reflect the physical origin of a particular GRB. In the collapsar scenario, a GRB can be observed when the jet breaks out from the envelope and circumstellar medium successfully. The observed GRB duration reflects only the time that the engine operates after the jet breaks out. This work studies the propagation of jets driven by the neutrino annihilation or Blandford-Znajek mechanism in massive collapsars. The signatures of the progenitors for producing LGRBs, SGRBs, and failed GRBs in the collapsar scenario are exhibited. The competition between the mass supply onto the BH hyperaccretion and jet propagation into the envelope are definitely dependent on the density profiles of the collapsars. We show that duration and isotropic energy $E_{\rm{\gamma,iso}}$ of GRBs can help constrain the density profiles of collapsars. Finally, we propose that a collapsar-origin SGRB, GRB 200826A, might originate from a neutrino-annihilation-dominated jet launched by a $\sim 10~M_\odot$ collapsar whose progenitor's envelope has been stripped.

To more thoroughly study the effects of radiative forces on the orbits of small, astronomical bodies, we introduce the Yarkovsky effect into REBOUNDx, an extensional library for the N-body integrator REBOUND. Two different versions of the Yarkovsky effect (the 'Full Version' and the 'Simple Version') are available for use, depending on the needs of the user. We provide demonstrations for both versions of the effect and compare their computational efficiency with another previously implemented radiative force. In addition, we show how this effect can be used in tandem with other features in REBOUNDx by simulating the orbits of asteroids during the asymptotic giant branch (AGB) phase of a 2 $M_{\odot}$ star. This effect is made freely available for use with the latest release of REBOUNDx.

The 2017 Event Horizon Telescope (EHT) observations of M87* detected a ring-shaped feature $\sim40\mu$as in diameter, consistent with the event horizon scale of a black hole of the expected mass. The thickness of this ring, however, proved difficult to measure, despite being an important parameter for constraining the observational appearance. In the first paper of this series we asked whether the width of the ring was sensitive to the choice of likelihood function used to compare observed closure phases and closure amplitudes to model predictions. In this paper we investigate whether the ring width is robust to changes in the model itself. We construct a more realistic geometric model with two new features: an adjustable radial falloff in brightness, and a secondary "photon ring" component in addition to the primary annulus. This thin, secondary ring is predicted by gravitational lensing for any black hole with an optically thin accretion flow. Analyzing the data using the new model, we find that the primary annulus remains narrow (fractional width $\leq0.25$) even with the added model freedom. This provides further evidence in favor of a narrow ring for the true sky appearance of M87*, a surprising feature that, if confirmed, would demand theoretical explanation. Comparing the Bayesian evidence for models with and without a secondary ring, we find no evidence for the presence of a lensed photon ring in the 2017 observations. However, the techniques we introduce may prove useful for future observations with a larger and more sensitive array.

We present the analysis of the absorption troughs of six outflows observed in quasar SDSS J1130+0411 ($z \approx 3.98$) with radial velocities ranging from $-2400$ to $-15,400$ km s$^{-1}$. These spectra were taken with the Very Large Telescope/Ultraviolet and Visual Echelle Spectrograph over the rest frame wavelength range of $1135-1890$ \r{A}. In the main outflow system ($v \approx -3200$ km s$^{-1}$), we identify Fe II and several Fe II* absorption troughs as well as Si II and Si II* troughs, which we use to determine the electron number density $\log n_e = 2.6_{-0.7}^{+0.8}$ cm$^{-3}$. Using the column densities of these and other ions, we determine a photoionization solution with hydrogen column density $\log N_H = 21.44_{-0.33}^{+0.24}$ cm$^{-2}$ and ionization parameter $\log U_H = -1.75_{-0.45}^{+0.28}$. From these values we derive the distance $R = 16_{-11}^{+23}$ kpc, the average mass flow rate $\dot{M} = 4100_{-2400}^{+6600}$ $M_{\odot}$ yr$^{-1}$, and the kinetic luminosity $\log \dot{E}_k = 46.13_{-0.37}^{+0.41}$ erg s$^{-1}$. This $\dot{E}_k$ is $1.4_{-0.8}^{+2.2}$% of the quasar's Eddington luminosity, and therefore contributes significantly to AGN feedback.

To investigate the effects of environment in the quenching phase, we study the empirical relations for green valley (GV) galaxies between overdensity and other physical properties (i.e., effective radius $r_{\rm e}$, S\'{e}rsic indices $n$, and specific star formation rate sSFR). Based on five 3D-{\it HST}/CANDELS fields, we construct a large sample of 2126 massive ($M_{\star} > 10^{10} M_{\sun}$) GV galaxies at $0.5<z<2.5$ and split it into the higher overdensity quarter and the lower overdensity quarter. The results shows that GV galaxies in denser environment have higher $n$ values and lower sSFR at $0.5< z <1$, while there is no discernible distinction at $1 < z < 2.5$. No significant enlarging or shrinking is found for GV galaxies in different environments within the same redshift bin. It suggests that a dense environment would promote the growth of bulge and suppress star formation activity of GV galaxies at $0.5< z <1.5$, but would not affect the galaxy size. We also study the dependence of the fraction of three populations (Blue Cloud, Green Valley, and Red Sequence) on both environments and $M_{\star}$. At a given $M_{\star}$, blue cloud fraction goes down with increasing environment density, while red sequence fraction is opposite. For the most massive GV galaxies, a sharp drop appears in the denser environment. Coupled with the mass dependence of three fractions in different redshift bins, our result implies that stellar mass and environments jointly promote the quenching process. Such dual effect is also confirmed by re-calculating the new effective GV fraction as the number of GV galaxies over the number of non-quiescent galaxies.

Among the ~2157 unassociated sources in the third data release (DR3) of the fourth Fermi catalog, ~1200 were observed with the Neil Gehrels Swift Observatory pointed instruments. These observations yielded 238 high S/N X-ray sources within the 95% Fermi uncertainty regions. Recently, Kerby et al. employed neural networks to find blazar candidates among these 238 X-ray counterparts to the 4FGL unassociated sources and found 112 likely blazar counterpart sources. A complete sample of blazars, along with their sub-classification, is a necessary step to help understand the puzzle of the blazar sequence and for the overall completeness of the gamma-ray emitting blazar class in the Fermi catalog. We employed a multi-perceptron neural network classifier to identify FSRQs and BL Lacs among these 112 blazar candidates using the gamma-ray, X-ray, UV/optical, and IR properties. This classifier provided probability estimates for each source to be associated with one or the other category, such that P_fsrq represents the probability for a source to be associated with the FSRQ subclass. Using this approach, 4 FSRQs and 50 BL Lacs are classified as such with >99% confidence, while the remaining 58 blazars could not be unambiguously classified as either BL Lac or FSRQ.

We present a detailed analysis of a reflecting intensity perturbation in a large coronal loop that appeared as sloshing oscillation and lasted for at least one and a half periods. The perturbation is initiated by a microflare at one footpoint of the loop, propagates along the loop and is eventually reflected at the remote footpoint where significant brightenings are observed in all the AIA extreme-ultraviolet (EUV) channels. This unique observation provides us with the opportunity to better understand not only the thermal properties and damping mechanisms of the sloshing oscillation, but also the energy transfer at the remote footpoint. Based on differential emission measures (DEM) analysis and the technique of coronal seismology, we find that 1) the calculated local sound speed is consistent with the observed propagation speed of the perturbation during the oscillation, which is suggestive of a slow magnetoacoustic wave; 2) thermal conduction is the major damping mechanism of the wave but additional damping mechanism such as anomalous enhancement of compressive viscosity or wave leakage is also required to account for the rapid decay of the observed waves; 3) the wave produced a nanoflare at the remote footpoint, with a peak thermal energy of $\thicksim10^{24}-10^{25}$ erg. This work provides a consistent picture of the magnetoacoustic wave propagation and reflection in a coronal loop, and reports the first solid evidence of a wave-induced nanoflare. The results reveal new clues for further simulation studies and may help solving the coronal heating problem.

The correlation between host star iron abundance and the exoplanet occurrence rate is well-established and arrived at in several studies. Similar correlations may be present for the most abundant elements, such as carbon and oxygen, which also control the dust chemistry of the protoplanetary disk. In this paper, using a large number of stars in the Kepler field observed by the LAMOST survey, it has been possible to estimate the planet occurrence rate with respect to the host star carbon abundance. Carbon abundances are derived using synthetic spectra fit of the CH G-band region in the LAMOST spectra. The carbon abundance trend with metallicity is consistent with the previous studies and follows the Galactic chemical evolution (GCE). Similar to [Fe/H], we find that the [C/H] values are higher among giant planet hosts. The trend between [C/Fe] and [Fe/H] in planet hosts and single stars is similar; however, there is a preference for giant planets around host stars with a sub-solar [C/Fe] ratio and higher [Fe/H]. Higher metallicity and sub-solar [C/Fe] values are found among younger stars as a result of GCE. Hence, based on the current sample, it is difficult to interpret the results as a consequence of GCE or due to planet formation.

We present new wide field survey strategies for Chilean Large Aperture Telescopes (LAT) measuring the Cosmic Microwave Background (CMB), which we call Sinusoidal Modulated High Cadence Survey Strategies. The strategies were developed during the process of optimizing LAT measurements for the CMB-S4, Simons Observatory, and CCAT-prime collaborations. Observing more than $f_{sky} \sim 0.5$, the telescope consistently achieves high observation efficiency, even with Sun-avoidance enabled. Classical azimuthal scan survey strategies observing fields of equal size suffer from problems of observation depth non-uniformity relative to declination and lack of crosslinking. The new survey strategies described here significantly improve both uniformity and crosslinking while also enabling higher cadence observations for time-domain astrophysics. Uniformity and crosslinking are improved by modulation of azimuthal angular velocity and sinusoidal elevation nods, respectively. In particular, there is nearly uniform observation depth and crosslinking is improved from total lack of crosslinking near -40 degree declination to clearing the strictest thresholds for crosslinking across the entire field. The simulated strategies are compared to the strategies used for the Atacama Cosmology Telescope and previously studied Simons Observatory survey strategies.

A sample of 279 massive red spirals was selected optically by Guo et al. (2020), among which 166 galaxies have been observed by the ALFALFA survey. In this work, we observe HI content of the rest 113 massive red spiral galaxies using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). 75 of the 113 galaxies have HI detection with a signal-to-noise ratio (S/N) greater than 4.7. Compared with the red spirals in the same sample that have been observed by the ALFALFA survey, galaxies observed by FAST have on average a higher S/N, and reach to a lower HI mass. To investigate why many red spirals contain a significant amount of HI mass, we check color profiles of the massive red spirals using images observed by the DESI Legacy Imaging Surveys. We find that galaxies with HI detection have bluer outer disks than the galaxies without HI detection, for both ALFALFA and FAST samples. For galaxies with HI detection, there exists a clear correlation between galaxy HI mass and g-r color at outer radius: galaxies with higher HI masses have bluer outer disks. The results indicate that optically selected massive red spirals are not fully quenched, and the HI gas observed in many of the galaxies may exist in their outer blue disks.

It is known that the blazar jet emissions are dominated by non-thermal radiation while the accretion disk jets are normally dominated by thermal emission. In this work, our aim is to study the connection between the two types of emission by investigating the correlation between the blazar emission line intensity property, which embodies the nature of accretion disk, and the $\gamma$-ray flux property, which is the representative of jet emission. We compiled a sample of 656 blazars with available emission line equivalent widths ($EW$), the GeV $\gamma$-ray flux, and the SED information from the literature. In this work, we found 55 previous BCUs are now identified as FSRQs, and found 52 Changing-look blazars based on their $EW$ and 45 of them are newly confirmed. These Changing-look blazars have a larger accretion ratio (${\dot M}/{\dot M}_{\rm Edd}$) than BL Lac objects. In addition, we suggest that the lower synchrotron peak blazars (LSPs) could be the source of Changing-look blazars because 90.7\% of the Changing-look blazars in this work are confirmed as LSPs. An anti-correlation between $EW$ and continuum intensity, the so-called global Baldwin effect (BEff) has been confirmed. We suggest the steeper global BEff observed for blazar than for radio-quiet active galactic nuclei (RQ-AGNs) is caused by the inverse Compton scattering of broad-emission-line photons. This interpretation is further supported by the positive correlation between the emission line $EW$ and intrinsic inverse Compton luminosity.

Buckminsterfullerene, C60 , is the largest molecule observed to date in interstellar and circumstellar environments. The mechanism of formation of this molecule is actively debated. Despite targeted searches in primitive carbonaceous chondrites, no unambiguous detection of C60 in a meteorite has been reported to date. Here we report the first firm detection of fullerenes, from C30 to at least C100 , in the Almahata Sitta (AhS) polymict ureilite meteorite. This detection was achieved using highly sensitive laser desorption laser ionization mass spectrometry. Fullerenes have been unambiguously detected in seven clasts of AhS ureilites. Molecular family analysis shows that fullerenes are from a different reservoir compared to the polycyclic aromatic hydrocarbons detected in the same samples. The fullerene family correlates best with carbon clusters, some of which may have been formed by the destruction of solid carbon phases by the impacting laser. We show that the detected fullerenes are not formed in this way. We suggest that fullerenes are an intrinsic component of a specific carbon phase that has yet to be identified. The nondetection of fullerenes in the Murchison and Allende bulk samples, while using the same experimental conditions, suggests that this phase is absent or less abundant in these primitive chondrites. The former case would support the formation of fullerenes by shock-wave processing of carbonaceous phases in the ureilite parent body. However, there are no experimental data to support this scenario. This leaves open the possibility that fullerenes are an interstellar heritage and a messenger of interstellar processes.

Solar radio zebras are used in the determination of the plasma density and magnetic field in solar flare plasmas. Analyzing observed zebra stripes and assuming their generation by the double-plasma resonance (DPR) instability, high values of the gyro-harmonic number are found. In some cases, they exceed one hundred, in disagreement with the DPR growth rates computed up to now, which decrease with increasing gyro-harmonic number. We address the question of how the zebras with high values of the gyro-harmonic numbers $s$ are generated. For this purpose, we compute growth rates of the DPR instability in a very broad range of $s$, considering a loss-cone $\kappa$-distribution of superthermal electrons and varying the loss-cone angle, electron energies, and background plasma temperature. We numerically calculated dispersion relations and growth rates of the upper-hybrid waves and found that the growth rates increase with increasing gyro-harmonic numbers if the loss-cone angles are $\sim80^\circ$. The highest growth rates for these loss-cone angles are obtained for the velocity $v_\kappa = 0.15\,c$. The growth rates as function of the gyro-harmonic number still show well distinct peaks, which correspond to zebra-stripe frequencies. The contrast of the growth rate peaks to surrounding growth rate levels increases as the $\kappa$ index increases and the background temperature decreases. Zebras with high values of $s$ can be generated in regions where loss-cone distributions of superthermal electrons with large loss-cone angles ($\sim80^\circ$) are present. Furthermore, owing to the high values of $s$, the magnetic field is relatively weak and has a small spatial gradient in such regions.

The spherical harmonic transform is a powerful tool in the analysis of spherical data sets, such as the cosmic microwave background data. In this work, we present a new scheme for the spherical harmonic transforms that supports both CPU and GPU computations, which is especially efficient on a large number of sky maps. By comparing our implementation with the standard \textsl{Libsharp-HEALPix} code, which is one of the top-speed programs, we demonstrate 4--10 times speedup for the CPU implementation, and up to 60 times speedup when a state-of-the-art GPU is employed. This new scheme is available at the following GitHub repository: \href{https://github.com/liuhao-cn/fastSHT}{\textsl{fastSHT}}.

The Minkowski tensors (MTs) can be used to probe anisotropic signals in a field, and are well suited for measuring the redshift space distortion (RSD) signal in large scale structure catalogs. We consider how the linear RSD signal can be extracted from a field without resorting to the plane parallel approximation. A spherically redshift space distorted field is both anisotropic and inhomogeneous. We derive expressions for the two point correlation functions that elucidate the inhomogeneity, and then explain how the breakdown of homogeneity impacts the volume and ensemble averages of the tensor Minkowski functionals. We construct the ensemble average of these quantities in curvilinear coordinates and show that the ensemble and volume averages can be approximately equated, but this depends on our choice of definition of the volume average of a tensor and the radial distance between the observer and field. We then extract the tensor Minkowski functionals from spherically redshift space distorted, Gaussian random fields and gravitationally evolved dark matter density fields at $z=0$ to test if we can successfully measure the Kaiser RSD signal. For the dark matter field we find a significant, $\sim 10\%$ anomalous signal in the MT component parallel to the line of sight that is present even on large scales $R_{\rm G} \gtrsim 15 \, {\rm Mpc}$, in addition to the Kaiser effect. This is due to the line of sight component of the MT being significantly contaminated by the Finger of God effect, which can be approximately modelled by an additional damping term in the cumulants.

Strong gravitational lensing effect is a powerful tool to probe cosmological models and gravity theories. Recently, the time-delay cosmography from strong lensing and the stellar kinematics of the deflector, which encode the Hubble constant and the post-Newtonian parameter via two distance ratios reflecting the lensing mass and dynamical mass respectively, have been proposed to investigate these two parameters simultaneously. Among strong lensing systems with different sources, strongly lensed fast radio bursts (FRBs) have been proposed as precision probes of the universe since the time delay $\sim$ 10 days between images could be measured extremely precisely because of their short duration of a few milliseconds. In this work, we investigate the ability of strongly lensed FRBs on simultaneously estimating these two parameters via simulations. Take the expected FRB detection rate of upcoming facilities and lensing probability into consideration, it is likely to accumulate 10 lensed FRBs in several years and we find that $H_0$ could be determined to a $\sim1.5\%$ precision and $\gamma_{\rm PPN}$ could be constrained to a $\sim8.7\%$ precision simultaneously from them. These simultaneous estimations will be helpful for properly reflecting the possible correlation between these two fundamental parameters.

The defining trait of magnetars, the most strongly magnetized neutron stars (NSs), is their transient activity in the X/$\gamma$-bands. In particular, many of them undergo phases of enhanced emission, the so-called outbursts, during which the luminosity rises by a factor $\sim$10$-$1000 in a few hours to then decay over months/years. Outbursts often exhibit a thermal spectrum, associated with the appearance of hotter regions on the surface of the star, which subsequently change in shape and cool down. Here we simulate the unfolding of a sudden, localized heat injection in the external crust of a NS with a 3D magneto-thermal evolution code, finding that this can reproduce the main features of magnetar outbursts. A full 3D treatment allows us to study for the first time the inherently asymmetric hot-spots which appear on the surface of the star as the result of the injection and to follow the evolution of their temperature and shape. We investigate the effects produced by different physical conditions in the heated region, highlighting in particular how the geometry of the magnetic field plays a key role in determining the properties of the event.

Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magnetohydrodynamic (MHD) simulations of SNRs expanding into a weak, regular magnetic field. We present the intrinsic FDFs of the shocked region of SNRs for different viewing angles. We find that the FDFs are generally Faraday complex, which implies that conventional rotation measure study is not sufficient to obtain the information on the magnetic fields in the shocked region and Faraday tomography is necessary. We also show that the FDF allows to derive the physical-depth distribution of polarization intensity when the line of sight is parallel to the initial magnetic field orientation. Furthermore, we demonstrate that the location of contact discontinuity can be identified from the radial profile of the width of the FDF with the accuracy of 0.1-0.2 pc.

Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate for the Large Magellanic Cloud (LMC). The binary models are calculated through Roche-lobe overflow and mass transfer, until the initially more massive star exhausts helium in its core. We distinguish models of WR and helium stars based on the estimated stellar wind optical depth. We use these models to build a synthetic WR population, assuming constant star formation. Our models can reproduce the WR population of the LMC to significant detail, including the number and luminosity functions of the main WR subtypes. We find that for binary fractions of 100% (50%), all LMC WR stars below $10^6\,L_{\odot}$ ($10^{5.7}\,L_{\odot}$) are stripped binary mass donors. We also identify several insightful mismatches. With a single star fraction of 50\%, our models produce too many yellow supergiants, calling either for a larger initial binary fraction, or for enhanced mass-loss near the Humphreys-Davidson limit. Our models predict more long-period WR binaries than observed, arguably due to an observational bias towards short periods. Our models also underpredict the shortest-period WR binaries, which may have implications for understanding the progenitors of double black hole mergers. The fraction of binary produced WR stars may be larger than often assumed, and outline the risk to mis-calibrate stellar physics when only single star models are used to reproduce the observed WR stars.

Planet 9 is an hypothetical object in the outer Solar system, which is as yet undiscovered. It has been speculated that it may be a terrestrial planet or gas/ice giant, or perhaps even a primordial black hole (or dark matter condensate). State-of-the-art models indicate that the semimajor axis of Planet 9 is $\sim 400$ AU. If the location of Planet 9 were to be confirmed and pinpointed in the future, this object constitutes an interesting target for a future space mission to characterize it further. In this paper, we describe various mission architectures for reaching Planet 9 based on a combination of chemical propulsion and flyby maneuvers, as well as more advanced options (with a $\sim 100$ kg spacecraft payload) such as nuclear thermal propulsion (NTP) and laser sails. The ensuing mission duration for solid chemical propellant ranges from 45 years to 75 years, depending on the distance from the Sun for the Solar Oberth maneuver. NTP can achieve flight times of about 40 years with only a Jupiter Oberth maneuver whereas, in contrast, laser sails might engender timescales as little as 7 years. We conclude that Planet 9 is close to the transition point where chemical propulsion approaches its performance limits, and alternative advanced propulsion systems (e.g., NTP and laser sails) apparently become more attractive.

$\require{mediawiki-texvc}$ Supermassive black holes at the centre of active galactic nuclei (AGN) produce relativistic jets that can affect the star formation characteristics of the AGN hosts. Observations in the ultraviolet (UV) band can provide an excellent view of the effect of AGN jets on star formation. Here, we present a census of star formation properties in the Northern Star-forming Region (NSR) that spans about 20 kpc of the large radio source Centaurus A hosted by the giant elliptical galaxy NGC 5128. In this region, we identified 352 UV sources associated with Cen A using new observations at an angular resolution of $<$1.5 arcseconds observed with the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat. These observations were carried out in one far-ultraviolet (FUV; $\lambda_{\text{mean}}$ = 1481 $\AA$) and three near-ultraviolet (NUV; with $\lambda_{\text{mean}}$ of 2196 $\AA$, 2447 $\AA$, and 2792 $\AA$, respectively) bands. The star-forming sources identified in UV tend to lie in the direction of the jet of Cen A, thereby suggesting jet triggering of star formation. Separating the NSR into Outer and Inner regions, we found the stars in the Inner region to have a relatively younger age than the Outer region, suggesting that the two regions may have different star formation histories. We also provide the UVIT source catalogue in the NSR.

The Mingantu Spectral Radioheliograph (MUSER), a new generation of solar dedicated radio imaging-spectroscopic telescope, has realized high-time, high-angular, and high-frequency resolution imaging of the sun over an ultra-broadband frequency range. Each pair of MUSER antennas measures the complex visibility in the aperture plane for each integration time and frequency channel. The corresponding radio image for each integration time and frequency channel is then obtained by inverse Fourier transformation of the visibility data. In general, the phase of the complex visibility is severely corrupted by instrumental and propagation effects. Therefore, robust calibration procedures are vital in order to obtain high-fidelity radio images. While there are many calibration techniques available -- e.g., using redundant baselines, observing standard cosmic sources, or fitting the solar disk -- to correct the visibility data for the above-mentioned phase errors, MUSER is configured with non-redundant baselines and the solar disk structure cannot always be exploited. Therefore it is desirable to develop alternative calibration methods in addition to these available techniques whenever appropriate for MUSER to obtain reliable radio images. In the case that a point-like calibration source containing an unknown position error, we have for the first time derived a mathematical model to describe the problem and proposed an optimization method to calibrate this unknown error by studying the offset of the positions of radio images over a certain period of the time interval. Simulation experiments and actual observational data analyses indicate that this method is valid and feasible. For MUSER's practical data the calibrated position errors are within the spatial angular resolution of the instrument. This calibration method can also be used in other situations for radio aperture synthesis observations.

KASCADE and its extension array of KASCADE-Grande were devoted to measure individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The experiment has substantially contributed to investigate the energy spectrum and mass composition of cosmic rays in the transition region from galactic to extragalactic origin of cosmic rays as well as to quantify the characteristics of hadronic interaction models in the air shower development through validity tests using the multi-detector information from KASCADE-Grande. Although the data accumulation was completed in 2013, data analysis is still continuing. Recently, we investigated the reliability of the new hadronic interactions models of the SIBYLL version 2.3d with the data from KASCADE-Grande. The evolution of the muon content of high energy air showers in the atmosphere is studied as well. In this talk, recent results from KASCADE-Grande and the update of the KASCADE Cosmic Ray Data Centre (KCDC) will be discussed.

We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts $1.0 \lesssim z \lesssim 1.7$ ($z_\mathrm{median} = 1.4$) and Sunyaev Zel'dovich (SZ) detection significance $\xi > 6.0$ from the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F606W and F814W images and used additional observations from HST/WFC3 in F110W and VLT/FORS2 in $U_\mathrm{HIGH}$ to preferentially select background galaxies at $z\gtrsim 1.8$, achieving a high purity. We combined recent redshift estimates from the CANDELS/3D-HST and HUDF fields to infer an improved estimate of the source redshift distribution. We measured weak lensing masses by fitting the tangential reduced shear profiles with spherical Navarro-Frenk-White (NFW) models. We obtained the largest lensing mass in our sample for the cluster SPT-CLJ2040$-$4451, thereby confirming earlier results that suggest a high lensing mass of this cluster compared to X-ray and SZ mass measurements. Combining our weak lensing mass constraints with results obtained by previous studies for lower redshift clusters, we extended the calibration of the scaling relation between the unbiased SZ detection significance $\zeta$ and the cluster mass for the SPT-SZ survey out to higher redshifts. We found that the mass scale inferred from our highest redshift bin ($1.2 < z < 1.7$) is consistent with an extrapolation of constraints derived from lower redshifts, albeit with large statistical uncertainties. Thus, our results show a similar tendency as found in previous studies, where the cluster mass scale derived from the weak lensing data is lower than the mass scale expected in a Planck $\nu\Lambda$CDM (i.e. $\nu$ $\Lambda$ Cold Dark Matter) cosmology given the SPT-SZ cluster number counts.

[Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive corotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, deplete their corotation regions, loose the support of the corotation torque and the trapping efficiency then becomes uncertain. We study what it takes to trap such gap-opening planets at the inner disk rim. We perform 2D locally isothermal and non-isothermal hydrodynamic simulations of planet migration. A viscosity transition is introduced in the disk to (i) create a density drop and (ii) mimic the viscosity increase as the planet migrates from a dead zone towards a region with active magneto-rotational instability (MRI). We choose TOI-216b as a Neptune-like upper-limit test case but we also explore different planetary masses, both on fixed and evolving orbits. For planet-to-star mass ratios $q\simeq(4$-$8)\times10^{-5}$, the density drop at the disk rim becomes reshaped due to gap opening and is often replaced with a small density bump centered on the planet's corotation. Trapping is possible only if the bump retains enough gas mass and if the corotation region becomes azimuthally asymmetric, with an island of librating streamlines that accumulate a gas overdensity ahead of the planet. The overdensity exerts a positive torque that can counteract the negative torque of spiral arms. In our model, efficient trapping depends on the $\alpha$-viscosity and its contrast across the viscosity transition. In order to trap TOI-216b, $\alpha_{\mathrm{DZ}}=10^{-3}$ in the dead zone requires $\alpha_{\mathrm{MRI}}\gtrsim5\times10^{-2}$ in the MRI-active zone. If $\alpha_{\mathrm{DZ}}=5\times10^{-4}$, $\alpha_{\mathrm{MRI}}\gtrsim7.5\times10^{-2}$ is needed.

The goal of this thesis was to develop a procedure to optimize the mirror suppression, which reduces the dynamic range of the Fast Fourier Transform spectrometers (FFTS), and implement this procedure in the Field Programmable Gate Array (FPGA) of the FFTS. This is achieved by applying a calibration on the complex amplitude spectrum. Modern multi-beam heterodyne receivers are equipped with a large number of independent signal chains, which leads to a complex system that consumes space as well as power and is expensive. To reduce the complexity of the receiver array the intermediate frequency (IF) band is 4-8 GHz, which requires the FFTS to sample this band directly. The input bandwidth of 4 GHz of the current generation of FFTS cannot be sampled with a single analog-digital converter (ADC). To reach that bandwidth, multiple ADCs sample the same signal at different points in time. The method is called \"time interleaved sampling\". The characteristics of the ADCs differ from each other, which leads to mirror signals in the spectrum. These mirror signals can reduce the dynamic range of the spectrometers. Therefore the mismatches must be actively corrected. The current generation of FFTS use a frequency independent calibration, which optimizes the mirror suppression at one fixed frequency. This it is not optimal as some of the mismatches are frequency dependent and their impact increases with frequency. The procedure developed during this thesis work involves three steps. First the mismatches are measured over the frequency band. Then filter coefficients are calculated from these mismatches using a model of a time interleaved ADC. These filter coefficients are then applied to the complex amplitude spectra of the ADCs. The newly implemented calibration improves the mirror suppression up to 20 dB compared to the frequency independent calibration.

Many astronomical systems display complex, rapid variability across the electromagnetic spectrum. Studying this variability is the key to understanding the underlying processes at play. However, a combination of limited or patchy telescope availability, viewing constraints, and the unpredictable and faint nature of many sources mean that obtaining data well-suited to this task can be tricky, especially when it comes to simultaneous multiwavelength observations. As a result, researchers can often find themselves tuning observational parameters in real-time by trial-and-error, or may realise later that their observation was not long enough to achieve their goals. Here, we present a code to mitigate some of these issues and aid planning of multiwavelength coordinated observations. The code takes a temporal and Fourier model of a system (i.e. Power Spectra, Coherence, and Lags) as input, and returns a simulated multiwavelength observation, including effects of noise, telescope parameters, and finite sampling. The goals of this are: (i) To simulate a potential observation (to inform decisions about an observation's feasibility); (ii) To investigate how different Fourier models affect the variability of the system (e.g. how altering the frequency-dependent coherence or Fourier lags between bands can affect data products like cross-correlation functions); and (iii) To simulate existing data and investigate its trustworthiness. We outline the methodology behind the code, show how a variety of parameters (e.g. noise sources, observation length, and telescope choice) can affect data, and present examples of the software in action.

We describe advances on a method designed to derive accurate parameters of M dwarfs. Our analysis consists in comparing high-resolution infrared spectra acquired with the near-infrared spectro-polarimeter SPIRou to synthetic spectra computed from MARCS model atmospheres, in order to derive the effective temperature ($T_{\rm eff}$), surface gravity ($\rm \log{g}$), metallicity ([M/H]) and alpha-enhancement ($\rm [\alpha/Fe]$) of 44 M dwarfs monitored within the SPIRou Legacy Survey (SLS). Relying on 12 of these stars, we calibrated our method by refining our selection of well modelled stellar lines, and adjusted the line list parameters to improve the fit when necessary. Our retrieved $T_{\rm eff}$, $\rm \log{g}$ and [M/H] are in good agreement with literature values, with dispersions of the order of 50 K in $T_{\rm eff}$ and 0.1 dex in $\rm \log{g}$ and [M/H]. We report that fitting $\rm [\alpha/Fe]$ has an impact on the derivation of the other stellar parameters, motivating us to extend our fitting procedure to this additional parameter. We find that our retrieved $\rm [\alpha/Fe]$ are compatible with those expected from empirical relations derived in other studies.

We report the discovery and characterization of two small transiting planets orbiting the bright M3.0V star TOI-1468 (LSPM J0106+1913), whose transit signals were detected in the photometric time series in three sectors of the TESS mission. We confirm the e planetary nature of both of them using precise radial velocity measurements from the CARMENES and MAROON-X spectrographs, and supplement them with ground-based transit photometry. A joint analysis of all these data reveals that the shorter-period planet, TOI-1468 b ($P_{\rm b}$ = 1.88 d), has a planetary mass of $M_{\rm b} = 3.21\pm0.24$ $M_{\oplus}$ and a radius of $R_{\rm b} =1.280^{+0.038}_{-0.039} R_{\oplus}$, resulting in a density of $\rho_{\rm b} = 8.39^{+ 1.05}_{- 0.92}$ g cm$^{-3}$, which is consistent with a mostly rocky composition. For the outer planet, TOI-1468 c ($P_{\rm c} = 15.53$ d), we derive a mass of $M_{\rm c} = 6.64^{+ 0.67}_{- 0.68}$ $M_{\oplus}$, a radius of $R_{\rm c} = 2.06\pm0.04\,R_{\oplus}$, and a bulk density of $\rho_{c} = 2.00^{+ 0.21}_{- 0.19}$ g cm$^{-3}$, which corresponds to a rocky core composition with a H/He gas envelope. These planets are located on opposite sides of the radius valley, making our system an interesting discovery as there are only a handful of other systems with the same properties. This discovery can further help determine a more precise location of the radius valley for small planets around M dwarfs and, therefore, shed more light on planet formation and evolution scenarios.

The petrologic and geochemical diversity of meteorites is a function of the bulk composition of their parent bodies, but also the result of how and when internal differentiation took place. Here we focus on this second aspect considering the two principal parameters involved: size and accretion time of the body. We discuss the interplay of the various time scales related to heating, cooling and drainage of silicate liquids. Based on two phase flow modelling in 1-D spherical geometry, we show that drainage time is proportional to two independent parameters: $\mu_m/R^2$, the ratio of the matrix viscosity to the square of the body radius and $\mu_f/a^2$, the ratio of the liquid viscosity to the square of the matrix grain size. We review the dependence of these properties on temperature, thermal history and degree of melting, demonstrating that they vary by several orders of magnitude during thermal evolution. These variations call into question the results of two phase flow modelling of small body differentiation that assume constant properties.For example, the idea that liquid migration was efficient enough to remove $^{26}$Al heat sources from the interior of bodies and dampen their melting (e.g. Moskovitz and Gaidos, 2011; Neumann et al., 2012) relies on percolation rates of silicate liquids overestimated by six to eight orders of magnitude. In bodies accreted during the first few million years of solar-system history, we conclude that drainage cannot prevent the occurrence of a global magma ocean. These conditions seem ideal to explain the generation of the parent-bodies of iron meteorites. A map of the different evolutionary scenarios of small bodies as a function of size and accretion time is proposed.

The missing baryons in the universe are assumed to be hidden in the whole space as a warm-hot intergalactic medium (WHIM). Finding them is one of the important subjects in modern cosmology. In this paper, we point out that the very high energy electron/positron rays may light up the WHIM due to the anomalous bremsstrahlung according to the improved Bethe-Heitler formula. The resulting excess of the extragalactic background light (EBL) can be observed by the direct measurement method. A possible explanation on the difference between the direct and indirect measurements of EBL is also proposed. Thus, we open a new window to probe the WHIM properties via the EBL.

We use optical spectroscopy from VLT/MUSE, as well as photometric observations from HST and VLT/HAWK-I, to study the morpho-kinematics of 17 low mass MACS J0416.1-2403 cluster galaxies at R200 and 5 field galaxies with a redshift of z~0.4. By measuring fluxes of emission lines, we have recovered the star formation rates, gas-phase metallicities and the spatially resolved gas kinematics. We have analysed the structure and morphology of the galaxies from the photometric data, performing a multi-component decomposition into a bulge and a disk. The spatially resolved gas velocity fields were modelled using a 3D approach, which allowed us to retrieve their intrinsic gas kinematics, including the maximum rotation velocity and velocity dispersion. The cluster and field population can be classified as star forming main-sequence galaxies, with only a sub-sample of 4 quenched systems. The lowest mass systems deviate from the predictions of the fundamental metallicity relation, showing higher metallicities, whereas the higher mass ones are in accordance with the model predictions. This might hint at the cut-off of pristine gas inflow and/or the removal of the hot halo gas as the mechanisms driving these offsets. Our morpho-kinematic analysis reveals a sub-sample of dwarfs with maximum velocities vmax <50 km/s, which depart from the Tully-Fisher relation. This might indicate that their interstellar medium is affected by ram pressure stripping. However, ~30% of the cluster galaxies have rotation-dominated gas disks and follow the Tully-Fisher relation. In the stellar mass-S0.5 plane, the systems follow a tight sequence, with only a sub-population of 5 galaxies strongly departing from this relation. Both the morphology and kinematics of the outlier galaxies hint at a combination of pre-processing and cluster-specific interactions affecting their stellar and gas disks.

We present an analysis of an unusual pattern of water vapor emission from the $\sim$2 Myr-old low-mass binary system VV CrA, as observed in infrared spectra obtained with VLT-CRIRES, VLT-VISIR, and Spitzer-IRS. Each component of the binary shows emission from water vapor in both the L ($\sim3\,\mu$m) and N ($\sim 12\,\mu$m) bands. The N-band and Spitzer spectra are similar to those previously observed from young stars with disks, and are consistent with emission from an extended protoplanetary disk. Conversely, the CRIRES L-band data of VV CrA A show an unusual spectrum, which requires the presence of a water reservoir with high temperature ($T\gtrsim1500$ K), column density ($N_\mathrm{H2O}\sim 3\times10^{20}\ \mathrm{cm}^{-2}$), and turbulent broadening ($v\sim 10$ km s$^{-1}$), but very small emitting area ($A\lesssim0.005$ AU$^2$). Similarity with previously observed water emission from V1331 Cyg (Doppmann et al. 2011) and SVS 13 (Carr et al. 2004) suggests that the presence of such a reservoir may be linked to evolutionary state, perhaps related to the presence of high accretion rates or winds. While the inner disk may harbor such a reservoir, simple Keplerian models do not match well with emitting line shapes, and alternative velocity fields must be considered. We also present a new idea, that the unusual emission could arise in a circumplanetary disk, embedded within the larger VV CrA A protoplanetary disk. Additional data are likely required to determine the true physical origin of this unusual spectral pattern.

Massive binaries play an important role in fields ranging from gravitational wave astronomy to stellar evolution. We provide several lines of evidence that classical OBe stars in the Small Magellanic Cloud (SMC) obtain their rapid rotation from mass and angular momentum transfer in massive binaries, which predicts that the subsequent supernovae should often eject OBe stars into the field. We find that (1) OBe stars have a higher field frequency than OB stars; (2) our cumulative distribution function (CDF) of stellar distances from O stars shows that OBe stars are indeed much more isolated than ordinary OB stars of corresponding spectral types; (3) the CDFs of OBe stars approach that of high-mass X-ray binaries (HMXBs), which are confirmed post-supernova objects; and (4) Oe stars are as isolated from clusters as Be stars, implying that their final masses are relatively independent of their initial masses, consistent with major mass transfer. Lastly, we also find that the spatial distribution of supergiant OBe stars differs from that of classical OBe stars, consistent with the different mechanism responsible for their emission-line spectra.

We propose a novel method to measure the $E_G$ statistic from clustering alone. The $E_G$ statistic provides an elegant way of testing the consistency of General Relativity by comparing the geometry of the Universe, probed through gravitational lensing, with the motion of galaxies in that geometry. Current $E_G$ estimators combine galaxy clustering with gravitational lensing, measured either from cosmic shear or from CMB lensing. In this paper, we construct a novel estimator for $E_G$, using only clustering information obtained from two tracers of the large-scale structure: intensity mapping and galaxy clustering. In this estimator, both the velocity of galaxies and gravitational lensing are measured through their impact on clustering. We show that with this estimator, we can suppress the contaminations that affect other $E_G$ estimators and consequently test the validity of General Relativity robustly. We forecast that with the coming generation of surveys like HIRAX and $\textit{Euclid}$, we will measure $E_G$ with a precision of up to 7% (3.9% for the more futuristic SKA2).

The detection of Lyman-$\alpha$ emitting galaxies (LAEs) puts severe constraints on the reionization history. In this Paper we derive the properties of very high-$z$ LAEs predicted in the only two reionization scenarios shown in a previous Paper to be consistent with current data on 15 independent evolving global (or averaged) cosmic properties regarding luminous objects and the intergalactic medium and the optical depth to electron scattering of ionized hydrogen to CMB photons: one with a monotonic behavior, which is completed by $z=6$, as commonly considered, and another one with a non-monotonic behavior with two full ionization events at $z=6$ and $z=10$. We find that the \Lya luminosity functions of very high-$z$ LAEs are very distinct in those two scenarios. Thus, comparing these predictions to the observations that will soon be available thanks to new instruments such as the James Webb Space Telescope, it should be possible to unveil the right reionization scenario. In the meantime, we can compare the predicted redshift distribution and UV (or Lyman-$\alpha$) luminosities of very high-$z$ LAEs to those of the few objects already observed at $z>7.5$. By doing that we find that such data are in tension with the single reionization scenario, while they are fully compatible with the double reionization scenario.

TESS (Transiting Exoplanet Survey Satellite) was launched in 2018 with the purpose of observing bright stars in the solar neighbourhood to search for transiting exoplanets. After the completion of the two year nominal mission, TESS has provided 2\,minute cadence photometry of over 200\,000 stars. This large collection of light curves opens the possibility to study the statistical and temporal properties of this ensemble of stars. Most of the currently available data pipelines are designed to work on single sector at a time. We present a new TESS data pipeline called {\tt Taranga}, with the purpose of merging multi-sector light curves, whilst performing a period search for all the observed stars, and stores the statistical results in a database. {\tt Taranga} pipeline has three components which 1) processes the PDCSAP fluxes of each sector and creates merged PDCSAP light curve, 2) performs a similar operation on the SAP fluxes, and 3) generates the periodograms of the merged SAP and PDCSAP light curves while performing peak identification. For all the 232\,122 stars observed in short cadence in the nominal TESS mission, we provide the merged PDCSAP and SAP light-curves along with their periodograms. We provide a database that has the statistics of all the results produced from {\tt Taranga} of these stars.

We present the clump populations detected in 18 lensed galaxies at redshifts 1 to 8 within the lensing cluster field SMACS0723. The recent JWST Early Release Observations of this poorly known region of the sky, have revealed numerous point-like sources within and surrounding their host galaxies, undetected in the shallower HST images. We use JWST multiband photometry and the lensing model of this galaxy cluster to estimate the intrinsic sizes and magnitudes of these stellar clumps. We derive optical restframe effective radii from $<$10 to 100s pc and masses ranging from $\sim10^5$ to $10^8$ M$_{\odot}$, overlapping with massive star clusters in the local universe. The ages range from 1 to 100s Myr. We compare the crossing time to the age of the clumps and determine that between 30 and 50 \% of the detected clumps are consistent with being gravitationally bound. The lack of Gyr old clumps suggest that the dissolution time scales are shorter than 1 Gyr. We see a significant increase in the luminosity (mass) surface density of the clumps with redshift. Clumps in galaxies at the reionisation era have stellar densities higher than massive clusters in the local universe. We zoom-in into single galaxies at redshift $<6$ and find for two galaxies, the Sparkler and the Firework, that their star clusters/clumps show distinctive color distributions and location within and surrounding their host galaxy that are compatible with recent (100-500 Myr) merger events. Our study, conducted on a small sample of galaxies, shows the potential of JWST observations for understanding the conditions under which star clusters form in rapidly evolving galaxies.

The distribution of dark matter halo masses can be accurately predicted in the $\Lambda$CDM cosmology. The presence of a single massive halo or galaxy at a particular redshift, assuming some baryon and stellar fraction for the latter, can therefore be used to test the underlying cosmological model. A number of recent measurements of very large galaxy stellar masses at high redshift ($z > 8$) motivate an investigation into whether any of these objects are in tension with $\Lambda$CDM. We use extreme value statistics to generate confidence regions in the mass-redshift plane for the most extreme mass haloes and galaxies. Tests against numerical models show no tension, neither in their dark matter halo masses nor their galaxy stellar masses. However, we find tentative $> 3\sigma$ tension with recent observational determinations of galaxy masses at high redshift from both HST & JWST, despite using conservative estimates for the stellar fraction ($f_{\star} \sim 1$). Either these galaxies are in tension with $\Lambda$CDM, or there are unaccounted for uncertainties in their stellar mass or redshift estimates.

Current and future B-mode polarization data are the most powerful observables to constrain gravitational waves from the Early Universe. We set conservative constraints on tensor modes when relaxing the inflationary consistency condition $n_t=-r/8$ between the tensor tilt $n_t$ and the tensor-to-scalar ratio r. By adding a power-law spectrum of tensor perturbations to $\Lambda$CDM, and parameterizing this tensor contribution by two independent primordial tensor-to-scalar ratios $(r_1,r_2)$ at $k_1 = 0.005$ Mpc$^{-1}$ and $k_2 = 0.02$ Mpc$^{-1}$, Planck and BICEP/Keck Array 2018 data (BK18) lead to constraints $r_{0.005} < 0.030$ and $r_{0.02} < 0.098$ at 95% CL. The corresponding upper bound $r_{0.01} < 0.039$ is by a factor 2 tighter than the one obtained with Planck 2018 and the older BK15 data. We then study the perspectives for future CMB experiments that will measure both the reionization bump and recombination peak of the B-mode polarization angular power spectrum, such as LiteBIRD. We test the robustness of the results to the choice of the scales for $(r_1,r_2)$ in these future perspectives. Whereas distinguishing $n_t=-r/8$ from exact scale invariance is impossible as expected, we show how radical, theoretically motivated departures from $n_t=-r/8$, which are consistent with the current data, could be distinguished with LiteBIRD. Moreover, LiteBIRD will be able to shrink the allowed parameter space area in the $(r_{0.005},r_{0.02})$ plane to less than one hundredth of the currently allowed area by Planck 2018 and BK18.

Binary neutron star mergers, which can lead to less massive black holes relative to other known astrophysical channels, have the potential to probe modifications to general relativity that arise at smaller curvature scales compared to more massive compact object binaries. As a representative example of this, here we study binary neutron star mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet gravity using evolutions of the full, non-perturbative evolution equations. We find that the impact on the inspiral is small, even at large values of the modified gravity coupling (as expected, as neutron stars do not have scalar charge in this theory). However, post-merger there can be strong scalar effects, including radiation. When a black hole forms, it develops scalar charge, impacting the ringdown gravitational wave signal. In cases where a longer-lived remnant star persists post-merger, we find that the oscillations of the star source levels of scalar radiation similar to the black hole formation cases. In remnant stars, we further find that at coupling values comparable to the maximum value for which black hole solutions of the same mass exist, there is significant nonlinear enhancement in the scalar field, which if sufficiently large lead to a breakdown in the evolution, seemingly due to loss of hyperbolicity of the underlying equations.

If mirror matter exists but cannot comprise all of the dark matter in the universe, we can expect that the additional dark matter component may only interact with the other sectors gravitationally. A natural candidate of a gravitationally interacting component is a primordial black hole (PBH). We show constraints on PBH with the mirror dark matter. Particularly, the initial PBH mass is estimated to be $10^{17} \ {\rm g} \lesssim M_{\rm PBH} \lesssim 10^{25} \ {\rm g}$, if the dark matter comprises mirror baryons and PBHs.

Asymptotic normalization coefficients (ANC) determine the overall normalization of cross sections of peripheral radiative capture reactions. In the present paper, we treat the ANC $C$ for the virtual decay $^{16}$O$(0^+; 6.05$ MeV)$\to \alpha+^{12}$C(g.s.), the known values of which are characterized by a large spread $(0.29-1.65)\times 10^3$ fm$^{-1/2}$. The ANC $C$ is found by analytic continuation in the energy of the $\alpha^{12}$C $s$-wave scattering amplitude, known from the phase-shift analysis of experimental data, to the pole corresponding to the $^{16}$O bound state and lying in the unphysical region of negative energies. To determine $C$, two different methods of analytic continuation are used. In the first method, the scattering data are approximated by the sum of polynomials in energy in the physical region and then extrapolated to the pole. The best way of extrapolation is chosen on the basis of the exactly solvable model. Within the second approach, the ANC $C$ is found by solving the Schr\"odinger equation for the two-body $\alpha^{12}$C potential, the parameters of which are selected from the requirement of the best description of the phase-shift analysis data at a fixed experimental binding energy of $^{16}$O$(0^+; 6.05$ MeV) in the $\alpha+^{12}$C channel. The values of the ANC $C$ obtained within these two methods lie in the interval (886--1139) fm$^{-1/2}$.

We investigate a static, cylindrically symmetric cosmic string on the brane without a perturbative approximation. We find there could be a (large) enhancement of the (effective) string tension when the energy density at the center of the string is (much) larger than twice the brane tension. We also point out a new way to evade the cosmic string problem when the energy density at the center of the string approaches twice the brane tension. These findings could have experimental and theoretical implications for searching for cosmic strings on the brane, in particular for cosmic strings generated after inflation (such as D-term inflation) on the brane.

To understand turbulent convection at very high Rayleigh numbers typical of natural phenomena, computational studies in slender cells are an option if the needed resources have to be optimized within available limits. However, the accompanying horizontal confinement affects some properties of the flow. Here, we explore the characteristics of turbulent fluctuations in the velocity and temperature fields in a cylindrical convection cell of aspect ratio 0.1 by varying the Prandtl number $Pr$ between 0.1 and 200 at a fixed Rayleigh number $Ra = 3 \times 10^{10}$, and find that the fluctuations weaken with increasing $Pr$, quantitatively as in aspect ratio 25. The probability density function (PDF) of temperature fluctuations in the bulk region of the slender cell remains mostly Gaussian, but increasing departures occur as $Pr$ increases beyond unity. We assess the intermittency of the velocity field by computing the PDFs of velocity derivatives and of the kinetic energy dissipation rate, and find increasing intermittency as $Pr$ decreases. In the bulk region of convection, a common result applicable to the slender cell, large aspect ratio cells, as well as in 2D convection, is that the turbulent Prandtl number decreases as $Pr^{-1/3}$.

We compute the expected response of detector arms of gravitational wave observatories to polymerized gravitational waves. The mathematical and theoretical features of these waves were discussed in our previous work. In the present manuscript, we find both perturbative analytical, and full nonperturbative numerical solutions to the equations of motion of the detector arms using the method of geodesic deviations. These results show the modifications to both frequency and amplitude of the signal measured by the detector. Furthermore, we study the detectability of these signals in LISA by analyzing the modes in the frequency space.

We show that viable electroweak baryogenesis can be realized without a first-order phase transition if plasma is heated inhomogeneously by non-gravitational structure formation in some particle species. Yukawa interactions can mediate relatively long-range attractive forces in the early Universe. This creates an instability and leads to growth of structure in some species even during the radiation dominated era. At temperatures below the elecroweak scale, the collapsing and annihilating halos can heat up plasma in fireballs that expand and create the out-of-equilibrium high-temperature environment suitable for generating the baryon asymmetry. The plasma temperature at the time of baryogenesis can be as low as 10 MeV, making it consistent with both standard and low-reheat cosmologies.

We derive the equations of motion for an $N$-body system in the Einstein-Cartan gravity theory at the first post-Newtonian order by exploiting the Weyssenhoff fluid as the spin model. Our approach consists in performing the point-particle limit of the continuous description of the gravitational source. The final equations provide a hint for the validity of the effacing principle at 1PN level in Einstein-Cartan model. The analogies with the general relativistic dynamics involving the macroscopic angular momentum are also discussed.

A long-lived scalar field ($\Phi$) which couples weakly to the right-handed (RH) neutrinos ($N_{Ri}$), generates small RH neutrino masses ($M_i$) in Low-Scale-Leptogenesis (LSL) mechanisms, despite having a large vacuum expectation value $v_\Phi$. In this case, the correlation shared by the $M_i$s and the duration of the non-standard cosmic history driven by the $\Phi$ provides an excellent opportunity to study LSL signatures on primordial gravitational waves (GWs). We find it engaging, specifically for the gravitational waves that originate due to the inflationary blue-tilted tensor power spectrum and propagate through the non-standard cosmic epoch. Depending on $M_i$, broadly, the scenario has two significant consequences. First, if LSL is at play, GWs with a sizeable blue tilt do not contradict the Big-Bang-Nucleosynthesis (BBN) bound even for the post-inflationary models with very high-scale reheating. Second, it opens up a possibility to probe LSLs via a low-frequency and a complementary high-frequency measurement of GW-spectral shapes which are typically double-peaked. For a case study, we consider the recent results on GWs from the Pulsar-Timing-Arrays (PTAs) as a `measurement' at the low frequencies and forecast the signatures of LSL mechanisms at the higher frequencies.

In this paper we discuss on the phenomenological footprints of theories where the gravitational effects are due not only to spacetime curvature, but also to nonmetricity. These theories are characterized by gauge invariance. Due to their simplicity, here we focus in theories with vectorial nonmetricity. We make special emphasis in gradient nonmetricity theories which are based in Weyl integrable geometry (WIG) spaces. While arbitrary and vectorial nonmetricities may have played a role in the quantum epoch, gradient nonmetricity can be important for the description of gravitational phenomena in our classical world instead. This would entail that gauge symmetry may be an actual symmetry of our past, present and future universe, without conflict with the standard model of particles (SMP). We show that, in a gauge invariant world modeled by WIG spacetime, the vacuum energy density is a dynamical quantity, so that the cosmological constant problem (CCP) may be avoided. Besides, due to gauge invariance, and to the fact that photons and radiation do not interact with nonmetricity, the accelerated pace of cosmic expansion can be explained without the need for the dark energy. We also discuss on the ``many-worlds'' interpretation of the resulting gauge invariant framework, where general relativity (GR) is just a specific gauge of the theory. The unavoidable discrepancy between the present value of the Hubble parameter computed on the GR basis and its value according to the gauge invariant theory, may explain the Hubble tension issue. It will be shown also that, due to gauge freedom, inflation is not required in order to explain the flatness, horizon and relict particles abundance problems within the present framework.

We consider a supergravity model that has the modular $A_4$ symmetry and discuss the interplay between the neutrino mixing and inflation. The model contains right-handed neutrinos that have the Majorana masses and additional Yukawa couplings to the waterfall field. In the model an active neutrino is massless and we find that only the inverted hierarchy is allowed and the Majorana phase is predicted to be around $\pm (120\text{--}180)^\circ$ from the observed neutrino mixing data. In the early universe, one of right-handed sneutrinos plays the role of the inflaton field. Focusing on the subcritical regime of the hybrid inflation that is consistent with the cosmic microwave background data, we analyze the dynamics of the scalar sector and derive an upper bound $\mathcal{O}(10^{10})~{\rm GeV}$ on the scale of the Majorana mass.

For many classes of astronomical and astrophysical binary systems, long observational records of their radial velocity $V$, which is their directly observable quantity, are available. For exoplanets close to their parent stars, they cover several full orbital revolutions, while for wide binaries like, e.g., the Proxima/$\alpha$ Centauri AB system, only relatively short orbital arcs are sampled by existing radial velocity measurements. Here, the changes $\Delta V$ induced on a binary's radial velocity by some long-range modified models of gravity are analytically calculated. In particular, extra-potentials proportional to $r^{-N},\,N=2,\,3$ and $r^2$ are considered; the Cosmological Constant $\Lambda$ belongs to the latter group. Both the net shift per orbit and the instantaneous one are explicitly calculated for each model. The Cosmological Constant induces a shift in the radial velocity of the Proxima/$\alpha$ Centauri AB binary as little as $\left|\Delta V\right|\lesssim 10^{-7}\,\mathrm{m\,s}^{-1}$, while the present-day accuracy in measuring its radial velocity is $\sigma_V\simeq 30\,\mathrm{m\,s}^{-1}$. The calculational scheme presented here is quite general, and can be straightforwardly extended to any other modified gravity.

We argue that the Big Bang can be understood as a type of mirror. We show how reflecting boundary conditions for spinors and higher spin fields are fixed by local Lorentz and gauge symmetry, and how a temporal mirror (like the Bang) differs from a spatial mirror (like the AdS boundary), providing a possible explanation for the observed pattern of left- and right-handed fermions. By regarding the Standard Model as the limit of a minimal left-right symmetric theory, we obtain a new, cosmological solution of the strong $CP$ problem, without an axion.

We study the direct radiation excited by oscillating axion (or axion-like particle) dark matter in a homogenous magnetic field and its detection scheme. We concretely derive the analytical expression of the axion-induced radiated power for a cylindrical uniform magnetic field. In the long wave limit, the radiation power is proportional to the square of the B-field volume and the axion mass $m_a$, whereas it oscillate as approaching the short wave limit and the peak powers are proportional to the side area of the cylindrical magnetic field and $m_a^{-2}$. The maximum power locates at mass $m_a\sim\frac{3\pi}{4R}$ for fixed radius $R$. Based on this characteristic of the power, we discuss a scheme to detect the axions in the mass range $1-10^4$\,neV, where four detectors of different bandwidths surround the B-field. The expected sensitivity for $m_a\lesssim1\,\mu$eV under typical-parameter values can far exceed the existing constraints.

The relation between the fourth-order symmetry energy $E_{\rm{sym,4}}(\rho_0)$ of nuclear matter at saturation density $\rho_0$ and its counterpart $a_{\rm{sym,4}}(A)$ of finite nuclei in a semi-empirical nuclear mass formula is revisited by considering the high-order isospin-dependent surface tension contribution to the latter. We derive the full expression of $a_{\rm{sym,4}}(A)$ which includes explicitly the high-order isospin-dependent surface tension effects, and find that the value of $E_{\rm{sym,4}}(\rho_0)$ cannot be extracted from the measured $a_{\rm{sym,4}}(A)$ before the high-order surface tension is well constrained. Our results imply that a large $a_{\rm{sym,4}}(A)$ value of several MeVs obtained from analyzing nuclear masses can nicely agree with the empirical constraint of $E_{\rm{sym,4}}(\rho_0)\lesssim 2$ MeV from mean-field models and does not necessarily lead to a large $E_{\rm{sym,4}}(\rho_0)$ value of $\sim 20$ MeV obtained previously without considering the high-order surface tension. Furthermore, we also give the expression for the sixth-order symmetry energy $a_{\rm{sym,6}}(A)$ of finite nuclei, which involves more nuclear matter bulk parameters and the higher-order isospin-dependent surface tension.

For over a decade, the SpacePy project has contributed open-source solutions for the production and analysis of heliophysics data and simulation results. Here we review SpacePy's design principles and functionality, before examining recent advances and the future of SpacePy in the broader scientific Python ecosystem. We conclude with some of the work that has used SpacePy.