pith. sign in

arxiv: 2604.18908 · v1 · submitted 2026-04-20 · ⚛️ physics.plasm-ph

Intense tunable terahertz radiation from phase-matched difference frequency generation in strongly magnetized plasmas

Pith reviewed 2026-05-10 02:44 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords terahertz radiationmagnetized plasmadifference frequency generationphase matchinglaser-plasma interactionparticle-in-cell simulation
0
0 comments X

The pith

Two-color laser pulses in strongly magnetized plasma generate tunable terahertz radiation with fields exceeding 500 GV/m by using extraordinary modes to minimize phase mismatch.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that high-efficiency, tunable terahertz pulses with extreme field strengths can be produced by sending two-color laser pulses through a strongly magnetized plasma. This works through difference frequency generation where phase mismatch is reduced by selecting two extraordinary-mode branches of the plasma waves. A sympathetic reader would care because nonlinear crystals suffer from low conversion efficiency and optical damage, while prior plasma approaches lack comparable intensity and tunability. The authors derive the required phase-matching conditions analytically, characterize the nonlinear coupling, and confirm the predictions with particle-in-cell simulations.

Core claim

Propagating two-color laser pulses through a strongly magnetized plasma enables phase-matched difference frequency generation of terahertz radiation with tunable frequency and field strengths exceeding 500 GV/m. This is substantially enhanced by utilizing two extraordinary-mode branches to minimize the phase mismatch. The phase-matching conditions are derived and the nonlinear coupling is characterized analytically, with the predictions validated through particle-in-cell simulations.

What carries the argument

Difference frequency generation between two extraordinary-mode branches of plasma waves, which minimizes phase mismatch for efficient terahertz production.

If this is right

  • Field strengths exceed those of existing crystal-based and plasma-based terahertz sources.
  • Frequency tunability arises from adjustments to laser wavelengths or plasma density.
  • Analytical expressions for phase-matching and coupling efficiency guide experimental design.
  • Performance is limited by plasma maintenance rather than material damage thresholds.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Laboratory tests with high-power lasers and strong magnets would directly test whether the predicted field strengths are reached.
  • The approach might combine with plasma wakefield acceleration for compact integrated THz sources.
  • Similar phase-matching strategies could be explored in other wave regimes or plasma geometries.

Load-bearing premise

The plasma remains in a strongly magnetized, low-collision state under the intense two-color laser drive without instabilities or nonlinear effects disrupting the phase-matching conditions.

What would settle it

A particle-in-cell simulation or experiment in which the terahertz field strength stays well below 500 GV/m despite satisfying the derived phase-matching conditions for the two extraordinary modes.

Figures

Figures reproduced from arXiv: 2604.18908 by Matthew R. Edwards, Sida Cao.

