Discrete Symmetry Research Articles

Stabilization of symmetry-protected long-range entanglement in stochastic quantum circuits

Long-range entangled states are vital for quantum information processing and quantum metrology. Preparing such states by combining measurements with unitary gates opened new possibilities for efficient protocols with finite-depth quantum circuits. The complexity of these algorithms is crucial for the resource requirements on a large-scale noisy quantum device, while their stability to perturbations decides the fate of their implementation. In this work, we consider stochastic quantum circuits in one and two dimensions comprising randomly applied unitary gates and local measurements. These operations preserve a class of discrete local symmetries, which are broken due to the stochasticity arising from timing and gate imperfections. In the absence of randomness, the protocol generates a symmetry-protected long-range entangled state in a finite-depth circuit. In the general case, by studying the time evolution under this hybrid circuit, we analyze the time to reach the target entangled state. We find two important time scales that we associate with the emergence of certain symmetry generators. The quantum trajectories embody the local symmetry with a time scaling logarithmically with system size, while global symmetries require exponentially long times. We devise error-mitigation protocols that significantly lower both time scales and investigate the stability of the algorithm to perturbations that naturally arise in experiments. We also generalize the protocol to realize toric code and Xu-Moore states in two dimensions, opening avenues for future studies of anyonic excitations. Our results unveil a fundamental relationship between symmetries and dynamics across a range of lattice geometries, which contributes to a broad understanding of the stability of preparation algorithms in terms of phase transitions. Our work paves the way for efficient error correction for quantum state preparation.

Observation of a phase transition from a continuous to a discrete time crystal

Discrete (DTCs) and continuous time crystals (CTCs) are novel dynamical many-body states, that are characterized by robust self-sustained oscillations, emerging via spontaneous breaking of discrete or continuous time translation symmetry. DTCs are periodically driven systems that oscillate with a subharmonic of the external drive, while CTCs are continuously driven and oscillate with a frequency intrinsic to the system. Here, we explore a phase transition from a continuous time crystal to a discrete time crystal. A CTC with a characteristic oscillation frequency ωCTC is prepared in a continuously pumped atom-cavity system. Modulating the pump intensity of the CTC with a frequency ωdr close to 2ωCTC leads to robust locking of ωCTC to ωdr/2 , and hence a DTC arises. This phase transition in a quantum many-body system is related to subharmonic injection locking of non-linear mechanical and electronic oscillators or lasers.

The lion, the witch, and the wormhole: ensemble averaging the symmetric product orbifold

We consider the ensemble average of two dimensional symmetric product orbifold CFTs SymN(\U0001d54bD) over the Narain moduli space. We argue for a bulk dual given by N copies of an abelian Chern-Simons theory coupled to topological gravity, endowed with a discrete gauge symmetry exchanging the N copies. As a check of this proposal, we calculate the ensemble average of various partition and correlation functions of the symmetric product orbifold theory and compare the resulting expressions to gauge theory quantities in the bulk. We comment on the ensemble average of the tensionless string partition function on AdS3 × S3 × \U0001d54b4 by considering the specific case of D = 4 with the addition of supersymmetry.

Symmetry-Enforced Many-Body Separability Transitions

We study quantum many-body mixed states with a symmetry from the perspective of , i.e., whether a mixed state can be expressed as an ensemble of short-range-entangled symmetric pure states. We provide evidence for “symmetry-enforced separability transitions” in a variety of states, where in one regime the mixed state is expressible as a convex sum of symmetric short-range-entangled pure states, while in the other regime, such a representation is not feasible. We first discuss the Gibbs state of Hamiltonians that exhibit spontaneous breaking of a discrete symmetry, and argue that the associated thermal phase transition can be thought of as a symmetry-enforced separability transition. Next we study cluster states in various dimensions subjected to local decoherence, and identify several distinct mixed-state phases and associated separability phase transitions, which also provides an alternative perspective on recently discussed “average symmetry-protected topological order.” We also study decohered p+ip superconductors, and find that if the decoherence breaks the fermion parity explicitly, then the resulting mixed state can be expressed as a convex sum of nonchiral states, while a fermion parity–preserving decoherence results in a phase transition at a nonzero threshold that corresponds to spontaneous breaking of fermion parity. Finally, we briefly discuss systems that satisfy the no low-energy trivial state property, such as the recently discovered good low-density parity-check codes, and argue that the Gibbs state of such systems exhibits a temperature-tuned separability transition. Published by the American Physical Society 2024

