We study the dynamics of the Peccei-Quinn (PQ) phase transition for the QCD axion. In weakly coupled models the transition is typically second order except in the region of parameters where the PQ symmetry is broken through the Coleman-Weinberg mechanism. In strongly coupled realizations the transition is often first order. We show examples where the phase transition leads to strong supercooling lowering the nucleation temperature and enhancing the stochastic gravitational wave signals. The models predict a frequency peak in the range 100-1000 Hz with an amplitude that is already within the sensitivity of LIGO and can be thoroughly tested with future gravitational wave interferometers.

Dark Matter can form bound states upon the emission of quanta of energy equal to the binding energy. The rate of this process is large for strongly-interacting Dark Matter, and further enhanced by long distance effects. The resulting monochromatic and diffuse $\gamma$-rays can be tested in indirect detection experiments. If Dark Matter has electroweak charge, indirect signals include multiple observable photon lines for masses in the TeV range. Else if it couples only via a dark photon portal, diffuse spectra from dwarf galaxies and CMB reionization set powerful limits for masses below a TeV. This mechanism provides a powerful means of probing Asymmetric Dark Matter today.

We study the formation of Dark Matter nuclei in scenarios where DM particles are baryons of a new confining gauge force. The dark nucleosynthesis is analogous to the formation of light elements in the SM and requires as a first step the formation of dark deuterium. We compute this process from first principles, using the formalism of pion-less effective theory for nucleon-nucleon interactions. This controlled effective field theory expansion allows us to systematically compute the cross sections for generic SM representations under the assumption of shallow bound states. In the context of vector-like confinement models we find that, for nucleon masses in the TeV range, baryonic DM made of electro-weak constituents can form a significant fraction of dark deuterium and a much smaller fraction of dark tritium. Formation of dark nuclei can also lead to monochromatic photon lines in indirect detection. Models with singlets do not undergo nucleosynthesis unless a dark photon is added to the theory.

We study a scenario in which the baryon asymmetry of the universe arises from a cosmological phase transition where lepton-number is spontaneously broken. If the phase transition is first order, a lepton-number asymmetry can arise at the bubble wall, through dynamics similar to electroweak baryogenesis, but involving right-handed neutrinos. In addition to the usual neutrinoless double beta decay in nuclear experiments, the model may be probed through a variety of "baryogenesis by-products," which include a stochastic background of gravitational waves created by the colliding bubbles. Depending on the model, other aspects may include a network of topological defects that produce their own gravitational waves, additional contribution to dark radiation, and a light pseudo-Goldstone boson (majoron) as dark matter candidate.

We propose a framework in which the QCD axion has an exponentially large coupling to photons, relying on the "clockwork" mechanism. We discuss the impact of present and future axion experiments on the parameter space of the model. In addition to the axion, the model predicts a large number of pseudo-scalars which can be light and observable at the LHC. In the most favorable scenario, axion Dark Matter will give a signal in multiple axion detection experiments and the pseudo-scalars will be discovered at the LHC, allowing to determine most of the parameters of the model.