Source author record

David J. Weir

David J. Weir appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2023arXiv

Droplet collapse during strongly supercooled transitions

We simulate the decay of isolated, spherically symmetric droplets in a cosmological phase transition. It has long been posited that such heated droplets of the metastable state could form, and they have recently been observed in 3D multi-bubble simulations. In those simulations, the droplets were associated with a reduction in the wall velocity and a decrease in the kinetic energy of the fluid, with a consequent suppression in the gravitational wave power spectrum. In the present work, we track the wall speed and kinetic energy production in isolated droplets and compare them to those found in multi-bubble collisions. The late-time wall velocities that we observe match those of the 3D simulations, though we find that the spherical simulations are a poor predictor of the kinetic energy production. This implies that spherically symmetric simulations could be used to refine baryogenesis predictions due to the formation of droplets, but not to estimate any accompanying suppression of the gravitational wave signal.

preprint2022arXiv

Generation of gravitational waves from freely decaying turbulence

We study the stochastic gravitational wave background (SGWB) produced by freely decaying vortical turbulence in the early Universe. We thoroughly investigate the time correlation of the velocity field, and hence of the anisotropic stresses producing the gravitational waves. With hydrodynamical simulations, we show that the unequal time correlation function (UETC) of the Fourier components of the velocity field is Gaussian in the time difference, as predicted by the "sweeping" decorrelation model. We introduce a decorrelation model that can be extended to wavelengths around the integral scale of the flow. Supplemented with the evolution laws of the kinetic energy and of the integral scale, this provides a new model UETC of the turbulent velocity field consistent with the simulations. We discuss the UETC as a positive definite kernel, and propose to use the Gibbs kernel for the velocity UETC as a natural way to ensure positive definiteness of the SGWB. The SGWB is given by a 4-dimensional integration of the resulting anisotropic stress UETC with the gravitational wave Green's function. We perform this integration using a Monte Carlo algorithm based on importance sampling, and find that the result matches that of the simulations. Furthermore, the SGWB obtained from the numerical integration and from the simulations show close agreement with a model in which the source is constant in time and abruptly turns off after a few eddy turnover times. Based on this assumption, we provide an approximate analytical form for the SGWB spectrum and its scaling with the initial kinetic energy and integral scale. Finally, we use our model and numerical integration algorithm to show that including an initial growth phase for the turbulent flow heavily influences the spectral shape of the SGWB. This highlights the importance of a complete understanding of the turbulence generation mechanism.

preprint2021arXiv

Detecting gravitational waves from cosmological phase transitions with LISA: an update

We investigate the potential for observing gravitational waves from cosmological phase transitions with LISA in light of recent theoretical and experimental developments. Our analysis is based on current state-of-the-art simulations of sound waves in the cosmic fluid after the phase transition completes. We discuss the various sources of gravitational radiation, the underlying parameters describing the phase transition and a variety of viable particle physics models in this context, clarifying common misconceptions that appear in the literature and identifying open questions requiring future study. We also present a web-based tool, PTPlot, that allows users to obtain up-to-date detection prospects for a given set of phase transition parameters at LISA.

preprint2021arXiv

Gravitational waves from vacuum first order phase transitions II: from thin to thick walls

In a vacuum first-order phase transition, gravitational waves are generated from collision of bubbles of the true vacuum. The spectrum from such collisions takes the form of a broken power law. We consider a toy model for such a phase transition, where the dynamics of the scalar field depends on a single parameter $\overlineλ$, which controls how thin the bubble wall is at nucleation and how close to degenerate the vacua are relative to the barrier. We extend on our previous work by performing a series of simulations with a range of $\overlineλ$. The peak of the gravitational-wave power spectrum varies by up to a factor of $1.3$, which is probably an unobservable effect. We find that the ultraviolet (UV) power law in the gravitational-wave spectrum becomes steeper as $\overlineλ \rightarrow 0$, varying between $k^{-1.4}$ and $k^{-2.2}$ for the $\overlineλ$ considered. This provides some evidence that the form of the underlying effective potential of a vacuum first-order phase transition could be determined from the gravitational-wave spectrum it produces.

preprint2020arXiv

Shape of the acoustic gravitational wave power spectrum from a first order phase transition