Finite modular majoron

We point out that the accidental U(1)B−L symmetry can arise from a finite modular symmetry ΓN in the type-I seesaw. The finite modular symmetry is spontaneously broken in such a way that the residual ℤNT\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\mathbb{Z}}_N^T $$\\end{document} discrete symmetry, associated with the T-transformation which shifts the modulus τ → τ + 1, remains unbroken. This discrete ℤNT\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\mathbb{Z}}_N^T $$\\end{document} symmetry mimics U(1)B−L, and hence the majoron appears as a pseudo Nambu-Goldstone boson of U(1)B−L. Without introducing additional interactions, the modulus τ can be stabilized by the Coleman-Weinberg (CW) potential given by the Majorana mass terms of the right-handed neutrinos. We study cosmological implications of the majoron, with particular interests in the dark matter and dark radiation, where the latter may alleviate the Hubble tension. We also find that the CW potential can have a wide range of nearly exponential shape which prevents τ from overshooting, and makes the amount of dark radiation not too large.

Hamiltonian Cycles on Ammann-Beenker Tilings

We provide a simple algorithm for constructing Hamiltonian graph cycles (visiting every vertex exactly once) on a set of arbitrarily large finite subgraphs of aperiodic two-dimensional Ammann-Beenker (AB) tilings. Using this result, and the discrete scale symmetry of AB tilings, we find exact solutions to a range of other problems which lie in the complexity class NP-complete for general graphs. These include the equal-weight traveling salesperson problem, providing, for example, the most efficient route a scanning tunneling microscope tip could take to image the atoms of physical quasicrystals with AB symmetries; the longest path problem, whose solution demonstrates that collections of flexible molecules of any length can adsorb onto AB quasicrystal surfaces at density one, with possible applications to catalysis; and the three-coloring problem, giving ground states for the q-state Potts model (q≥3) of magnetic interactions defined on the planar dual to AB, which may provide useful models for protein folding. Published by the American Physical Society 2024

A realistic neutrino mixing scheme arising from A4 symmetry

We propose a unique lepton mixing scheme and its association with an exact hierarchy-philic neutrino mass matrix texture in the light of Type-I+Type-II seesaw mechanism under the framework of discrete flavour symmetry. The proposed model successfully predicts the normal ordering of neutrino masses and the two Majorana phases. Additionally, the analysis extends to the effective Majorana neutrino mass, in the context of neutrinoless double beta .

Total ankle arthroplasty improves discrete and continuous stance phase gait symmetry

BackgroundTotal ankle arthroplasty (TAA) is used to treat symptomatic end-stage ankle arthritis (AA). However, little is known about TAA's effects on gait symmetry. Research questionDetermine if symmetry changes from before surgery through two years following TAA utilizing the normalized symmetry index (NSI) and statistical parametric mapping (SPM). Methods141 patients with end-stage unilateral AA were evaluated from a previously collected prospective database, where each participant was tested within two weeks of surgery (Pre-Op), one year and two years following TAA. Walking speed, hip extension angle and moment, hip flexion angle, ankle plantarflexion angle and moment, ankle dorsiflexion angle, weight acceptance (GRF1), and propulsive (GRF2) vertical ground reaction forces were calculated for each limb. Gait symmetry was assessed using the NSI. A linear mixed effects model with a single response for each gait symmetry variable was used to examine the fixed effect of follow-up time (Pre-Op, Post-1 yr, Post-2 yr) and the random effect of participant with gait speed as a covariate in the model. A one-dimensional repeated measures analysis of variance (ANOVA) statistical parameter mapping (SPM) was completed to examine differences in the time-series NSI to determine regions of significant differences between follow-up times. ResultsRelative to Pre-Op values, GRF1, and GRF2 showed increased symmetry for discrete metrics and the time-series NSI across sessions. Hip extension moment had the largest symmetry improvement. Ankle plantarflexion angle was different between Pre-Op and Post-2 yr (p=0.010); and plantarflexion moment was different between Pre- Op and each post-operative session (p<0.001). The time-series Ankle Angle NSI was greater during the early stance phase in the Pre-Op session compared to Post-2 yr. SignificanceSymmetry across most of the stance phase improved following TAA indicating that TAA successfully improves gait symmetry and future work should determine if these improvements restore symmetry to levels equivalent with health age-match controls.

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus (T2) and associated systems on the plane (R2), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

