Search

Determining the bubble wall velocity in first-order phase transitions is a challenging task, requiring the solution of (coupled) equations of motion for the scalar field and Boltzmann equations for the particles in the plasma. The collision terms appearing in the Boltzmann equation present a prominent source of uncertainty as they are often known only at leading log accuracy. In this paper, we derive upper and lower bounds on the wall velocity, corresponding to the local thermal equilibrium and ballistic limits. These bounds are completely independent of the collision terms. For the ballistic approximation, we argue that the inhomogeneous plasma temperature and velocity distributions across the bubble wall should be taken into account. This way, the hydrodynamic obstruction previously observed in local thermal equilibrium is also present for the ballistic approximation. This is essential for the ballistic approximation to provide a lower bound on the wall velocity. We use a model-independent approach to study the behaviour of the limiting wall velocities as a function of a few generic parameters, and we test our developments in the singlet extended Standard Model., Comment: 47 pages, 8 figures

WallGo is an open source software for the computation of the bubble wall velocity in first-order cosmological phase transitions. It also computes the energy budget available for the generation of gravitational waves. The main part of WallGo, built in Python, determines the wall velocity by solving the scalar-field(s) equation of motion, the Boltzmann equations and energy-momentum conservation for the fluid velocity and temperature. WallGo also includes two auxiliary modules: WallGoMatrix, which computes matrix elements for out-of-equilibrium particles, and WallGoCollision, which performs higher-dimensional integrals for Boltzmann collision terms. Users can implement custom models by defining an effective potential and specifying a list of out-of-equilibrium particles and their interactions. As the first public software to compute the wall velocity including out-of-equilibrium contributions, WallGo improves the precision of the computation compared to common assumptions in earlier computations. It utilises a spectral method for the deviation from equilibrium and collision terms that provides exponential convergence in basis polynomials, and supports multiple out-of-equilibrium particles, allowing for Boltzmann mixing terms. WallGo is tailored for non-runaway wall scenarios where leading-order coupling effects dominate friction. While this work introduces the software and the underlying theory, a more detailed documentation can be found in https://wallgo.readthedocs.io., Comment: 69 pages, 8 figures

The production of three Higgs bosons at hadron colliders can be enhanced by a double-resonant effect in the $\mathbb{Z}_2$-symmetric two-real-singlet extension of the Standard Model, making it potentially observable in future LHC runs. The production rate is maximized for large scalar couplings, which prompts us to carefully reconsider the perturbativity constraints on the theory. This leads us to construct a new set of 140 benchmark points that have a triple Higgs boson production cross-section at least 100 times larger than the SM value. Furthermore, we study the dynamics of the electroweak phase transition, both analytically at leading order, and numerically without the high-temperature expansion. Both analyses indicate that a first-order phase transition is incompatible with the requirement that both singlets have a non-zero vev in the present-day vacuum, as required by doubly-enhanced triple Higgs boson production. Allowing instead one of the singlets to remain at zero field value opens up the possibility of a first-order phase transition, while di-Higgs boson production can still be enhanced by a (single) resonance., Comment: 34 pages, 5 figures. 140 new benchmark points for the TRSM are included in the ancillary file New_BM_points_arXiv.dat

Models with radiative symmetry breaking typically feature strongly supercooled first-order phase transitions, which result in an observable stochastic gravitational wave background. In this work, we analyse the role of higher order thermal corrections for these transitions, applying high-temperature dimensional reduction to a theory with dimensional transmutation. In particular, we study to what extent high-temperature effective field theories (3D EFT) can be used. We find that despite significant supercooling down from the critical temperature, the high-temperature expansion for the bubble nucleation rate can be applied using the 3D EFT framework, and we point out challenges in the EFT description. We compare our findings to previous studies, and find that the next-to-leading order corrections obtained in this work have a significant effect on the predictions for GW observables, motivating a further exploration of higher order thermal effects., Comment: Revised version, published in JHEP. 45 pages, 9 figures

