Hubble Time Research Articles

The Potential for Long-lived Intermediate-mass Black Hole Binaries in the Lowest Density Dwarf Galaxies

Intermediate-mass black hole (IMBH) mergers with masses 104–106 M ⊙ are expected to produce gravitational waves detectable by the Laser Interferometer Space Antenna (LISA) with high signal-to-noise ratios from the present day to cosmic dawn. IMBH mergers are expected to take place within dwarf galaxies; however, the dynamics, timescales, and effect on their hosts are largely unexplored. In a previous study, we examined how IMBHs would pair and merge within nucleated dwarf galaxies. IMBHs in nucleated hosts evolve very efficiently, forming a binary system and coalescing within a few hundred million years. Although the fraction of dwarf galaxies (107 M ⊙ ≤ M ⋆ ≤ 1010 M ⊙) hosting nuclear star clusters is between 60% and 100%, this fraction drops to 20%–70% for lower-mass dwarfs (M ⋆ ≈ 107 M ⊙), with the largest drop in low-density environments. Here, we extend our previous study by performing direct N-body simulations to explore the dynamics and evolution of IMBHs within nonnucleated dwarf galaxies, under the assumption that IMBHs exist within these dwarfs. To our surprise, none of the IMBHs in our simulation suite merge within a Hubble time, despite many attaining high eccentricities e ∼ 0.7–0.95. We conclude that extremely low stellar density environments in the centers of nonnucleated dwarfs do not provide an ample supply of stars to interact with an IMBH binary, resulting in its stalling, in spite of triaxiality and high eccentricity, common means to drive a binary to coalescence. Our findings underline the importance of considering all detailed host properties to predict IMBH merger rates for LISA.

Hidden Companions to Intermediate-mass Stars. Upgraded Multiplicity 2 → 3. XXIII. Discovery of a 0.94M ⊙, 9.2 au Companion to the Ap Star Gamma Piscis Austrini A* *Based on observations collected at the European Southern Observatory, Chile, Program ID 113.26GM.001.

Gamma Piscis Austrini is a nearby ρ = 4.″1 ↔ 284 au visual binary consisting of an A0Vp primary (A) and an F5V secondary (B). Here we report that the primary is actually a closer binary based on a VLTI/GRAVITY interferometric observation. The newly discovered companion Ab has a K band flux ratio 4.8% ± 0.2% and a projected separation ρ = 143.50 mas ↔ 9.2 au relative to Aa. We estimate isochrone masses M Aa = 2.65M ⊙, M Ab = 0.94M ⊙ and M B = 1.20M ⊙. γ PsA will become a triple white dwarf system within a Hubble time, with a low enough hierarchy for possibly interesting dynamical effects with a timescale P KL ∼ 200 kyr.

Explaining the oblate morphology of dwarf spheroidals with wave dark matter perturbations

ABSTRACT We investigate whether the oblate, spheroidal morphology of common dwarf spheroidal galaxies (dSph) may result from the slow relaxation of stellar orbits within a halo of wave dark matter ($\psi$DM) when starting from an initial disc of stars. Stellar orbits randomly walk over a Hubble time, perturbed by the pervasive ‘granular’ interference pattern of $\psi$DM, that fully modulates the dark matter density on the de Broglie scale. Our simulations quantify the level of stellar disc thickening over the Hubble time, showing that distribution of stars is predicted to become an oblate spheroid of increasing radius, that plausibly accounts for the morphology of dSph galaxies. We predict a low level of residual rotation remains after a Hubble time at the 1–3 km/s level, depending on orientation, that compares with recent claims of rotation for some well-studied local dSph galaxies. This steady internal dynamical evolution may be witnessed directly with JWST for well-resolved dwarf galaxies, appearing more oblate with look back time and tending to small discs of young stars at high redshift.

Primordial black holes from supercooled phase transitions

