Semi-Dirac Quantum Critical Region

• The semi-Dirac quantum critical region is defined by a unique band-touching with linear dispersion in one momentum direction and quadratic in the other, reflecting its mixed kinematics.
• This regime underpins transitions such as the Dirac semimetal–band insulator merger, magnetic-field-induced quantum Hall effects, and excitonic pairing that leads to non-Fermi-liquid behavior.
• Key signatures include anisotropic thermodynamic scaling, direction-dependent optical conductivity, and quantized Hall conductivity despite broken rotational symmetry.

The semi-Dirac quantum critical region is the regime controlled by a semi-Dirac fixed point, where the low-energy spectrum is linear in one momentum direction and quadratic in the orthogonal direction. In two dimensions this structure arises most directly at the critical point where two Dirac cones merge and annihilate, separating a Dirac semimetal from a band insulator, but the same anisotropic kinematics also organizes magnetic-field-induced quantum Hall response, interaction-driven excitonic transitions, and Yukawa criticality toward density-wave or superconducting order (Hoyos et al., 2020, Roy et al., 2017, Wang et al., 2016, Uryszek et al., 2019). The defining feature throughout is anisotropic scaling, commonly summarized as $E\sim k_{\mathrm{lin}}\sim k_{\mathrm{quad}}^{2}$, together with a vanishing but enhanced low-energy density of states relative to an isotropic Dirac semimetal.

1. Definition, microscopic origin, and scope

A canonical semi-Dirac Hamiltonian at the Dirac-semimetal–band-insulator boundary is

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$

with spectrum

$E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$

At $\Delta=0$ the node is gapless and semi-Dirac; for $\Delta<0$ the system is a Dirac semimetal with two Dirac points at $k_x=0$, $k_y=\pm\sqrt{-\Delta/b}$; for $\Delta>0$ it is a trivial band insulator (Roy et al., 2017). An equivalent formulation used in quantum Hall effective-action work is

$$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$

which gives

$$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$

and reduces at $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$0 to a spectrum linear along one axis and quadratic along the other (Hoyos et al., 2020).

A notational difference across the literature is that the linear and quadratic axes are sometimes interchanged by convention. The underlying content is unchanged: semi-Dirac kinematics is the mixed linear–quadratic band touching itself. In the generic family $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$1, the physical two-dimensional semi-Dirac case corresponds to $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$2, for which the low-energy density of states obeys $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$3 (Roy et al., 2017).

The literature uses the semi-Dirac quantum critical region in several closely related senses. The common structure is summarized below.

Setting	Control parameter	Characteristic regime
DSM–BI merging point	$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$4	Gapless semi-Dirac node at $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$5
Quantum Hall semi-Dirac regime	$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$6	Universal scaling for $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$7
Excitonic semimetal–insulator QCP	$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$8	Finite-temperature NFL fan above $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$9
Yukawa ordering transitions	$E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$0 or $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$1	Anisotropic criticality with $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$2

Because there is no symmetry distinction between the Dirac semimetal and the symmetry-preserving band insulator in the merging problem, the direct DSM–BI transition is topological rather than symmetry-breaking; the semi-Dirac point is the associated quantum critical point (Roy et al., 2017).

2. Anisotropic scaling and the quantum critical fan

The hallmark of semi-Dirac criticality is anisotropic scaling. In the covariant quantum Hall formulation, with $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$3 transverse and $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$4 along the anisotropy vector, the semi-Dirac point $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$5 is invariant under

$E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$6

so that $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$7, $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$8, and $E_{\mathbf k} = \pm\sqrt{v^2 k_x^2 + b^2 k_y^4}.$9 (Hoyos et al., 2020). In the merging-transition literature the same physics is written as

$\Delta=0$0

which for $\Delta=0$1 again implies $\Delta=0$2 (Roy et al., 2017). The two formulations differ only by choice of reference scaling variable.

This anisotropy controls the thermodynamic and transport scaling of the semi-Dirac fan. For the physical $\Delta=0$3 case, the quantum critical fan appears for $\Delta=0$4 with $\Delta=0$5, and similarly for $\Delta=0$6 or $\Delta=0$7. In that regime the density of states obeys $\Delta=0$8, the compressibility scales as $\Delta=0$9, the specific heat at charge neutrality as $\Delta<0$0, the optical conductivity as $\Delta<0$1 and $\Delta<0$2, and the diamagnetic susceptibility as $\Delta<0$3 or $\Delta<0$4 in the weak-field finite-$\Delta<0$5 regime (Roy et al., 2017).

A second, interaction-driven usage of the term refers to the finite-temperature fan above the excitonic semimetal–insulator quantum critical point. There the control parameter is the Coulomb strength $\Delta<0$6, the zero-temperature critical point sits at $\Delta<0$7, and the quantum critical region is entered when $\Delta<0$8 exceeds the zero-temperature gap or correlation scale associated with $\Delta<0$9; the schematic criterion is $k_x=0$0 (Wang et al., 2016). This fan is labeled NFL in the phase diagram because the critical excitonic fluctuations destroy Fermi-liquid behavior.