Neutrino Phenomenology and keV Dark Matter in the Two-Higgs Doublet Model with A4 Symmetry

Abstract We propose a minimal extended seesaw scheme based on the discrete symmetry $A_4\times Z_4\times Z_2\times Z_8$, which can successfully address neutrino phenomenology and keV sterile neutrino dark matter. The lepton mass hierarchy is naturally achieved. Active neutrino mixing angles can reached the best-fit points with the predictive Dirac CP violation phase. The active–sterile mixing matrix elements are small enough to access the observed cosmological dark matter abundance constraint with keV sterile neutrino dark matter. The effective neutrino masses are predicted to be in the ranges of the recent experimental limits.

Application of discrete symmetry to natural convection in vertical porous microchannels

Abstract This work focuses on the study of natural convection in a flat porous microchannel with asymmetric heating. The novelty of the work lies in the fact that for the first time the method of discrete symmetries was used to analyze the complete system of Navier–Stokes and energy equations in a two-dimensional approximation. Analytical solutions for velocity and temperature profiles have been derived based on symmetry analysis, taking into account boundary conditions such as slip and temperature jump at the channel walls. The effect of Grashof, Knudsen, Darcy, and Prandtl numbers on the flow characteristics in the microchannel and heat transfer coefficients was elucidated. At high Grashof numbers, an ascending flow near the hot wall and a descending flow near the cold wall arise. Increasing the Knudsen number leads to an increase in the velocity, temperature jump at the walls and a decrease in heat transfer coefficients. As the Darcy number increases, velocities amplify in both ascending and descending flows. The temperature jump at the hot wall grows up, while it remains unchanged at the cold wall. In the same time, the heat transfer coefficient at the hot wall decreases.

Disconnected gauge groups in the infrared

Gauging a discrete 0-form symmetry of a QFT is a procedure that changes the global form of the gauge group but not its perturbative dynamics. In this work, we study the Seiberg-Witten solution of theories resulting from the gauging of charge conjugation in 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 theories with SU(N) gauge group and fundamental hypermultiplets. The basic idea of our procedure is to identify the ℤ2 action at the level of the SW curve and perform the quotient, and it should also be applicable to non-lagrangian theories. We study dynamical aspects of these theories such as their moduli space singularities and the corresponding physics; in particular, we explore the complex structure singularity arising from the quotient procedure. We also discuss some implications of our work in regards to three problems: the geometric classification of 4d SCFTs, the study of non-invertible symmetries from the SW geometry, and the String Theory engineering of theories with disconnected gauge groups.

The Schwinger-Keldysh coset construction

The coset construction is a tool for systematically building low energy effective actions for Nambu-Goldstone modes. This technique is typically used to compute time-ordered correlators appropriate for S-matrix computations for systems in their ground state. In this paper, we extend this technique to the Schwinger-Keldysh formalism, which enables one to calculate a wider variety of correlators and applies also to systems in a mixed state. We focus our attention on internal symmetries and demonstrate that, after identifying the appropriate symmetry breaking pattern, Schwinger-Keldysh effective actions for Nambu-Goldstone modes can be constructed using the standard rules of the coset construction. Particular emphasis is placed on the thermal state and ensuring that correlators satisfy the KMS relation. We also discuss explicitly the power counting scheme underlying our effective actions. We comment on the similarities and differences between our approach and others that have previously appeared in the literature. In particular, our prescription does not require the introduction of additional “diffusive” symmetries and retains the full non-linear structure generated by the coset construction. We conclude with a series of explicit examples, including a computation of the finite-temperature two-point functions of conserved spin currents in non-relativistic paramagnets, antiferromagnets, and ferromagnets. Along the way, we also clarify the discrete symmetries that set antiferromagnets apart from ferromagnets, and point out that the dynamical KMS symmetry must be implemented in different ways in these two systems.

Finding origins of CMB anomalies in the inflationary quantum fluctuations