Cosmological first-order phase transitions (1stOPTs) are said to be strongly supercooled when the nucleation temperature is much smaller than the critical temperature. These are often encountered in theories that admit a nearly scale-invariant potential, for which the bounce action decreases only logarithmically with temperature. During supercooled 1stOPTs the equation of state of the universe undergoes a rapid and drastic change, transitioning from vacuum domination to radiation domination. The statistical variations in bubble nucleation histories imply that distinct causal patches percolate at slightly different times. Patches which percolate the latest undergo the longest vacuum-domination stage and as a consequence develop large overdensities triggering their collapse into primordial black holes (PBHs). We derive an analytical approximation for the probability of a patch to collapse into a PBH as a function of the 1stOPT duration, β−1, and deduce the expected PBH abundance. We find that 1stOPTs which take more than 15% of a Hubble time to complete (β/H≲7) produce observable PBHs. Their abundance is independent of the duration of the supercooling phase, in agreement with the de Sitter no hair conjecture. Published by the American Physical Society 2024

Cosmic purity lost: perturbative and resummed late-time inflationary decoherence

We compute the rate with which unobserved fields decohere other fields to which they couple, both in flat space and in de Sitter space, for spectator scalar fields prepared in their standard adiabatic vacuum. The process is very efficient in de Sitter space once the modes in question pass outside the Hubble scale, displaying the tell-tale phenomenon of secular growth that indicates the breakdown of perturbative methods on a time scale parameterically long compared with the Hubble time. We show how to match the perturbative evolution valid at early times onto a late-time Lindblad evolution whose domain of validity extends to much later times, thereby allowing a reliable resummation of the perturbative result beyond the perturbative regime. Super-Hubble modes turn out to be dominantly decohered by unobserved modes that are themselves also super-Hubble. If applied to curvature perturbations during inflation our observations here could close a potential loophole in recent calculations of the late-time purity of the observable primordial fluctuations.

Evolution and final fate of massive post-common-envelope binaries

Mergers of neutron stars (NSs) and black holes (BHs) are nowadays observed routinely thanks to gravitational-wave (GW) astronomy. In the isolated binary-evolution channel, a common-envelope (CE) phase of a red supergiant (RSG) and a compact object is crucial to sufficiently shrink the orbit and thereby enable a merger via GW emission. Here, we use the outcomes of two three-dimensional (3D) magneto-hydrodynamic CE simulations of an initially 10.0 solar-mass RSG with a 5.0 solar-mass BH and a 1.4 solar-mass NS, respectively, to explore the further evolution and final fate of the remnant binaries (post-CE binaries). Notably, the 3D simulations reveal that the post-CE binaries are likely surrounded by circumbinary disks (CBDs), which contain substantial mass and angular momentum to influence the subsequent evolution. The binary systems in MESA modelling undergo another phase of mass transfer and we find that most donor stars do not explode in ultra-stripped supernovae (SNe), but rather in Type Ib/c SNe. Without NS kicks, the final orbital configurations of our models with the BH companion are too wide to allow for a compact object merger within a Hubble time. NS kicks are actually required to sufficiently perturb the orbit and thus facilitate a merger via GW emission. Moreover, we explore the influence of CBDs observed in 3D CE simulations on the evolution and final fate of the post-CE binaries. We find that mass accretion from the disk widens the binary orbit, while resonant interactions between the CBD and the binary can shrink the separation and increase the eccentricity of the binary depending on the disk mass and lifetime. Efficient resonant contractions may even enable a BH or NS to merge with the remnant He stars before a second SN explosion, which may be observed as gamma-ray burst-like transients, luminous fast blue optical transients, and Thorne-Żytkow objects. For the surviving post-CE binaries, the CBD-binary interactions may significantly increase the GW-induced double compact merger fraction. We conclude that accounting for CBD may be crucial to better understand observed GW mergers.

Observational constraints on FLRW, Bianchi type I and V brane models