3. Field-theoretic formulations

Semi-Dirac criticality admits both geometric effective-action and fermion-boson Yukawa descriptions. In the geometric construction for Hall transport, anisotropy is encoded by a unit spacelike vector $k_x=0$1 and the projector $k_x=0$2. In flat space the fermionic action is

$k_x=0$3

with

$k_x=0$4

This explicitly breaks $k_x=0$5 down to $k_x=0$6, leaving relativistic dynamics only in the transverse $k_x=0$7-dimensional plane and quadratic dynamics along the longitudinal direction (Hoyos et al., 2020).

In curved space the same theory becomes

$k_x=0$8

with background fields $k_x=0$9, vielbeins, $k_y=\pm\sqrt{-\Delta/b}$0, and $k_y=\pm\sqrt{-\Delta/b}$1. The relevant covariant building blocks are the electromagnetic field strength $k_y=\pm\sqrt{-\Delta/b}$2, the reduced curvature $k_y=\pm\sqrt{-\Delta/b}$3, and the anisotropy “field strength” $k_y=\pm\sqrt{-\Delta/b}$4, all arranged in a derivative expansion consistent with diffeomorphism invariance, local $k_y=\pm\sqrt{-\Delta/b}$5, and imposed $k_y=\pm\sqrt{-\Delta/b}$6 symmetry (Hoyos et al., 2020).

For interaction-driven transitions, the natural description is a Yukawa theory of semi-Dirac fermions and an order-parameter field. For CDW, SDW, and $k_y=\pm\sqrt{-\Delta/b}$7-wave SC, the generic action is

$k_y=\pm\sqrt{-\Delta/b}$8

with semi-Dirac fermion kinetic term

$k_y=\pm\sqrt{-\Delta/b}$9

and Yukawa coupling $\Delta>0$0 (Uryszek et al., 2019). A decisive feature in $\Delta>0$1 dimensions is that the bosonic infrared propagator is not relativistic. After integrating over fermions one obtains

$\Delta>0$2

so that at criticality

$\Delta>0$3

This non-analytic kernel regularizes the loop integrals and transmits the fermionic anisotropy directly to the order-parameter sector (Uryszek et al., 2019).

The excitonic semimetal–insulator transition has an analogous structure. There the critical real scalar $\Delta>0$4 couples through $\Delta>0$5, while the bosonic propagator is dominated by the fermionic polarization

$\Delta>0$6

Because the boson dynamics is governed by this non-analytic fermion-induced term, the Hertz-Millis procedure of integrating out fermions is invalid; fermions and bosons must be kept on equal footing (Wang et al., 2016).

4. Magnetic field, effective action, and Hall response

In the quantum Hall setting the semi-Dirac quantum critical region is the regime $\Delta>0$7, where the Landau scale and the anisotropic fixed-point symmetry dominate the effective action (Hoyos et al., 2020). The associated Landau-level scaling is

$\Delta>0$8

and the low-energy derivative expansion is organized so that transverse derivatives are $\Delta>0$9, longitudinal derivatives are $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$0, the magnetic field is $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$1, and electric fields are small (Hoyos et al., 2020).

At leading order, the effective Lagrangian contains an energy-density term and an ordinary Chern-Simons term,

$$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$2

Consequently, the Hall conductivity remains

$$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$3

exactly as in isotropic relativistic Dirac systems. The semi-Dirac anisotropy therefore does not alter the quantization of $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$4 (Hoyos et al., 2020).

The anisotropic structure appears instead in subleading parity-odd, $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$5-allowed terms and in the viscosity sector. At $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$6 the effective action contains six invariants built from $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$7, $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$8, and longitudinal derivatives of the anisotropic magnetic scalars; among their coefficients, $$\mathcal H = \mathbf d(\mathbf p)\cdot\boldsymbol\sigma,\qquad \mathbf d(\mathbf p)=\left(p_x,\frac{p_y^2}{2m}-\Delta,0\right),$$9 is especially important. At $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$0 the theory contains an Euler-current coupling with coefficient $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$1 and torsional-like terms involving $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$2, with coefficient $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$3 entering the final viscosity formulas (Hoyos et al., 2020).

The odd viscosity tensor in a $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$4-invariant anisotropic state decomposes into isotropic, nematic, and vorticity-related pieces. The corresponding coefficients are

$$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$5

$$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$6

$$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$7

At the semi-Dirac point the coefficients scale as

$$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$8

so both $$E_\pm(\mathbf p)=\pm \sqrt{p_x^2+\left(\frac{p_y^2}{2m}-\Delta\right)^2},$$9 and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$00 scale linearly with $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$01, up to the topological contribution $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$02 in $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$03 (Hoyos et al., 2020). A central consequence is that the Hall conductivity remains topologically quantized while the Hall viscosities split into anisotropy-sensitive components.