The bubble expansion velocity is an important parameter in the prediction of gravitational waves from first order phase transitions. This parameter is difficult to compute, especially in phase transitions in strongly coupled theories. In this work, we present a method to estimate the wall velocity for phase transitions with a large enthalpy jump, valid for weakly and strongly coupled theories. We find that detonations are disfavored in this limit, but wall velocities are not necessarily small. We also investigate the effect of two other features in the equation of state: non-conformal sound speeds and a limited range of temperatures for which the phases exist. We find that the former can increase the wall velocity for a given nucleation temperature, and the latter can restrict the wall velocities to small values. To test our approach, we use holographic phase transitions, which typically display these three features. We find excellent agreement with numerically obtained values of the wall velocity. We also demonstrate that the implications for gravitational waves can be significant., Comment: 12 pages, 5 figures

Accurately determining bubble wall velocities in first-order phase transitions is of great importance for the prediction of gravitational wave signals and the matter-antimatter asymmetry. However, it is a challenging task which typically depends on the underlying particle physics model. Recently, it has been shown that assuming local thermal equilibrium can provide a good approximation when calculating the bubble wall velocity. In this paper, we provide a model-independent determination of bubble wall velocities in local thermal equilibrium. Our results show that, under the reasonable assumption that the sound speeds in the plasma are approximately uniform, the hydrodynamics can be fully characterized by four quantities: the phase strength $\alpha_n$, the ratio of the enthalpies in the broken and symmetric phases, $\Psi_n$, and the sound speeds in both phases, $c_s$ and $c_b$. We provide a code snippet that allows for a determination of the wall velocity and energy fraction in local thermal equilibrium in any model. In addition, we present a fit function for the wall velocity in the case $c_s = c_b = 1/\sqrt 3$., Comment: 34 pages, 6 figures; v2 matches published version

First order phase transitions are well-motivated and extensively studied sources of gravitational waves (GWs) from the early Universe. The vacuum energy released during such transitions is assumed to be transferred primarily either to the expanding walls of bubbles of true vacuum, whose collisions source GWs, or to the surrounding plasma, producing sound waves and turbulence, which act as GW sources. In this Letter, we study an alternative possibility that has so far not been considered: the released energy gets transferred primarily to feebly interacting particles that do not admit a fluid description but simply free-stream individually. We develop the formalism to study the production of GWs from such configurations, and demonstrate that such GW signals have qualitatively distinct characteristics compared to conventional sources and are potentially observable with near-future GW detectors., Comment: 5 pages+appendices, 8 figures

The energy budget for gravitational waves of a cosmological first order phase transitions depends on the speed of sound in the thermal plasma in both phases around the bubble wall. Working in the real-singlet augmented Standard Model, which admits a strong two-step electroweak phase transition, we compute higher order corrections to the pressure and sound speed. We compare our result to lower-order approximations to the sound speed and the energy budget and investigate the impact on the gravitational wave signal. We find that deviations in the speed of sound from $c_s^2 = 1/3$ are enhanced up to $\mathcal O(5\%)$ in our higher-order computation. This results in a suppression in the energy budget of up to $\mathcal O (50\%)$ compared to approximations assuming $c_s^2 = 1/3$. The effect is most significant for hybrid and detonation solutions. We generalise our discussion to the case of multiple inert scalars and the case of a reduced number of fermion families in order to mimic hypothetical dark sector phase transitions. In this sector with modified field content, the sound speed can receive significant suppression, with potential order-of-magnitude impact on the gravitational wave amplitude., Comment: 43 pages, 10 figures

We re-derive the vev-insertion approximation (VIA) source in electroweak baryogenesis. In contrast to the original derivation, we rely solely on 1-particle-irreducible self-energy diagrams. We solve the Green's function equations both perturbatively and resummed over all vev-insertions. The VIA source corresponds to the leading order contribution in the gradient expansion of the Kadanoff-Baym (KB) equations. We find that it vanishes both for bosons and fermions, both in the perturbative and in the resummed approach. Interestingly, the non-existence of the source is a result of a cancellation between different terms in the KB equations, and not of a pathology in the vev-insertion approximation itself., Comment: 29 Pages

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