This study explores the compatibility of Covariant Extrinsic Gravity (CEG), a braneworld scenario with an arbitrary number of non-compact extra dimensions, with current cosmological observations. We employ the chi-square statistic and Markov Chain Monte Carlo (MCMC) methods to fit the Friedmann–Lemaître–Robertson–Walker (FLRW) and Bianchi type-I and V brane models to the latest datasets, including Hubble, Pantheon+ Supernova samples, Big Bang Nucleosynthesis (BBN), Baryon Acoustic Oscillations (BAO), and the structure growth rate, fσ8(z). Parameters for FLRW universe consist Ω0(b),Ω0(cd),Ω0(k),H0,γ,σ8, while for the Bianchi model are Ω0(b),Ω0(cd),Ω0(β),H0,γ,Ω0(θ),σ8. By comparing our models to observational data, we determine the best values for cosmological parameters. For the FLRW model, these values depend on the sign of γ (which gives the time variation of gravitational constant in Hubble time unit): γ>0 yields γ=0.00008−0.00011+0.00015, and Ω0(k)=0.014−0.022+0.024 and γ<0 leads to γ=−0.0226−0.0062+0.0054, and Ω0(k)=0.023−0.041+0.039. It should be noted that in both cases Ω0(k)>0, which represents a closed universe. Similarly, for the Bianchi type-V brane model, the parameter values vary with the sign of γ, resulting in γ=0.00084−0.00021+0.00019, Ω0(β)=0.0258−0.0063+0.0052, and Ω0θ(×10−5)=4.19−0.75+0.67 (as with the density parameter of stiff matter) for γ>0, and γ=−0.00107−0.00020+0.00019, Ω0(β)=0.0259−0.0062+0.0050, and Ω0θ(×10−5)=4.17−0.98+0.91 for γ<0. In both cases Ω0(β)>0, which represents the Bianchi type-V, because in the Bianchi type-I, β=0. Subsequently, utilizing these obtained best values, we analyze the behavior of key cosmological parameters such as the Hubble parameter, deceleration parameter, distance modulus, equation of state, and density parameters that characterize both matter and the geometric component of dark energy, as functions of redshift. Our results notably show that the FLRW model with γ<0 is more compatible with observational data than the Bianchi model, based on various statistical criteria.

The torsion of stellar streams and the overall shape of galactic gravity’s source

Context. Flat rotation curves, v(r), are naturally explained by elongated (prolate) dark matter (DM) distributions, and we have provided competitive fits to the SPARC database. To further probe the geometry of the halo, or the equivalent source of gravity in other formulations, one needs observables outside the galactic plane. Stellar streams, poetically analogous to airplane contrails, but caused by tidal dispersion of massive substructures such as satellite dwarf galaxies, would lie on their own plane (consistently with angular momentum conservation) should the DM-halo gravitational field be spherically symmetric. Tracks resembling entire orbits are seldom available because their periods are commensurable with Hubble time, with streams often presenting themselves as short segments. Aims. Therefore, we aim to establish stellar stream torsion, a local observable that measures the deviation from planarity in differential curve geometry, as a diagnostic providing sensitivity to aspherical DM distributions and ensuring the use of even relatively short streams. Methods. We performed small-scale simulations of tidally distorted star clusters to check that indeed a central force center produces negligible torsion, while distorted halos can generate it. Turning to observational data, we identified among the known streams those that are at the largest distance from the Galactic center, and that are likely not affected by the Magellanic clouds, as the most promising for the study, and by means of polynomial fits we extracted their differential torsion. Results. We find that the torsion of the few known streams that should be sensitive to most of the Milky Way’s DM halo is much larger than expected for a central spherical bulb alone. This is consistent with the nonsphericity of the halo. Conclusions. Future studies of stellar stream torsion with larger samples and further out of the galactic plane should be able to extract the ellipticity of the halo to see whether it is just a slight distortion of a spherical shape or whether it rather resembles a more elongated cigar.

Fluctuation-dissipation relation in cosmic microwave background

We study the fluctuation-dissipation relation for sound waves in the cosmic microwave background (CMB), employing effective field theory (EFT) for fluctuating hydrodynamics.Treating sound waves as the linear response to thermal radiation, we establish the fluctuation-dissipation relation within a cosmological framework.While dissipation is elucidated in established linear cosmological perturbation theory, the standard Boltzmann theory overlooks the associated noise, possibly contributing to inconsistencies in Lambda Cold Dark Matter (ΛCDM) cosmology. This paper employs EFT for fluctuating hydrodynamics in cosmological perturbation theory, deriving sound wave noise.Notably, the long-time limit of the noise spectrum is independent of viscosity details, resembling a Brownian motion bounded in a harmonic potential.The net energy transfer between the sound wave system and the radiation environment reaches a balance within Hubble time, suggesting the thermal equilibrium of the sound waves themselves.The induced density power spectrum is characterized as white noise dependent on the inverse of the entropy density, which is negligibly small on the CMB scale. The energy density of the entire sound wave system scales as a -4, akin to radiation. While the numerical factor is not determined in the present calculation, the back reaction of the sound wave system to the background radiation may not be negligible, serving as a potential source for various fitting issues in ΛCDM cosmology.

