Papers for Wednesday, Jun 16 2021

The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 microns and KS band at 2.1 microns using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 earth radii or 4500 km. We also detected a medium-depth transit which can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5 to 1.0 microns, depending on the composition of the debris. This leads to a high mass loss rate of about 3e8 kg/s that is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 earth masses predicted from hydrodynamical models.

Globular clusters (GCs) are bright objects that span a wide range of galactocentric distances, and are thus probes of the structure of dark matter (DM) haloes. In this work, we explore whether the projected radial profiles of GCs can be used to infer the structural properties of their host DM haloes. We use the simulated GC populations in a sample of 166 central galaxies from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. We find that more massive galaxies host stellar and GC populations with shallower density profiles that are more radially extended. In addition, the metal-poor GC subpopulations tend to have shallower and more extended profiles than the metal-rich subsamples, which we relate to the preferentially accreted origin of the metal-poor GCs. We find strong correlations between the slopes and effective radii of the radial profiles of the GC populations and the structural properties of the DM haloes, such as their power-law slopes, scale radii, and concentration parameters. Accounting for a dependence on the galaxy stellar mass decreases the scatter of the two-dimensional relations. This suggests that the projected number counts of GCs, combined with their galaxy mass, trace the density profile of the DM halo of their host galaxy. When applied to extragalactic GC systems, we recover the scale radii and the extent of the DM haloes of a sample of ETGs with uncertainties smaller than $0.2~\rm dex$. Thus, extragalactic GC systems provide a novel avenue to explore the structure of DM haloes beyond the Local Group.

Characterizing the fundamental parameters of stars from observations is crucial for studying the stars themselves, their planets, and the galaxy as a whole. Stellar evolution theory predicting the properties of stars as a function of stellar age and mass enables translating observables into physical stellar parameters by fitting the observed data to synthetic isochrones. However, the complexity of overlapping evolutionary tracks often makes this task numerically challenging, and with a precision that can be highly variable, depending on the area of the parameter space the observation lies in. This work presents StelNet, a Deep Neural Network trained on stellar evolutionary tracks that quickly and accurately predicts mass and age from absolute luminosity and effective temperature for stars with close to solar metallicity. The underlying model makes no assumption on the evolutionary stage and includes the pre-main sequence phase. We use bootstrapping and train many models to quantify the uncertainty of the model. To break the model's intrinsic degeneracy resulting from overlapping evolutionary paths, we also built a hierarchical model that retrieves realistic posterior probability distributions of the stellar mass and age. We further test and train StelNet using a sample of stars with well-determined masses and ages from the literature.

Accurate stellar properties are essential for precise stellar astrophysics and exoplanetary science. In the M-dwarf regime, much effort has gone into defining empirical relations that can use readily accessible observables to assess physical stellar properties. Often, these relations for the quantity of interest are cast as a non-linear function of available data; however, in Bayesian modeling the reverse is needed. In this article, we introduce a new Bayesian framework to self-consistently and simultaneously apply multiple empirical calibrations to fully characterize the mass, luminosity, radius, and effective temperature of a field age M-dwarf. This framework includes a new M-dwarf mass-radius relation with a scatter of 3.1% at fixed mass. We further introduce the M-dwarf Ultraviolet Spectroscopic Sample (MUSS), and apply our methodology to provide consistent stellar parameters for these nearby low-mass stars, selected as having available spectroscopic data in the ultraviolet. These targets are of interest largely as either exoplanet hosts or benchmarks in multi-wavelength stellar activity. We use the field MUSS stars to define a low-mass main sequence in the solar neighborhood through Gaussian Process (GP) regression. These results enable us to empirically measure a feature in the GP derivative at $\mathcal{M}_{b} = 0.337 \pm ^{0.013}_{0.026}$ $M_{\odot}$ that indicates where the M-dwarf Ultraviolet Spectroscopic Sample transitions from fully to partly convective interiors.

We report the discovery of a new low-mass X-ray binary near the center of the unassociated Fermi GeV gamma-ray source 4FGL J0540.0-7552. The source shows the persistent presence of an optical accretion disk and exhibits extreme X-ray and optical variability. It also has an X-ray spectrum well-fit by a hard power law with a Gamma = 1.8 and a high ratio of X-ray to gamma-ray flux. Together, these properties are consistent with the classification of the binary as a transitional millisecond pulsar (tMSP) in the sub-luminous disk state. Uniquely among the candidate tMSPs, 4FGL J0540.0-7552 shows consistent optical, X-ray, and gamma-ray evidence for having undergone a state change, becoming substantially brighter in the optical and X-rays and fainter in GeV gamma-rays sometime in mid-2013. In its current sub-luminous disk state, and like one other candidate tMSP in the Galactic field, 4FGL J0540.0-7552 appears to always be in an X-ray "flare mode", indicating that this could be common phenomenology for tMSPs.

We identify hierarchical structures in the Vela OB2 complex and the cluster pair Collinder 135 and UBC 7 with Gaia EDR3 using the neural network machine learning algorithm StarGO. Five second-level substructures are disentangled in Vela OB2, which are referred to as Huluwa 1 (Gamma Velorum), Huluwa 2, Huluwa 3, Huluwa 4 and Huluwa 5. Huluwa 1-5 may have originated through sequential star formation. The clusters Huluwa 1-3 are the older generation, with ages of 10-20 Myr, and generated stellar feedback that caused turbulence which fostered the formation of the younger-generation clusters, Huluwa 4-5 (10 Myr). The cluster pair Collinder 135 and UBC 7 was likely formed from the same molecular cloud, are coeval (40 Myr). The 3D morphology of Huluwa 1-5 resembles a shell-like structure, right along the rim of the Vela IRAS shell. A supernova explosion located inside the Vela IRAS shell quenched star formation in the younger generation Huluwa 4-5 and rapidly expelled the remaining gas from the clusters. This resulted in mass stratification across the shell, with more low-mass stars located along the inner rim of the shell, and more massive stars in the outer region of the shell. Mass segregation is observed only in the lowest-mass cluster Huluwa 5. Significant expansion is observed in Vela OB2 with a 1D expansion rate of $(6.9-7.9)\times10^{-2}$~km $\rm s^{-1} pc^{-1}$. Expansion in the cluster pair is moderate. The velocity dispersions suggest that both Vela OB2 and the cluster pair are supervirial and are undergoing disruption. $N$-body simulations predict that Huluwa 1-5 in Vela OB2 and the cluster pair will continue to expand in the future 100 Myr and eventually dissolve.are supervirial and are undergoing disruption. $N$-body simulations predict that Huluwa 1-5 in Vela OB2 and the cluster pair will continue to expand in the future 100 Myr and eventually dissolve.