We present results from large-scale numerical simulations of a first order thermal phase transition in the early universe, in order to explore the shape of the acoustic gravitational wave and the velocity power spectra. We compare the results with the predictions of the recently proposed sound shell model. For the gravitational wave power spectrum, we find that the predicted $k^{-3}$ behaviour, where $k$ is the wavenumber, emerges clearly for detonations. The power spectra from deflagrations show similar features, but exhibit a steeper high-$k$ decay and an extra feature not accounted for in the model. There are two independent length scales: the mean bubble separation and the thickness of the sound shell around the expanding bubble of the low temperature phase. It is the sound shell thickness which sets the position of the peak of the power spectrum. The low wavenumber behaviour of the velocity power spectrum is consistent with a causal $k^{3}$, except for the thinnest sound shell, where it is steeper. We present parameters for a simple broken power law fit to the gravitational wave power spectrum for wall speeds well away from the speed of sound where this form can be usefully applied. We examine the prospects for the detection, showing that a LISA-like mission has the sensitivity to detect a gravitational wave signal from sound waves with an RMS fluid velocity of about $0.05c$, produced from bubbles with a mean separation of about $10^{-2}$ of the Hubble radius. The shape of the gravitational wave power spectrum depends on the bubble wall speed, and it may be possible to estimate the wall speed, and constrain other phase transition parameters, with an accurate measurement of a stochastic gravitational wave background.

preprint2020arXiv

Vorticity, kinetic energy, and suppressed gravitational wave production in strong first order phase transitions

We have performed the first 3-dimensional simulations of strong first-order thermal phase transitions in the early Universe. For deflagrations, we find that the rotational component of the fluid velocity increases as the transition strength is increased. For detonations, however, the rotational velocity component remains constant and small. We also find that the efficiency with which kinetic energy is transferred to the fluid falls below theoretical expectations as we increase the transition strength. The probable origin of the kinetic energy deficit is the formation of reheated droplets of the metastable phase during the collision, slowing the bubble walls. The rate of increase in the gravitational wave energy density for deflagrations in strong transitions is suppressed compared to that predicted in earlier work. This is largely accounted for by the reduction in kinetic energy. Current modelling therefore substantially overestimates the gravitational wave signal for strong transitions with deflagrations, in the most extreme case by a factor of $10^{3}$. Detonations are less affected.

preprint2016arXiv

Improving cosmic string network simulations

In real-time lattice simulations of cosmic strings in the Abelian Higgs model, the broken translational invariance introduces lattice artefacts; relativistic strings therefore decelerate and radiate. We introduce two different methods to construct a moving string on the lattice, and study in detail the lattice effects on moving strings. We find that there are two types of lattice artefact: there is an effective maximum speed with which a moving string can be placed on the lattice, and a moving string also slows down, with the deceleration approximately proportional to the exponential of the velocity. To mitigate this, we introduce and study an improved discretisation, based on the tree-level Lüscher-Weisz action, which is found to reduce the deceleration by an order of magnitude, and to increase the string speed limit by an amount equivalent to halving the lattice spacing. The improved algorithm is expected to be very useful for 3D simulations of cosmic strings in the early universe, where one wishes to simulate as large a volume as possible.

preprint2016arXiv

Lattice Calculation of the Decay of Primordial Higgs Condensate

We study the resonant decay of the primordial Standard Model Higgs condensate after inflation into $SU(2)$ gauge bosons on the lattice. We find that the non-Abelian interactions between the gauge bosons quickly extend the momentum distribution towards high values, efficiently destroying the condensate after the onset of backreaction. For the inflationary scale $H = 10^8$ GeV, we find that 90% of the Higgs condensate has decayed after $n \sim 10$ oscillation cycles. This differs significantly from the Abelian case where, given the same coupling strengths, most of the condensate would persist after the resonance.

preprint2016arXiv

New solutions for non-Abelian cosmic strings

We study the properties of classical vortex solutions in a non-Abelian gauge theory. A system of two adjoint Higgs fields breaks the SU(2) gauge symmetry to $Z_2$, producing 't Hooft-Polyakov monopoles trapped on cosmic strings, termed beads; there are two charges of monopole and two degenerate string solutions. The strings break an accidental discrete $Z_2$ symmetry of the theory, explaining the degeneracy of the ground state. Further symmetries of the model, not previously appreciated, emerge when the masses of the two adjoint Higgs fields are degenerate. The breaking of the enlarged discrete symmetry gives rise to additional string solutions and splits the monopoles into four types of `semipole': kink solutions that interpolate between the string solutions, classified by a complex gauge invariant magnetic flux and a $Z_4$ charge. At special values of the Higgs self-couplings, the accidental symmetry broken by the string is continuous, giving rise to supercurrents on the strings. The SU(2) theory can be embedded in a wide class of Grand Unified Theories, including SO(10). We argue that semipoles and supercurrents are generic on GUT strings.

preprint2016arXiv

Numerical simulations of acoustically generated gravitational waves at a first order phase transition