The baryon asymmetry of the Universe is one of the central motivations to expect physics beyond the Standard Model. In this Snowmass white paper, we review the challenges and opportunities in testing some of the central paradigms that predict physics at scales low enough to expect new experimental data in the next decade. Focusing on theoretical ideas and some of their experimental implications, in particular, we discuss neutron-antineutron transformations, flavor observables, next generation colliders, future neutron facilities, gravitational waves, searches for permanent electric dipole moments, $0\nu \beta \beta $ decay and some future large underground experiments as methods to test post-sphaleron baryogenesis, electroweak baryogenesis, mesogenesis and low scale leptogenesis. Finally, we comment on the cases where high scale physics can be probed through some of these same mechanisms., Comment: White paper to be submitted to the Snowmass Process' Rare Processes and Precision Measurements Frontier, Baryon and Lepton Number Violation Topical Group (This v2 version corrects a missing author and an author's name which was misspelled. Apologies to these authors for this oversight on the part of the submitter!)

Theoretical investigations into the evolution of the early universe are an essential part of particle physics that allow us to identify viable extensions to the Standard Model as well as motivated parameter space that can be probed by various experiments and observations. In this white paper, we review particle physics models of the early universe. First, we outline various models that explain two essential ingredients of the early universe (dark matter and baryon asymmetry) and those that seek to address current observational anomalies. We then discuss dynamics of the early universe in models of neutrino masses, axions, and several solutions to the electroweak hierarchy problem. Finally, we review solutions to naturalness problems of the Standard Model that employ cosmological dynamics., Comment: contribution to Snowmass 2021

The Standard Model of Particle Physics cannot explain the observed baryon asymmetry of the Universe. This observation is a clear sign of new physics beyond the Standard Model. There have been many recent theoretical developments to address this question. Critically, many new physics models that generate the baryon asymmetry have a wide range of repercussions for many areas of theoretical and experimental particle physics. This white paper provides an overview of such recent theoretical developments with an emphasis on experimental testability., Comment: Contribution to Snowmass 2021. Solicited white paper from TF08. v2: includes additional contributions

We study the effect of density perturbations on the process of first-order phase transitions and gravitational wave production in the early Universe. We are mainly interested in how the distribution of nucleated bubbles is affected by fluctuations in the local temperature. We find that large-scale density fluctuations ($H_* < k_* < \beta$) result in a larger effective bubble size at the time of collision, enhancing the produced amplitude of gravitational waves. The amplitude of the density fluctuations necessary for this enhancement is ${\cal P}_\zeta (k_*) \gtrsim (\beta / H_*)^{-2}$, and therefore the gravitational wave signal from first-order phase transitions with relatively large $\beta / H_*$ can be significantly enhanced by this mechanism even for fluctuations with moderate amplitudes., Comment: 22 pages, 14 figures, published version

We study to what extend LISA can observe features of gravitational wave spectra originating from cosmological first-order phase transitions. We focus on spectra which are of the form of double-broken power laws. These spectra are predicted by hydrodynamic simulations and also analytical models such as the sound shell model. We argue that the ratio of the two break frequencies is an interesting observable since it can be related to the wall velocity while overall amplitude and frequency range are often degenerate for the numerous characteristics of the phase transition. Our analysis uses mock data obtained from the power spectra predicted by the simplified simulations and the sound shell model and analyzes the detection prospects using $\chi^2$-minimization and likelihood sampling. We point out that the prospects of observing two break frequencies from the electroweak phase transition is hindered by a shift of the spectrum to smaller frequencies for strong phase transitions. Finally, we also highlight some differences between signals from the sound shell model compared to simulations., Comment: 32 pages, 12 figures