The classic classification scheme for Active Galactic Nuclei (AGNs) was recently challenged by the discovery of the so-called changing-state (changing-look) AGNs (CSAGNs). The physical mechanism behind this phenomenon is still a matter of open debate and the samples are too small and of serendipitous nature to provide robust answers. In order to tackle this problem, we need to design methods that are able to detect AGN right in the act of changing-state. Here we present an anomaly detection (AD) technique designed to identify AGN light curves with anomalous behaviors in massive datasets. The main aim of this technique is to identify CSAGN at different stages of the transition, but it can also be used for more general purposes, such as cleaning massive datasets for AGN variability analyses. We used light curves from the Zwicky Transient Facility data release 5 (ZTF DR5), containing a sample of 230,451 AGNs of different classes. The ZTF DR5 light curves were modeled with a Variational Recurrent Autoencoder (VRAE) architecture, that allowed us to obtain a set of attributes from the VRAE latent space that describes the general behaviour of our sample. These attributes were then used as features for an Isolation Forest (IF) algorithm, that is an anomaly detector for a ''one class'' kind of problem. We used the VRAE reconstruction errors and the IF anomaly score to select a sample of 8,809 anomalies. These anomalies are dominated by bogus candidates, but we were able to identify 75 promising CSAGN candidates.

The LIGO-Virgo collaboration reported in their third run the coalescence event GW190814 involving a 2.6 $M_{\odot}$ object with a 23 $M_{\odot}$ black hole. In this letter we study the conditions under which Thorne-$\dot{\textrm{Z}}$ytkow objects (T$\dot{\textrm{Z}}$Os) can be connected to that type of events. We evaluate first the rate of appearance of T$\dot{\textrm{Z}}$Os in the local Universe. Under the assumption that T$\dot{\textrm{Z}}$Os eventually become low mass gap black holes we evaluate how those black holes end up in binaries with other stellar mass black holes and compare to the reported rate for GW190814-type of events (1-23 $\textrm{Gpc}^{-3} \textrm{yr}^{-1}$). We find that T$\dot{\textrm{Z}}$Os in dense stellar clusters can not explain the LIGO-Virgo rate without a T$\dot{\textrm{Z}}$Os population in the field providing a dominant contribution. We also find that T$\dot{\textrm{Z}}$Os formed within hierarchical triple systems in the field with the third more distant star being the progenitor of a stellar mass black hole, may be able to give a rate comparable to that of GW190814-type events. In that case, future observations should discover mergers between stellar mass and low mass gap black holes, with the lower mass spanning the entire low mass gap range.

The largest fluctuation in the CMB sky is the CMB dipole, which is believed to be caused by the motion of our observation frame with respect to the CMB rest frame. This motion accounts for the known motion of the Solar System barycentre with a best-fit amplitude of 369 km/s, in the direction ($\ell= 264^\circ$, $b=48^\circ$) in galactic coordinates. Along with the CMB dipole signal, this motion also causes an inevitable signature of statistical anisotropy in the higher multipoles due to the modulation and aberration of the CMB temperature and polarization fields. This leads to a correlation between adjacent CMB multipoles causing a non-zero value of the off-diagonal terms in the covariance matrix which can be captured in terms of the dipolar spectra of the bipolar spherical harmonics (BipoSH). In our work, we jointly infer the CMB power spectrum and the BipoSH spectrum in a Bayesian framework using the $\textit{Planck}$-2018 $\texttt{SMICA}$ temperature map. We detect amplitude and direction of the local motion consistent with the canonical value $v=369$ km/s inferred from CMB dipole with a statistical significance of $4.54\sigma$, $4.97\sigma$ and $5.23\sigma$ respectively from the masked temperature map with the available sky fraction $40.1\%$, $59.1\%$, and $72.2\%$, confirming the common origin of both the signals. The Bayes factor in favor of the canonical value is between $7$ to $8$ depending on the choice of mask. But it strongly disagrees with the value inferred from quasar distribution from the Wide-field Infrared Survey Explorer data set with a value of the Bayes factor about $10^{-11}$.

We present a catalogue of white dwarf candidates selected from Gaia early data release three (EDR3). We applied several selection criteria in absolute magnitude, colour, and Gaia quality flags to remove objects with unreliable measurements while preserving most stars compatible with the white dwarf locus in the Hertzsprung-Russell diagram. We then used a sample of over 30 000 spectroscopically confirmed white dwarfs and contaminants from the Sloan Digital Sky Survey (SDSS) to map the distribution of these objects in the Gaia absolute magnitude-colour space. Finally, we adopt the same method presented in our previous Gaia DR2 work to calculate a probability of being a white dwarf (Pwd) for $\simeq$1.3 million sources which passed our quality selection. The Pwd values can be used to select a sample of $\simeq$359 000 high-confidence white dwarf candidates in the magnitude range 8< G <21. We calculated stellar parameters (effective temperature, surface gravity, and mass) for all these stars by fitting Gaia astrometry and photometry with synthetic models. We estimate an upper limit of 93 per cent for the overall completeness of our catalogue for white dwarfs with G $\leq$20 mag and effective temperature (Teff)>7000K, at high Galactic latitudes (|b|>20{\deg}). Alongside the main catalogue we include a reduced-proper-motion extension containing $\simeq$10 200 white dwarf candidates with unreliable parallax measurements which could, however be identified on the basis of their proper motion. We also performed a cross-match of our catalogues with SDSS DR16 spectroscopy and provide spectral classification based on visual inspection for all resulting matches.

A strong super-Alfv\'{e}nic drift of energetic particles (or cosmic rays, CRs) in a magnetized plasma can amplify the magnetic field significantly through non-resonant streaming instability (NRSI). While the traditional analysis is done for an ion current, here we use kinetic particle-in-cell simulations to study how the NRSI behaves when it is driven by electrons or by a mixture of electrons and positrons. In particular, we characterize growth rate, spectrum, and helicity of the unstable modes, as well the level of magnetic field at saturation. Our results are potentially relevant for several space/astrophysical environments (e.g, electron strahl in the solar wind, at oblique non-relativistic shocks, around pulsar wind nebulae) and also in laboratory experiments.