5. Interaction-driven criticality and non-Fermi-liquid behavior

Long-range Coulomb interaction is substantially more effective in a semi-Dirac semimetal than in graphene-like Dirac systems. In the model

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$04

the spectrum

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$05

and the unusual dynamical screening make the system more susceptible to excitonic pairing (Wang et al., 2016). Solving the Dyson-Schwinger gap equation with three standard kernel approximations yields, for $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$06 and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$07, critical couplings

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$08

in the instantaneous approximation,

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$09

in the Khveshchenko approximation, and

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$10

in the GGG approximation. The corresponding values for a two-dimensional Dirac semimetal are much larger: approximately $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$11, $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$12, and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$13, respectively (Wang et al., 2016). Additional short-range four-fermion coupling in the same $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$14 channel enhances excitonic pairing and lowers $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$15, with $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$16 when $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$17 (Wang et al., 2016).

At the excitonic quantum critical point, the Yukawa coupling to the critical Ising-like field $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$18 produces a genuine non-Fermi-liquid regime. To leading order in $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$19,

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$20

so $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$21 as $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$22, and

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$23

The decay rate is therefore parametrically larger than the Fermi-liquid form $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$24, and the quasiparticle peak becomes ill-defined (Wang et al., 2016). This is the microscopic basis of the NFL fan above $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$25.

For ordering transitions out of a semi-Dirac semimetal driven by short-range interactions, the large-$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$26 Yukawa theory in $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$27 dimensions yields a controlled anisotropic critical point. With $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$28, the one-loop results are

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$29

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$30

together with distinct spatial correlation-length exponents

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$31

In the $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$32 limit one obtains

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$33

Thus order-parameter correlations inherit the electronic anisotropy of the semi-Dirac fermions, and the two spatial correlation lengths diverge with different powers already at mean-field level (Uryszek et al., 2019).

6. Multicriticality, emergent phases, and materials

Weak short-range interactions do not immediately destabilize the semi-Dirac point. In the generalized anisotropic semimetal with $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$34, every quartic coupling has engineering dimension

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$35

so for the physical $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$36 case the semi-Dirac anisotropic semimetal is stable against weak local interactions (Roy et al., 2017). At stronger coupling, however, the direct DSM–BI transition can be preempted either by a first-order transition or by an intervening broken-symmetry phase. The dominant ordered states identified in the spinful problem are charge density wave, antiferromagnet, and singlet $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$37-wave superconductor, and the corresponding interacting quantum critical points have

$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$38

hence $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$39 for $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$40 (Roy et al., 2017).

The fixed-point structure is organized through a controlled deformation away from the one-dimensional limit. In the spinful case the RG yields six fixed points: the trivial noninteracting ASM fixed point, three interacting QCPs with one unstable direction in interaction space, and two bicritical points. FP$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$41 describes a pseudospin-SU(2)-symmetric transition toward coexisting CDW and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$42-wave SC; FP$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$43 controls an ASM–CDW transition; FP$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$44 controls an ASM–AFM transition (Roy et al., 2017). This suggests that the semi-Dirac quantum critical region is not tied to a single universality class, but to a family of anisotropic critical theories built on the same semi-Dirac kinematics.

The mean-field and large-$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$45 analyses of Yukawa criticality also indicate unusually soft spatial gradients. In the order-parameter free energy one finds non-analytic gradient terms, and this led to the conjecture that proximity to the critical point may stabilize novel modulated order phases (Uryszek et al., 2019). That conjecture is specific, not generic: it follows from the anisotropic stiffness structure rather than from symmetry alone.

Material platforms proposed across the literature include uniaxially strained honeycomb lattices, optical honeycomb lattices for ultracold fermions, the organic material $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$46-($$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$47, black phosphorus or few-layer phosphorene, oxide interfaces such as TiO$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$48/VO$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$49, and artificial or molecular graphene (Roy et al., 2017). In the excitonic context, TiO$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$50/VO$$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$51 is singled out because first-principles work predicts intrinsic two-dimensional semi-Dirac bands, and the estimated parameters $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$52 m/s, $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$53, and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$54 imply $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$55, close to the calculated critical couplings (Wang et al., 2016).

Two common confusions are resolved by the existing literature. First, semi-Dirac criticality is not merely isotropic Dirac criticality with rescaled velocities: it generically involves non-analytic bosonic dynamics and distinct spatial exponents $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$56 and $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$57 (Uryszek et al., 2019). Second, anisotropy does not automatically spoil topological Hall response: in the quantum Hall regime $$H(\mathbf k,\Delta)= \sigma_0 \left[ v k_x \,\tau_1 + (b k_y^2 + \Delta)\,\tau_2\right],$$58 remains quantized, while the viscosity sector develops genuinely anisotropic, nematic, and vorticity-related structures (Hoyos et al., 2020). Taken together, these results define the semi-Dirac quantum critical region as a broadly applicable but sharply structured framework for anisotropic quantum criticality in two-dimensional systems with merged Dirac cones.