Figure 1
Figure 1. Figure 1: FIG. 1. A schematic of THz generation using an ultrashort [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. PIC simulation results of generated terahertz pulses [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Conversion efficiency scan of THz generation [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Three-dimensional PIC simulation of THz generation [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) The conversion efficiency to THz radiation for [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. (a) Phase-matching conditions as a function of [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
read the original abstract

High-energy terahertz pulses are challenging to produce due to the low conversion efficiency and limited optical damage threshold of nonlinear crystals. Here, we describe the high-efficiency generation of terahertz radiation pulses with tunable frequency and field strengths exceeding 500 GV/m by propagating two-color laser pulses through a strongly magnetized plasma. The field strength is substantially enhanced by utilizing two extraordinary-mode branches to minimize the phase mismatch. We derive the phase-matching conditions and characterize the nonlinear coupling analytically, and validate these predictions with particle-in-cell simulations. These results establish a new pathway toward next-generation intense terahertz sources with performance well beyond the limits of existing plasma mechanisms and conventional crystal-based approaches.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims that propagating two-color laser pulses through a strongly magnetized plasma enables high-efficiency, tunable terahertz generation via difference-frequency mixing, with field strengths exceeding 500 GV/m achieved by phase-matching two extraordinary-mode branches. It derives the phase-matching conditions and nonlinear coupling coefficients analytically from cold-plasma dispersion relations and validates the predictions with particle-in-cell simulations.

Significance. If the central result holds, the work would establish a plasma-based route to intense tunable THz sources with performance well beyond crystal limits and prior plasma mechanisms, supported by the analytical phase-matching derivation and PIC validation. These elements provide a concrete, falsifiable pathway that could be tested experimentally.

major comments (2)
  1. [Analytical derivation and PIC validation sections] The phase-matching conditions and coupling analysis rest on the assumption of a uniform, unperturbed background B-field and density throughout the interaction length. At the laser intensities needed to reach >500 GV/m THz fields, ponderomotive expulsion, relativistic mass increase, and possible filamentation or two-stream instabilities can locally modify the extraordinary-mode dispersion; the manuscript must demonstrate (via additional PIC diagnostics or analytic estimates) that these effects remain negligible over the quoted propagation distance.
  2. [PIC simulations] The quoted field strength of 500 GV/m is load-bearing for the central claim, yet the support is moderate because the PIC runs appear to use idealized fixed external B and uniform initial density. A quantitative comparison showing the actual THz field amplitude extracted from the simulation versus the analytic prediction, including any degradation due to self-consistent plasma evolution, is required.
minor comments (2)
  1. [Abstract] The abstract and introduction should explicitly state the representative values of B0 and ne used to reach the 500 GV/m figure so that readers can assess the magnetization parameter ωce/ωpe.
  2. [Throughout] Notation for the two extraordinary branches (e.g., labeling of the dispersion branches) should be introduced once and used consistently in equations and figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and positive evaluation of the significance of our work. We address each major comment below, providing clarifications and committing to revisions that strengthen the manuscript's analytical and numerical support.

read point-by-point responses
  1. Referee: [Analytical derivation and PIC validation sections] The phase-matching conditions and coupling analysis rest on the assumption of a uniform, unperturbed background B-field and density throughout the interaction length. At the laser intensities needed to reach >500 GV/m THz fields, ponderomotive expulsion, relativistic mass increase, and possible filamentation or two-stream instabilities can locally modify the extraordinary-mode dispersion; the manuscript must demonstrate (via additional PIC diagnostics or analytic estimates) that these effects remain negligible over the quoted propagation distance.

    Authors: We agree that the high intensities required for >500 GV/m THz fields necessitate explicit verification that nonlinear effects do not invalidate the uniform-background assumption. Our original PIC runs already evolve the plasma self-consistently, but we did not present supporting diagnostics. In the revision we will add (i) analytic estimates of the ponderomotive scale length, relativistic mass correction, and instability growth rates (filamentation and two-stream) for the quoted parameters, showing that perturbations remain <5% over the several-millimeter propagation distance, and (ii) new PIC diagnostics (line-outs of density and B-field at multiple times) confirming uniformity. These additions will be placed in a new subsection of the PIC validation section. revision: yes

  2. Referee: [PIC simulations] The quoted field strength of 500 GV/m is load-bearing for the central claim, yet the support is moderate because the PIC runs appear to use idealized fixed external B and uniform initial density. A quantitative comparison showing the actual THz field amplitude extracted from the simulation versus the analytic prediction, including any degradation due to self-consistent plasma evolution, is required.

    Authors: We acknowledge that a direct, quantitative comparison between the simulated THz amplitude and the analytic prediction was not included. Although the simulations use self-consistent evolution (not fixed external B or frozen density), we will revise the manuscript to extract the peak THz electric-field amplitude directly from the simulation data at the end of the interaction region. This value will be compared side-by-side with the analytic result obtained from the derived coupling coefficient and phase-matching condition. Any amplitude reduction due to self-consistent effects will be quantified and discussed, including results from an ensemble of runs to indicate statistical variation. revision: yes

Circularity Check

0 steps flagged

Derivation relies on standard cold-plasma dispersion relations with no reduction to fitted inputs or self-citations

full rationale

The phase-matching conditions for the two extraordinary-mode branches are obtained directly from the cold-plasma dispersion relation in a uniform magnetized plasma, which is an external, well-established result independent of any quantities fitted or defined inside this paper. The nonlinear coupling coefficient is likewise derived analytically from first-principles fluid or kinetic equations without circular substitution of the target THz field strength. PIC simulations are used only for validation, not for parameter fitting that would force the analytic predictions. No self-citation chain is load-bearing for the central claim, and no ansatz or uniqueness theorem is smuggled in via prior work by the same authors. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The claim rests on standard plasma dispersion relations and nonlinear wave coupling under the assumption of a uniform strong magnetic field; no new entities are introduced.

free parameters (2)
  • magnetic field strength
    Selected to satisfy phase-matching condition between the two extraordinary modes; value not derived from first principles but chosen for the target THz frequency.
  • plasma density
    Tuned to place the operating point on the desired dispersion branches.
axioms (2)
  • domain assumption Plasma is cold, collisionless, and uniformly magnetized with no thermal effects or density gradients affecting wave propagation.
    Invoked to derive the extraordinary-mode dispersion and phase-matching conditions.
  • domain assumption Laser intensities remain below the threshold for relativistic or other higher-order nonlinearities that would invalidate the analytic coupling description.
    Required for the nonlinear coupling characterization to hold as stated.

pith-pipeline@v0.9.0 · 5405 in / 1521 out tokens · 38701 ms · 2026-05-10T02:44:36.510083+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages

  1. [1]

    1 Bc η 0 50 100 150 200 10−4 10−3 10−2 10−1 L/λ 0 η ∆ k ̸= 0 ∆ k = 0 −50 0 50−200 0 200 t [fs] E3 [GV/ m] −50 0 50−200 0 200 t [fs] (a) (b) FIG. 3. (a) Conversion efficiency scan of THz generation withα= 0.1at varied plasma densityN 0 and magnetic field strengthB c. The plasma length isL= 100λ 0. The solid line marks the phase-matching condition and the i...

  2. [2]

    G. A. Mourou, T. Tajima, and S. V. Bulanov, Optics in the relativistic regime, Rev. Mod. Phys.78, 309 (2006)

  3. [3]

    Liao and Y

    G. Liao and Y. Li, Perspectives on ultraintense laser- driven terahertz radiation from plasmas, Phys. Plasmas 30(2023)

  4. [4]

    Pálfalvi, J

    L. Pálfalvi, J. A. Fülöp, G. Tóth, and J. Hebling, Evanescent-wave proton postaccelerator driven by in- tense thz pulse, Phys. Rev. Accel. Beams17, 031301 (2014)

  5. [5]

    Sharma, Z

    A. Sharma, Z. Tibai, and J. Hebling, Intense tera-hertz laser driven proton acceleration in plasmas, Phys. Plas- mas23(2016)

  6. [6]

    Sharma, Z

    A. Sharma, Z. Tibai, J. Hebling, and J. A. Fülöp, Terahertz-driven wakefield electron acceleration, J. Phys. B: At. Mol. Opt. Phys.51, 204001 (2018)

  7. [7]

    Hernández-García, J

    C. Hernández-García, J. A. Pérez-Hernández, T. Pop- mintchev, M. M. Murnane, H. C. Kapteyn, A. Jaron- Becker, A. Becker, and L. Plaja, Zeptosecond high har- monic keV x-ray waveforms driven by midinfrared laser pulses, Phys. Rev. Lett.111, 033002 (2013)

  8. [8]

    Liang, P

    H. Liang, P. Krogen, Z. Wang, H. Park, T. Kroh, K. Za- wilski, P. Schunemann, J. Moses, L. F. DiMauro, F. X. Kärtner,et al., High-energy mid-infrared sub-cycle pulse synthesis from a parametric amplifier, Nat. Commun.8, 141 (2017)

  9. [9]

    I.Thiele, E.Siminos,andT.Fülöp,Electronbeamdriven generation of frequency-tunable isolated relativistic sub- cycle pulses, Phys. Rev. Lett.122, 104803 (2019)

  10. [10]

    Qi, Y.-H

    T. Qi, Y.-H. Shin, K.-L. Yeh, K. A. Nelson, and A. M. Rappe, Collective coherent control: Synchronization of polarization in ferroelectricPbTiO 3 by shaped THz fields, Phys. Rev. Lett.102, 247603 (2009)

  11. [11]

    Serrat and J

    C. Serrat and J. Biegert, All-regions tunable high har- monic enhancement by a periodic static electric field, Phys. Rev. Lett.104, 073901 (2010)

  12. [12]

    Balogh, K

    E. Balogh, K. Kovacs, P. Dombi, J. A. Fulop, G. Farkas, J. Hebling, V. Tosa, and K. Varju, Single attosecond pulse from terahertz-assisted high-order harmonic gen- eration, Phys. Rev. A84, 023806 (2011)

  13. [13]

    Kovács, E

    K. Kovács, E. Balogh, J. Hebling, V. Toşa, and K. Varjú, Quasi-phase-matching high-harmonic radiation using chirped THz pulses, Phys. Rev. Lett.108, 193903 (2012)

  14. [14]

    Rumiantsev, E

    B. Rumiantsev, E. Migal, A. Pushkin, and F. Potemkin, Observation of terahertz-field-induced coherent control of high-order harmonic generation in a noble gas, Phys. Rev. A111, 023117 (2025)

  15. [15]

    Cook and R

    D. Cook and R. Hochstrasser, Intense terahertz pulses by four-wave rectification in air, Opt. Lett.25, 1210 (2000)

  16. [16]

    Bartel, P

    T. Bartel, P. Gaal, K. Reimann, M. Woerner, and T. El- saesser, Generation of single-cycle THz transients with high electric-field amplitudes, Opt. Lett.30, 2805 (2005)

  17. [17]

    A. Sell, A. Leitenstorfer, and R. Huber, Phase-locked 6 generation and field-resolved detection of widely tunable terahertz pulses with amplitudes exceeding 100 MV/cm, Opt. Lett.33, 2767 (2008)

  18. [18]

    K.-Y. Kim, A. J. Taylor, J. Glownia, and G. Rodriguez, Coherent control of terahertz supercontinuum generation in ultrafast laser–gas interactions, Nat. Photon.2, 605 (2008)

  19. [19]

    Junginger, A

    F. Junginger, A. Sell, O. Schubert, B. Mayer, D. Brida, M. Marangoni, G. Cerullo, A. Leitenstorfer, and R. Hu- ber, Single-cycle multiterahertz transients with peak fields above 10 MV/cm, Opt. Lett.35, 2645 (2010)

  20. [20]

    Vicario, B

    C. Vicario, B. Monoszlai, and C. P. Hauri, Gv/m single- cycle terahertz fields from a laser-driven large-size par- titioned organic crystal, Phys. Rev. Lett.112, 213901 (2014)

  21. [21]

    Knorr, J

    M. Knorr, J. Raab, M. Tauer, P. Merkl, D. Peller, E. Wittmann, E. Riedle, C. Lange, and R. Huber, Phase-locked multi-terahertz electric fields exceeding 13 MV/cm at a 190 kHz repetition rate, Opt. Lett.42, 4367 (2017)

  22. [22]

    Liang, S

    Y. Liang, S. Antipov, Q. Tian, L. Yan, Y. Du, C. Cheng, R. Li, W. Huang, and C. Tang, Generation of a chirp- controlled terahertz pulse by a tapered corrugated wake- field structure, Phys. Rev. Appl.19, 054085 (2023)

  23. [23]

    H. Kim, C. Kang, D. Jang, Y. Roh, S. H. Lee, J. W. Lee, J. H. Sung, S. K. Lee, and K.-Y. Kim, Ionizing terahertz waves with 260 MV/cm from scalable optical rectifica- tion, Light Sci. Appl.13, 118 (2024)

  24. [24]

    Hamster, A

    H. Hamster, A. Sullivan, S. Gordon, W. White, and R. Falcone, Subpicosecond, electromagnetic pulses from intense laser-plasma interaction, Phys. Rev. Lett.71, 2725 (1993)

  25. [25]

    Liao, Y.-T

    G.-Q. Liao, Y.-T. Li, Y.-H. Zhang, H. Liu, X.-L. Ge, S. Yang, W.-Q. Wei, X.-H. Yuan, Y.-Q. Deng, B.-J. Zhu,et al., Demonstration of coherent terahertz tran- sition radiation from relativistic laser-solid interactions, Phys. Rev. Lett.116, 205003 (2016)

  26. [26]

    T. Pak, M. Rezaei-Pandari, S. B. Kim, G. Lee, D. H. Wi, C. I. Hojbota, M. Mirzaie, H. Kim, J. H. Sung, S. K. Lee,et al., Multi-millijoule terahertz emission from laser- wakefield-accelerated electrons, Light Sci. Appl.12, 37 (2023)

  27. [27]

    T. T. Simpson, J. J. Pigeon, K. G. Miller, D. Ramsey, D. H. Froula, and J. P. Palastro, Dephasingless two-color terahertz generation, Sci. Rep.14, 26587 (2024)

  28. [28]

    Maity and G

    S. Maity and G. Arora, Enhanced terahertz emission from the wakefield ofCO 2-laser-created plasma, Phys. Rev. E111, 045205 (2025)

  29. [29]

    Sheng, K

    Z.-M. Sheng, K. Mima, J. Zhang, and H. Sanuki, Emission of electromagnetic pulses from laser wakefields through linear mode conversion, Phys. Rev. Lett.94, 095003 (2005)

  30. [30]

    G. Liao, Y. Li, C. Li, L. Su, Y. Zheng, M. Liu, W. Wang, Z. Hu, W. Yan, J. Dunn,et al., Bursts of terahertz ra- diation from large-scale plasmas irradiated by relativis- tic picosecond laser pulses, Phys. Rev. Lett.114, 255001 (2015)

  31. [31]

    Pukhov, A

    A. Pukhov, A. Golovanov, and I. Kostyukov, Efficient narrow-band terahertz radiation from electrostatic wake- fields in nonuniform plasmas, Phys. Rev. Lett.127, 175001 (2021)

  32. [32]

    L. Wang, Z. Zhang, S. Chen, Y. Chen, X. Hu, M. Zhu, W. Yan, H. Xu, L. Sun, M. Chen,et al., Millijoule tera- hertz radiation from laser wakefields in nonuniform plas- mas, Phys. Rev. Lett.132, 165002 (2024)

  33. [33]

    K. B. Kwon, T. Kang, H. S. Song, Y.-K. Kim, B. Ersfeld, D. A. Jaroszynski, and M. S. Hur, High-energy, short- duration bursts of coherent terahertz radiation from an embedded plasma dipole, Sci. Rep.8, 145 (2018)

  34. [34]

    J. Lee, H. S. Song, D. Park, M. Kumar, B. Ersfeld, S. R. Yoffe, D. A. Jaroszynski, and M. S. Hur, Intense nar- rowband terahertz pulses produced by obliquely colliding laser pulses in helium gas, Phys. Plasmas30(2023)

  35. [35]

    Kumar, B

    M. Kumar, B. Ersfeld, J. Lee, D. Park, S. Kim, I. Nam, M. Kim, S. Jeon, D. A. Jaroszynski, H. Suk,et al., Narrowband terahertz emission from a plasma oscillator imbedded in a plasma density gradient, Phys. Rev. Lett. 134, 015001 (2025)

  36. [36]

    X. Xie, J. Dai, and X.-C. Zhang, Coherent control of thz wave generation in ambient air, Phys. Rev. Lett.96, 075005 (2006)

  37. [37]

    Andreeva, O

    V. Andreeva, O. Kosareva, N. Panov, D. Ship- ilo, P. Solyankin, M. Esaulkov, P. González de Alaiza Martínez, A. Shkurinov, V. Makarov, L. Bergé, et al., Ultrabroad terahertz spectrum generation from an air-based filament plasma, Phys. Rev. Lett.116, 063902 (2016)

  38. [38]

    A. D. Koulouklidis, C. Gollner, V. Shumakova, V. Y. Fe- dorov, A. Pugžlys, A. Baltuška, and S. Tzortzakis, Ob- servation of extremely efficient terahertz generation from mid-infraredtwo-colorlaserfilaments,Nat.Commun.11, 292 (2020)

  39. [39]

    Leemans, C

    W. Leemans, C. Geddes, J. Faure, C. Tóth, J. Van Tilborg, C. Schroeder, E. Esarey, G. Fubiani, D. Auerbach, B. Marcelis,et al., Observation of tera- hertz emission from a laser-plasma accelerated electron bunch crossing a plasma-vacuum boundary, Phys. Rev. Lett.91, 074802 (2003)

  40. [40]

    A. J. Pearson, J. Palastro, and T. M. Antonsen, Simula- tion of terahertz generation in corrugated plasma waveg- uides, Phys. Rev. E83, 056403 (2011)

  41. [41]

    Clerici, M

    M. Clerici, M. Peccianti, B. E. Schmidt, L. Caspani, M. Shalaby, M. Giguère, A. Lotti, A. Couairon, F. Lé- garé, T. Ozaki,et al., Wavelength scaling of terahertz generation by gas ionization, Phys. Rev. Lett.110, 253901 (2013)

  42. [42]

    C. Miao, J. P. Palastro, and T. M. Antonsen, Laser pulse driven terahertz generation via resonant transition ra- diation in inhomogeneous plasmas, Phys. Plasmas23 (2016)

  43. [43]

    C. Miao, J. P. Palastro, and T. M. Antonsen, High-power tunablelaserdrivenTHzgenerationincorrugatedplasma waveguides, Phys. Plasmas24(2017)

  44. [44]

    Déchard, A

    J. Déchard, A. Debayle, X. Davoine, L. Gremillet, and L. Bergé, Terahertz pulse generation in underdense rela- tivistic plasmas: From photoionization-induced radiation to coherent transition radiation, Phys. Rev. Lett.120, 144801 (2018)

  45. [45]

    Déchard, X

    J. Déchard, X. Davoine, and L. Bergé, THz generation from relativistic plasmas driven by near-to far-infrared laser pulses, Phys. Rev. Lett.123, 264801 (2019)

  46. [46]

    Yoshii, C

    J. Yoshii, C. Lai, T. Katsouleas, C. Joshi, and W. Mori, Radiation from cerenkov wakes in a magnetized plasma, Phys. Rev. Lett.79, 4194 (1997)

  47. [47]

    Spence, T

    N. Spence, T. Katsouleas, P. Muggli, W. Mori, and R. Hemker, Simulations of cerenkov wake radiation sources, Phys. Plasmas8, 4995 (2001)

  48. [48]

    Yugami, T

    N. Yugami, T. Higashiguchi, H. Gao, S. Sakai, K. Taka- 7 hashi, H. Ito, Y. Nishida, and T. Katsouleas, Experimen- tal observation of radiation from cherenkov wakes in a magnetized plasma, Phys. Rev. Lett.89, 065003 (2002)

  49. [49]

    W.-M.Wang, P.Gibbon, Z.-M.Sheng,andY.-T.Li,Tun- able circularly polarized terahertz radiation from magne- tized gas plasma, Phys. Rev. Lett.114, 253901 (2015)

  50. [50]

    Bogatskaya, N

    A. Bogatskaya, N. Gnezdovskaia, and A. Popov, Circu- larly polarized terahertz pulse generation in a plasma channel created by a uv high-intense laser pulse in the presence of a static magnetic field, Phys. Rev. E102, 043202 (2020)

  51. [51]

    Kumar, T

    M. Kumar, T. Kang, S. Kylychbekov, H. S. Song, and M. S. Hur, Simulation study of phase-matched THz emis- sionfromanaxiallymodulatedmagnetizedplasma,Phys. Plasmas28(2021)

  52. [52]

    Tailliez, X

    C. Tailliez, X. Davoine, A. Debayle, L. Gremillet, and L. Bergé, Terahertz pulse generation by strongly mag- netized, laser-created plasmas, Phys. Rev. Lett.128, 174802 (2022)

  53. [53]

    Tailliez, X

    C. Tailliez, X. Davoine, L. Gremillet, and L. Bergé, Ter- ahertz pulse generation from relativistic laser wakes in axially magnetized plasmas, Phys. Rev. Res.5, 023143 (2023)

  54. [54]

    J. Cai, Y. Shou, Z. Gong, H. Wen, L. Han, J. Yu, and X. Yan, Dynamic spin-polarization control of terahertz waves in magnetized plasmas, Matter Radiat. Extrem. 10(2025)

  55. [55]

    J. Cai, Y. Shou, H. Wen, L. Han, Z. Gong, T. Tajima, J. Yu, and X. Yan, Generation of strong terahertz pulse with topologically tunable polarization features, Ultra- fast Sci.5, 0116 (2025)

  56. [56]

    See Supplemental Material for the detailed derivation of the electron velocities, nonlinear current, and the nonlin- earities

  57. [57]

    T. D. Arber, K. Bennett, C. S. Brady, A. Lawrence- Douglas, M. G. Ramsay, N. J. Sircombe, P. Gillies, R. G. Evans, H. Schmitz, A. R. Bell, and C. P. Ridgers, Con- temporary particle-in-cell approach to laser-plasma mod- elling, Plasma Phys. Contr. F.57, 113001 (2015). 8 SUPPLEMENT AL MA TERIAL Derivation of the nonlinear interaction strength To understand...

  58. [58]

    These relations are justified by Fig. 6. As shown in Fig. 6(a), for different output THz frequency, the phase-matching condition requires tenuous plasmaN≪1. Therefore, we simplify Eq. (17) using ω2 1 −ω 2 p ≈ω 2 1 andω 2 2 −ω 2 p ≈ω 2

  59. [59]

    As the output frequency varies from 0.02ω0 to0.2ω 0, the indices of the three waves are all around 1 and have a maximum fluctuation about3%

    Figure 6(b) shows the refractive indices of the three waves at a fixedBc but different plasma densities due to phase-matching at different output THz frequencies. As the output frequency varies from 0.02ω0 to0.2ω 0, the indices of the three waves are all around 1 and have a maximum fluctuation about3%. Therefore, we approximatedn X,ω1 ≈n X,ω2 ≈n X,ω3 = 1t...