(Abridged:) It has been hypothesized that the location of Herbig Ae/Be stars (HAeBes) within the empirical relation between the inner disk radius (r$_{in}$), inferred from K-band interferometry, and the stellar luminosity (L$_*$), is related to the presence of the innermost gas, the disk-to-star accretion mechanism, the dust disk properties inferred from the spectral energy distributions (SEDs), or a combination of these effects. This work aims to test whether the previously proposed hypotheses do, in fact, serve as a general explanation for the distribution of HAeBes in the size-luminosity diagram. GRAVITY/VLTI spectro-interferometric observations at $\sim $2.2 $ \mu$m have been obtained for five HBes representing two extreme cases concerning the presence of innermost gas and accretion modes. V590 Mon, PDS 281, and HD 94509 show no excess in the near-ultraviolet, Balmer region of the spectra ($\Delta$D$_B$), indicative of a negligible amount of inner gas and disk-to-star accretion, whereas DG Cir and HD 141926 show such strong $\Delta$D$_B$ values that cannot be reproduced from magnetospheric accretion, but probably come from the alternative boundary layer mechanism. Additional data for these and all HAeBes resolved through K-band interferometry have been compiled from the literature and updated using Gaia EDR3 distances, almost doubling previous samples used to analyze the size-luminosity relation. We find no general trend linking the presence of gas inside the dust destruction radius or the accretion mechanism with the location of HAeBes in the size-luminosity diagram. Underlying trends are present and must be taken into account when interpreting the size-luminosity correlation. Still, it is argued that the size-luminosity correlation is most likely to be physically relevant in spite of the previous statistical warning concerning dependencies on distance.

Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies.

In this paper, we revisit the governing equations for linear magnetohydrodynamic (MHD) waves and instabilities existing within a magnetized, plane-parallel, self-gravitating slab. Our approach allows for fully non-uniformly magnetized slabs, which deviate from isothermal conditions, such that the well-known Alfv\'en and slow continuous spectra enter the description. We generalize modern MHD textbook treatments, by showing how self-gravity enters the MHD wave equation, beyond the frequently adopted Cowling approximation. This clarifies how Jeans' instability generalizes from hydro to magnetohydrodynamic conditions without assuming the usual Jeans' swindle approach. Our main contribution lies in reformulating the completely general governing wave equations in a number of mathematically equivalent forms, ranging from a coupled Sturm-Liouville formulation, to a Hamiltonian formulation linked to coupled harmonic oscillators, up to a convenient matrix differential form. The latter allows us to derive analytically the eigenfunctions of a magnetized, self-gravitating thin slab. In addition, as an example we give the exact closed form dispersion relations for the hydrodynamical p- and Jeans-unstable modes, with the latter demonstrating how the Cowling approximation modifies due to a proper treatment of self-gravity. The various reformulations of the MHD wave equation open up new avenues for future MHD spectral studies of instabilities as relevant for cosmic filament formation, which can e.g. use modern formal solution strategies tailored to solve coupled Sturm-Liouville or harmonic oscillator problems.

The search for signs of life through the detection of exoplanet atmosphere biosignature gases is gaining momentum. Yet, only a handful of rocky exoplanet atmospheres are suitable for observation with planned next-generation telescopes. To broaden prospects, we describe the possibilities for an aerial, liquid water cloud-based biosphere in the atmospheres of sub Neptune-sized temperate exoplanets, those receiving Earth-like irradiation from their host stars. One such planet is known (K2-18b) and other candidates are being followed up. Sub Neptunes are common and easier to study observationally than rocky exoplanets because of their larger sizes, lower densities, and extended atmospheres or envelopes. Yet, sub Neptunes lack any solid surface as we know it, so it is worthwhile considering whether their atmospheres can support an aerial biosphere. We review, synthesize, and build upon existing research. Passive microbial-like life particles must persist aloft in a region with liquid water clouds for long enough to metabolize, reproduce, and spread before downward transport to lower altitudes that may be too hot for life of any kind to survive. Dynamical studies are needed to flesh out quantitative details of life particle residence times. A sub Neptune would need to be a part of a planetary system with an unstable asteroid belt in order for meteoritic material to provide nutrients, though life would also need to efficiently reuse and recycle metals. The origin of life may be the most severe limiting challenge. Regardless of the uncertainties, we can keep an open mind to the search for biosignature gases as a part of general observational studies of sub Neptune exoplanets.

Interstellar travel in the Milky Way is commonly thought to be a long and dangerous enterprise, but are all galaxies so hazardous? I introduce the concept of galactic traversability to address this question. Stellar populations are one factor in traversability, with higher stellar densities and velocity dispersions aiding rapid spread across a galaxy. The interstellar medium (ISM) is another factor, as gas, dust grains, and cosmic rays (CRs) all pose hazards to starfarers. I review the current understanding of these components in different types of galaxies, and conclude that red quiescent galaxies without star formation have favorable traversability. Compact elliptical galaxies and globular clusters could be "super-traversable", because stars are packed tightly together and there are minimal ISM hazards. Overall, if the ISM is the major hindrance to interstellar travel, galactic traversability increases with cosmic time as gas fractions and star formation decline. Traversability is a consideration in extragalactic surveys for the search for extraterrestrial intelligence (SETI).

A parsec-scale dusty torus is thought to be the cause of Active Galactic Nuclei (AGN) dichotomy in the 1/2 types, narrow/broad emission lines. In a previous work, on the basis of parsec-scale resolution infrared / optical dust maps it was found that dust filaments, few parsecs wide, several hundred parsecs long, were ubiquitous features crossing the centre of type 2 AGN, their optical thickness being sufficient to fully obscure the optical nucleus. This work presents the complementary view for type 1 and intermediate-type AGN. The same type of narrow, collimated, dust filaments are equally found at the centre of these AGN. The difference now resides in their location with respect to the nucleus, next to it but not crossing it, as it is the case in type 2, and their reduced optical thickness towards the centre, $A_V \lesssim 1.5 \rm{mag}$, insufficient to obscure at UV nucleus wavelengths. It is concluded that large scale, hundred pc to kpc long, dust filaments and lanes, reminiscent of those seen in the Milky Way, are a common ingredient to the central parsec of galaxies. Their optical thickness changes along their structure, in type 2 reaching optical depths high enough to obscure the nucleus in full. Their location with respect to the nucleus and increasing gradient in optical depth towards the centre could naturally lead to the canonical type 1/2 AGN classification, making these filaments to play the role of the torus. Dust filaments and lanes show equivalent morphologies in molecular gas. Available gas kinematic indicates mass inflows at rates $ < ~ 1 M\odot~ yr^{-1}$.