We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions ($3\%

We present here far-infrared photometry of galaxies in a sample that is relatively unexplored at these wavelengths: low-metallicity dwarf galaxies with moderate star formation rates. Four dwarf irregular galaxies from the $\mathrm{L{\small{ITTLE}}}$ $\mathrm{T{\small{HINGS}}}$ survey are considered, with deep $\textit{Herschel}$ PACS and SPIRE observations at 100 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m. Results from modified-blackbody fits indicate that these galaxies have low dust masses and cooler dust temperatures than more actively star-forming dwarfs, occupying the lowest $L_\mathrm{TIR}$ and $M_\mathrm{dust}$ regimes seen among these samples. Dust-to-gas mass ratios of $\sim$10$^{-5}$ are lower, overall, than in more massive and active galaxies, but are roughly consistent with the broken power law relation between the dust-to-gas ratio and metallicity found for other low-metallicity systems. Chemical evolution modeling suggests that these dwarf galaxies are likely forming very little dust via stars or grain growth, and have very high dust destruction rates., Comment: 26 pages, 19 figures. To be published in The Astronomical Journal

We study particle production and unitarity violation caused by a curved target space right after inflation. We use the inflaton field value instead of cosmic time as the time variable, and derive a semiclassical formula for the spectrum of produced particles. We then derive a simple condition for unitarity violation during preheating, which we confirm by our semiclassical method and numerical solution. This condition depends not only on the target space curvature but also on the height of the inflaton potential at the end of inflation. This condition tells us, for instance, that unitarity is violated in running kinetic inflation and Higgs inflation, while unitarity is conserved in $\alpha$-attractor inflation and Higgs-Palatini inflation., Comment: 25 pages, 6 figures; v2: clarifications and references added

We provide an easy method to obtain the kinetic energy fraction in gravitational waves, generated during a cosmological first-order phase transition, as a function of only the wall velocity and quantities that can be determined from the particle physics model at the nucleation temperature. This generalizes recent work that achieved this goal for detonations. Here we present the corresponding results for deflagrations and hybrids. Unlike for detonations, the sound speed in the symmetric phase also enters the analysis. We perform a detailed comparison between our model-independent approach and other approaches in the literature. We provide a Python code snippet to determine the kinetic energy fraction $K$ as a function of the wall velocity, the two speeds of sound and the strength parameter of the phase transition. We also assess how realistic sizable deviations in speed of sound are close to the phase transition temperature in a specific model., Comment: 31 pages, 5 figures, v2 matches published version

We derive one- and two-loop renormalization group equations (RGEs) of Higgs-$R^2$ inflation. This model has a non-minimal coupling between the Higgs and the Ricci scalar and a Ricci scalar squared term on top of the standard model. The RGEs derived in this paper are valid as long as the energy scale of interest (in the Einstein frame) is below the Planck scale. We also discuss implications to the inflationary predictions and the electroweak vacuum metastability., Comment: 35 pages, 1 figure; v2: references added, version accepted by JHEP

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium., Comment: 56 pages, 30 figures, accepted for publication in MNRAS

Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems., Comment: 18 pages, 11 figures, accepted for publication in A&A

We have conducted extensive lattice simulations to study the post-inflation dynamics of multifield models involving nonminimal couplings. We explore the parameter dependence of preheating in these models and describe the various time-scales that control such nonlinear processes as energy transfer, re-scattering, and the approach to radiation-domination and thermalization. In the limit of large nonminimal couplings ($\xi_I \sim 100$), we find that efficient transfer of energy from the inflaton condensate to radiative degrees of freedom, emergence of a radiation-dominated equation of state, and the onset of thermalization each consistently occur within $N_{\rm reh} \lesssim 3$ $e$-folds after the end of inflation, largely independent of the values of the other couplings in the models. The exception is the case of negative ellipticity, in which there is a misalignment between the dominant direction in field-space along which the system evolves and the larger of the nonminimal couplings $\xi_I$. In those cases, the field-space-driven parametric resonance is effectively shut off. More generally, the competition between the scalar fields' potential and the field-space manifold structure can yield interesting phenomena such as two-stage resonances. Despite the explosive particle production, which can lead to a quick depletion of the background energy density, the nonlinear processes do not induce any super-horizon correlations after the end of inflation in these models, which keeps predictions for CMB observables unaffected by the late-time amplification of isocurvature fluctuations. Hence the excellent agreement between primordial observables and recent observations is preserved for this class of models, even when we consider post-inflation dynamics., Comment: 32 pages (plus appendices), 17 figures. References added and minor edits to match published version

We study the energy budget of a first-order cosmological phase transition, which is an important factor in the prediction of the resulting gravitational wave spectrum. Formerly, this analysis was based mostly on simplified models as for example the bag equation of state. Here, we present a model-independent approach that is exact up to the temperature dependence of the speed of sound in the broken phase. We find that the only relevant quantities that enter in the hydrodynamic analysis are the speed of sound in the broken phase and a linear combination of the energy and pressure differences between the two phases which we call pseudotrace (normalized to the enthalpy in the broken phase). The pseudotrace quantifies the strength of the phase transition and yields the conventional trace of the energy-momentum tensor for a relativistic plasma (with speed of sound squared of one third). We study this approach in several realistic models of the phase transition and also provide a code snippet that can be used to determine the efficiency coefficient for a given phase transition strength and speed of sound. It turns out that our approach is accurate to the percent level for moderately strong phase transitions, while former approaches give at best the right order of magnitude., Comment: 24 pages, 2 figures

Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged), Comment: 23 pages, 16 figures, accepted for publication in A&A

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

