Aims: We investigate electron temperature (T_e) and gas-phase oxygen abundance (Z_Te) measurements for galaxies in the local Universe (z < 0.25). Our sample comprises spectra from a total of 264 emission-line systems, ranging from individual HII regions to whole galaxies, including 23 composite HII regions from star-forming main sequence galaxies in the MaNGA survey.
Methods: We utilise 130 of these systems with directly measurable T_e(OII) to calibrate a new metallicity-dependent T_e(OIII)-T_e(OII) relation that provides a better representation of our varied dataset than existing relations from the literature. We also provide an alternative T_e(OIII)-T_e(NII) calibration. This new T_e method is then used to obtain accurate Z_Te estimates and form the mass - metallicity relation (MZR) for a sample of 118 local galaxies.
Results: We find that all the T_e(OIII)-T_e(OII) relations considered here systematically under-estimate Z_Te for low-ionisation systems by up to 0.6 dex. We determine that this is due to such systems having an intrinsically higher O⁺ abundance than O⁺⁺ abundance, rendering Z_Te estimates based only on [OIII] lines inaccurate. We therefore provide an empirical correction based on strong emission lines to account for this bias when using our new T_e(OIII)-T_e(OIII) and T_e(OIII)-T_e(NII) relations. This allows for accurate metallicities (1σ = 0.08 dex) to be derived for any low-redshift system with an [OIII]λ4363 detection, regardless of its physical size or ionisation state. The MZR formed from our dataset is in very good agreement with those formed from direct measurements of metal recombination lines and blue supergiant absorption lines, in contrast to most other T_e-based and strong-line-based MZRs. Our new T_e method therefore provides an accurate and precise way of obtaining Z_Te for a large and diverse range of star-forming systems in the local Universe.

We present a real-space renormalization group transformation with continuous scale change to calculate the continuous renormalization group β function in nonperturbative lattice simulations. Our method is motivated by the connection between Wilsonian renormalization group and the gradient flow transformation. It does not rely on the perturbative definition of the renormalized coupling and is also valid at nonperturbative fixed points. Although our method requires an additional extrapolation compared to traditional step scaling calculations, it has several advantages which compensates for this extra step even when applied in the vicinity of the perturbative fixed point. We illustrate our approach by calculating the β function of 2-flavor QCD and show that lattice predictions from individual lattice ensembles, even without the required continuum and finite volume extrapolations, can be very close to the result of the full analysis. Thus our method provides a nonperturbative framework and intuitive understanding into the structure of strongly coupled systems, in addition to being complementary to existing lattice determinations.

