Flavor-violating Higgs decays and stellar cooling anomalies in axion models

We study a class of DFSZ-like models for the QCD axion that can address observed anomalies in stellar cooling. Stringent constraints from SN1987A and neutron stars are avoided by suppressed couplings to nucleons, while axion couplings to electrons and photons are sizable. All axion couplings depend on few parameters that also control the extended Higgs sector, in particular lepton flavor-violating couplings of the Standard Model-like Higgs boson h. This allows us to correlate axion and Higgs phenomenology, and we find that BR(h → τe) can be as large as the current experimental bound of 0.22%, while BR(h → μμ) can be larger than in the Standard Model by up to 70%. Large parts of the parameter space will be tested by the next generation of axion helioscopes such as the IAXO experiment.