The double low-mass white dwarf eclipsing binary system J2102–4145 and its possible evolution

In recent years, about 150 low-mass white dwarfs (WDs), typically with masses below 0.4 M⊙, have been discovered. The majority of these low-mass WDs are observed in binary systems as they cannot be formed through single-star evolution within Hubble time. In this work, we present a comprehensive analysis of the double low-mass WD eclipsing binary system J2102−4145. Our investigation encompasses an extensive observational campaign, resulting in the acquisition of approximately 28 h of high-speed photometric data across multiple nights using NTT/ULTRACAM, SOAR/Goodman, and SMARTS-1m telescopes. These observations have provided critical insights into the orbital characteristics of this system, including parameters such as inclination and orbital period. To disentangle the binary components of J2102−4145, we employed the XTGRID spectral fitting method with GMOS/Gemini-South and X-shooter data. Additionally, we used the PHOEBE package for light curve analysis on NTT/ULTRACAM high-speed time-series photometry data to constrain the binary star properties. Our analysis unveils remarkable similarities between the two components of this binary system. For the primary star, we determine Teff,1 = 13 688−72+65 K, log g1 = 7.36 ± 0.01, R1 = 0.0211 ± 0.0002 R⊙, and M1 = 0.375 ± 0.003 M⊙, while, the secondary star is characterised by Teff,2 = 12952−66+53 K, log g2 = 7.32 ± 0.01, R2 = 0.0203−0.0003+0.0002 R⊙, and M2 = 0.314 ± 0.003 M⊙. Furthermore, we found a notable discrepancy between Teff and R of the less massive WD, compared to evolutionary sequences for WDs from the literature, which has significant implications for our understanding of WD evolution. We discuss a potential formation scenario for this system which might explain this discrepancy and explore its future evolution. We predict that this system will merge in ∼800 Myr, evolving into a helium-rich hot subdwarf star and later into a hybrid He/CO WD.

Investigating the Chemically Homogeneous Evolution Channel and Its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980

Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf–Rayet binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries, then use a Bayesian approach to compare them with the stars’ observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ∼3 days, widening to the present-day ∼20 days orbit as a result of mass loss—correspondingly, they have very high initial stellar masses (≳150 M ⊙). We explore variations in stellar-wind mass loss and internal mixing efficiency, and find that models assuming enhanced mass loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favored over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ∼19–37 M ⊙. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source.

White Dwarf Merger Remnants: The DAQ Subclass

Four years after the discovery of a unique DAQ white dwarf with a hydrogen-dominated and carbon-rich atmosphere, we report the discovery of four new DAQ white dwarfs, including two that were not recognized properly in the literature. We find all five DAQs in a relatively narrow mass and temperature range of M = 1.14–1.19 M ⊙ and T eff = 13,000–17,000 K. In addition, at least two show photometric variations due to rapid rotation with ≈10 minute periods. All five are also kinematically old, but appear photometrically young, with estimated cooling ages of about 1 Gyr based on standard cooling tracks, and their masses are roughly twice the mass of the most common white dwarfs in the solar neighborhood. These characteristics are smoking gun signatures of white dwarf merger remnants. Comparing the DAQ sample with warm DQ white dwarfs, we demonstrate that there is a range of hydrogen abundances among the warm DQ population and that the distinction between DAQ and warm DQ white dwarfs is superficial. We discuss the potential evolutionary channels for the emergence of the DAQ subclass, suggesting that DAQ white dwarfs are trapped on the crystallization sequence and may remain there for a significant fraction of the Hubble time.

Dynamical Analysis of the Redshift Drift in FLRW Universes

Redshift drift is the phenomenon whereby the observed redshift between an emitter and observer comoving with the Hubble flow in an expanding FLRW universe will slowly evolve—on a timescale comparable to the Hubble time. In a previous article, three of the current authors performed a cosmographic analysis of the redshift drift in an FLRW universe, temporarily putting aside the issue of dynamics (the Friedmann equations). In the current article, we add dynamics while still remaining within the framework of an exact FLRW universe. We developed a suitable generic matter model and applied it to both standard FLRW and various dark energy models. Furthermore, we present an analysis of the utility of alternative cosmographic variables to describe the redshift drift data.

The abundance of core-collapsed subhaloes in SIDM: insights from structure formation in ΛCDM