We point out that a model with scalar fields with a large nonminimal coupling to the Ricci scalar, such as Higgs inflation, can be regarded as a nonlinear sigma model (NLSM). With the inclusion of not only the scalar fields but also the conformal mode of the metric, our definition of the target space of the NLSM is invariant under the frame transformation. We show that the $\sigma$-meson that linearizes this NLSM to be a linear sigma model (LSM) corresponds to the scalaron, the degree of freedom associated to the $R^2$ term in the Jordan frame. We demonstrate that quantum corrections inevitably induce this $\sigma$-meson in the large-$N$ limit, thus providing a frame independent picture for the emergence of the scalaron. The resultant LSM only involves renormalizable interactions and hence its perturbative unitarity holds up to the Planck scale unless it hits a Landau pole, which is in agreement with the renormalizability of quadratic gravity., Comment: 18 pages + appendices, 1 figures; v2: references added; v3: explanation clarified, published version; v4: one reference added

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies., Comment: Accepted for publication in A&A, 21 pages, 14 Figures, 7 Tables, 2 Appendices (with 3 Tables). Published version, unchanged results with respect to the previously posted version

In electroweak baryogenesis the baryon asymmetry of the universe is created during the electroweak phase transition. The quantum transport equations governing the dynamics of the plasma particles can be derived in the vev-insertion approximation, which treats the vev-dependent part of the particle masses as a perturbation. We calculate the next-to-leading order (NLO) contribution to the CP-violating source term and CP-conserving relaxation rate, corresponding to Feynman diagrams for the self-energies with four mass insertions. We consider both a pair of Weyl fermions and a pair of complex scalars, that scatter off the bubble wall. We find: (i) The NLO correction becomes large for $\mathcal O(1)$ couplings. If only the Standard Model (SM) Higgs obtains a vev during the phase transition, this implies the vev-insertion approximation breaks down for top quarks. (ii) The resonant enhancement of the source term and relaxation rate, that exists at leading order in the limit of degenerate thermal masses for the fermions/scalars, persists at NLO., Comment: 28 pages