We present details of numerical simulations of the gravitational radiation produced by a first order thermal phase transition in the early universe. We confirm that the dominant source of gravitational waves is sound waves generated by the expanding bubbles of the low-temperature phase. We demonstrate that the sound waves have a power spectrum with a power-law form between the scales set by the average bubble separation (which sets the length scale of the fluid flow $L_\text{f}$) and the bubble wall width. The sound waves generate gravitational waves whose power spectrum also has a power-law form, at a rate proportional to $L_\text{f}$ and the square of the fluid kinetic energy density. We identify a dimensionless parameter $\tildeΩ_\text{GW}$ characterising the efficiency of this "acoustic" gravitational wave production whose value is $8π\tildeΩ_\text{GW} \simeq 0.8 \pm 0.1$ across all our simulations. We compare the acoustic gravitational waves with the standard prediction from the envelope approximation. Not only is the power spectrum steeper (apart from an initial transient) but the gravitational wave energy density is generically larger by the ratio of the Hubble time to the phase transition duration, which can be 2 orders of magnitude or more in a typical first order electroweak phase transition.

preprint2016arXiv

Revisiting the envelope approximation: gravitational waves from bubble collisions

We study the envelope approximation and its applicability to first-order phase transitions in the early universe. We demonstrate that the power laws seen in previous studies exist independently of the nucleation rate. We also compare the envelope approximation prediction to results from large-scale phase transition simulations. For phase transitions where the contribution to gravitational waves from scalar fields dominates over that from the coupled plasma of light particles, the envelope approximation is in agreement, giving a power spectrum of the same form and order of magnitude. In all other cases the form and amplitude of the gravitational wave power spectrum is markedly different and new techniques are required.

preprint2016arXiv

Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions

We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational waves. The projected sensitivity of eLISA to cosmological phase transitions is computed in a model-independent way for various detector designs and configurations. By applying these results to several specific models, we demonstrate that eLISA is able to probe many well-motivated scenarios beyond the Standard Model of particle physics predicting strong first-order cosmological phase transitions in the early Universe.

preprint2014arXiv

Gravitational waves from the sound of a first order phase transition

We report on the first 3-dimensional numerical simulations of first-order phase transitions in the early universe to include the cosmic fluid as well as the scalar field order parameter. We calculate the gravitational wave (GW) spectrum resulting from the nucleation, expansion and collision of bubbles of the low-temperature phase, for phase transition strengths and bubble wall velocities covering many cases of interest. We find that the compression waves in the fluid continue to be a source of GWs long after the bubbles have merged, a new effect not taken properly into account in previous modelling of the GW source. For a wide range of models the main source of the GWs produced by a phase transition is therefore the sound the bubbles make.

preprint2014arXiv

On the quantum stability of Q-balls

We consider the evolution and decay of Q-balls under the influence of quantum fluctuations. We argue that the most important effect resulting from these fluctuations is the modification of the effective potential in which the Q-ball evolves. This is in addition to spontaneous decay into elementary particle excitations and fission into smaller Q-balls previously considered in the literature, which -- like most tunnelling processes -- are likely to be strongly suppressed. We illustrate the effect of quantum fluctuations in a particular model $ϕ^6$ potential, for which we implement the inhomogeneous Hartree approximation to quantum dynamics and solve for the evolution of Q-balls in 3+1 dimensions. We find that the stability range as a function of (field space) angular velocity $ω$ is modified significantly compared to the classical case, so that small-$ω$ Q-balls are less stable than in the classical limit, and large-$ω$ Q-balls are more stable. This can be understood qualitatively in a simple way.

preprint2013arXiv

Form factor and width of a quantum string

In the Yang-Mills theory, the apparent thickness of the confining string is known to grow logarithmically when its length increases. The same logarithmic broadening also happens to strings in other quantum field theories and domain walls in statistical physics models. Even in quantum field theories, the correlators used to measure and characterise this phenomenon are analogous to those in statistical physics. In this paper we describe it using the string form factor which is a meaningful quantum observable, obtainable in principle from scattering experiments. We show how the form factor can be obtained from field correlation functions calculated in lattice Monte Carlo simulations. We apply this method to 2+1-dimensional scalar theory in the strong coupling limit, where it is equivalent to the 3D Ising model, and through duality also to 2+1-dimensional $\mathbb{Z}_2$ gauge theory. We measure the string form factor by simulating the Ising model, and demonstrate that it displays the same logarithmic broadening as observed by other quantities.

preprint2013arXiv

Investigating the Sharpe-Singleton scenario on the lattice by direct eigenvalue computation

We investigate the phase structure of lattice QCD with dynamical Wilson fermions. Wilson chiral perturbation theory predicts that the Aoki phase and the Sharpe-Singleton scenario manifest themselves in very distinct behavior of the Wilson Dirac eigenvalue spectrum. To test this prediction we perform a direct calculation of the eigenvalues of the non-Hermitian Wilson Dirac operator in dynamical lattice simulations. Moreover, we demonstrate an unexpected quark mass dependence on the shape of the eigenvalue distribution in the positive quark mass side.