ABSTRACT Dark matter haloes enter a phase of gravothermal core-collapse in the presence of self-interactions. This phase is expected to be subdominant given the long time-scales involved. However, it has been shown that collapse is accelerated in subhaloes in strong tidal environments. Cosmological simulations in ΛCDM give complete distributions of subhalo orbits and their profiles. We use these to estimate the fraction of the subhaloes in different host halo environments that are in the core-collapse phase. Our hosts range from the Large Magellanic cloud (LMC)-like hosts to clusters. We use fluid simulations of self-interacting dark matter (SIDM) to evolve subhaloes in their hosts, including the effect of tidal truncation at the time of pericentre crossing. We find, for parameters that allow the interaction cross-section to be high at dwarf scales, at least $10\%$ of all subhaloes are expected to have intrinsically collapsed within Hubble time up to the group mass host scales. This fraction increases significantly, becoming at least 20% when tidal interactions are considered. To identify these objects in observations, we find that we need to measure their densities at very small radial scales, where the subhaloes show a bimodal distribution of densities or we need to measure the slopes of their inner density profiles near the scale radius, which are much steeper than NFW slopes. Current measurements of central slopes of classical dwarfs do not show preference for being completely collapsed, however they are consistent with being cored or in a phase of transition between expansion and collapse.

The X-ray binaries in M83: Will any of them form gravitational wave sources for LIGO-VIRGO-KAGRA?

There are 214 X-ray point sources (LX &gt; 1035 erg s−1) identified as X-ray binaries (XRBs) in the nearby spiral galaxy M83. Since XRBs are powered by accretion onto a neutron star (NS) or a black hole (BH) from a companion or donor star, these systems are promising progenitors of merging double compact objects (DCOs): BH-BH, BH-NS, or NS-NS systems. The connection (i.e., XRBs evolving into DCOs) may provide some hints to the as-yet-unanswered question: what is the origin of the LIGO, Virgo, and KAGRA mergers? Available observations do not allow us to determine what the final fate of the XRBs observed in M83 will be. However, we can use an evolutionary model of isolated binaries to reproduce the population of XRBs in M83 by matching model XRB numbers, types, and luminosities to observations. Knowing the detailed properties of M83 model XRBs (donor and accretor masses, and their evolutionary ages and orbits), we follow their evolution to the deaths of donor stars to check whether any merging DCOs are formed. Although all merging DCOs in our isolated binary evolution model go through the XRB phase (defined as reaching X-ray luminosity from RLOF or wind accretion onto NSs or BHs above 1035 erg s−1), only very few XRBs evolve to form merging (in Hubble time) DCOs. For M83, with its solar-like metallicity stars and continuous star formation, we find that only ∼1 − 2% of model XRBs evolve into merging DCOs depending on the adopted evolutionary physics. This is caused by (i) the merger of the donor star with a compact object during the common envelope phase, (ii) a binary disruption at the supernova explosion of a donor star, (iii) the formation of a DCO on a wide orbit (merger time longer than Hubble time).

Are there any extragalactic high speed dark matter particles in the Solar neighborhood?

We use the APOSTLE suite of cosmological hydrodynamical simulations of the Local Group to examine the high speed tail of the local dark matter velocity distribution in simulated Milky Way analogues. The velocity distribution in the Solar neighborhood is well approximated by a generalized Maxwellian distribution sharply truncated at a well-defined maximum “escape” speed. The truncated generalized Maxwellian distribution accurately models the local dark matter velocity distribution of all our Milky Way analogues, with no evidence for any separate extragalactic high-speed components. The local maximum speed is well approximated by the terminal velocity expected for particles able to reach the Solar neighborhood in a Hubble time from the farthest confines of the Local Group. This timing constraint means that the local dark matter velocity distribution is unlikely to contain any high-speed particles contributed by the Virgo Supercluster “envelope”, as argued in recent work. Particles in the Solar neighborhood with speeds close to the local maximum speed can reach well outside the virial radius of the Galaxy, and, in that sense, belong to the Local Group envelope posited in earlier work. The local manifestation of such envelope is thus not a distinct high-speed component, but rather simply the high-speed tail of the truncated Maxwellian distribution.

Non-perturbative wavefunction of the universe in inflation with (resonant) features