This paper reports on a new analysis of archival ALMA $870\,\mu$m dust continuum observations. Along with the previously observed bright inner ring ($r \sim 20-40\,$au), two addition substructures are evident in the new continuum image: a wide dust gap, $r \sim 40-150\,$au, and a faint outer ring ranging from $r \sim 150\,$au to $r \sim 250\,$au and whose presence was formerly postulated in low-angular-resolution ALMA cycle 0 observations but never before observed. Notably, the dust emission of the outer ring is not homogeneous, and it shows two prominent azimuthal asymmetries that resemble an eccentric ring with eccentricity $e = 0.07 $. The characteristic double-ring dust structure of HD 100546 is likely produced by the interaction of the disk with multiple giant protoplanets. This paper includes new smoothed-particle-hydrodynamic simulations with two giant protoplanets, one inside of the inner dust cavity and one in the dust gap. The simulations qualitatively reproduce the observations, and the final masses and orbital distances of the two planets in the simulations are 3.1 $M_{J}$ at 15 au and 8.5 $M_{J}$ at 110 au, respectively. The massive outer protoplanet substantially perturbs the disk surface density distribution and gas dynamics, producing multiple spiral arms both inward and outward of its orbit. This can explain the observed perturbed gas dynamics inward of 100 au as revealed by ALMA observations of CO. Finally, the reduced dust surface density in the $\sim 40-150\,$au dust gap can nicely clarify the origin of the previously detected H$_2$O gas and ice emission.

We discuss a probe of the contribution of wind-related shocks to the radio emission in otherwise radio-quiet quasars. Given 1) the non-linear correlation between UV and X-ray luminosity in quasars, 2) that such correlation leads to higher likelihood of radiation-line-driven winds in more luminous quasars, and 3) that luminous quasars are more abundant at high redshift, deep radio observations of high-redshift quasars are needed to probe potential contributions from accretion disk winds. We target a sample of 50 $z\simeq 1.65$ color-selected quasars that span the range of expected accretion disk wind properties as traced by broad CIV emission. 3-GHz observations with the Very Large Array to an rms of $\approx10\mu$Jy beam$^{-1}$ probe to star formation rates of $\approx400\,M_{\rm Sun}\,{\rm yr}^{-1}$, leading to 22 detections. Supplementing these pointed observations are survey data of 388 sources from the LOFAR Two-metre Sky Survey Data Release 1 that reach comparable depth (for a typical radio spectral index), where 123 sources are detected. These combined observations reveal a radio detection fraction that is a non-linear function of \civ\ emission-line properties and suggest that the data may require multiple origins of radio emission in radio-quiet quasars. We find evidence for radio emission from weak jets or coronae in radio-quiet quasars with low Eddingtion ratios, with either (or both) star formation and accretion disk winds playing an important role in optically luminous quasars and correlated with increasing Eddington ratio. Additional pointed radio observations are needed to fully establish the nature of radio emission in radio-quiet quasars.

We take a broad look at the problem of identifying the magnetic solar causes of space weather. With the lackluster performance of extrapolations based upon magnetic field measurements in the photosphere, we identify a region in the near UV part of the spectrum as optimal for studying the development of magnetic free energy over active regions. Using data from SORCE, Hubble Space Telescope, and SKYLAB, along with 1D computations of the near-UV (NUV) spectrum and numerical experiments based on the MURaM radiation-MHD and HanleRT radiative transfer codes, we address multiple challenges. These challenges are best met through a combination of near UV lines of bright \ion{Mg}{2}, and lines of \ion{Fe}{2} and \ion{Fe}{1} (mostly within the $4s-4p$ transition array) which form in the chromosphere up to $2\times10^4$ K. Both Hanle and Zeeman effects can in principle be used to derive vector magnetic fields. However, for any given spectral line the $\tau=1$ surfaces are generally geometrically corrugated owing to fine structure such as fibrils and spicules. By using multiple spectral lines spanning different optical depths, magnetic fields across nearly-horizontal surfaces can be inferred in regions of low plasma $\beta$, from which free energies, magnetic topology and other quantities can be derived. Based upon the recently-reported successful suborbital space measurements of magnetic fields with the CLASP2 instrument, we argue that a modest space-borne telescope will be able to make significant advances in the attempts to predict solar eruptions. Difficulties associated with blended lines are shown to be minor in an Appendix.

The $24^{\circ}$ N jet borders the North Tropical Belt and North Tropical Zone, and is the fastest prograde jet on Jupiter, reaching speeds above $170$ m/s. In this region, observations have shown several periodic convective plumes, likely from latent heat release from water condensation, which affect the cloud and zonal wind structure of the jet. We model this region with the Explicit Planetary hybrid-Isentropic Coordinate model using its active microphysics scheme to study the phenomenology of water and ammonia clouds within the jet region. On perturbing the atmosphere, we find that an upper tropospheric wave develops that directly influences the cloud structure within the jet. This wave travels at $\sim75$ m/s in our model, and leads to periodic chevron-shaped features in the ammonia cloud deck. These features travel with the wave speed, and are subsequently much slower than the zonal wind at the cloud deck. The cloud structure, and the slower drift rate, were both observed following the convective outbreak in this region in 2016 and 2020. We find that an upper level circulation is responsible for these cloud features in the aftermath of the convective outbursts. The comparatively slower observed drift rates of these features, relative to the wind speed of the jet, provides constraints on the vertical wind shear above the cloud tops, and we suggest that wind velocities determined from cloud tracking should correspond to a different altitude compared to the $680$ hPa pressure level. We also diagnose the convective potential of the atmosphere due to water condensation, and find that it is strongly coupled to the wave.

The energy cascade rate of turbulence can be measured with the structure function. In practice, the 3D velocity of the gas in molecular cloud is hard to measure, which makes the measurement of structure function difficult. In the case of thin molecular clouds perpendicular to the line of sight, the structure function $S^2_{ tt}$ can be measured with core velocity dispersion (CVD), ${\rm CVD}^2=\frac{1}{2}S^2_{ tt}$. This method was extended to the case when the thin molecular cloud is not perpendicular to the line of sight, with intersection angle $\theta$, ${\rm CVD}^2=\frac{1}{2}S^2_{ tt}\left(1-\frac{1}{8}\cos^2\theta\right)R^{2/3}$, where $R$ can be expressed with elliptic integrals of the second kind $E(k,\varphi)$ as $R=\frac{2}{\pi}E(\cos\theta,\frac{\pi}{2})$.