In this paper, we present compelling evidence for the parity asymmetry (a discrete symmetry separate from isotropy) in the Cosmic Microwave Background (CMB) map, measured through two-point temperature correlations. Any asymmetry associated with discrete symmetries, such as parity, challenges our understanding of quantum physics associated with primordial physics rather than LCDM (Λ Cold-Dark-Matter) itself. We commence by conducting a comprehensive analysis of the Planck CMB, focusing on the distribution of power in low-multipoles and temperature anticorrelations at parity conjugate points in position space. We find tension with the near scale-invariant power-law power spectrum of Standard Inflation (SI), with p-values of the order 𝒪(10-4 - 10-3). Alternatively, we explore the framework of direct-sum inflation (DSI), where a quantum fluctuation arises as a direct sum of two components evolving forward and backward in time at parity conjugate points in physical space. This mechanism results in a parity-asymmetric scale-dependent power spectrum, particularly prominent at low-multipoles, without any additional free model parameters. Our findings indicate that DSI is consistent with data on parity asymmetry, the absence of power at θ > 60°, and power suppression at low-even-multipoles which are major data anomalies in the SI model. Furthermore, we discover that the parameters characterizing the hemispherical power asymmetry anomaly become statistically insignificant when the large SI quadrupole amplitude is reduced to align with the data. DSI explains this low quadrupole with a p-value of 3.5%, 39 times higher than SI. Combining statistics from parameters measuring parity and low-ℓ angular power spectrum, we find that DSI is 50-650 times more probable than SI. In summary, our investigation suggests that while CMB temperature fluctuations exhibit homogeneity and isotropy, they also display parity-asymmetric behavior consistent with predictions of DSI. This observation provides a tantalizing evidence for the quantum mechanical nature of gravity.

Leptogenesis, dark matter and gravitational waves from discrete symmetry breaking

We analyse a model that connects the neutrino sector and the dark sector of the universe via a mediator Φ, stabilised by a discrete 𝒵 4 symmetry that breaks to a remnant 𝒵 2 upon Φ acquiring a non-zero vacuum expectation value (vϕ ). The model accounts for the observed baryon asymmetry of the universe via additional contributions to the canonical Type-I leptogenesis. The 𝒵 4 symmetry breaking scale (vϕ ) in the model not only establishes a connection between the neutrino sector and the dark sector, but could also lead to gravitational wave signals that are within the reach of current and future experimental sensitivities.

Discrete global symmetries: gauging and twisted compactification

Discrete global symmetries of 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 SCFTs are studied via two operations: gauging and twisted compactification. We consider gauging of discrete symmetries in several well-known 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 SCFTs, including SU(n) SQCD with 2n flavors, theories of class S\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{S} $$\\end{document} of type A2n−1, and Argyres-Douglas theories of type (AN, AN), as well as propose new 4d SCFTs as a result. The wreathing technique, which involves gauging a subgroup of the automorphism group of the quiver diagram of the corresponding 3d mirror theory, is exploited. This allows us to understand several properties of discretely gauged theories, including moduli spaces and how discrete gauging affects the mixed ’t Hooft anomaly between the 1-form symmetry and the 0-form flavor symmetry. Many examples are viewed through the lens of the Argyres-Seiberg duality and its generalization. We also examine discrete gauging of SU(2) SQCD with 4 flavors by various ℤ2 and ℤ2 × ℤ2 subgroups of the permutation group S4 using the superconformal index. Regarding compactification, we propose a magnetic quiver for 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 SU(n) SQCD with 2n flavors compactified on a circle with a ℤ2 twist. The twisted compactification by non-invertible symmetries of the 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM theory with gauge group SU(N) is revisited. The non-invertible symmetry naturally gives rise to a ℤk action on the scalar fields parametrizing the moduli space. Upon examining the ℤk invariant chiral ring of the Higgs branch, we find that, in addition to the largest branch of the moduli space that is expected to be captured by the ABJ(M) theory, there exist in general nilpotent operators that lead to a branch of the moduli space which is a radical ideal.

Reading between the rational sections: Global structures of 4d $\mathcal{N}=2$ KK theories

We study how the global structure of rank-one 4d \mathcal{N}=2𝒩=2 supersymmetric field theories is encoded into global aspects of the Seiberg-Witten elliptic fibration. Starting with the prototypical example of the \mathfrak{su}(2)𝔰𝔲(2) gauge theory, we distinguish between relative and absolute Seiberg-Witten curves. For instance, we discuss in detail the three distinct absolute curves for the SU(2)SU(2) and SO(3)_±SO(3)± 4d \mathcal{N}=2𝒩=2 gauge theories. We propose that the 11-form symmetry of an absolute theory is isomorphic to a torsion subgroup of the Mordell-Weil group of sections of the absolute curve, while the full defect group of the theory is encoded in the torsion sections of a so-called relative curve. We explicitly show that the relative and absolute curves are related by isogenies (that is, homomorphisms of elliptic curves) generated by torsion sections - hence, gauging a one-form symmetry corresponds to composing isogenies between Seiberg-Witten curves. We apply this approach to Kaluza-Klein (KK) 4d \mathcal{N}=2𝒩=2 theories that arise from toroidal compactifications of 5d and 6d SCFTs to four dimensions, uncovering an intricate pattern of 4d global structures obtained by gauging discrete 00-form and/or 11-form symmetries. Incidentally, we propose a 6d BPS quiver for the 6d M-string theory on \mathbb{R}^4× T^2ℝ4×T2.