The Herschel-ATLAS unbiased survey of cold dust in the local Universe is dominated by a surprising population of very blue (FUV-K < 3.5), dust-rich galaxies with high gas fractions (f_HI = M_HI/(M*+M_HI)>0.5)). Dubbed `Blue and Dusty Gas Rich Sources' (BADGRS) they have cold diffuse dust temperatures, and the highest dust-to-stellar mass ratios of any galaxies in the local Universe. Here, we explore the molecular ISM in a representative sample of BADGRS, using very deep CO(J_up=1,2,3) observations across the central and outer disk regions. We find very low CO brightnesses (Tp=15-30 mK), despite the bright far-infrared emission and metallicities in the range 0.5

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made., Comment: 13 pages, 6 figures, 2 tables, A&A accepted

In this paper, we use the Herschel Reference Survey to make a direct test of the hypothesis that the growth of a stellar bulge leads to a reduction in the star-formation efficiency of a galaxy (or conversely a growth in the gas-depletion timescale) as a result of the stabilisation of the gaseous disk by the gravitational field of the bulge. We find a strong correlation between star-formation efficiency and specific star-formation rate in galaxies without prominent bulges and in galaxies of the same morphological type, showing that there must be some other process besides the growth of a bulge that reduces the star-formation efficiency in galaxies. However, we also find that galaxies with more prominent bulges (Hubble types E to Sab) do have significantly lower star-formation efficiencies than galaxies with later morphological types, which is at least consistent with the hypothesis that the growth of a bulge leads to the reduction in the star-formation efficiency. The answer to the question in the title is therefore, yes and no: bulges may reduce the star-formation efficiency in galaxies but there must also be some other process at work. We also find that there is a significant but small difference in the star-formation efficiencies of galaxies with and without bars, in the sense that galaxies with bars have slightly higher star-formation efficiencies., Comment: Accepted for publication in MNRAS

The dust mass absorption coefficient, $\kappa_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $\kappa_{d}$ in galaxies. We determine $\kappa_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $\kappa_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it., Comment: Corrected typographical error in Equation A1, as per erratum accepted by MNRAS on 20th April 2022. No results or conclusions effected by the error, or by this correction

We study the post-inflation dynamics of multifield models involving nonminimal couplings using lattice simulations to capture significant nonlinear effects like backreaction and rescattering. We measure the effective equation of state and typical time-scales for the onset of thermalization, which could affect the usual mapping between predictions for primordial perturbation spectra and measurements of anisotropies in the cosmic microwave background radiation. For large values of the nonminimal coupling constants, we find efficient particle production that gives rise to nearly instantaneous preheating. Moreover, the strong single-field attractor behavior that was previously identified persists until the end of preheating, thereby suppressing typical signatures of multifield models. We therefore find that predictions for primordial observables in this class of models retain a close match to the latest observations., Comment: 6 pages, 7 figures. References added and minor edits made to match version to be published in Physical Review Letters

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment., Comment: Accepted for publication in Astronomy and Astrophysics

We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics., Comment: White paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics

Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate., Comment: 22 pages, 18 figures

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass., Comment: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

Galaxies continuously reprocess their interstellar material. One can therefore expect changing dust grain properties in galaxies which have followed different evolutionary pathways. Determining the intrinsic dust grain mix of a galaxy helps in reconstructing its evolutionary history. Early-type galaxies occasionally display regular dust lanes in their central regions. Due to the relatively simple geometry and composition of their stellar bodies, these galaxies are ideal to disentangle dust mix variations from geometric effects. We therefore model the various components of such a galaxy (FCC 167). We reconstruct its recent history, and investigate the possible fate of the dust lane. MUSE and ALMA observations reveal a nested ISM structure. An ionised-gas disk pervades the central regions of FCC 167, including those occupied by the main dust lane. Inward of the dust lane, we also find a disk/ring of cold molecular gas where stars are forming and HII regions contribute to the ionised-gas emission. Further in, the gas ionisation points towards an active galactic nucleus and the fuelling of a central supermassive black hole from its surrounding ionised and molecular reservoir. Observational constraints and radiative transfer models suggest the dust and gas are distributed in a ring-like geometry and the dust mix lacks small grains. The derived dust destruction timescales from sputtering in hot gas are short and we conclude that the dust must be strongly self-shielding and clumpy, or will quickly be eroded and disappear. Our findings show how detailed analysis of individual systems can complement statistical studies of dust-lane ETGs., Comment: 14 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

We investigate the role of leptons in electroweak baryogenesis by studying a relatively simple framework inspired by effective field theory that satisfies all Sakharov conditions. In particular, we study the effectiveness of CP-violating source terms induced by dimension-six Yukawa interactions for quarks and charged leptons. Despite the relatively small Yukawa coupling, CP-violating source terms involving taus are quite effective and can account for the observed matter-antimatter asymmetry. We obtain analytical and numerical expressions for the total baryon asymmetry, the former providing important insight into what makes lepton CP violation relatively effective compared to quark CP violation. Leptons also play an important role if the CP-violating source involves top quarks. While the tau Yukawa coupling in the Standard Model is small, it significantly enhances the baryon asymmetry by transferring the chiral asymmetry in quarks, which is washed out by strong sphalerons, to a chiral asymmetry in leptons. We conclude that leptons should not be ignored even if CP violation is limited to the quark sector. The role of leptons can be further increased in scenarios of new physics with additional chiral-symmetry-breaking interactions between quarks and leptons, as can happen in models with additional Higgs bosons or leptoquarks. Finally, we study CP-violating dimension-six Yukawa interactions for lighter quarks and leptons but conclude that these lead to too small baryon asymmetries., Comment: 36 pages

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500{\mu}m, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single S\'ersic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the S\'ersic index $n$, which characterises the shape of the S\'ersic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250{\mu}m is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.751.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed., Comment: 24 pages, 17 figures, Accepted for publication in Astronomy & Astrophysics

We study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of f_abs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of f_abs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light., Comment: 21 pages, 14 figures, A&A accepted. Template SEDs available at the DustPedia Archive: http://dustpedia.astro.noa.gr/

We perform an extensive analysis of linear fluctuations during preheating in Higgs inflation in the Einstein frame, where the fields are minimally coupled to gravity, but the field-space metric is nontrivial. The self-resonance of the Higgs and the Higgsed gauge bosons are governed by effective masses that scale differently with the nonminimal couplings and evolve differently in time. Coupled metric perturbations enhance Higgs self-resonance and make it possible for Higgs inflation to preheat solely through this channel. For $\xi\gtrsim 100$ the total energy of the Higgs-inflaton condensate can be transferred to Higgs particles within $3$ $e$-folds after the end of inflation. For smaller values of the nonminimal coupling preheating takes longer, completely shutting off at around $\xi\simeq 30$. The production of gauge bosons is dominated by the gauge boson mass and the field space curvature. For large values of the nonminimal coupling $\xi \gtrsim 1000$, it is possible for the Higgs condensate to transfer the entirety of its energy into gauge fields within one oscillation. For smaller values of the nonminimal coupling gauge bosons decay very quickly into fermions, thereby shutting off Bose enhancement. Estimates of non-Abelian interactions indicate that they will not suppress preheating into gauge bosons for $\xi \gtrsim 1000$., Comment: 57 pages, 15 figures

The galaxies found in optical surveys fall in two distinct regions of a diagram of optical colour versus absolute magnitude: the red sequence and the blue cloud with the green valley in between. We show that the galaxies found in a submillimetre survey have almost the opposite distribution in this diagram, forming a `green mountain'. We show that these distinctive distributions follow naturally from a single, continuous, curved Galaxy Sequence in a diagram of specific star-formation rate versus stellar mass without there being the need for a separate star-forming galaxy Main Sequence and region of passive galaxies. The cause of the red sequence and the blue cloud is the geometric mapping between stellar mass/specific star-formation rate and absolute magnitude/colour, which distorts a continuous Galaxy Sequence in the diagram of intrinsic properties into a bimodal distribution in the diagram of observed properties. The cause of the green mountain is Malmquist bias in the submillimetre waveband, with submillimetre surveys tending to select galaxies on the curve of the Galaxy Sequence, which have the highest ratios of submillimetre-to-optical luminosity. This effect, working in reverse, causes galaxies on the curve of the Galaxy Sequence to be underrepresented in optical samples, deepening the green valley. The green valley is therefore not evidence (1) for there being two distinct populations of galaxies, (2) for galaxies in this region evolving more quickly than galaxies in the blue cloud and the red sequence, (c) for rapid quenching processes in the galaxy population., Comment: Accepted by Monthly Notices of the Royal Astronomical Society