We study the statistics of scalar perturbations in models of inflation with small and rapid oscillations in the inflaton potential (resonant non-Gaussianity). We do so by deriving the wavefunction Ψ[ζ(x)] non-perturbatively in ζ, but at first order in the amplitude of the oscillations. The expression of the wavefunction of the universe (WFU) is explicit and does not require solving partial differential equations. One finds qualitative deviations from perturbation theory for |ζ| ≳ α−2, where α ≫ 1 is the number of oscillations per Hubble time. Notably, the WFU exhibits distinct behaviours for negative and positive values of ζ (troughs and peaks respectively). While corrections for ζ < 0 remain relatively small, of the order of the oscillation amplitude, positive ζ yields substantial effects, growing exponentially as eπα/2 in the limit of large ζ. This indicates that even minute oscillations give large effects on the tail of the distribution.

Formation of PSR J1012+5307 with an Extremely Low-mass White Dwarf: Testing Magnetic Braking Models

PSR J1012+5307 is a millisecond pulsar with an extremely low-mass white dwarf (WD) companion in an orbit of 14.5 hr. Magnetic braking (MB) plays an important role in influencing the orbital evolution of binary systems with a low-mass (≲1–2 M ⊙) donor star. At present, there exist several different MB descriptions. In this paper, we investigate the formation of PSR J1012+5307 as a probe to test plausible MB models. Employing a detailed stellar evolution model with the MESA code, we find that the convection- and rotation-boosted MB and “Intermediate” MB models can reproduce the WD mass, WD radius, WD surface gravity, neutron star mass, and orbital period observed in PSR J1012+5307. However, our simulated WD has a higher effective temperature than the observations. Three other MB mechanisms, including the standard MB model, are too weak to account for the observed orbital period in the Hubble time. A long cooling timescale caused by H-shell flashes of the WD may alleviate the discrepancy between the simulated effective temperature and the observed value.

Chaotic Diffusion in a Triaxial Galactic Model: an Example of Global Stable Chaos

In this work we focus on the chaotic diffusion in the phase space of a triaxial potential resembling an elliptical galaxy. The transport process is studied in two different action-like starting spaces in order to cope with circulating and noncirculating orbits. Estimates of the diffusion rate obtained by means of the variance approach are discussed in detail and their limitations are exposed. After revisiting the Shannon-entropy-based method from a conceptual point of view in the framework of simple arguments taken from the information theory, we apply it to measure changes in the unperturbed actions or integrals of motion of the system for different sets of small ensembles of random initial conditions. For such sets of ensembles, estimates of the Lyapunov times are also provided. The results show that, within the chaotic component of the phase space, the Lyapunov times are shorter than any physical time scale as the Hubble time, but the diffusion times are much larger than the latter. Thus, we conclude that stable chaos dominates the dynamics of realistic galactic models.

Eternal inflation and collapse theories

The eternal inflation problem continues to be considered one of standard's cosmology most serious shortcomings. This arises when one considers the effects of “quantum fluctuations” on the zero mode of inflaton field during a Hubble time in the inflationary epoch.In the slow-roll regime it is quite clear that such quantum fluctuations could dwarf the classical rolling down of the inflaton, and with overwhelming probability this prevents inflation from ever ending. When one recognizes that quantum fluctuations can not be taken as synonymous of stochastic fluctuations, but rather intrinsic levels of indefiniteness in the quantities in question, one concludes that the eternal inflation problem simply does not exist. However, the same argument would serve to invalidate the account for the generation of the primordial seeds of cosmic structure as has been amply discussed elsewhere [1,2,3]. In order to do address that issue, one must explain the breaking of homogeneity and isotropy of the situation prevailing during the early inflationary epoch (at both the quantum and classical levels of the description). For that one needs to rely on some additional element, beyond those present in the traditional treatments. The so called spontaneous collapse theories offer a viable candidate for that element, namely the stochastic and spontaneous state reduction characteristic of those proposals possesses the basic features to break those symmetries. In fact, a version of the CSL theory adapted to the cosmological context has been shown to offer a satisfactory account for the origin the seeds of cosmic structure with an adequate power spectrum [2], and will serve as the basis of our analysis. However, once such stochastic collapse is introduced into the theoretical framework the eternal inflation problem has the potential reappear.In this manuscript we explore those issues in detail and discuss an avenue that seems to allow for a satisfactory account for the generation of the primordial inhomogeneities and anisotropies while freeing the theory from the eternal inflation problem.