We analyze the prospects of probing the C P -odd i κ ∼t ¯γ ⁵th interaction at the LHC and its projected upgrades, the high-luminosity and high-energy LHC, directly using associated on-shell Higgs boson and top quark or top quark pair production. To this end we first construct a C P -odd observable based on top quark polarization in Wb → th scattering with optimal linear sensitivity to κ ∼. For the corresponding hadronic process pp → thj we present a method of extracting the phase-space dependent weight function that allows to retain close to optimal sensitivity to κ ∼. We project future sensitivity to the signal in pp → t(→ ℓ

In this paper, we extend the collinear superspace formalism to include the full range of N = 1 supersymmetric interactions. Building on the effective field theory rules developed in a companion paper — Navigating Collinear Superspace [1] — we construct collinear superspace Lagrangians for theories with non-trivial F- and D-term auxiliary fields. For (massless) Wess-Zumino models, the key ingredient is a novel type of Grassmann-valued supermultiplet whose lowest component is a (non-propagating) fermionic degree of freedom. For gauge theories coupled to charged chiral matter, the key ingredient is a novel type of vector superfield whose lowest component is a non-propagating gauge potential. This unique vector superfield is used to construct a gauge-covariant derivative; while such an object does not appear in the standard full superspace formalism, it is crucial for modeling gauge interactions when the theory is expres sed on a collinear slice. This brings us full circle, by showing that all types of N = 1 theories in four dimensions can beconstructed in collinear superspace from purely infrared considerations. We speculate that supersymmetric theories with N > 1 could also be implemented using similar collinear superspace constructions.

We recently introduced in [9] a boundary-to-bound dictionary between gravitational scattering data and observables for bound states of non-spinning bodies. In this paper, we elaborate further on this holographic map. We start by deriving the following — remarkably simple — formula relating the periastron advance to the scattering angle: ΔΦ (" separators=",J E )=χ (" separators=",J E )+χ (" separators=",-J E ), via analytic continuation in angular momentum and binding energy. Using explicit expressions from [9], we confirm its validity to all orders in the Post-Minkowskian (PM) expansion. Furthermore, we reconstruct the radial action for the bound state directly from the knowledge of the scattering angle. The radial action enables us to write compact expressions for dynamical invariants in terms of the deflection angle to all PM orders, which can also be written as a function of the PM-expanded amplitude. As an example, we reproduce our result in [9] for the periastron advance, and compute the radial and azimuthal frequencies and redshift variable to two-loops. Agreement is found in the overlap between PM and Post-Newtonian (PN) schemes. Last but not least, we initiate the study of our dictionary including spin. We demonstrate that the same relation between deflection angle and periastron advance applies for aligned-spin contributions, with J the (canonical) total angular momentum. Explicit checks are performed to display perfect agreement using state-of-the-art PN results in the literature. Using the map between test- and two-body dynamics, we also compute the periastron advance up to quadratic order in spin, to one-loop and to all orders in velocity. We conclude with a discussion on the generalized `impetus formula' for spinning bodies and black holes as `elementary particles'. Our findings here and in [9] imply that the deflection angle already encodes vast amount of physical information for bound orbits, encouraging independent derivations using numerical and/or self-force methodologies.

Galaxy assembly bias, the correlation between galaxy properties and halo properties at fixed halo mass, could be an important ingredient in halo-based modelling of galaxy clustering. We investigate the central galaxy assembly bias by studying the relation between various galaxy and halo properties in the Illustris hydrodynamic galaxy formation simulation. Galaxy stellar mass M_* is found to have a tighter correlation with peak maximum halo circular velocity V_peak than with halo mass M_h. Once the correlation with V_peak is accounted for, M_* has nearly no dependence on any other halo assembly variables. The correlations between galaxy properties related to star formation history and halo assembly properties also show a cleaner form as a function of V_peak than as a function of M_h, with the main correlation being with halo formation time and to a less extent halo concentration. Based on the galaxy-halo relation, we present a simple model to relate the bias factors of a central galaxy sample and the corresponding halo sample, both selected based on assembly-related properties. It is found that they are connected by the correlation coefficient of the galaxy and halo properties used to define the two samples, which provides a reasonable description for the samples in the simulation and suggests a simple prescription to incorporate galaxy assembly bias into the halo model. By applying the model to the local galaxy clustering measurements in Lin et al., we infer that the correlation between star formation history or specific star formation rate and halo formation time is consistent with being weak.

In light of prospects for measurements of B_c → D^(∗)lv decays in the upcoming Upgrade II of the LHC, we show that by using calculated B_c → D^(∗) form factors a competitive extraction of the |V_ub| CKM matrix element from the B_c→Dμ v_¯μ decay might be possible. To minimize experimental and theoretical uncertainties we provide the ratio |V_ub|/|V_cb| by normalizing the B_c→D^(∗)μ v_¯μ to B_c→J /ψμ v_¯μ decay. We also briefly examine the suggestion to extract |V_ub|/|V_cs| from the theoretically interesting ratio of B_c→D⁰e v_¯e and B_c→B_se v_¯e decay rates in the zero-recoil limit. With the present average value of |V_ub|, the predicted branching ratios are estimated to be BR (B_c→D⁰μ ν_¯μ )=(2.4 ±0.4 ).10^-5 and BR (B_c→D^∗μ ν_¯μ )=(7 ±3 ).10^-5, and the semileptonic ratios for testing the lepton flavour universality in these B_c decays are R_c(D⁰) = 0.64 ± 0.05 and R_c(D^∗) = 0.55 ± 0.05. We also provide q² distributions and various angular observables of B_c→ D^(∗)lν decays.

Effective field theory (EFT) approaches are widely used at the Large Hadron Collider (LHC), such that it is important to study their validity and ease of matching to specific new physics models. In this paper, we consider an extension of the Standard Model (SM) in which a top quark couples to a new heavy scalar. We find the dimension six operators generated by this theory at low energy and match the EFT to the full theory up to the next-to-leading order (NLO) precision in the simplified model coupling. We then examine the range of validity of the EFT description in top pair production, finding excellent validity even if the scalar mass is only slightly above LHC energies, provided NLO corrections are included. In the absence of the latter, the LO EFT overestimates kinematic distributions, such that overoptimistic constraints on beyond the Standard Model (BSM) contributions are obtained. We next examine the constraints on the EFT and full models that are expected to be obtained from both top pair and four top production at the LHC, finding for low scalar masses that both processes show similar exclusion power. However, for larger masses, estimated LHC uncertainties push constraints into the nonperturbative regime, where the full model is difficult to analyze, and thus is not perturbatively matchable to the EFT. This highlights the necessity to improve uncertainties of SM hypotheses in top final states.

Three hidden-charm pentaquark P_c states, P_c(4312 ), P_c(4440 ), and P_c(4457 ) were revealed in the Λ_b⁰→J /ψ p K^- process measured by LHCb using both run I and run II data. Their nature is under lively discussion, and their quantum numbers have not been determined. We analyze the J /ψ p invariant mass distributions under the assumption that the crossed-channel effects provide a smooth background. For the first time, such an analysis is performed employing a coupled-channel formalism with the scattering potential involving both one-pion exchange as well as short-range operators constrained by heavy quark spin symmetry. We find that the data can be well described in the hadronic molecular picture, which predicts seven Σ_c^(*)D^¯(*) molecular states in two spin multiplets, such that the P_c(4312 ) is mainly a Σ_cD ¯ bound state with J^P=1 /2^-, while P_c(4440 ) and P_c(4457 ) are Σ_cD^¯* bound states with quantum numbers 3 /2^- and 1 /2^-, respectively. We also show that there is evidence for a narrow Σ_c^*D ¯ bound state in the data which we call P_c(4380 ), different from the broad one reported by LHCb in 2015. With this state included, all predicted Σ_cD ¯, Σ_c^*D ¯, and Σ_cD^¯* hadronic molecules are seen in the data, while the missing three Σ_c^*D^¯* states are expected to be found in future runs of the LHC or in photoproduction experiments.

We explored the possibility that Higgs coupling to new physics violates flavor universality. In particular, we parametrize such models with dimension-six effective operators which modify the coupling between the first generation quarks, Higgs boson, and Z boson. Through the use of boosted Higgsstrahlung events at both the HL-LHC and potential future hadron colliders, as well as existing ATLAS data for background estimates, we projected constraints on the scale of new physics as a function of the Wilson coefficient. The high energy Z h process will provide unique information about these class of operators, and the sensitivity is competitive with the LEP electroweak precision measurements. We include different scenarios of the overall systematic uncertainties and the PDF uncertainties when presenting the projected sensitivities. We also discuss the constraints from flavor changing neutral currents to these flavor-violating models and the complementarity of the exotic Higgs decay to the Z h process.

For direct detection of sub-MeV dark matter, a promising strategy is to search for individual phonon excitations in a crystal. We perform an analytic calculation of the rate for light dark matter (keV <m_DM<MeV ) to produce two acoustic phonons through scattering in cubic crystals such as GaAs, Ge, Si, and diamond. The multiphonon rate is always smaller than the rate to produce a single optical phonon, whenever the latter is kinematically accessible. In Si and diamond, there is a dark matter mass range for which multiphonon production can be the most promising process, depending on the experimental threshold.

We combine the NLTE spectral analysis of the detached O-type eclipsing binary OGLE-LMC-ECL-06782 with the analysis of the radial velocity curve and light curve to measure an independent distance to the Large Magellanic Cloud (LMC). In our spectral analysis we study composite spectra of the system at quadrature and use the information from radial velocity and light curve about stellar gravities, radii, and component flux ratio to derive effective temperature, reddening, extinction, and intrinsic surface brightness. We obtain a distance modulus to the LMC of m - M = 18.53 ± 0.04 mag. This value is 0.05 mag larger than the precision distance obtained recently from the analysis of a large sample of detached, long period late spectral type eclipsing binaries but agrees within the margin of the uncertainties. We also determine the surface brightnesses of the system components and find good agreement with the published surface brightness-color relationship. A comparison of the observed stellar parameters with the prediction of stellar evolution based on the MESA stellar evolution code shows reasonable agreement, but requires a reduction of the internal angular momentum transport to match the observed rotational velocities.

Origins of contemporary $B$-physics. Mesons with beauty and charm. Stable tetraquarks? Flavor and the problem of identity. Top matters. Electroweak symmetry breaking and the Higgs sector. Future instruments.

We show that for a wide range of stellar masses, from 0.3 to 20 M_⊙, and for evolutionary phases from the main sequence to the beginning of the red giant stage, the stellar flux-weighted gravity, g_F ≅ g/ ${T}_{\mathrm{eff}}^{4}$ , is tightly correlated with absolute bolometric magnitude ${M}_{\mathrm{bol}}$ . Such a correlation is predicted by stellar evolution theory. We confirm this relation observationally, using a sample of 445 stars with precise stellar parameters. It holds over 17 stellar magnitudes from ${M}_{\mathrm{bol}}$ = 9.0 to -8.0 mag with a scatter of 0.17 mag above ${M}_{\mathrm{bol}}$ = -3.0 and 0.29 mag below this value. We then test the relation with 2.2 million stars with 6.5 mag ≥ ${M}_{\mathrm{bol}}$ ≥ 0.5 mag, where "mass-produced" but robust $\mathrm{log}\,g$ , ${T}_{{\rm{e}}{\rm{f}}{\rm{f}}},$ and ${M}_{\mathrm{bol}}$ from LAMOST DR5 and Gaia DR2 are available. We find that the same relation holds with a scatter of ∼0.2 mag for single stars offering a simple spectroscopic distance estimate good to ∼10%.

We carry out an analysis of the full set of ten B¯→D(∗) form factors within the framework of the Heavy-Quark Expansion (HQE) to order (αs,1/mb,1/m2c), both with and without the use of experimental data. This becomes possible due to a recent calculation of these form factors at and beyond the maximal physical recoil using QCD light-cone sum rules, in combination with constraints from lattice QCD, QCD three-point sum rules and unitarity. We find good agreement amongst the various theoretical results, as well as between the theoretical results and the kinematical distributions in B¯→D(∗){e−,μ−}ν¯ measurements. The coefficients entering at the 1/m2c level are found to be of (1), indicating convergence of the HQE. The phenomenological implications of our study include an updated exclusive determination of |Vcb| in the HQE, which is compatible with both the exclusive determination using the BGL parametrization and with the inclusive determination. We also revisit predictions for the lepton-flavour universality ratios RD(∗), the τ polarization observables PD(∗)τ, and the longitudinal polarization fraction FL. Posterior samples for the HQE parameters are provided as ancillary files, allowing for their use in subsequent studies.

The axion provides a compelling solution to the strong CP problem as well as a candidate for the dark matter of the universe. However, the axion solution relies on the spontaneous breaking of a global U(1)_PQ symmetry, which is also explicitly violated by quantum gravity. To preserve the axion solution, gravitational violations of the U(1)_PQ symmetry must be suppressed to sufficiently high order. We present a simple, geometric solution of the axion quality problem by modelling the axion with a bulk complex scalar field in a slice of AdS₅, where the U(1)_PQ symmetry is spontaneously broken in the bulk but explicitly broken on the UV brane. By localising the axion field towards the IR brane, gravitational violations of the PQ symmetry on the UV brane are sufficiently sequestered. This geometric solution is holographically dual to 4D strong dynamics where the global U(1)_PQ is an accidental symmetry to sufficiently high order.

We study UV-complete Froggatt-Nielsen-like models for the generation of mass and mixing hierarchies, assuming that the integrated heavy fields are chiral with respect to an abelian Froggatt-Nielsen symmetry. It modifies the mixed anomalies with respect to the Standard Model gauge group, which opens up the possibility to gauge the Froggatt-Nielsen symmetry without the need to introduce additional spectator fermions, while keeping mass matrices usually associated to anomalous flavour symmetries. We give specific examples where this happens, and we study the flavourful axion which arises from an accidental Peccei-Quinn symmetry in some of those models. Such an axion is typically more coupled to matter than in models with spectator fermions.

This document is the final report of the ATLAS-CMS Dark Matter Forum, a forum organized by the ATLAS and CMS collaborations with the participation of experts on theories of Dark Matter, to select a minimal basis set of dark matter simplified models that should support the design of the early LHC Run-2 searches. A prioritized, compact set of benchmark models is proposed, accompanied by studies of the parameter space of these models and a repository of generator implementations. This report also addresses how to apply the Effective Field Theory formalism for collider searches and present the results of such interpretations.

We present a detailed analysis of the Magellanic Bridge Cepheid sample constructed using the Optical Gravitational Lensing Experiment Collection of Variable Stars. Our updated Bridge sample contains 10 classical and 13 anomalous Cepheids. We calculate their individual distances using optical period-Wesenheit relations and construct three-dimensional maps. Classical Cepheid (CC) on-sky locations match very well neutral hydrogen and young stars distributions; thus, they add to the overall young Bridge population. In three dimensions, 8 out of 10 CCs form a bridge-like connection between the Magellanic Clouds. The other two are located slightly farther away and may constitute the Counter Bridge. We estimate ages of our Cepheids to be less than 300 Myr for from 5 up to 8 out of 10, depending on whether the rotation is included. This is in agreement with a scenario where these stars were formed in situ after the last encounter of the Magellanic Clouds. Cepheids' proper motions reveal that they are moving away from both Large and Small Magellanic Clouds. Anomalous Cepheids are more spread than CCs in both two and three dimensions, even though they form a rather smooth connection between the Magellanic Clouds. However, this connection does not seem to be bridge-like, as there are many outliers around both Magellanic Clouds.

We use the extended and updated Optical Gravitational Lensing Experiment (OGLE) Collection of Variable Stars to thoroughly analyze the distribution of RR Lyrae stars in the Magellanic Bridge. We use photometric metallicities to derive the absolute Wesenheit magnitude and individual distance of each RR Lyrae star. We confirm results from our earlier study showing that RR Lyrae stars are present in between the Magellanic Clouds, though their three-dimensional distribution more resembles two extended overlapping structures than a strict bridge-like connection. The contours do connect in the southern parts of the Bridge, albeit on a level too low to state that an evident connection exists. To test the sample numerically, we use multi-Gaussian fitting and conclude that there is no additional population or overdensity located in the Bridge. We also try to reproduce results on the putative RR Lyrae Magellanic Bridge stream by selecting RR Lyrae candidates from Gaia Data Release 1. We show that we are not able to obtain the evident connection of the Clouds without many spurious sources in the sample, as the cuts are not able to remove artifacts without eliminating the evident connection at the same time. Moreover, for the first time, we present the Gaia Data Release 2 RR Lyrae stars in the Magellanic Bridge area and show that their distribution matches our results.

We show that the leading coupling between a shift symmetric inflaton and the standard model fermions leads to an induced electroweak symmetry breaking due to particle production during inflation, and as a result, a unique oscillating feature in non-Gaussianities. In this one parameter model, the enhanced production of standard model fermions dynamically generates a new electroweak symmetry breaking minimum, where the Higgs field classically rolls to. The production of fermions stops when the Higgs expectation value and hence the fermion masses become too large, suppressing fermion production. The balance between the above-mentioned effects gives the standard model fermions masses that are uniquely determined by their couplings to the inflaton. In particular, the heaviest standard model fermion, the top quark, can produce a distinct cosmological collider physics signature characterized by a one-to-one relation between amplitude and frequency of the oscillating signal, which is observable at future 21-cm surveys.

We explore the possibility that lepton family numbers and baryon number are such good symmetries of Nature because they are the global remnant of a spontaneously broken gauge symmetry. An almost arbitrary linear combination of these symmetries (together with a component of global hypercharge) can be consistently gauged, if the Standard Model (SM) fermion content is augmented by three chiral SM singlet states. Within this framework of U (1 ) extensions of the SM one generically expects flavor nonuniversality to emerge in the charged leptons, in such a way that naturally prevents lepton flavor violation, by aligning the mass and weak eigenbases. For quarks, all the SM Yukawa couplings responsible for their observed masses and mixings arise at the renormalizable level. We perform fits to show that models in this class can explain R_K^(*) and the other neutral current B anomaly data if we introduce a heavy vectorlike quark to mediate the required quark flavor violation, while simultaneously satisfying other constraints from direct Z^' searches at the LHC, B_s meson mixing, a number of electroweak precision observables, and neutrino trident production. Within this symmetry-motivated framework of models, we find interesting implications for the flavor anomalies; notably, any axial couplings of the Z^' to electrons and muons must be flavor universal, with the flavor universality violation arising solely from the vectorlike couplings. We also comment on the generation of neutrino masses in these models.

Context. Planets in accretion disks can excite spiral shocks and if these planets are massive enough, they can even open gaps in their vicinity. Both of these effects can influence the overall thermal structure of the disk.
Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and to highlight the mutable, non-axisymmetric nature of ice lines in systems with massive planets.
Methods: We conducted a parameter study using numerical hydrodynamics simulations where we treated viscous heating, thermal cooling, and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters.
Results: Both gap opening and shock heating can displace the ice line, with the effects amplified for massive planets in optically thick disks. The gap region can split an initially hot (T > 170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric ice line into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context.
Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and, in general, they can move or reshape the water ice line. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones which lead to "islands" of vapor and ice around spirals that could affect the accretion or growth of icy aggregates.

We introduce a — somewhat holographic — dictionary between gravitational observables for scattering processes (measured at the boundary) and adiabatic invariants for bound orbits (in the bulk), to all orders in the Post-Minkowskian (PM) expansion. Our map relies on remarkable connections between the relative momentum of the twobody problem, the classical limit of the scattering amplitude and the deflection angle in hyperbolic motion. These relationships allow us to compute observables for generic orbits (such as the periastron advance ∆Φ) through analytic continuation, via a radial action depending only on boundary data. A simplified (more geometrical) map can be obtained for circular orbits, enabling us to extract the orbital frequency as a function of the (conserved) binding energy, Ω(E), directly from scattering information. As an example, using the results in Bernet al. [36, 37], we readily derive Ω(E) and ∆Φ(J, E) to two-loop orders. We also provide closed-form expressions for the orbital frequency and periastron advance at tree-level and one-loop order, respectively, which capture a series of exact terms in the Post-Newtonian expansion. We then perform a partial PM resummation, using a no-recoil approximation for the amplitude. This limit is behind the map between the scattering angle for a test-particle and the two-body dynamics to 2PM. We show that it also captures a subset of higher order terms beyond the test-particle limit. While a (rather lengthy) Hamiltonian may be derived as an intermediate step, our map applies directly between gauge invariant quantities. Our findings provide a starting point for an alternative approach to the binary problem. We conclude with future directions and some speculations on the classical double copy.

In models with extended scalars and C P violation, resonance searches in double Higgs final states stand in competition with related searches in top quark final states as optimal channels for the discovery of beyond the Standard Model (BSM) physics. This complementarity is particularly relevant for benchmark scenarios that aim to highlight multi-Higgs production as a standard candle for the study of BSM phenomena. In this paper, we compare interference effects in t t ¯ final states with correlated phenomena in double Higgs production in the complex singlet and the complex two-Higgs-doublet models. Our results indicate that the BSM discovery potential in di-Higgs searches can be underestimated in comparison to t t ¯ resonance searches. Top pair final states are typically suppressed due to destructive signal-background interference, while h h final states can be enhanced due to signal-signal interference. For parameter choices where the two heavy Higgs resonances are well separated in mass, top final states are suppressed relative to the naive signal expectation, while estimates of the production cross section times branching ratio remain accurate at the O (10 %) level for double Higgs final states.

Searching for new physics in large data sets needs a balance between two competing effects---signal identification vs background distortion. In this work, we perform a systematic study of both single variable and multivariate jet tagging methods that aim for this balance. The methods preserve the shape of the background distribution by either augmenting the training procedure or the data itself. Multiple quantitative metrics to compare the methods are considered, for tagging 2-, 3-, or 4-prong jets from the QCD background. This is the first study to show that the data augmentation techniques of Planing and PCA based scaling deliver similar performance as the augmented training techniques of Adversarial NN and uBoost, but are both easier to implement and computationally cheaper.

The detection of an oscillating pattern in the bispectrum of density perturbations could suggest the existence of a high-energy second minimum in the Higgs potential. If the Higgs field resided in this new minimum during inflation and was brought back to the electroweak vacuum by thermal corrections during reheating, the coupling of Standard Model particles to the inflaton would leave its imprint on the bispectrum. We focus on the fermions, whose dispersion relation can be modified by the coupling to the inflaton, leading to an enhanced particle production during inflation even if their mass during inflation is larger than the Hubble scale. This results in a large non-analytic contribution to non-Gaussianities, with an amplitude f_NL as large as 100 in the squeezed limit, potentially detectable in future 21-cm surveys. Measuring the contributions from two fermions would allow us to compute the ratio of their masses, and to ascribe the origin of the signal to a new Higgs minimum. Such a discovery would be a tremendous step towards understanding the vacuum instability of the Higgs potential, and could have fascinating implications for anthropic considerations.

Deformed relativistic kinematics, expected to emerge in a flat-spacetime limit of quantum gravity, predicts violation of discrete symmetries at energy scale in the vicinity of the Planck mass. Momentum-dependent deformations of the C, P and T invariance are derived from the \k{appa}-deformed Poincaré algebra. Deformation of the CPT symmetry leads to a subtle violation of Lorentz symmetry. This entails some small but measurable phenomenological consequences, as corrections to characteristics of time evolution: particle lifetimes or frequency of flavour oscillations in two-particle states at high energy. We argue here that using current experimental precisions on the muon lifetime one can bound the deformation parameter \k{appa} > 10^14 GeV at LHC energy and move this limit even to 10^16 GeV at Future Circular Collider, planned at CERN. Weaker limits on deformation can be also obtained from interference of neutral mesons. In case of B0s from {\Upsilon} decay it amounts to \k{appa} > 10^8 GeV at confidence level 99%.

The (tree) amplituhedron A(n,k,m) is the image in the Grassmannian Gr(k,k+m) of the totally nonnegative part of Gr(k,n), under a (map induced by a) linear map which is totally positive. It was introduced by Arkani-Hamed and Trnka in 2013 in order to give a geometric basis for the computation of scattering amplitudes in N=4 supersymmetric Yang-Mills theory. In the case relevant to physics (m=4), there is a collection of recursively-defined 4k-dimensional BCFW cells in the totally nonnegative part of Gr(k,n), whose images conjecturally "triangulate" the amplituhedron--that is, their images are disjoint and cover a dense subset of A(n,k,4). In this paper, we approach this problem by first giving an explicit (as opposed to recursive) description of the BCFW cells. We then develop sign-variational tools which we use to prove that when k=2, the images of these cells are disjoint in A(n,k,4). We also conjecture that for arbitrary even m, there is a decomposition of the amplituhedron A(n,k,m) involving precisely M(k, n-k-m, m/2) top-dimensional cells (of dimension km), where M(a,b,c) is the number of plane partitions contained in an a x b x c box. This agrees with the fact that when m=4, the number of BCFW cells is the Narayana number N(n-3, k+1).

We report the discovery of 10 kpc [C II] 158 μm halos surrounding star-forming galaxies in the early universe. We choose deep Atacama Large Millimeter/submillimeter Array data for 18 galaxies, each with a star formation rate of ≃10-70 M _⊙ with no signature of an active galactic nucleus whose [C II] lines are individually detected at z = 5.153-7.142, and we conduct stacking of the [C II] lines and dust continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10 kpc scale [C II] halo at the 9.2σ level, significantly more extended than the Hubble Space Telescope stellar continuum data by a factor of ∼5 on the exponential-profile basis, as well as the dust continuum. We compare the radial profiles of [C II] and Lyα halos universally found in star-forming galaxies at this epoch, and we find that the scale lengths agree within the 1σ level. While two independent hydrodynamic zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [C II] line emission. The existence of the extended [C II] halo is evidence of outflow remnants in the early galaxies and suggests that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function slopes, with the goal of connecting the statistical properties of neutral hydrogen (H I) with the turbulent drivers. The new methods were applied to the H I observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties of H I, with no evidence for local enhancements of turbulence due to stellar feedback. These properties could be caused by a significant turbulent driving on large scales. Alternatively, the significant line-of-sight depth of the SMC could be masking out localized regions with a steeper SPS slope caused by stellar feedback. In contrast to the SMC, the LMC H I shows a large diversity in terms of its turbulent properties. Across most of the LMC, the small-scale SPS slope is steeper than the large-scale slope due to the presence of the H I disk. On small spatial scales, we find several areas of localized steepening of the SPS slope around major H II regions, with the 30 Doradus region being the most prominent. This is in agreement with predictions from numerical simulations, which suggest a steepening of the SPS slope due to stellar feedback that erodes and destroys interstellar clouds. We also find a localized steepening of the large-scale SPS slope in the outskirts of the LMC. This is likely caused by the flaring of the H I disk, or alternatively, by ram-pressure stripping of the LMC disk due to the interactions with the surrounding halo gas.

We investigate the [C II] line intensity mapping (IM) signal from galaxies in the Epoch of Reionization (EoR) to assess its detectability, the possibility to constrain the L_{C II}-SFR relation, and to recover the [C II] luminosity function (LF) from future experiments. By empirically assuming that log L_{C II}=log A+γ SFR± σ _ L, we derive the [C II] LF from the observed UV LF, and the [C II] IM power spectrum. We study the shot noise and the full power spectrum separately. Although, in general, the shot-noise component has a much higher signal-to-noise ratio than the clustering one, it cannot be used to put independent constraints on log A and γ. Full power spectrum measurements are crucial to break such degeneracy and reconstruct the [C II] LF. In our fiducial survey S1 (inspired by CCAT-p/1000 h) at z ∼ 6, the shot-noise (clustering) signal is detectable for two (one) of the five considered L_{C II}-SFR relations. The shot noise is generally dominated by galaxies with L_{C II}≳ 10⁸-10⁹ L_⊙ (M_UV ∼ -20 to -22), already at reach of ALMA pointed observations. However, given the small field of view of such telescope, an IM experiment would provide unique information on the bright end of the LF. The detection depth of an IM experiment crucially depends on the (poorly constrained) L_{C II}-SFR relation in the EoR. If the L_{C II}-SFR relation varies in a wide log A-γ range, but still consistent with ALMA [C II] LF upper limits, even the signal from galaxies with L_{C II} as faint as ∼10⁷ L_⊙ could be detectable. Finally, we consider the contamination by continuum foregrounds (cosmic infrared background, dust, cosmic microwave background) and CO interloping lines, and derive the requirements on the residual contamination level to reliably extract the [C II] signal.

Dark matter substructure can contribute significantly to local dark matter searches and may provide a large uncertainty in the interpretation of those experiments. For direct detection experiments, sub-halos give rise to an additional dark matter component on top of the smooth dark matter distribution of the host halo. In the case of dark matter capture in the Sun, sub-halo encounters temporarily increase the number of captured particles. Even if the encounter happened in the past, the number of dark matter particles captured by the Sun can still be enhanced today compared to expectations from the host halo as those enhancements decay over time. Using results from an analytical model of the sub-halo population of a Milky Way-like galaxy, valid for sub-halo masses between 10^-5 M_solar and 10¹¹ M_solar, we assess the impact of sub-halos on direct dark matter searches in a probabilistic way. We find that the impact on direct detection can be sizable, with a probability of ~ 10^-3 to find an Script O(1) enhancement of the recoil rate. In the case of the capture rate in the Sun, we find that Script O(1) enhancements are very unlikely, with probability lesssim 10^-5, and are even impossible for some dark matter masses.

We present two determinations of the strong coupling α_s. The first one is from the static energy at three-loop accuracy, and may be considered an update of earlier determinations by some of us. The new analysis includes new lattice data at smaller lattice spacings, and reaches distances as short as 0.0237 fm. We present a comprehensive and detailed estimate of the error sources that contribute to the uncertainty of the final result, α_s(M_Z)=0.1166 0_-0.00056^+0.00110. The second determination is based on lattice data for the singlet free energy at finite temperature up to distances as small as 0.0081 fm, from which we obtain α_s(M_Z)=0.1163 8_-0.00087^{+0.0009 5}.

We calculate the step scaling function, the lattice analog of the renormalization group β -function, for an SU(3) gauge theory with twelve flavors. The gauge coupling of this system runs very slowly, which is reflected in a small step scaling function, making numerical simulations particularly challenging. We present a detailed analysis including the study of systematic effects of our extensive data set generated with twelve dynamical flavors using the Symanzik gauge action and three times stout smeared Möbius domain wall fermions. Using up to 32⁴ volumes, we calculate renormalized couplings for different gradient flow schemes and determine the step-scaling β function for a scale change s =2 on up to five different lattice volume pairs. Our preferred analysis is fully O (a²) Symanzik improved and uses Zeuthen flow combined with the Symanzik operator. We find an infrared fixed point within the range 5.2 ≤g_c²≤6.4 in the c =0.250 finite volume gradient flow scheme. We account for systematic effects by calculating the step-scaling function based on alternative flows (Wilson or Symanzik) as well as operators (Wilson plaquette, clover) and also explore the effects of the perturbative tree-level improvement.

The detection of the high-energy neutrino event, IceCube-170922A, demonstrated that multimessenger particle astrophysics triggered by neutrino alerts is feasible. We consider time delay signatures caused by secret neutrino interactions with the cosmic neutrino background and dark matter and suggest that these can be used as a novel probe of neutrino interactions beyond the standard model (BSM). The tests with BSM-induced neutrino echoes are distinct from existing constraints from the spectral modification and will be enabled by multimessenger observations of bright neutrino transients with future experiments such as IceCube-Gen2, KM3Net, and Hyper-Kamiokande. The constraints are complementary to those from accelerator and laboratory experiments and powerful for testing various particle models that explain tensions prevailing in the cosmological data.

The existence of millicharged dark matter (mDM) can leave a measurable imprint on 21-cm cosmology through mDM-baryon scattering. However, the minimal scenario is severely constrained by existing cosmological bounds on both the fraction of dark matter that can be millicharged and the mass of mDM particles. We point out that introducing a long-range force between a millicharged subcomponent of dark matter and the dominant cold dark matter (CDM) component leads to efficient cooling of baryons in the early Universe, while also significantly extending the range of viable mDM masses. Such a scenario can explain the anomalous absorption signal in the sky-averaged 21-cm spectrum observed by EDGES and leads to a number of testable predictions for the properties of the dark sector. The mDM mass can then lie between 10 MeV and a few hundreds of GeVs, and its scattering cross section with baryons lies within an unconstrained window of parameter space above direct detection limits and below current bounds from colliders. In this allowed region, mDM can make up as little as 10^-8 of the total dark matter energy density. The CDM mass ranges from 10 MeV to a few GeVs and has an interaction cross section with the Standard Model that is induced by a loop of mDM particles. This cross section is generically within reach of near-future low-threshold direct detection experiments.

Physics beyond the Standard Model can manifest itself as both new light states and heavy degrees of freedom. In this paper, we assume that the former comprise only a sterile neutrino, N . Therefore, the most agnostic description of the new physics is given by an effective field theory built upon the Standard Model fields as well as N . We show that Higgs phenomenology provides a sensitive and potentially crucial tool to constrain effective gauge interactions of sterile neutrinos, not yet probed by current experiments. In parallel, this motivates a range of new Higgs decay channels with clean signatures as candidates for the next LHC runs, including h →γ +p_T^miss and h →γ γ +p_T^miss .

We study charged Dirac quasinormal modes (QNMs) on Reissner-Nordström-anti-de Sitter (RN-AdS) black holes with generic Robin boundary conditions, by extending our earlier work of neutral Dirac QNMs on Schwarzschild-AdS black holes. We first derive the equations of motion for charged Dirac fields on a RN-AdS background. To solve these equations we impose a requirement on the Dirac field: that its energy flux should vanish at asymptotic infinity. A set of two Robin boundary conditions compatible with QNMs is consequently found. By employing both analytic and numeric methods, we then obtain the quasinormal spectrum for charged Dirac fields and analyze the impact of various parameters, in particular of electric charges. An analytic calculation shows explicitly that the charge coupling between the black hole and the Dirac field does not trigger super-radiant instabilities in the small black hole and low frequency limit. Numeric calculations, on the other hand, show quantitatively that Dirac QNMs may change substantially due to the electric charge. Our results illustrate how vanishing energy flux boundary conditions, as a generic principle, are applicable not only to neutral but also to electrically charged fields.

Consistency relations for the large scale structure are exact equalities between correlation functions of different order. These relations descend from the equivalence principle and hold for primordial perturbations generated by single-field models of inflation. They are not affected by nonlinearities and hold also for biased tracers and in redshift space. We show that baryonic acoustic oscillations in the bispectrum (BS) in the squeezed limit are suppressed with respect to those in the power spectrum by a coefficient that depends on the BS configuration and on the bias parameter (and, in redshift space, also on the growth rate). We test these relations using large volume N -body simulations and show that they provide a novel way to measure large scale halo bias and, potentially, the growth rate. Since bias is obtained by comparing two directly observable quantities, the method is free from theoretical uncertainties both on the computational scheme and on the underlying cosmological model.

By means of 3D hydrodynamical simulations, we evaluate here the impact that supernova (SN) explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both the dust generated within the ejecta and the dust initially present in the ambient gas and later locked up in the surrounding wind-driven shell (WDS). The collision of the SN blast wave with the WDS leads to a transmitted shock that moves into the shell and a reflected shock that moves into the ejecta. The transmitted shock is capable of destroying large amounts of the dust locked in the shell, but only if the mass of the WDS is small, less than a few tens the ejected mass. Conversely, massive WDSs, with several times the ejected mass, lead upon the interaction to strong radiative cooling, which inhibits the Sedov-Taylor phase and weakens the transmitted shock, making it unable to traverse the WDS. In such a case, the destruction/disruption of the ambient dust is largely inhibited. On the other hand, the SN remnants grow rapidly in the very tenuous region excavated by the stellar winds, and thus a large fraction of the dust generated within the ejecta is not efficiently destroyed by the SN reverse shock, nor by the reflected shock. Our calculations favor a scenario in which core-collapse SNe within sufficiently massive WDSs supply more dust to the interstellar medium than they are able to destroy.

We introduce k-evolution, a relativistic N-body code based on gevolution, which includes clustering dark energy among its cosmological components. To describe dark energy, we use the effective field theory approach. In particular, we focus on k-essence with a speed of sound much smaller than unity but we lay down the basis to extend the code to other dark energy and modified gravity models. We develop the formalism including dark energy non-linearities but, as a first step, we implement the equations in the code after dropping non-linear self-coupling in the k-essence field. In this simplified setup, we compare k-evolution simulations with those of CLASS and gevolution 1.2, showing the effect of dark matter and gravitational non-linearities on the power spectrum of dark matter, of dark energy and of the gravitational potential. Moreover, we compare k-evolution to Newtonian N-body simulations with back-scaled initial conditions and study how dark energy clustering affects massive halos.

We introduce a new technique to estimate the comet nuclear size frequency distribution (SFD) that combines a cometary activity model with a survey simulation and apply it to 150 long period comets (LPC) detected by the Pan-STARRS1 near-Earth object survey. The debiased LPC size-frequency distribution is in agreement with previous estimates for large comets with nuclear diameter ≳1 km but we measure a significant drop in the SFD slope for small objects with diameters <1 km and approaching only 100 m diameter. Large objects have a slope α_big = 0.72 ± 0.09(stat.) ± 0.15(sys.) while small objects behave as α_small = 0.07 ± 0.03(stat.) ± 0.09(sys.) where the SFD is ∝ 1^{0 αH_N} and H_N represents the cometary nuclear absolute magnitude. The total number of LPCs that are >1 km diameter and have perihelia q < 10 au is 0.46 ± 0.15 × 10⁹ while there are only 2.4 ± 0.5(stat.) ± 2(sys.) × 10⁹ objects with diameters >100 m due to the shallow slope of the SFD for diameters <1 km. We estimate that the total number of 'potentially active' objects with diameters ≥1 km in the Oort cloud, objects that would be defined as LPCs if their perihelia evolved to <10 au, is (1.5 ± 1) × 10¹² with a combined mass of 1.3 ± 0.9 M_⊕. The debiased LPC orbit distribution is broadly in agreement with expectations from contemporary dynamical models but there are discrepancies that could point towards a future ability to disentangle the relative importance of stellar perturbations and galactic tides in producing the LPC population.

We study the role of anisotropic escape in generating the elliptic flow of bottomonia produced in ultrarelativistic heavy-ion collisions. We implement temperature-dependent decay widths for the various bottomonium states to calculate their survival probability when traversing through the anisotropic hot medium formed in noncentral collisions. We employ the recently developed 3 +1 -dimensional quasiparticle anisotropic hydrodynamic simulation to model the space-time evolution of the quark-gluon plasma. We provide a quantitative prediction for the transverse momentum dependence of bottomonium elliptic flow and the nuclear modification factor for Pb +Pb collisions in √{s_NN}=2.76 TeV at the CERN Large Hadron Collider.

Confining hidden sectors are an attractive possibility for physics beyond the Standard Model (SM). They are especially motivated by neutral naturalness theories, which reconcile the lightness of the Higgs with the strong constraints on colored top partners. We study hidden QCD with one light quark flavor, coupled to the SM via effective operators suppressed by the mass M of new electroweak-charged particles. This effective field theory is inspired by a new tripled top model of supersymmetric neutral naturalness. The hidden sector is accessed primarily via the Z and Higgs portals, which also mediate the decays of the hidden mesons back to SM particles. We find that exotic Z decays at the LHC and future Z factories provide the strongest sensitivity to this scenario, and we outline a wide array of searches. For a larger hidden confinement scale Λ ∼ O (10) GeV, the exotic Z decays dominantly produce final states with two hidden mesons. ATLAS and CMS can probe their prompt decays up to M ∼ 3 TeV at the high luminosity phase, while a TeraZ factory would extend the reach up to M ∼ 20 TeV through a combination of searches for prompt and displaced signals. For smaller Λ ∼ O (1) GeV, the Z decays to the hidden sector produce jets of hidden mesons, which are long-lived. LHCb will be a powerful probe of these emerging jets. Furthermore, the light hidden vector meson could be detected by proposed dark photon searches.

We discuss how LHC di-muon data collected to study $B_q \to \mu\mu$ can be used to constrain light particles with flavour-violating couplings to $b$-quarks. Focussing on the case of a flavoured QCD axion, $a$, we compute the decay rates for $B_q \to \mu \mu a$ and the SM background process $B_q \to \mu \mu \gamma$ near the kinematic endpoint. These rates depend on non-perturbative $B_q \to \gamma^{(*)}$ form factors with on- or off-shell photons. The off-shell form factors -- relevant for generic searches for beyond-the-SM particles -- are discussed in full generality and computed with QCD sum rules for the first time. With these results, we analyse available LHCb data to obtain the sensitivity on $B_q \to \mu \mu a$ at present and future runs. We find that the full LHCb dataset alone will allow to probe axion-coupling scales of the order of $10^6$ GeV for both $b\to d$ and $b \to s$ transitions.

The gradient flow transformation can be interpreted as continuous real-space renormalization group transformation if a coarse-graining step is incorporated as part of calculating expectation values. The method allows to predict critical properties of strongly coupled systems including the renormalization group $\beta$ function and anomalous dimensions at nonperturbative fixed points. In this contribution we discuss a new analysis of the continuous renormalization group $\beta$ function for $N_f=2$ and $N_f=12$ fundamental flavors in SU(3) gauge theories based on this method. We follow the approach developed and tested for the $N_f=2$ system in arXiv:1910.06408. Here we present further information on the analysis, emphasizing the robustness and intuitive features of the continuous $\beta$ function calculation. We also discuss the applicability of the continuous $\beta$ function calculation in conformal systems, extending the possible phase diagram to include a 4-fermion interaction. The numerical analysis for $N_f=12$ uses the same set of ensembles that was generated and analyzed for the step scaling function in arXiv:1909.05842. The new analysis uses volumes with $L \ge 20$ and determines the $\beta$ function in the $c=0$ gradient flow renormalization scheme. The continuous $\beta$ function predicts the existence of a conformal fixed point and is consistent between different operators. Although determinations of the step scaling and continuous $\beta$ function use different renormalization schemes, they both predict the existence of a conformal fixed point around $g^2\sim 6$.

We present a Markov-Chain Monte-Carlo (MCMC) forecast for the precision of neutrino mass and cosmological parameter measurements with a Euclid-like galaxy clustering survey. We use a complete perturbation theory model for the galaxy one-loop power spectrum and tree-level bispectrum, which includes bias, redshift space distortions, IR resummation for baryon acoustic oscillations and UV counterterms. The latter encapsulate various effects of short-scale dynamics which cannot be modeled within perturbation theory. Our MCMC procedure consistently computes the non-linear power spectra and bispectra as we scan over different cosmologies. The second ingredient of our approach is the theoretical error covariance which captures uncertainties due to higher-order non-linear corrections omitted in our model. Having specified characteristics of a Euclid-like spectroscopic survey, we generate and fit mock galaxy power spectrum and bispectrum likelihoods. Our results suggest that even under very agnostic assumptions about non-linearities and short-scale physics a future Euclid-like survey will be able to measure the sum of neutrino masses with a standard deviation of 28 meV . When combined with the Planck cosmic microwave background likelihood, this uncertainty decreases to 13 meV . Over-optimistically reducing the theoretical error on the bispectrum down to the two-loop level marginally tightens this bound to 11 meV . Moreover, we show that the future large-scale structure (LSS) spectroscopic data will greatly improve constraints on the other cosmological parameters, e.g. reaching a percent (per mille) error on the Hubble constant with LSS alone (LSS + Planck).

We propose a mechanism to solve the Higgs naturalness problem through a cosmological selection process. The discharging of excited field configurations through membrane nucleation leads to discrete jumps of the cosmological constant and the Higgs mass, which vary in a correlated way. The resulting multitude of universes are all empty, except for those in which the cosmological constant and the Higgs mass are both nearly vanishing. Only under these critical conditions can inflation be activated and create a non-empty universe.

A tight relation between the [C II] 158 μm line luminosity and star formation rate is measured in local galaxies. At high redshift (z > 5), though, a much larger scatter is observed, with a considerable (15-20 per cent) fraction of the outliers being [C II]-deficient. Moreover, the [C II] surface brightness (Σ_[C II]) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [C II]-deficiency, we have developed an analytical model that fits local [C II] data and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of Σ_[C II] with Σ_SFR. However, for Σ_SFR {≳} 1 M_⊙ yr^{-1} kpc^{-2}, Σ_[C II] saturates. We conclude that underluminous [C II] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation (κ_s ≫ 1), parametrized by the `burstiness' parameter κ_s; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [C II] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [O III]88 μm, C III]1909 Å, CO lines). Simple formulae are given to interpret available data for low- and high-z galaxies.

Two-fluid (electron-positron) plasma modeling has shown that inductive acceleration can convert Poynting flux directly into bulk kinetic energy in the relativistic flows driven by rotating magnetized neutron stars and black holes. Here, we generalize this approach by adding an ion fluid. Solutions are presented in which all particles are accelerated as the flow expands, with comparable power channeled into each of the plasma components. In an ion-dominated flow, each species reaches the limiting rigidity, according to Hillas’ criterion, in a distance significantly shorter than in a lepton-dominated flow. These solutions support the hypothesis that newly born magnetars and pulsars are potential sources of ultrahigh energy cosmic rays. The competing process of Poynting flux dissipation by magnetic reconnection is shown to be ineffective in low-density flows in which the conventionally defined electron multiplicity satisfies {κ }_{{e}}≲ {10}⁵{≤ft(4π {L}₃₈/{{Ω }}\right)}^1/4/{{Max}}≤ft({η }_ion}^1/2,1\right), where L ₃₈ × 10³⁸ erg s^-1 is the power carried by the flow in a solid angle Ω, and {η }_ion} is the ratio of the ion to lepton power at launch.

Recent high-resolution interferometric observations of protoplanetary disks at (sub)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process. We investigated if ongoing planetesimal formation could explain the “fine-tuned” optical depths in the DSHARP rings by removing dust and transforming it into “invisible” planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations. We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.

We present Period-Luminosity and Period-Luminosity-Color relations at maximum light for Mira variables in the Magellanic Clouds using time-series data from the Optical Gravitational Lensing Experiment (OGLE-III) and Gaia data release 2. The maximum-light relations exhibit a scatter typically up to ∼30% smaller than their mean-light counterparts. The apparent magnitudes of oxygen-rich Miras at maximum light display significantly smaller cycle-to-cycle variations than at minimum light. High-precision photometric data for Kepler Mira candidates also exhibit stable magnitude variations at the brightest epochs, while their multi-epoch spectra display strong Balmer emission lines and weak molecular absorption at maximum light. The stability of maximum-light magnitudes for Miras possibly occurs due to the decrease in the sensitivity to molecular bands at their warmest phase. At near-infrared wavelengths, the period-luminosity relations (PLRs) of Miras display similar dispersion at mean and maximum light with limited time-series data in the Magellanic Clouds. A kink in the oxygen-rich Mira PLRs is found at 300 days in the VI-bands, which shifts to longer periods (∼350 days) at near-infrared wavelengths. Oxygen-rich Mira PLRs at maximum light provide a relative distance modulus, Δμ = 0.48 ± 0.08 mag, between the Magellanic Clouds with a smaller statistical uncertainty than the mean-light relations. The maximum-light properties of Miras can be very useful for stellar atmosphere modeling and distance scale studies provided their stability and the universality can be established in other stellar environments in the era of extremely large telescopes.

Neutron star mergers can form a hypermassive neutron star (HMNS) remnant, which may be the engine of a short gamma-ray burst (SGRB) before it collapses to a black hole, possibly several hundred milliseconds after the merger. During the lifetime of an HMNS, numerical relativity simulations indicate that it will undergo strong oscillations and emit gravitational waves with frequencies of a few kilohertz, which are unfortunately too high for detection to be probable with the Advanced Laser Interferometer Gravitational-Wave Observatory. Here we discuss the current and future prospects for detecting these frequencies as modulation of the SGRB. The understanding of the physical mechanism responsible for the HMNS oscillations will provide information on the equation of state of the hot HMNS, and the observation of these frequencies in the SGRB data would give us insight into the emission mechanism of the SGRB.

We present results from a comparative study of light curves of Cepheid and RR Lyrae stars in the Galaxy and the Magellanic Clouds with their theoretical models generated from the stellar pulsation codes. Fourier decomposition method is used to analyse the theoretical and the observed light curves at multiple wavelengths. In case of RR Lyrae stars, the amplitude and Fourier parameters from the models are consistent with observations in most period bins except for low metal-abundances (Z < 0:004). In case of Cepheid variables, we observe a greater offset between models and observations for both the amplitude and Fourier parameters. The theoretical amplitude parameters are typically larger than those from observations, except close to the period of 10 days. We find that these discrepancies between models and observations can be reduced if a higher convective efficiency is adopted in the pulsation codes. Our results suggest that a quantitative comparison of light curve structure is very useful to provide constraints for the input physics to the stellar pulsation models.

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimetre broad-band flux to total infrared luminosity (L_IR), and hence star formation rate, in high-redshift galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g. redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (M_star > 10^{10} M_{\odot }) z = 2-6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. At z ∼ 2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature, which is required to accurately estimate the total dust mass, does not strongly evolve with redshift over z = 2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. We also analyse an `equivalent' dust temperature for converting (sub)millimetre flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature (`warmer dust') do not necessarily have higher mass-weighted temperatures. A `two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures.

We study for the first time the collider reach on the derivative Higgs portal, the leading effective interaction that couples a pseudo Nambu-Goldstone boson (pNGB) scalar Dark Matter to the Standard Model. We focus on Dark Matter pair production through an off-shell Higgs boson, which is analyzed in the vector boson fusion channel. A variety of future high-energy lepton colliders as well as hadron colliders are considered, including CLIC, a muon collider, the High-Luminosity and High-Energy versions of the LHC, and FCC-hh. Implications on the parameter space of pNGB Dark Matter are discussed. In addition, we give improved and extended results for the collider reach on the marginal Higgs portal, under the assumption that the new scalars escape the detector, as motivated by a variety of beyond the Standard Model scenarios.

Atmospheric neutrinos produced by cosmic-ray interactions around the globe provide a beam for the study of neutrino properties. They are also a background in searches for neutrinos of astrophysical origin. Both aspects are addressed in this chapter, which begins with a brief introduction on neutrino oscillations in relation to the spectrum of atmospheric neutrinos. Section 2 describes the cascade equation for hadrons in the atmosphere and the main features of atmospheric leptons from their decays. Next, uncertainties in the fluxes that arise from limited knowledge of the primary spectrum and of particle production are discussed. The final section covers aspects specific to neutrino telescopes.

We describe the minimal space of polylogarithmic functions that is required to express the six-particle amplitude in planar N = 4 super-Yang-Mills theory through six and seven loops, in the NMHV and MHV sectors respectively. This space respects a set of extended Steinmann relations that restrict the iterated discontinuity structure of the amplitude, as well as a cosmic Galois coaction principle that constrains the functions and the transcendental numbers that can appear in the amplitude at special kinematic points. To put the amplitude into this space, we must divide it by the BDS-like ansatz and by an additional zeta-valued constant ρ. For this normalization, we conjecture that the extended Steinmann relations and the coaction principle hold to all orders in the coupling. We describe an iterative algorithm for constructing the space of hexagon functions that respects both constraints. We highlight further simplifications that begin to occur in this space of functions at weight eight, and distill the implications of imposing the coaction principle to all orders. Finally, we explore the restricted spaces of transcendental functions and constants that appear in special kinematic configurations, which include polylogarithms involving square, cube, fourth and sixth roots of unity.

Using the latest LHC data, we analyse and compare the lower limits on the masses of gluinos and the lightest stop in two natural supersymmetric motivated scenarios: one with a neutralino being the lightest supersymmetric particle (LSP) and the other one with gravitino as the LSP and neutralino as the next-to-lightest supersymmetric particle. In the second case our analysis applies to neutralinos promptly decaying to very light gravitinos, which are of cosmological interest, and are generic for low, of order O (100) TeV, messenger scale in gauge mediation models. We find that the lower bounds on the gluino and the lightest stop masses are stronger for the gravitino LSP scenarios due to the extra handle from the decay products of neutralinos. Generally, in contrast to the neutralino LSP case the limits now extend to a region of compressed spectrum. In bino scenarios the highest excluded stop mass increases from 1000 GeV to almost 1400 GeV. Additionally, in the higgsino-like NLSP scenario the higgsinos below 650 GeV are universally excluded and the stop mass limit is {m}_{\tilde{t}} > 1150 GeV, whereas there is no limit on stops in the higgsino LSP model for {m}_{\tilde{h}} = 650 GeV. Nevertheless, we find that the low messenger scale still ameliorates the fine tuning in the electroweak potential.

Inverse Compton-pair cascades are initiated when gamma-rays are absorbed on an ambient soft photon field to produce relativistic pairs, which in turn up-scatter the same soft photons to produce more gamma-rays. If the Compton scatterings take place in the deep Klein-Nishina regime, then triplet pair production (e{γ }_b\to {{ee}}⁺{e}^-) becomes relevant and may even regulate the development of the cascade. We investigate the properties of pair-Compton cascades with triplet pair production in accelerating gaps, i.e., regions with an unscreened electric field. Using the method of transport equations for the particle evolution, we compute the growth rate of the pair cascade as a function of the accelerating electric field in the presence of blackbody and power-law ambient photon fields. Informed by the numerical results, we derive simple analytical expressions for the peak growth rate and the corresponding electric field. We show that for certain parameters, which can be realized in the vicinity of accreting supermassive black holes at the centers of active galactic nuclei, the pair cascade may well be regulated by inverse Compton scattering in the deep Klein-Nishina regime and triplet pair production. We present indicative examples of the escaping gamma-ray radiation from the gap, and discuss our results in application to the TeV observations of radio galaxy M87.

Rare B-decays induced by flavour-changing neutral currents (FCNC) is one of the promising candidates for probing physics beyond the Standard model. However, for identifying potential new physics from the data, reliable control over QCD contributions is necessary. We focus on one of such QCD contributions - the charming loops - that potentially can lead to difficulties in disentangling new physics effects from the observable and discuss the possibility to gain control over theoretical predictions for charming loops.

Dark matter evolution during the process of cosmological structure formation can be described in terms of a one-particle irreducible effective action at a characteristic scale k_m and a loop expansion below this scale, based on the effective propagators and vertices. We calculate the form of the effective vertices and compute the bispectrum of density perturbations within a one-loop approximation. We find that the effective vertices play a subdominant role as compared to the effective viscosity and sound velocity that modify the (inverse) propagators. For the bispectrum we reproduce the results of standard perturbation theory in the range where it is applicable, and find a slightly improved agreement with N-body simulations at larger wavenumbers.

In the simple Higgs-portal dark matter model with a conserved dark matter number, we show that there exists a non-topological soliton state of dark matter. This state has smaller energy per dark matter number than a free particle state and has its interior in the electroweak symmetric vacuum. It could be produced in the early universe from first-order electroweak phase transition and contribute most of dark matter. This electroweak symmetric dark matter ball is a novel macroscopic dark matter candidate with an energy density of the electroweak scale and a mass of 1 gram or above. Because of its electroweak-symmetric interior, the dark matter ball has a large geometric scattering cross section off a nucleon or a nucleus. Dark matter and neutrino experiments with a large-size detector like Xenon1T, BOREXINO and JUNO have great potential to discover electroweak symmetric dark matter balls. We also discuss the formation of bound states of a dark matter ball and ordinary matter.

We study the soft limit of a recently proposed generalization of the biadjoint scalar amplitudes $m^{(k)}_{n}$, which have been conjectured to have a relation to the tropical Grassmannian $\text{Tr G}(k,n)$. Using the CHY formulation along with the Global Residue Theorem, we prove the soft factorization for $m^{(k)}_{n}$ amplitudes for arbitrary $k$ and $n$. We find that the soft factors are in direct correspondence to vertices of the associahedron $\mathcal{A}_{k-1}$, and hence take the form of $m^{(2)}_{n}$ amplitudes. This entails that all scattering amplitudes of the ordinary biadjoint scalar theory can be interpreted as an infinite family of soft factors. Additionally, Grassmannian duality reveals that generalized amplitudes $m^{(k)}_{n}$ with $k>2$ satisfy not only a soft theorem, but also a non-trivial "hard" theorem. We perform numerical checks of our theorems against previous results for $\text{Tr G}(4,7)$ and $\text{Tr G}(5,8)$, thereby providing strong evidence of their relation with the CHY formulation.

We present a method to flexibly and self-consistently determine individual galaxies' star formation rates (SFRs) from their host haloes' potential well depths, assembly histories, and redshifts. The method is constrained by galaxies' observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, ultraviolet (UV) luminosity functions, UV-stellar mass relations, IRX-UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0 < z < 10. Key findings include the following: galaxy assembly correlates strongly with halo assembly; quenching correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies' star formation histories at fixed mass differ most at z < 0.5; satellites have large scatter in quenching time-scales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar - halo mass ratios at early times; and, none the less, stellar mass-halo mass ratios rise at z > 5. Also presented are revised stellar mass - halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching time-scale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>10⁵ cores) implementation (the UNIVERSEMACHINE), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.

We present a simple and well defined prescription to compare absorption lines in supernova (SN) spectra with lists of transitions drawn from the National Institute of Standards and Technology database. The method is designed to be applicable to simple spectra where the photosphere can be mostly described by absorptions from single transitions with a single photospheric velocity. These conditions are plausible for SN spectra obtained shortly after explosion. Here we show that the method also works well for spectra of hydrogen-poor (Type I) superluminous supernovae (SLSNe-I) around peak. Analysis of high signal to noise spectra leads to clear identification of numerous spectroscopic features arising from ions of carbon and oxygen, which account for the majority of absorption features detected in the optical range, suggesting the outer envelope of SLSN-I progenitors is dominated by these elements. We find that the prominent absorption features seen in the blue are dominated by numerous lines of O II, as previously suggested, and that the apparent absorption feature widths are dominated by line density and not by Doppler broadening. In fact, we find that while the expansion velocities of SLSNe-I around peak are similar to those of normal SNe, the apparent velocity distribution (manifested as the width of single transition features) is much lower (∼1500 km s^-1) indicating emission from a very narrow photosphere in velocity space that is nevertheless expanding rapidly. We inspect the controversial case of ASASSN-15lh, and find that the early spectrum of this object is not consistent with those of SLSNe-I. We also show that SLSNe that initially lack hydrogen features but develop these at late phases, such as iPTF15esb and iPTF16bad, also differ in their early spectra from standard SLSNe-I.

This is the first installment of a series of three papers in which we describe a method to determine higher-point correlation functions in one-loop open-superstring amplitudes from first principles. In this first part, we exploit the synergy between the co-homological features of pure-spinor superspace and the pure-spinor zero-mode integration rules of the one-loop amplitude prescription. This leads to the study of a rich variety of multiparticle superfields which are local, have covariant BRST variations, and are compatible with the particularities of the pure-spinor amplitude prescription. Several objects related to these superfields, such as their non-local counterparts and the so-called BRST pseudo-invariants, are thoroughly reviewed and put into new light. Their properties will turn out to be mysteriously connected to products of one-loop worldsheet functions in packages dubbed "generalized elliptic integrands", whose prominence will be seen in the later parts of this series of papers.

This is the second installment of a series of three papers in which we describe a method to determine higher-point correlation functions in one-loop open-superstring amplitudes from first principles. In this second part, we study worldsheet functions defined on a genus-one surface built from the coefficient functions of the Kronecker-Einsenstein series. We construct two classes of worldsheet functions whose properties lead to several simplifying features within our description of one-loop correlators with the pure-spinor formalism. The first class is described by functions with prescribed monodromies, whose characteristic shuffle-symmetry property leads to a Lie-polynomial structure when multiplied by the local superfields from part I of this series. The second class is given by so-called generalized elliptic integrands (GEIs) that are constructed using the same combinatorial patterns of the BRST pseudo-invariant superfields from part I. Both of them lead to compact and combinatorially rich expressions for the correlators in part III. The identities obeyed by the two classes of worldsheet functions exhibit striking parallels with those of the superfield kinematics. We will refer to this phenomenon as a duality between worldsheet functions and kinematics.

In this final part of a series of three papers, we will assemble supersymmetric expressions for one-loop correlators in pure-spinor superspace that are BRST invariant, local, and single valued. A key driving force in this construction is the generalization of a so far unnoticed property at tree-level; the correlators have the symmetry structure akin to Lie polynomials. One-loop correlators up to seven points are presented in a variety of representations manifesting different subsets of their defining properties. These expressions are related via identities obeyed by the kinematic superfields and worldsheet functions spelled out in the first two parts of this series and reflecting a duality between the two kinds of ingredients. Interestingly, the expression for the eight-point correlator following from our method seems to capture correctly all the dependence on the worldsheet punctures but leaves undetermined the coefficient of the holomorphic Eisenstein series G₄. By virtue of chiral splitting, closed-string correlators follow from the double copy of the open-string results.

We study the formation and evolution of a sample of Lyman break galaxies in the epoch of reionization by using high-resolution (∼10 pc), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium thermochemical non-equilibrium evolution and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([C II], [N II], and [O III]). At z = 8, the most massive galaxy, `Freesia', has an age t_\star ∼eq 409 Myr, stellar mass M_⋆ ≃ 4.2 × 10⁹M_⊙, and a star formation rate (SFR), SFR∼eq 11.5 M_{⊙ } yr^{-1}, due to a recent burst. Freesia has two stellar components (A and B) separated by ≃ 2.5 kpc; other 11 galaxies are found within 56.9 ± 21.6 kpc. The mean ISRF in the Habing band is G = 7.9 G_0 and is spatially uniform; in contrast, the ionization parameter is U = 2^{+20}_{-2} × 10^{-3}, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionizing escape fraction from Freesia is f_esc∼eq 2{{ per cent}}. While [C II] emission is extended (radius 1.54 kpc), [O III] is concentrated in Freesia-A (0.85 kpc), where the ratio Σ _[O III]/Σ _[C II]≃ 10. As many high-z galaxies, Freesia lies below the local [C II]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) that disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [C II]-SFR deviations.

We study the kinematical properties of galaxies in the Epoch of Reionization via the [C II]158 μm line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disc and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disc inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use `Althæa', a highly resolved (30 pc) simulated prototypical Lyman-break galaxy, in the redshift range z = 6-7, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: (I) merger, (II) spiral disc, and (III) disturbed disc. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in (I), (II), and (III), respectively. We derive a generalized dynamical mass versus [C II]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor 2 (4). A Tully-Fisher relation is found for the observed high-z galaxies, i.e. L_{[C II]} ∝ (FWHM)^{1.80 ± 0.35} for which we provide a simple, physically based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [C II]. When seen face-on, Althæa is always detected at >5σ; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects.

We compute the six-particle maximally-helicity-violating (MHV) and next-to-MHV (NMHV) amplitudes in planar maximally supersymmetric Yang-Mills theory through seven loops and six loops, respectively, as an application of the extended Steinmann relations and using the cosmic Galois coaction principle. Starting from a minimal space of functions constructed using these principles, we identify the amplitude by matching its symmetries and predicted behavior in various kinematic limits. Through five loops, the MHV and NMHV amplitudes are uniquely determined using only the multi-Regge and leading collinear limits. Beyond five loops, the MHV amplitude requires additional data from the kinematic expansion around the collinear limit, which we obtain from the Pentagon Operator Product Expansion, and in particular from its single-gluon bound state contribution. We study the MHV amplitude in the self-crossing limit, where its singular terms agree with previous predictions. Analyzing and plotting the amplitudes along various kinematical lines, we continue to find remarkable stability between loop orders.

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that the object is made up of a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at 2.5 m s^-1. The expansion of the ring defines a cone with a 40° half-opening. We propose that the P/2016 G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.

We report the spectroscopic confirmation and modeling of the quadruply imaged quasar GRAL 113100-441959, the first gravitational lens (GL) to be discovered from a machine learning technique that only relies on the relative positions and fluxes of the observed images without considering colour informations. Follow-up spectra obtained with Keck/LRIS reveal the lensing nature of this quadruply imaged quasar with redshift z_s = 1.090 ± 0.002, but show no evidence of the central lens galaxy. Using the image positions and G-band flux ratios provided by Gaia Data Release 2 as constraints, we modeled the system with a singular power-law elliptical mass distribution (SPEMD) plus external shear, to different levels of complexity. We show that relaxing the isothermal constraint of the SPEMD does not lead to statistically significant different results in terms of fitting the lensing data. We thus simplified the SPEMD to a singular isothermal ellipsoid to estimate the Einstein radius of the main lens galaxy θ_E = 0.″851, the intensity and position angle of the external shear (γ,θ_γ) = (0.044, 11.°5), and we predict the lensing galaxy position to be (θ_gal,1, θ_gal,2) = (-0.″424, -0.″744) with respect to image A. We provide time delay predictions for pairs of images, assuming a plausible range of lens redshift values z_l between 0.5 and 0.9. Finally, we examine the impact on time delays of the so-called source position transformation, a family of degeneracies existing between different mass density profiles that reproduce most of the lensing observables equally well. We show that this effect contributes significantly to the time delay error budget and cannot be ignored during the modeling. This has implications for robust cosmography applications of lensed systems. GRAL 113100-441959 is the first in a series of seven new spectroscopically confirmed GLs discovered from Gaia Data Release 2.

After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete onto the star+disk system, while other material may fly by in a hyperbolic orbit. In contrast to the original cloud collapse event, this process will be a "cloudlet flyby" and/or "cloudlet capture" event: A Bondi-Hoyle-Lyttleton type accretion event, driven by the relative velocity between the star and the cloudlet. As we will show in this paper, if the cloudlet is small enough and has an impact parameter similar or less than GM_*/v_∞² (with v_∞ being the approach velocity), such a flyby and/or capture event would lead to arc-shaped or tail-shaped reflection nebulosity near the star. Those shapes of reflection nebulosity can be seen around several transitional disks and FU Orionis stars. Although the masses in the those arcs appears to be much less than the disk masses in these sources, we speculate that higher-mass cloudlet capture events may also happen occasionally. If so, they may lead to the tilting of the outer disk, because the newly infalling matter will have an angular momentum orientation entirely unrelated to that of the disk. This may be one possible explanation for the highly warped/tilted inner/outer disk geometries found in several transitional disks. We also speculate that such events, if massive enough, may lead to FU Orionis outbursts.

We study the photoevaporation of Jeans-unstable molecular clumps by isotropic FUV (6 eV < hν < 13.6 eV) radiation, through 3D radiative transfer hydrodynamical simulations implementing a non-equilibrium chemical network that includes the formation and dissociation of H₂. We run a set of simulations considering different clump masses (M=10 - 200 M_{\odot }) and impinging fluxes (G₀ = 2 × 10³ to 8 × 10⁴ in Habing units). In the initial phase, the radiation sweeps the clump as an R-type dissociation front, reducing the H₂ mass by a factor 40 - 90{{ per cent}}. Then, a weak (M∼eq 2) shock develops and travels towards the centre of the clump, which collapses while losing mass from its surface. All considered clumps remain gravitationally unstable even if radiation rips off most of the clump mass, showing that external FUV radiation is not able to stop clump collapse. However, the FUV intensity regulates the final H₂ mass available for star formation: for example, for G₀ < 10⁴ more than 10 per cent of the initial clump mass survives. Finally, for massive clumps ({≳ } 100 M_{\odot }) the H₂ mass increases by 25 - 50{{ per cent}} during the collapse, mostly because of the rapid density growth that implies a more efficient H₂ self-shielding.

We discuss the impact of the recent untagged analysis of B⁰ →D^* lν_barl decays by the Belle Collaboration on the extraction of the CKM element |V_cb | and provide updated SM predictions for the b → cτν observables R (D^*), P_τ, and F_L^D* . The value of |V_cb | that we find is about 2σ from the one from inclusive semileptonic B decays, and is very sensitive to the slope of the form factor at zero recoil which should soon become available from lattice calculations.

We present a sample of 21 hydrogen-free superluminous supernovae (SLSNe-I) and one hydrogen-rich SLSN (SLSN-II) detected during the five-year Dark Energy Survey (DES). These SNe, located in the redshift range 0.220 < z < 1.998, represent the largest homogeneously selected sample of SLSN events at high redshift. We present the observed g, r, i, z light curves for these SNe, which we interpolate using Gaussian processes. The resulting light curves are analysed to determine the luminosity function of SLSNe-I, and their evolutionary time-scales. The DES SLSN-I sample significantly broadens the distribution of SLSN-I light-curve properties when combined with existing samples from the literature. We fit a magnetar model to our SLSNe, and find that this model alone is unable to replicate the behaviour of many of the bolometric light curves. We search the DES SLSN-I light curves for the presence of initial peaks prior to the main light-curve peak. Using a shock breakout model, our Monte Carlo search finds that 3 of our 14 events with pre-max data display such initial peaks. However, 10 events show no evidence for such peaks, in some cases down to an absolute magnitude of <-16, suggesting that such features are not ubiquitous to all SLSN-I events. We also identify a red pre-peak feature within the light curve of one SLSN, which is comparable to that observed within SN2018bsz.

Protoplanetary disc surveys conducted with Atacama Large Millimetre Array (ALMA) are measuring disc radii in multiple star-forming regions. The disc radius is a fundamental quantity to diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous protoplanetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust mass increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly on to the star, the residual dust retained on Myr time-scales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc flux (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time.

We report the very bright detection of cold molecular gas with the IRAM NOEMA interferometer of the strongly lensed source WISE J132934.18+224327.3 at z = 2.04, the so-called Cosmic Eyebrow. This source has a similar spectral energy distribution from optical-mid/IR to submillimeter/radio but significantly higher fluxes than the well-known lensed SMG SMMJ 2135, the Cosmic Eyelash at z = 2.3. The interferometric observations unambiguously identify the location of the molecular line emission in two components, component CO32-A with {I}_CO(3-2)}=52.2+/- 0.9 Jy km s^-1 and component CO32-B with {I}_CO(3-2)}=15.7+/- 0.7 Jy km s^-1. Thus, our NOEMA observations of the CO(3-2) transition confirm the SMG-nature of WISE J132934.18+224327.3, resulting in the brightest CO(3-2) detection ever of an SMG. In addition, we present follow-up observations of the brighter component with the Green Bank Telescope (CO(1-0) transition) and IRAM 30 m telescope (CO(4-3) and [C I](1-0) transitions). The star formation efficiency of ∼100 L _⊙/(K km s^-1 pc²) is at the overlap region between merger-triggered and disk-like star formation activity and the lowest seen for lensed dusty star-forming galaxies. The determined gas depletion time ∼60 Myr, intrinsic infrared star formation SFR_IR ≈ 2000 M _⊙ yr^-1, and gas fraction M _mol/M _* = 0.44 indicate a starburst/merger-triggered star formation. The obtained data of the cold ISM—from CO(1-0) and dust continuum—indicates a gas mass μM _mol ∼ 15 × 10¹¹ M _⊙ for component CO32-A. Its unseen brightness offers us the opportunity to establish the Cosmic Eyebrow as a new reference source at z = 2 for galaxy evolution.

The upcoming radio interferometer Square Kilometre Array (SKA) is expected to directly detect the redshifted 21-cm signal from the neutral hydrogen present during the Cosmic Dawn. Temperature fluctuations from X-ray heating of the neutral intergalactic medium can dominate the fluctuations in the 21-cm signal from this time. This heating depends on the abundance, clustering, and properties of the X-ray sources present, which remain highly uncertain. We present a suite of three new large-volume, 349 Mpc a side, fully numerical radiative transfer simulations including QSO-like sources, extending the work previously presented in Ross et al. (2017). The results show that our QSOs have a modest contribution to the heating budget, yet significantly impact the 21-cm signal. Initially, the power spectrum is boosted on large scales by heating from the biased QSO-like sources, before decreasing on all scales. Fluctuations from images of the 21-cm signal with resolutions corresponding to SKA1-Low at the appropriate redshifts are well above the expected noise for deep integrations, indicating that imaging could be feasible for all the X-ray source models considered. The most notable contribution of the QSOs is a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. However, in the case of late Lyman-α saturation, this non-Gaussianity could be dramatically decreased particularly when heating occurs earlier. We conclude that increased non-Gaussianity is a promising signature of rare X-ray sources at this time, provided that Lyman-α saturation occurs before heating dominates the 21-cm signal.

Context. Pre-main-sequence variability characteristics can be used to probe the physical processes leading to the formation and initial evolution of both stars and planets.
Aims: The photometric variability of pre-main-sequence stars is studied at optical wavelengths to explore star-disk interactions, accretion, spots, and other physical mechanisms associated with young stellar objects.
Methods: We observed a field of 16' × 16' in the star-forming region Pelican Nebula (IC 5070) at BVRI wavelengths for 90 nights spread over one year in 2012-2013. More than 250 epochs in the VRI bands are used to identify and classify variables up to V ∼ 21 mag. Their physical association with the cluster IC 5070 is established based on the parallaxes and proper motions from the Gaia second data release (DR2). Multiwavelength photometric data are used to estimate physical parameters based on the isochrone fitting and spectral energy distributions.
Results: We present a catalog of optical time-series photometry with periods, mean magnitudes, and classifications for 95 variable stars including 67 pre-main-sequence variables towards star-forming region IC 5070. The pre-main-sequence variables are further classified as candidate classical T Tauri and weak-line T Tauri stars based on their light curve variations and the locations on the color-color and color-magnitude diagrams using optical and infrared data together with Gaia DR2 astrometry. Classical T Tauri stars display variability amplitudes up to three times the maximum fluctuation in disk-free weak-line T Tauri stars, which show strong periodic variations. Short-term variability is missed in our photometry within single nights. Several classical T Tauri stars display long-lasting (≥10 days) single or multiple fading and brightening events of up to two magnitudes at optical wavelengths. The typical mass and age of the pre-main-sequence variables from the isochrone fitting and spectral energy distributions are estimated to be ≤1 M_⊙ and ∼2 Myr, respectively. We do not find any correlation between the optical amplitudes or periods with the physical parameters (mass and age) of pre-main-sequence stars.
Conclusions: The low-mass pre-main-sequence stars in the Pelican Nebula region display distinct variability and color trends and nearly 30% of the variables exhibit strong periodic signatures attributed to cold spot modulations. In the case of accretion bursts and extinction events, the average amplitudes are larger than one magnitude at optical wavelengths. These optical magnitude fluctuations are stable on a timescale of one year.

Full Tables 1 and 2 are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/627/A135

We present a detailed analysis of the effect of an observationally determined dark matter (DM) velocity distribution function (VDF) of the Milky Way (MW) on DM direct detection rates. We go beyond local kinematic tracers and use rotation curve data up to 200 kpc to construct a MW mass model and self-consistently determine the local phase-space distribution of DM. This approach mitigates any incomplete understanding of local dark matter-visible matter degeneracies that can affect the determination of the VDF. Comparing with the oft used Standard Halo Model (SHM), which assumes an isothermal VDF, we look at how the tail of the empirically determined VDF alters our interpretation of the present direct detection WIMP DM cross section exclusion limits. While previous studies have suggested a very large difference (of more than an order of magnitude) in the bounds at low DM masses, we show that accounting for the detector response at low threshold energies, the difference is still significant although less extreme. The change in the number of signal events, when using the empirically determined DM VDF in contrast to the SHM VDF, is most prominent for low DM masses for which the shape of the recoil energy spectrum depends sensitively on the detector threshold energy as well as detector response near the threshold. We demonstrate that these trends carry over to the respective DM exclusion limits, modulo detailed understanding of the experimental backgrounds. With the unprecedented precision of astrometric data in the GAIA era, use of observationally determined DM phase space will become a critical and necessary ingredient for DM searches. We provide an accurate fit to the current best observationally determined DM VDF (and self-consistent local DM density) for use in analyzing current DM direct detection data by the experimental community.

This work explores the detailed chemistry of the Milky Way bulge using the HERMES spectrograph on the Anglo-Australian Telescope. Here, we present the abundance ratios of 13 elements for 832 red giant branch and clump stars along the minor bulge axis at latitudes b = -10^○, - 7.5^○, and -5^○. Our results show that none of the abundance ratios vary significantly with latitude. We also observe disc-like [Na/Fe] abundance ratios, which indicate that the bulge does not contain helium-enhanced populations as observed in some globular clusters. Helium enhancement is therefore not the likely explanation for the double red-clump observed in the bulge. We confirm that bulge stars mostly follow abundance trends observed in the disc. However, this similarity is not confirmed across all elements and metallicity regimes. The more metal-poor bulge population at [Fe/H] ≲ - 0.8 is enhanced in the elements associated with core collapse supernovae (SNeII). In addition, the [La/Eu] abundance ratio suggests higher r-process contribution, and likely higher star formation in the bulge compared to the disc. This highlights the complex evolution in the bulge, which should be investigated further, both in terms of modelling; and with additional observations of the inner Galaxy.

Asteroids and other Small Solar System Bodies (SSSBs) are of high general and scientific interest in many aspects. The origin, formation, and evolution of our Solar System (and other planetary systems) can be better understood by analysing the constitution and physical properties of small bodies in the Solar System. Currently, two space missions (Hayabusa2, OSIRIS-REx) have recently arrived at their respective targets and will bring a sample of the asteroids back to Earth. Other small body missions have also been selected by, or proposed to, space agencies. The threat posed to our planet by near-Earth objects (NEOs) is also considered at the international level, and this has prompted dedicated research on possible mitigation techniques. The DART mission, for example, will test the kinetic impact technique. Even ideas for industrial exploitation have risen during the last years. Lastly, the origin of water and life on Earth appears to be connected to asteroids. Hence, future space mission projects will undoubtedly target some asteroids or other SSSBs. In all these cases and research topics, specific knowledge of the structure and mechanical behaviour of the surface as well as the bulk of those celestial bodies is crucial. In contrast to large telluric planets and dwarf planets, a large proportion of such small bodies is believed to consist of gravitational aggregates (`rubble piles') with no—or low—internal cohesion, with varying macro-porosity and surface properties (from smooth regolith covered terrain, to very rough collection of boulders), and varying topography (craters, depressions, ridges). Bodies with such structure can sustain some plastic deformation without being disrupted in contrast to the classical visco-elastic models that are generally valid for planets, dwarf planets, and large satellites. These SSSBs are hence better described through granular mechanics theories, which have been a subject of intense theoretical, experimental, and numerical research over the last four decades. This being the case, it has been necessary to use the theoretical, numerical and experimental tools developed within soil mechanics, granular dynamics, celestial mechanics, chemistry, condensed matter physics, planetary and computer sciences, to name the main ones, in order to understand the data collected and analysed by observational astronomy (visible, thermal, and radio), and different space missions. In this paper, we present a review of the multi-disciplinary research carried out by these different scientific communities in an effort to study SSSBs.

A correlation between proto-planetary disc radii and sub-mm fluxes has been recently reported. In this letter, we show that the correlation is a sensitive probe of grain growth processes. Using models of grain growth and drift, we have shown in a companion paper that the observed disc radii trace where the dust grains are large enough to have a significant sub-mm opacity. We show that the observed correlation emerges naturally if the maximum grain size is set by radial drift, implying relatively low values of the viscous α parameter ≲0.001. In this case, the relation has an almost universal normalization, while if the grain size is set by fragmentation the flux at a given radius depends on the dust-to-gas ratio. We highlight two observational consequences of the fact that radial drift limits the grain size. The first is that the dust masses measured from the sub-mm could be overestimated by a factor of a few. The second is that the correlation should be present also at longer wavelengths (e.g. 3mm), with a normalization factor that scales as the square of the observing frequency as in the optically thick case.

We review open questions and prospects for progress in ultrahigh-energy cosmic ray (UHECR) research, based on a series of discussions that took place during the `The High-Energy Universe: Gamma-Ray, Neutrino, and Cosmic-ray Astronomy' MIAPP workshop in 2018. Specifically, we overview open questions on the origin of the bulk of UHECRs, the UHECR mass composition, the origin of the end of the cosmic-ray spectrum, the transition from Galactic to extragalactic cosmic rays, the effect of magnetic fields on the trajectories of UHECRs, anisotropy expectations for specific astrophysical scenarios, hadronic interactions, and prospects for discovering neutral particles as well as new physics at ultrahigh energies. We also briefly present upcoming and proposed UHECR experiments and discuss their projected science reach.

Local, manifestly dual-conformally invariant loop integrands are now known for all finite quantities associated with observables in planar, maximally supersymmetric Yang-Mills theory through three loops. These representations, however, are not infrared-finite term by term and therefore require regularization; and even using a regulator consistent with dual-conformal invariance, ordinary methods of loop integration would naïvely obscure this symmetry. In this work, we show how any planar loop integral through at least two loops can be systematically regulated and evaluated directly in terms of strictly finite, manifestly dual-conformal Feynman-parameter integrals. We apply these methods to the case of the two-loop ratio and remainder functions for six particles, reproducing the known results in terms of individually regulated local loop integrals, and we comment on some of the novelties that arise for this regularization scheme not previously seen at one loop.

We present an improved determination of the Hubble constant from Hubble Space Telescope (HST) observations of 70 long-period Cepheids in the Large Magellanic Cloud (LMC). These were obtained with the same WFC3 photometric system used to measure extragalactic Cepheids in the hosts of SNe Ia. Gyroscopic control of HST was employed to reduce overheads while collecting a large sample of widely separated Cepheids. The Cepheid period-luminosity relation provides a zero-point-independent link with 0.4% precision between the new 1.2% geometric distance to the LMC from detached eclipsing binaries (DEBs) measured by Pietrzyński et al. and the luminosity of SNe Ia. Measurements and analysis of the LMC Cepheids were completed prior to knowledge of the new DEB LMC distance. Combined with a refined calibration of the count-rate linearity of WFC3-IR with 0.1% precision, these three improved elements together reduce the overall uncertainty in the geometric calibration of the Cepheid distance ladder based on the LMC from 2.5% to 1.3%. Using only the LMC DEBs to calibrate the ladder, we find H ₀ = 74.22 ± 1.82 km s^-1 Mpc^-1 including systematic uncertainties, 3% higher than before for this particular anchor. Combining the LMC DEBs, masers in NGC 4258, and Milky Way parallaxes yields our best estimate: H ₀ = 74.03 ± 1.42 km s^-1 Mpc^-1, including systematics, an uncertainty of 1.91%-15% lower than our best previous result. Removing any one of these anchors changes H ₀ by less than 0.7%. The difference between H ₀ measured locally and the value inferred from Planck CMB and ΛCDM is 6.6 ± 1.5 km s^-1 Mpc^-1 or 4.4σ (P = 99.999% for Gaussian errors) in significance, raising the discrepancy beyond a plausible level of chance. We summarize independent tests showing that this discrepancy is not attributable to an error in any one source or measurement, increasing the odds that it results from a cosmological feature beyond ΛCDM.

We calculate the contributions to the two-loop scattering amplitudes h → gg, h → ggg and h\to q\overline{q}g that arise from a modified trilinear Higgs coupling λ. Analytic expressions are obtained by performing an asymptotic expansion near the limit of infinitely heavy top quark. The calculated amplitudes are necessary to study the impact of the O(λ ) corrections to the transverse momentum distributions ( p _T,h ) in single-Higgs production at hadron colliders for low and moderate values of p _{T, h} .

Relativistic heavy ion collisions produce nuclei-sized droplets of quark-gluon plasma whose expansion is well described by viscous hydrodynamic calculations. Over the past half decade, this formalism was also found to apply to smaller droplets closer to the size of individual nucleons, as produced in p +p and p +A collisions. The hydrodynamic paradigm was further tested with a variety of collision species, including p +Au ,d +Au , and ³He+Au producing droplets with different geometries. Nevertheless, questions remain regarding the importance of pre-hydrodynamic evolution and the exact medium properties during the hydrodynamic evolution phase, as well as the applicability of alternative theories that argue the agreement with hydrodynamics is accidental. In this work we explore options for new collision geometries including p +O and O +O proposed for running at the Large Hadron Collider, as well as ⁴He+Au ,C +Au ,O +Au , and ^Be,9₇+Au at the Relativistic Heavy Ion Collider.

In a shallow near-infrared survey of the dwarf irregular galaxy, NGC 3109, near the periphery of the Local Group, we have found eight Mira variables, seven of which appear to be oxygen-rich (O-Miras). The periods range from about 430 d to almost 1500 d. Because of our relatively bright limiting magnitude, only 45 of the more than 400 known carbon stars were measured, but none was found to be a large amplitude variable. One of the Miras may be an unrecognized C star. Five of the O-Miras are probably hot-bottom burning stars considering that they are brighter than expected from the period-luminosity relation of Miras and that, by comparison with theoretical evolutionary tracks, they appear to have masses {≳}4 M_⊙. A census of very long period (P > 1000 d) Miras in the Galaxy and Magellanic Clouds is presented and discussed together with the newly discovered long-period, but relatively blue, variables in NGC 3109. New JHKL photometry is presented for three O-rich long-period Miras in the Small Magellanic Cloud (including a candidate super-AGB star).

We present a spatially resolved study of the relation between dust and metallicity in the nearby spiral galaxies M 101 (NGC 5457) and NGC 628 (M 74). We explore the relation between the chemical abundances of their gas and stars with their dust content and their chemical evolution. The empirical spatially resolved oxygen effective yield and the gas-to-dust mass ratio (GDR) across both disc galaxies are derived, sampling 1 dex in oxygen abundance. We find that the metal budget of the NGC 628 disc and most of the M 101 disc appears consistent with the predictions of the simple model of chemical evolution for an oxygen yield between half and one solar, whereas the outermost region (R ≥ 0.8 R_{25}) of M 101 presents deviations suggesting the presence of gas flows. The GDR-metallicity relation shows a two slopes behaviour, with a break at 12 + log(O/H) ≈ 8.4, a critical metallicity predicted by theoretical dust models when stardust production equals grain growth. A relation between GDR and the fraction of molecular to total gas, Σ _{H₂}/Σ _{gas} is also found. We suggest an empirical relationship between GDR and the combination of 12 + log(O/H), for metallicity, and Σ _{H₂}/Σ _{gas}, a proxy for the molecular clouds fraction. The GDR is closely related with metallicity at low abundance and with Σ _{H₂}/Σ _{gas} for higher metallicities suggesting interstellar medium dust growth. The ratio Σ _{dust}/Σ _{star} correlates well with 12 + log(O/H) and strongly with log(N/O) in both galaxies. For abundances below the critical one, the `stardust' production gives us a constant value suggesting a stellar dust yield similar to the oxygen yield.

We report on a superdense star-forming region with an effective radius (R_e) smaller than 13 pc identified at z = 6.143 and showing a star formation rate density Σ_SFR ∼ 1000 M_⊙ yr^-1 kpc^-2 (or conservatively >300 M_⊙ yr^-1 kpc^-2). Such a dense region is detected with S/N ≳ 40 hosted by a dwarf extending over 440 pc, dubbed D1. D1 is magnified by a factor 17.4(±5.0) behind the Hubble Frontier Field galaxy cluster MACS J0416 and elongated tangentially by a factor 13.2 ± 4.0 (including the systematic errors). The lens model accurately reproduces the positions of the confirmed multiple images with a rms of 0.35 arcsec. D1 is part of an interacting star-forming complex extending over 800 pc. The SED-fitting, the very blue ultraviolet slope (β ≃ -2.5, F_λ ∼ λ^β), and the prominent Lyα emission of the stellar complex imply that very young (<10-100 Myr), moderately dust-attenuated (E(B - V) < 0.15) stellar populations are present and organized in dense subcomponents. We argue that D1 (with a stellar mass of 2 × 10⁷ M_⊙) might contain a young massive star cluster of M ≲ 10⁶ M_⊙ and M_UV ≃ -15.6 (or m_UV = 31.1), confined within a region of 13 pc, and not dissimilar from some local super star clusters (SSCs). The ultraviolet appearance of D1 is also consistent with a simulated local dwarf hosting an SSC placed at z = 6 and lensed back to the observer. This compact system fits into some popular globular cluster formation scenarios. We show that future high spatial resolution imaging (e.g. E-ELT/MAORY-MICADO and VLT/MAVIS) will allow us to spatially resolve light profiles of 2-8 pc.

We present an implementation of the relativistic three-particle quantization condition including both s- and d-wave two-particle channels. For this, we develop a systematic expansion of the three-particle K matrix, K _df,3, about threshold, which is the generalization of the effective range expansion of the two-particle K matrix, K ₂. Relativistic invariance plays an important role in this expansion. We find that d-wave two-particle channels enter first at quadratic order. We explain how to implement the resulting multichannel quantization condition, and present several examples of its application. We derive the leading dependence of the threshold three-particle state on the two-particle d-wave scattering amplitude, and use this to test our implementation. We show how strong two-particle d-wave interactions can lead to significant effects on the finite-volume three-particle spectrum, including the possibility of a generalized three-particle Efimov-like bound state. We also explore the application to the 3 π ⁺ system, which is accessible to lattice QCD simulations, where we study the sensitivity of the spectrum to the components of K _df,3. Finally, we investigate the circumstances under which the quantization condition has unphysical solutions.

High-energy radiation from a planet host star can have strong influence on the final habitability of a system through several mechanisms. In this context we have constructed a catalogue containing the X-ray luminosities, as well as basic stellar and planetary properties of all known stars hosting giant planets (> 0.1 M_J) that have been observed by the Chandra X-ray Observatory, XMM-Newton, and/or ROSAT. Specifically in this paper we present a first application of this catalogue to search for a possible imprint of X-ray photoevaporation of planet-forming discs on the present-day orbital distribution of the observed giant planets. We found a suggestive void in the semimajor axis, a, versus X-ray luminosity, L_x, plane, roughly located between a ∼ 0.05-1 au and L_x ∼ 10²⁷-10^{29} erg s^{-1}, which would be expected if photoevaporation played a dominant role in the migration history of these systems. However, due to the small observational sample size, the statistical significance of this feature cannot be proven at this point.

We propose an expression for a local planetesimal formation rate proportional to the instantaneous radial pebble flux. The result—a radial planetesimal distribution—can be used as an initial condition to study the formation of planetary embryos. We follow the idea that one needs particle traps to locally enhance the dust-to-gas ratios sufficiently, such that particle gas interactions can no longer prevent planetesimal formation on small scales. The locations of these traps can emerge everywhere in the disk. Their occurrence and lifetime is subject to ongoing research; thus, here they are implemented via free parameters. This enables us to study the influence of the disk properties on the formation of planetesimals, predicting their time-dependent formation rates and the location of primary pebble accretion. We show that large α-values of 0.01 (strong turbulence) prevent the formation of planetesimals in the inner part of the disk, arguing for lower values of around 0.001 (moderate turbulence), at which planetesimals form quickly at all places where they are needed for proto-planets. Planetesimals form as soon as dust has grown to pebbles (mm to dm) and the pebble flux reaches a critical value, which is after a few thousand years at 2-3 au and after a few hundred thousand years at 20-30 au. Planetesimal formation lasts until the pebble supply has decreased below a critical value. The final spatial planetesimal distribution is steeper compared to the initial dust and gas distribution, which helps explain the discrepancy between the minimum mass solar nebula and viscous accretion disks.

In the presence of a relative suppression of the different quarkonium states and due to the feed downs, the statement on pg. 1 & 2 "one is entitled […] to square the measured suppression factor [12] in pPb collisions to extrapolate to PbPb collisions."

Context. Classical Cepheids (CCs) and RR Lyrae stars (RRLs) are important classes of variable stars used as standard candles to estimate galactic and extragalactic distances. Their multiplicity is imperfectly known, particularly for RRLs. Astoundingly, to date only one RRL has convincingly been demonstrated to be a binary, TU UMa, out of tens of thousands of known RRLs.
Aims: Our aim is to detect the binary and multiple stars present in a sample of Milky Way CCs and RRLs.
Methods: In the present article, we combine the HIPPARCOS and Gaia DR2 positions to determine the mean proper motion of the targets, and we search for proper motion anomalies (PMa) caused by close-in orbiting companions.
Results: We identify 57 CC binaries from PMa out of 254 tested stars and 75 additional candidates, confirming the high binary fraction of these massive stars. For 28 binary CCs, we determine the companion mass by combining their spectroscopic orbital parameters and astrometric PMa. We detect 13 RRLs showing a significant PMa out of 198 tested stars, and 61 additional candidates.
Conclusions: We determine that the binary fraction of CCs is likely above 80%, while that of RRLs is at least 7%. The newly detected systems will be useful to improve our understanding of their evolutionary states. The discovery of a significant number of RRLs in binary systems also resolves the long-standing mystery of their extremely low apparent binary fraction.

Full Tables A.1 and A.3 are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/623/A116