preprint2012arXiv

Nonperturbative study of the 't Hooft-Polyakov monopole form factors

The mass and interactions of a quantum 't Hooft-Polyakov monopole are measured nonperturbatively using correlation functions in lattice Monte Carlo simulations. A method of measuring the form factors for interactions between the monopole and fundamental particles, such as the photon, is demonstrated. These quantities are potentially of experimental relevance in searches for magnetic monopoles.

preprint2012arXiv

PT phase transition in multidimensional quantum systems

Non-Hermitian PT-symmetric quantum-mechanical Hamiltonians generally exhibit a phase transition that separates two parametric regions, (i) a region of unbroken PT symmetry in which the eigenvalues are all real, and (ii) a region of broken PT symmetry in which some of the eigenvalues are complex. This transition has recently been observed experimentally in a variety of physical systems. Until now, theoretical studies of the PT phase transition have generally been limited to one-dimensional models. Here, four nontrivial coupled PT-symmetric Hamiltonians, $H=p^2/2+x^2/2+q^2/2+y^2/2+igx^2y$, $H=p^2/2+x^2/2+q^2/2+y^2+igx^2y$, $H=p^2/2+x^2/2+q^2/2+y^2/2+r^2/2+z^2/2+igxyz$, and $H=p^2/2+x^2/2+q^2/2+y^2+r^2/2+3z^2/2+igxyz$ are examined. Based on extensive numerical studies, this paper conjectures that all four models exhibit a phase transition. The transitions are found to occur at $g\approx 0.1$, $g\approx 0.04$, $g\approx 0.1$, and $g\approx 0.05$. These results suggest that the PT phase transition is a robust phenomenon not limited to systems having one degree of freedom.

preprint2012arXiv

Variational matrix product ansatz for dispersion relations

A variational ansatz for momentum eigenstates of translation invariant quantum spin chains is formulated. The matrix product state ansatz works directly in the thermodynamic limit and allows for an efficient implementation (cubic scaling in the bond dimension) of the variational principle. Unlike previous approaches, the ansatz includes topologically non-trivial states (kinks, domain walls) for systems with symmetry breaking. The method is benchmarked using the spin-1/2 XXZ antiferromagnet and the spin-1 Heisenberg antiferromagnet and we obtain surprisingly accurate results.

preprint2011arXiv

Fluxoid formation: size effects and non-equilibrium universality

Simple causal arguments put forward by Kibble and Zurek suggest that the scaling behaviour of condensed matter at continuous transitions is related to the familiar universality classes of the systems at quasi-equilibrium. Although proposed 25 years ago or more, it is only in the last few years that it has been possible to devise experiments from which scaling exponents can be determined and in which this scenario can be tested. In previous work, an unusually high Kibble-Zurek scaling exponent was reported for spontaneous fluxoid production in a single isolated superconducting Nb loop, albeit with low density. Using analytic approximations backed up by Langevin simulations, we argue that densities as small as these are too low to be attributable to scaling, and are conditioned by the small size of the loop. We also reflect on the physical differences between slow quenches and small rings, and derive some criteria for these differences, noting that recent work on slow quenches does not adequately explain the anomalous behaviour seen here.

preprint2011arXiv

Soliton form factors from lattice simulations

The form factor provides a convenient way to describe properties of topological solitons in the full quantum theory, when semiclassical concepts are not applicable. It is demonstrated that the form factor can be calculated numerically using lattice Monte Carlo simulations. The approach is very general and can be applied to essentially any type of soliton. The technique is illustrated by calculating the kink form factor near the critical point in 1+1-dimensional scalar field theory. As expected from universality arguments, the result agrees with the exactly calculable scaling form factor of the two-dimensional Ising model.

preprint2010arXiv

Studying a relativistic field theory at finite chemical potential with the density matrix renormalization group

The density matrix renormalization group is applied to a relativistic complex scalar field at finite chemical potential. The two-point function and various bulk quantities are studied. It is seen that bulk quantities do not change with the chemical potential until it is larger than the minimum excitation energy. The technical limitations of the density matrix renormalization group for treating bosons in relativistic field theories are discussed. Applications to other relativistic models and to nontopological solitons are also suggested.

preprint2009arXiv

Quantum kink and its excitations

We show how detailed properties of a kink in quantum field theory can be extracted from field correlation functions. This makes it possible to study quantum kinks in a fully non-perturbative way using Monte Carlo simulations. We demonstrate this by calculating the kink mass as well as the spectrum and approximate wave functions of its excitations. This way of measuring the kink mass has clear advantages over the existing approaches based on creation and annihilation operators or the kink free energy. Our methods are straightforward to generalise to more realistic theories and other defect types.