Observationally, the interstellar gas-phase abundance of deuterium (D) is considerably depleted and the missing D atoms are often postulated to have been locked up into carbonaceous solids and polycyclic aromatic hydrocarbon (PAH) molecules. An accurate knowledge of the fractional amount of D (relative to H) tied up in carbon dust and PAHs has important cosmological implications since D originated exclusively from the Big Bang and the present-day D abundance, after accounting for the astration it has experienced during the Galactic evolution, provides essential clues to the primordial nucleosynthesis and the cosmological parameters. To quantitatively explore the extent to which PAHs could possibly accommodate the observed D depletion, we have previously quantum-chemically computed the infrared vibrational spectra of mono-deuterated PAHs and derived the mean intrinsic band strengths of the 3.3 $\mu$m C--H stretch (A$_{3.3}$) and the 4.4 $\mu$m C--D stretch (A$_{4.4}$). Here we extend our previous work to multi-deuterated PAH species of different deuterations, sizes and structures. We find that both the intrinsic band strengths A$_{3.3}$ and A$_{4.4}$ and their ratios A$_{4.4}$/A$_{3.3}$ not only show little variations among PAHs of different deuterations, sizes and structures, they are also closely similar to that of mono-deuterated PAHs. Therefore, a PAH deuteration level (i.e., the fraction of peripheral atoms attached to C atoms in the form of D) of ~2.4% previously estimated from the observed 4.4 $\mu$m to 3.3 $\mu$m band ratio based on the A$_{4.4}$/A$_{3.3}$ ratio of mono-deuterated PAHs is robust.

(abridged) We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. The softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