Spontaneous breaking of the SO(2N) symmetry in the Gross-Neveu model

The canonical Gross-Neveu model for N two-component Dirac fermions in 2+1 dimensions suffers a continuous phase transition at a critical interaction gc1∼1/N at large N, at which its continuous symmetry SO(2N) is preserved and a discrete (Ising) symmetry becomes spontaneously broken. A recent mean-field calculation, however, points to an additional transition at a different critical gc2∼−Ngc1, at which SO(2N)→SO(N)×SO(N). To study the latter phase transition we rewrite the Gross-Neveu interaction g(ψ¯ψ)2 in terms of three different quartic terms for the single (L=1) 4N-component real (Majorana) fermion, and then extend the theory to L&gt;1. This allows us to track the evolution of the fixed points of the renormalization group transformation starting from L≫1, where one can discern three distinct critical points which correspond to continuous phase transitions into (1) SO(2N)-singlet mass-order-parameter, (2) SO(2N)-symmetric-tensor mass-order-parameters, and (3) SO(2N)-adjoint nematic-order-parameters, down to L=1 value that is relevant to the standard Gross-Neveu model. Below the critical value of Lc(N)≈0.35N for N≫1 only the Gross-Neveu critical point (1) still implies a diverging susceptibility for its corresponding (SO(2N)-singlet) order parameter, whereas the two new critical points that existed at large L ultimately become equivalent to the Gaussian fixed point at L=1. We interpret this metamorphosis of the SO(2N)-symmetric-tensor fixed point from critical to spurious as an indication that the transition at gc2 in the original Gross-Neveu model is turned first-order by fluctuations. Published by the American Physical Society 2024

Thermal relic right-handed neutrino dark matter

It is known that two heavy Majorana right-handed neutrinos are sufficient to generate the baryon asymmetry in the present universe. Thus, it is interesting to identify the third right-handed neutrino N with the dark matter. We impose a new discrete symmetry Z2 on this dark matter neutrino to stabilize it. However, the U(1)B−L gauge boson A′ couples to the right-handed neutrino N. If the B−L breaking scale VB−L is sufficiently low, the dark matter neutrino N can be in the thermal bath. We find that the thermal relic N can explain the dark matter abundance for the B−L breaking scale VB−L∼O(10) TeV. After considering all the constraints from the existing experiments, a narrow mass region of the thermal produced right-handed neutrino dark matter N is still surviving.

A model of pseudo-Nambu-Goldstone dark matter with two complex scalars

Pseudo-Nambu-Goldstone (pNG) dark matter (DM) is a promising DM candidate and able to explain the measured DM abundance by the thermal freeze-out mechanism evading the stringent bound from DM direct detection experiments. We propose a new model providing a pNG DM by introducing two Standard-Model-singlet complex scalars with the same charges of a dark U(1) gauge symmetry. They are also charged under a U(1) global symmetry corresponding to their relative phase rotations, which is explicitly broken by a soft-breaking term in the scalar potential. The both U(1) symmetries are spontaneously broken by their vacuum expectation values, giving rise to one real pNG boson. We also introduce a discrete ℤ2 symmetry exchanging the two scalars to stabilize the pNG boson as DM. It is shown that this model reproduces the DM abundance consistently with the current bound from the direct detection experiments. The model has a gauge kinetic mixing between the dark and U(1)Y gauge fields, which allows the dark gauge boson to decay even with a relatively light mass and prevents it from being an additional DM component. The Landau pole is avoided thanks to the small gauge coupling constant. In addition, a DM pair dominantly annihilates into a pair of the dark gauge bosons if the gauge boson mass is lighter than the DM mass, and thus its cross section has significantly different parameter dependence from other pNG DM models. We also calculate the DM-nucleon scattering cross section at the loop level. It turns out that it is necessary to probe region covered by the neutrino fog in order to test this model.