We produce the first low to mid frequency radio simulation that incorporates both traditional extragalactic radio sources as well as synchrotron cosmic web emission. The FIlaments \& GAlactic RadiO (FIGARO) simulation includes ten unique \SI{4x4}{\degree} fields, incorporating active galactic nucleii (AGNs), star forming galaxies (SFGs) and synchrotron cosmic web emission out to a redshift of $z = 0.8$ and over the frequency range 100-1400 MHz. To do this, the simulation brings together a recent $100^3$ Mpc$^3$ magneto-hydrodynamic simulation (Vazza et al., 2019), calibrated to match observed radio relic population statistics, alongside updated `T-RECS' code for simulating extragalactic radio sources (Bonaldi et al., 2019). Uniquely, the AGNs and SFGs are populated and positioned in accordance with the underlying matter density of the cosmological simulation. In this way, the simulation provides an accurate understanding of the apparent morphology, angular scales, and brightness of the cosmic web as well as -- crucially -- the clustering properties of the cosmic web with respect to the embedded extragalactic radio population. We find that the synchrotron cosmic web does not closely trace the underlying mass distribution of the cosmic web, but is instead dominated by shocked shells of emission surrounding dark matter halos and resembles a large, undetected population of radio relics. We also show that, with accurate kernels, the cosmic web radio emission is clearly detectable by cross-correlation techniques and this signal is separable from the embedded extragalactic radio population. We offer the simulation as a public resource towards the development of techniques for detecting and measuring the synchrotron cosmic web.

Ten protostellar outflows in the Orion molecular clouds were mapped in the $^{12}$CO/$^{13}$CO ${J=6\rightarrow5}$ and $^{12}$CO ${J=7\rightarrow6}$ lines. The maps of these mid-$J$ CO lines have an angular resolution of about 10$''$ and a typical field size of about 100$''$. Physical parameters of the molecular outflows were derived, including mass transfer rates, kinetic luminosities, and outflow forces. The outflow sample was expanded by re-analyzing archival data of nearby low-luminosity protostars, to cover a wide range of bolometric luminosities. Outflow parameters derived from other transitions of CO were compared. The mid-$J$ ($J_{\rm up} \approx 6$) and low-$J$ ($J_{\rm up} \leq 3$) CO line wings trace essentially the same outflow component. By contrast, the high-$J$ (up to $J_{\rm up} \approx 50$) line-emission luminosity of CO shows little correlation with the kinetic luminosity from the ${J=6\rightarrow5}$ line, which suggests that they trace distinct components. The low/mid-$J$ CO line wings trace long-term outflow behaviors while the high-$J$ CO lines are sensitive to short-term activities. The correlations between the outflow parameters and protostellar properties are presented, which shows that the strengths of molecular outflows increase with bolometric luminosity and envelope mass.

We reported the gamma-ray observation towards the giant molecular cloud Polaris Flare. Together with the dust column density map, we derived the cosmic ray density and spectrum in this cloud. Compared with the CR measured locally, the CR density in Polaris Flare is significantly lower and the spectrum is softer. Such a different CR spectrum reveals either a rather large gradient of CR distribution in the direction perpendicular to the Galactic plane or a suppression of CR inside molecular clouds.

In 2013, we have published the first rotational analysis and detection of mono-deuterated dimethyl ether in the solar-type protostar IRAS 16293-2422 with the IRAM 30m telescope. Dimethyl ether is one of the most abundant complex organic molecules (COMs) in star-forming regions and their D-to-H (D/H) ratios is important to understand its chemistry and trace the source history. We present the first analysis of doubly-deuterated dimethyl ether (methoxy-d2-methane, 1,1-dideuteromethylether) in its ground-vibrational state, based on an effective Hamiltonian for an asymmetric rotor molecules with internal rotors. The analysis covers the frequency range 0.15-1.5THz. The laboratory rotational spectrum of this species was measured between 150 and 1500 GHz with the Lille's submillimeter spectrometer. For the astronomical detection, we used the Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Protostellar Interferometric Line Survey, PILS. New sets of spectroscopic parameters have been determined by a least squares fit with the ERHAM code for both symmetric and asymmetric conformers. As for the mono-deuterated species, these parameters have permitted the first identification in space of both conformers of a doubly-deuterated dimethyl ether via detection near the B component of the Class 0 protostar IRAS 16293-2422.

Thanks to modern understanding of stellar evolution, we can accurately measure the age of Open Clusters (OCs). Given their position, they are ideal tracers of the Galactic disc. Gaia data release 2, besides providing precise parallaxes, led to the detection of many new clusters, opening a new era for the study of the Galactic disc. However, detailed information on the chemical abundance for OCs is necessary to accurately date them and to efficiently use them to probe the evolution of the disc.Mapping and exploring the Milky Way structure %to combine accurate chemical information of OCs is the main aim of the Stellar Population Astrophysics (SPA) project. Part of this work involves the use of OCs and the derivation of their precise and accurate chemical composition.We analyze here a sample of OCs located within about 2 kpc from the Sun, with ages from about 50 Myr to a few Gyr.We used HARPS-N at the Telescopio Nazionale Gaileo and collected very high-resolution spectra (R = 115\,000) of 40 red giant/red clump stars in 18 OCs (16 never or scarcely studied plus two comparison clusters). We measured their radial velocities and derived the stellar parameters.We discussed the relationship between metallicity and Galactocentric distance, adding literature data to our results to enlarge the sample and taking also age into account. We compared the result of observational data with that from chemo-dynamical models. These models generally reproduce the metallicity gradient well. However, at young ages we found a large dispersion in metallicity, not reproduced by models. Several possible explanations are explored, including uncertainties in the derived metallicity. We confirm the difficulties in determining parameters for young stars (age < 200 Myr), due to a combination of intrinsic factors which atmospheric models can not easily reproduce and which affect the parameters uncertainty

We present a radio search for WIMP dark matter in the Large Magellanic Cloud (LMC). We make use of a recent deep image of the LMC obtained from observations of the Australian Square Kilometre Array Pathfinder (ASKAP), and processed as part of the Evolutionary Map of the Universe (EMU) survey. LMC is an extremely promising target for WIMP searches at radio frequencies because of the large J-factor and the presence of a substantial magnetic field. We detect no evidence for emission arising from WIMP annihilations and derive stringent bounds. This work excludes the thermal cross section for masses below 480 GeV and annihilation into quarks.

We study the propagation of a specific class of instrumental systematics to the reconstruction of the B-mode power spectrum of the cosmic microwave background (CMB). We focus on non-idealities of the half-wave plate (HWP), a polarization modulator that will be deployed by future CMB experiments, such as the phase-A satellite mission LiteBIRD. More in details, we study the effects of non-ideal HWP properties, such as transmittance, phase shift and cross-polarization. To this purpose, we develop a simple, yet stand-alone end-to-end simulation pipeline adapted to LiteBIRD. Through the latter, we analyze the effects of a possible mismatch between the measured frequency profiles of HWP properties (used in the mapmaking stage of the pipeline) and the actual profiles (used in the sky-scanning step). We simulate single-frequency, CMB-only observations to emphasize the effects of non-idealities on the BB power spectrum. We also consider multi-frequency observations to account for the frequency dependence of HWP properties and the contribution of foreground emission. We quantify the systematics effects in terms of a bias $\Delta r$ on the tensor-to-scalar ratio $r$ with respect to the ideal case of no-systematics. We derive the accuracy requirements on the measurements of HWP properties by requiring $\Delta r < 10^{-5}$ (1% of the expected LiteBIRD sensitivity on $r$). The analysis is introduced by a detailed presentation of the mathematical formalism employed in this work, including the use of the Jones and Mueller matrix representations.

The Calorimetric Electron Telescope (CALET), in operation on the International Space Station since 2015, collected a large sample of cosmic-ray iron over a wide energy interval. In this Letter a measurement of the iron spectrum is presented in the range of kinetic energy per nucleon from 10 GeV$/n$ to 2.0 TeV$/n$ allowing the inclusion of iron in the list of elements studied with unprecedented precision by space-borne instruments. The measurement is based on observations carried out from January 2016 to May 2020. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number $Z$ = 40). The energy is measured by a homogeneous calorimeter with a total equivalent thickness of 1.2 proton interaction lengths preceded by a thin (3 radiation lengths) imaging section providing tracking and energy sampling. The analysis of the data and the detailed assessment of systematic uncertainties are described and results are compared with the findings of previous experiments. The observed differential spectrum is consistent within the errors with previous experiments. In the region from 50 GeV$/n$ to 2 TeV$/n$ our present data are compatible with a single power law with spectral index -2.60 $\pm$ 0.03.

The electric field accelerating electrons in the Timokhin-Arons polar-cap model as applied to millisecond pulsars is so high that the electron Lorentz factors are limited either by radiation reaction or by the Breit-Wheeler process. In the former case, it is possible to obtain an upper limit for curvature radiation momentum components perpendicular to the local magnetic field which is independent of the flux-line radius of curvature. The threshold value for single-photon conversion to pairs is high but is possibly reached in J0030+0451. However, owing to the high polar-cap temperature reported, direct pair production by the Breit-Wheeler process is probably more important. If the existence of coherent radio emission in millisecond pulsars is assumed to need a high-multiplicity pair plasma, then it follows that the primary electrons also produce gamma-rays in the Fermi-LAT energy band for which the magnetosphere is completely transparent. The absence of these, in phase with the radio emission, would be an immediate indication that ultra-high energy electrons and an active Timokhin-Arons polar cap are not present in J0030+0451.

This article explains an optimization-based approach for counting and localizing stars within a small cluster, based on photon counts in a focal plane array. The array need not be arranged in any particular way, and relatively small numbers of photons are required in order to ensure convergence. The stars can be located close to one another, as the location and brightness errors were found to be low when the separation was larger than $0.2$ Rayleigh radii. To ensure generality of our approach, it was constructed as a special case of a general theory built upon topological signal processing using the mathematics of sheaves.

The second data release of \it Gaia \rm revealed a parallax zero point offset of $-0.029$~mas based on quasars. The value depended on the position on the sky, and also likely on magnitude and colour. The offset and its dependence on other parameters inhibited an improvement in the local distance scale using e.g. the Cepheid and RR Lyrae period-luminosity relations. Analysis of the recent \it Gaia \rm Early Data Release 3 (EDR3) reveals a mean parallax zero point offset of $-0.021$~mas based on quasars. The \it Gaia \rm team addresses the parallax zero point offset in detail and proposes a recipe to correct for it, based on ecliptic latitude, $G$-band magnitude, and colour information. This paper is a completely independent investigation into this issue focussing on the spatial dependence of the correction based on quasars and the magnitude dependence based on wide binaries. The spatial and magnitude corrections are connected to each other in the overlap region between $17 < G < 19$. The spatial correction is presented at several spatial resolutions based on the HEALPix formalism. The colour dependence of the parallax offset is unclear and in any case secondary to the spatial and magnitude dependence. The spatial and magnitude corrections are applied to two samples of brighter sources, namely a sample of $\sim$100 stars with independent trigonometric parallax measurements from \it HST \rm data, and a sample of 75 classical cepheids using photometric parallaxes. The mean offset between the observed GEDR3 parallax and the independent trigonometric parallax (excluding outliers) is about $-39$~\muas, and after applying the correction it is consistent with being zero. For the classical cepheid sample it is suggested that the photometric parallaxes may be underestimated by about 5\%.

The Zeeman effect and dust grain alignment are two major methods for probing magnetic fields (B-fields) in molecular clouds, largely motivated by the study of star formation, as the B-field may regulate gravitational contraction and channel turbulence velocity. This review summarizes our observations of B-fields over the past decade, along with our interpretation. Galactic B-fields anchor molecular clouds down to cloud cores with scales around 0.1 pc and densities of 10^4-5 H_2/cc. Within the cores, turbulence can be slightly super-Alfvenic, while the bulk volumes of pa-rental clouds are sub-Alfvenic. The consequences of these largely ordered cloud B-fields on fragmentation and star formation are observed. The above paradigm is very different from the generally accepted theory during the first decade of the century, when cloud turbulence was assumed to be highly super-Alfvenic. Thus, turbulence anisotropy and turbulence-induced ambipolar diffusion are also revisited.

We use ultraviolet imaging taken with the XMM-Newton Optical Monitor telescope (XMM-OM), covering 280 square arcminutes in the UVW1 band (effective wavelength 2910 Angstroms) to measure rest-frame ultraviolet (1500 Angstrom) luminosity functions of galaxies with redshifts z between 0.6 and 1.2. The XMM-OM data are supplemented by a large body of optical and infrared imaging to provide photometric redshifts. The XMM-OM data have a significantly narrower point-spread-function (resulting in less source confusion) and simpler K-correction than the GALEX data previously employed in this redshift range. Ultraviolet-bright active galactic nuclei are excluded to ensure that the luminosity functions relate directly to the star-forming galaxy population. Binned luminosity functions and parametric Schechter-function fits are derived in two redshift intervals: 0.6<z<0.8 and 0.8<z<1.2. We find that the luminosity function evolves such that the characteristic absolute magnitude M* is brighter for 0.8<z<1.2 than for 0.6<z<0.8.

The Fermi-LAT DR1 and DR2 4FGL catalogues feature more than 5000 gamma-ray sources of which about one fourth are not associated with already known objects, and approximately one third are associated with blazars of uncertain nature. We perform a three-category classification of the 4FGL DR1 and DR2 sources independently, using an ensemble of Artificial Neural Networks (ANNs) to characterise them based on the likelihood of being a Pulsar (PSR), a BL Lac type blazar (BLL) or a Flat Spectrum Radio Quasar (FSRQ). We identify candidate PSR, BLL and FSRQ among the unassociated sources with approximate equipartition among the three categories and select ten classification outliers as potentially interesting for follow up studies.

We examine the composition of barium stars in the context of mass transfer from an asymptotic giant branch (AGB) companion. We accrete between 0.01 and 0.5 M$_\odot$ of AGB ejecta on to low mass companions of [Fe/H] = -0.25 at the ages expected for the end of the lives of AGB stars of 2.5, 3 and 4M$_\odot$. In each case, we form a star of 2.5 M$_\odot$ which is thought to be a typical barium star mass. We discuss the extent of dilution of accreted material as the star evolves, and describe the impact on the surface abundances. For accretion from a 2.5\ms\ primary, if the secondary's initial mass is 2.45 M$_\odot$ or more, accretion takes place when the secondary is undergoing core helium burning. Using data from the sample of De Castro et al., we attempt to fit the observed properties of 74 barium giants using the models we have computed. We find that all but six of these objects are best fit using ejecta from 2.5 M$_\odot$ (32 objects) or 3 M$_\odot$ (36 objects) AGB stars. Higher accretion masses are typically required when accreting from a lower mass companion. We find accretion masses that are broadly consistent with recent hydrodynamical simulations of wind mass transfer, though the accretion efficiency is toward the upper limit found in these simulations. For the 18 stars with reported orbital periods, we find no strong correlations between period and accretion mass.

Asymptotic giant branch (AGB) stars are considered to be among the most significant contributors to the fluorine budget in our Galaxy. While at close-to-solar metallicity observations and theory agree, at lower metallicities stellar models overestimate the fluorine production with respect to heavy elements. We present ${}^{19}$F nucleosynthesis results for a set of AGB models with different masses and metallicities in which magnetic buoyancy acts as the driving process for the formation of the ${}^{13}$C neutron source (the so-called ${}^{13}$C pocket). We find that ${}^{19}$F is mainly produced as a result of nucleosynthesis involving secondary ${}^{14}$N during convective thermal pulses, with a negligible contribution from the ${}^{14}$N present in the ${}^{13}$C pocket region. A large ${}^{19}$F production is thus prevented, resulting in lower fluorine surface abundances. As a consequence, AGB stellar models with magnetic-buoyancy-induced mixing at the base of the convective envelope well agree with available fluorine spectroscopic measurements at both low and close-to-solar metallicity.

We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~\mu\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~\mu\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$.

Getman et al. (2021) reports the discovery, energetics, frequencies, and effects on environs of $>1000$ X-ray super-flares with X-ray energies $E_X \sim 10^{34}-10^{38}$~erg from pre-main sequence (PMS) stars identified in the Chandra MYStIX and SFiNCs surveys. Here we perform detailed plasma evolution modeling of 55 bright MYStIX/SFiNCs super-flares from these events. This is the largest sample of highly energetic flares analyzed in a uniform fashion. They are compared with published X-ray super-flares from young stars in the Orion Nebula Cluster, older active stars, and the Sun. Several results emerge. First, the properties of PMS super-flares are independent of the presence or absence of protoplanetary disks, supporting the solar-type model of PMS flaring magnetic loops with both footpoints anchored in the stellar surface. Second, most PMS super-flares resemble solar long-duration events associated with coronal mass ejections. Slow rise PMS super-flares are an interesting exception. Third, strong correlations of super-flare peak emission measure and plasma temperature with the stellar mass are similar to established correlations for the PMS X-ray emission likely composed of numerous smaller flares. Fourth, a new correlation of loop geometry is linked to stellar mass; more massive stars have thicker flaring loops. Finally, the slope of a long-standing relationship between the X-ray luminosity and magnetic flux of various solar-stellar magnetic elements appears steeper in PMS super-flares than for solar events.

Galaxy flybys are as common as mergers in low redshift universe and are important for galaxy evolution as they involve the exchange of significant amounts of mass and energy. In this study we investigate the effect of minor flybys on the bulges, disks, and spiral arms of Milky Way mass galaxies for two types of bulges - classical bulges and boxy/peanut pseudobulges. Our N-body simulations comprise of two disk galaxies of mass ratios 10:1 and 5:1, where the disks of the galaxies lie in their orbital plane and the pericenter distance is varied. We performed photometric and kinematic bulge-disk decomposition at regular time steps and traced the evolution of the disk size, spiral structure, bulge sersic index, bulge mass, and bulge angular momentum. Our results show that the main effect on the disks is disk thickening, which is seen as the increase in the ratio of disk scale height to scale radius. The strength of the spiral structure A2/A0 shows small oscillations about the mean time-varying amplitude in the pseudobulge host galaxies. The flyby has no significant effect on non-rotating classical bulge, which shows that these bulges are extremely stable in galaxy interactions. However, the pseudobulges become dynamically hotter in flybys indicating that flybys may play an important role in accelerating the rate of secular evolution in disk galaxies. This effect on pseudobulges is a result of their rotating nature as part of the bar. Also, flybys do not affect the time and strength of bar buckling.

The abundance, distribution and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross-section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos -- which may host orphan galaxies -- are taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structure of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak lensing profile and can be used to constrain $\sigma_T/m$ at the $\sim 1$ cm$^2$ g$^{-1}$ level at $v\sim 2000$ km s$^{-1}$.

We present a novel framework to solve simultaneously the electroweak hierarchy problem and the strong-CP problem. A small but finite Higgs vacuum expectation value and a small $\theta$-angle are selected after the QCD phase transition, without relying on the Peccei-Quinn mechanism or other traditional solutions. We predict a distinctive pattern of correlated signals at hadronic EDM, fuzzy dark matter and axion experiments.

This work presents the first steps to modelling synthetic rovibrational spectra for all molecules of astrophysical interest using the new code Prometheus. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is on molecules identified as signatures of life in planetary atmospheres (biosignatures). As a starting point, in this work we evaluate the accuracy of our method by studying the diatomics molecules H$_2$, O$_2$, N$_2$ and CO, all of which have well-known spectra. Prometheus uses the Transition-Optimised Shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental $\nu=0 \rightarrow \nu=1$ band, along with thermal profile modeling for the rotational transitions. We present a novel new application of the TOSH theory with regards to rotational constants. Our results show that this method can achieve results that are a better approximation than the ones produced through the basic harmonic method. We discuss the current limitations of our method. In particular, we compare our results with high-resolution HITRAN spectral data. We find that modelling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.

Demanding that charged Nariai black holes in (quasi-)de Sitter space decay without becoming super-extremal implies a lower bound on the masses of charged particles, known as the Festina Lente (FL) bound. In this paper we fix the $\mathcal{O}(1)$ constant in the bound and elucidate various aspects of it, as well as extensions to $d>4$ and to situations with scalar potentials and dilatonic couplings. We also discuss phenomenological implications of FL including an explanation of why the Higgs potential cannot have a local minimum at the origin, thus explaining why the weak force must be broken. For constructions of meta-stable dS involving anti-brane uplift scenarios, even though the throat region is consistent with FL, the bound implies that we cannot have any light charged matter fields coming from any far away region in the compactified geometry, contrary to the fact that they are typically expected to arise in these scenarios. This strongly suggests that introduction of warped anti-branes in the throat cannot be decoupled from the bulk dynamics as is commonly assumed. Finally, we provide some evidence that in certain situations the FL bound can have implications even with gravity decoupled and illustrate this in the context of non-compact throats.

In this paper, we study the spin transitions of neutrinos caused by the interaction with a gravitational field. We consider a model with a scalar field (describing screening effects) conformally coupled to matter and neutrinos. The presence of screening effects suppresses the neutrino spin-flip probability as compared with General Relativity predictions. Such a result could be used, combined with neutrino astronomy, for testing modified theories of gravity and, in turn, screening effects invoked to bypass the solar system and Lab tests. Such an analysis has been also extended to the case of the quintessence field surrounding a black hole. Here we investigate the flavor and spin transitions, showing that also in such a case exists a suppression of the effect compared to General Relativity prediction.

We study decoherence effects in neutrino flavor oscillations in curved spacetime with particular emphasis on the lensing in a Schwarzschild geometry. Assuming Gaussian wave packets for neutrinos, we argue that the decoherence length derived from the exponential suppression of the flavor transition probability depends on the proper time of the geodesic connecting the events of the production and detection in general gravitational setting. In the weak gravity limit, the proper time between two events of given proper distance is smaller than that in the flat spacetime. Therefore, in presence of a Schwarzschild object, the neutrino wave packets have to travel relatively more physical distance in space to lapse the same amount of proper time before they decoher. For non-radial propagation applicable to the lensing phenomena, we show that the decoherence, in general, is sensitive to the absolute values of neutrino masses as well as the classical trajectories taken by neutrinos between the source and detector along with the spatial widths of neutrino wave packets. At distances beyond the decoherence length, the probability of neutrino flavor transition due to lensing attains a value which depends only on the leptonic mixing parameters. Hence, the observability of neutrino lensing significantly depends on these parameters and in-turn the lensing can provide useful information about the latter.

Typicality arguments attempt to use the Copernican Principle to draw conclusions about the cosmos and presently unknown conscious beings within it. The most notorious is the Doomsday Argument, which purports to constrain humanity's future from its current lifespan alone. These arguments rest on a likelihood calculation that penalizes models in proportion to the number of distinguishable observers. I argue that such reasoning leads to solipsism, the belief that one is the only being in the world, and is therefore unacceptable. Using variants of the "Sleeping Beauty" thought experiment as a guide, I present a framework for evaluating observations in a large cosmos: Fine Graining with Auxiliary Indexicals (FGAI). FGAI requires the construction of specific models of physical outcomes and observations. Valid typicality arguments then emerge from the combinatorial properties of third-person physical microhypotheses. Indexical (observer-relative) facts do not directly constrain physical theories. Instead they serve to weight different provisional evaluations of credence. These weights define a probabilistic reference class of locations. As indexical knowledge changes, the weights shift. I show that the self-applied Doomsday Argument fails in FGAI, even though it can work for an external observer. I also discuss how FGAI could handle observations in large universes with Boltzmann brains.

In this master's thesis we study the production of axion dark matter through the so-called misalignment mechanism by considering that during that time, the universe was dominated by a new kind of fluid, different than radiation. We perform a very detailed analysis of the oscillation temperature and the relic density today, both analytically and numerically. Our findings show that on the one hand, the oscillation temperature is strongly influenced by the non-standard cosmology, affecting the relic density, and on the other hand, the energy density of the axion gets diluted, because the new fluid eventually decays, injecting entropy into the thermal bath. We find the predicted parameter space of axion dark matter for different non-standard cosmologies and we show its impact on the coupling of axions to two photons.

We present a novel universal relation for binary neutron star mergers with long-lived neutron star remnants: inspired by recent work based on numerical relativity simulations, we propose a novel approach using perturbative calculations that allow us to relate the pre-merger neutron star binary tidal deformability to the effective compactness of the post-merger remnant. Our results allow for the prediction of the stellar parameters of a long-lived remnant neutron star from the study of gravitational waves emitted during the pre-merger phase.

In this paper we study the possibility of having a wormhole (WH) as a candidate for the Sgr A$^\star$ central object and test this idea by constraining their geometry using the motion of S2 star and the reconstructed shadow images. In particular, we consider three WH models, including WHs in Einstein theory, brane-world gravity, and Einstein-Dirac-Maxwell theory. To this end, we have constrained the WH throat using the motion of S2 star and shown that the flare out condition is satisfied. We also consider the accretion of infalling gas model and study the accretion rate and the intensity of the electromagnetic radiation as well as the shadow images.

Light bosonic scalars (e.g. axions) may form clouds around black holes via superradiant instabilities, or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we study numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. We show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar the rapid force oscillations could in principle leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral.