Search NASA⌕ Search

SEARCH · Search NASA

Results for “TRANSIENT OSCILLATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theoretical analysis of the role of complex transition dipole phase in XUV transient-absorption probing of charge migration

We theoretically investigate the role of complex dipole phase in the attosecond probing of charge migration. The iodobromoacetylene ion (ICCBr + ) is considered as an example, in which one can probe charge migration by accessing both the iodine and bromine ends of the molecule with different spectral windows of an extreme-ultraviolet (XUV) pulse. The analytical expression for transient absorption shows that the site-specific information of charge migration is encoded in the complex phase of cross dipole products for XUV transitions between the I-4 d and Br-3 d spectral windows. Ab-initio quantum chemistry calculations on ICCBr + reveal that there is a constant π phase difference between the I-4 d and Br-3 d transient-absorption spectral windows, irrespective of the fine-structure energy splittings. Transient absorption spectra are simulated with a multistate model including the complex dipole phase, and the results correctly reconstruct the charge-migration dynamics via the quantum beats in the two element spectral windows, exhibiting out-of-phase oscillations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Stability enhancement module for grid-following converters: Hardware implementation and validation

Low-frequency oscillations limit power delivery of voltage-source converter interfaced resources, for example, wind and solar, operating in grid-following control mode. Our prior research identified the instability mechanism as the coupling between real power and voltage at point of common coupling. In turn, a stability enhancement control was proposed to weaken the coupling by introducing a simple voltage-based feedback loop. The stability enhancement control has been validated via electromagnetic transient simulation testbeds. The objectives of this paper are (a) to provide a more insightful explanation of the enhancement module by examining the impedance of the grid-following converter and (b) to verify the performance of the stability module using a hardware testbed. The theoretic examination demonstrates that stability enhancement module effectively reduces the negative resistance effect of $Z_{qq}$ and thus boosts stability. The experimental research leads to a hardware testbed of a grid-following converter with the ability of demonstrating low-frequency oscillations. Finally, the experimental results demonstrate the efficacy of the proposed stability enhancement control and the readiness of commercial realization of the stability enhancement unit.

42 ENGINEERING↗
DOE OSTI2025

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Online Optimization of NSLS-II Dynamic Aperture and Injection Transient (Facility Improvement Project)

The Facility Improvement Project (FIP) “Methods of online optimization of NSLS-II storage ring concurrent with user operations” proposed and approved in 2018 is now complete. The first goal is the online optimization of nonlinear beam dynamics to increase the beam lifetime and injection efficiency. We have developed a model-independent optimization technique using advanced algorithms. Using this technique, we increased the NSLS-II dynamic aperture by 20% and reduced the amplitude-dependent tune shift by a factor of two. We applied sextupole optimization to the new high-chromaticity lattice, which has been developed to improve the beam stability and to increase the single-bunch beam intensity. We were able to double the injection efficiency and exceed 90%. The second goal of this project is to provide the top-off injection with minimized perturbations of the beamline user operations. To achieve the transparent injection, we optimize the matching of four storage ring injection kickers reducing the perturbation of the stored beam orbit. The first application of the online optimization resulted in a reduction of the amplitude of residual beam oscillation by a factor of six, from 1300 μm down to 200 μm. Further improvement was limited by the timing jitter of the trigger boards. Then, we replaced all trigger boards with a new design and reduced the timing jitter by a factor of five, from 5 ns down to 1 ns. This upgrade resulted in the injection transient reduction from 200 μm to 120 μm. To optimize the full set of kicker parameters, including the trigger timing, amplitude, and pulse width, we upgraded all kicker power supplies with the capability of tunable waveform width. As a result, we have reduced the injection transient to the limit of 60 μm.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Dynamic magneto-chiral instability in photoexcited tellurium

In systems of charged chiral fermions out of equilibrium, an electric current parallel to a magnetic field can generate a dynamic instability that amplifies electromagnetic waves. Whether this mechanism also operates in chiral solid-state systems has remained uncertain. Here we observe signatures of a dynamic magneto-chiral instability in elemental tellurium, a structurally chiral crystal, using time-domain terahertz emission spectroscopy. Under transient photoexcitation in a moderate magnetic field, we observe terahertz radiation with coherent modes that grow in amplitude over time. We present a theoretical model that describes this behaviour based on a dynamic instability of electromagnetic waves interacting with infrared-active oscillators of acceptor states in tellurium, giving rise to an amplifying polariton. These results demonstrate that magneto-chiral instabilities can emerge in solid-state systems and establish a mechanism for terahertz-wave amplification in chiral materials.

42 ENGINEERING↗
DOE OSTI2024

Strain shift measured from stress-controlled oscillatory shear: Evidence for a continuous yielding transition and new techniques to determine recovery rheology measures

Understanding the yielding of complex fluids is an important rheological challenge that affects our ability to engineer and process materials for a wide variety of applications. Common theoretical understandings of yield stress fluids follow the Oldroyd–Prager formalism in which the material behavior below the yield stress is treated as solidlike, and above the yield stress as liquidlike, with an instantaneous transition between the two states. This formalism was built on a quasi-static approach to the yield stress, while most applications, ranging from material processing to end user applications, involve a transient approach to yielding over a finite timescale. Using stress-controlled oscillatory shear experiments, we show that yield stress fluids flow below their yield stresses. This is quantified through measuring the strain shift, which is the value about which the strain oscillates during a stress-controlled test and is a function of only the unrecoverable strain. Measurements of the strain shift are, therefore, measurements of flow having taken place. These experimental results are compared to the Herschel–Bulkley form of the Saramito model, which utilizes the Oldroyd–Prager formalism, and the recently published Kamani–Donley–Rogers (KDR) model, in which one constitutive equation represents the entire range of material responses. Scaling relationships are derived, which allow us to show why yield stress fluids will flow across all stresses, above and below their yield stress. Finally, derivations are presented that show strain shift can be used to determine average metrics previously attainable only through recovery rheology, and these are experimentally verified.

Mechanics↗
DOE OSTI2024

Adaptive Narrowband Damping for Improving Harmonic Stability of Modular Multilevel Converter: Preprint

Harmonic instability events between modular multilevel converter (MMC) and ac systems have been widely reported in recent years. To resolve it, this paper proposes an adaptive narrowband damping control that automatically programs, adds and adjusts damping around the oscillation frequency when detect the oscillation. The paper first presents a low-pass filter design for MMC control loops, which pushes all negative damping of MMC impedance down to the medium frequency range (< approximately 1000 Hz). This makes narrowband damping can be more targeted and easier to design because multiple harmonic oscillations are avoided. An adaptive damping control that uses online oscillation detection is then proposed, which can automatically configure the narrowband damper to provide damping to the MMC at the detected oscillation frequency. In contrast to existing narrowband damping methods, the proposed adaptive narrowband damper adjusts the damping gain (e.g., to zero when system resonance disappears) and the width of the damping range automatically based on continuous monitoring of system resonance conditions. Electromagnetic transient (EMT) simulation results validate the efficacy of the proposed method in two typical MMC-based power systems.

active damping↗
DOE OSTI2024

Adaptive Narrowband Damping for Improving Harmonic Stability of Modular Multilevel Converter

Harmonic instability events between modular multilevel converter (MMC) and ac systems have been widely reported in recent years. To damp harmonic resonance, this paper proposes an adaptive narrowband damping control that automatically programs, adds, and adjusts the damping around the oscillation frequency when an oscillation is detected. First, the paper presents a low-pass filter design for MMC control loops that pushes all negative damping of the MMC impedance down to the medium-frequency range (< ~ 1000 Hz). Then, an adaptive damping control that uses online oscillation detection is proposed, which can automatically configure the narrowband damper to provide positive damping to the MMC around the detected oscillation frequency. In contrast to existing narrowband damping methods, the proposed adaptive narrowband damper dynamically adjusts the damping gain and the width of the damping range based on continuous monitoring of system resonance conditions (e.g., adjust damping gain to zero when the system resonance disappears). Electromagnetic transient simulation results validate the efficacy of the proposed method in two typical MMC-based power systems.

active damping↗
DOE OSTI2025

Ultrafast low-temperature metal–insulator interface phonon dynamics and heat transport in a Pt/Gd 3 Fe 5 O 12 heterostructure

Interfacial thermal and acoustic phenomena have an important role in quantum science and technology, including in spintronic and spincaloritronic materials and devices. Simultaneous measurements of the low-temperature thermal and acoustic properties of a metal/insulator heterostructure reveal distinct dynamics in the characteristic phonon frequency ranges of acoustic and thermal transport. The measurements probed a heterostructure consisting of a thin film of Pt on the ferrimagnetic insulator gadolinium iron garnet (Gd 3 Fe 5 O 12 , GdIG) grown epitaxially on a gadolinium gallium garnet substrate. Ultrafast structural dynamics within the Pt layer were tracked using time-resolved ultrafast x-ray diffraction and analyzed to probe interfacial acoustic and thermal properties. The rapid heating of the Pt layer by a 400 nm wavelength femtosecond-duration optical pulse produced transient structural changes that provided the stimulus for these measurements. Rapid heating produced a broadband acoustic pulse that was partially reflected by the Pt/GdIG interface. Temporal frequencies up to 740 GHz, corresponding to angular frequencies of several THz, were detected in a wavelet analysis of the acoustic oscillations of the strain in the Pt layer. The structural results were analyzed to determine (i) the acoustic damping coefficient and phonon mean free path in Pt at frequencies of hundreds of GHz and (ii) the Grüneisen anharmonicity parameter. The thermal conductance of the Pt/GdIG interface was tracked using the slower, tens-of-picosecond-scale, dynamics of the initial cooling of the heated Pt layer. Analysis using a model based on the Boltzmann transport equation shows that the phonon transmission is lower at the phonon frequencies relevant to thermal transport than for subterahertz regime acoustics.

Acoustic phenomena↗
DOE OSTI2021

The large-scale environment of thermonuclear and core-collapse supernovae

ABSTRACT The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of ∼3.8 Mpc, constrained by the 2M+ + galaxy survey. We limit our analysis to redshift z < 0.036, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility’s Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments, and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe.

Astronomy & Astrophysics↗
DOE OSTI2021

Fault Location Using the Natural Frequency of Oscillation of Current Discharge in MTdc Networks

This paper discusses a novel offline-fault location approach using single ended measurements. The natural dissipation of the circuit parameters are considered for fault location. A relationship between the damped natural frequency of oscillation of the transmission line current and fault location is established. The hybrid dc circuit breaker (dcCB) interrupts the fault current, whereby the transmission line current attenuates under the absence of any driving voltage source. The line capacitance discharges into the fault at a specific frequency of oscillation and rate of attenuation. Utilizing this information, the fault location in a multi-terminal direct current (MTdc) network can be obtained. A three terminal radial model of a MTdc is used for performance evaluation of the proposed method using Power System Computer Aided Design (PSCAD)/Electromagnetic Transients including dc (EMTdc).

42 ENGINEERING↗
DOE OSTI2022

Transient quantum beatings of trions in hybrid organic tri-iodine perovskite single crystal

Utilizing the spin degree of freedom of photoexcitations in hybrid organic inorganic perovskites for quantum information science applications has been recently proposed and explored. However, it is still unclear whether the stable photoexcitations in these compounds correspond to excitons, free/trapped electron-hole pairs, or charged exciton complexes such as trions. Here we investigate quantum beating oscillations in the picosecond time-resolved circularly polarized photoinduced reflection of single crystal methyl-ammonium tri-iodine perovskite (MAPbI 3 ) measured at cryogenic temperatures. We observe two quantum beating oscillations (fast and slow) whose frequencies increase linearly with B with slopes that depend on the crystal orientation with respect to the applied magnetic field. We assign the quantum beatings to positive and negative trions whose Landé g-factors are determined by those of the electron and hole, respectively, or by the carriers left behind after trion recombination. These are $g^{e}_{[001]}$ = 2.52 and $g^{e}_{[1\bar10]}$ = 2.63 for electrons, whereas |$g^{h}_{[001]}$| = 0.28 and |$g^{h}_{[1\bar10]}$| = 0.57 for holes. The obtained g-values are in excellent agreement with an 8-band K.P calculation for orthorhombic MAPbI3. Using the technique of resonant spin amplification of the quantum beatings we measure a relatively long spin coherence time of ~ 11 (6) nanoseconds for electrons (holes) at 4 K.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dynamic signatures of spin-lattice coupling in the layered ferrimagnet Mn 3 Si 2 Te 6

Magnetic van der Waals (vdW) materials exhibit a profound interconnectedness between their various degrees of freedom, pointing to a wealth of potential applications in low-power and high-speed spintronic devices. Recently, light-matter interactions have been leveraged as robust, dynamic pathways to gain control over the properties of vdW magnets through the use of ultrafast pulses of light. Here, we utilize ultrafast photoexcitation to drive coherent lattice oscillations in the layered ferrimagnetic crystal Mn 3 Si 2 Te 6 , which significantly stiffen below the magnetic ordering temperature. We suggest that this is due to an exchange-mediated contraction of the lattice, stemming from strong magneto-structural coupling in this material. Furthermore, simulations of the transient incoherent response uncover the critical role of the spin-mediated electronic relaxation pathways. These results underscore the importance of spin-lattice coupling in vdW magnets and demonstrate a promising strategy for their dynamic optical control via their entangled degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Role of Coupled Feedbacks in the Decadal Variability of the Southern Hemisphere Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Additionally, Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Excited-State Dynamics and Nonlinear Optical Properties of Hyperpolarizable Chromophores Based on Conjugated Bis(terpyridyl)Ru(II) and Palladium and Platinum Porphyrinic Components: Impact of Heavy Metals upon Supermolecular Electro-Optic Properties

A new series of strongly coupled oscillators based upon (porphinato)Pd, (porphinato)Pt, and bis(terpyridyl)-ruthenium(II) building blocks is described. These RuPPd, RuPPt, RuPPdRu, and RuPPtRu chromophores feature bis- (terpyridyl)Ru(II) moieties connected to the (porphinato)metal unit via an ethyne linker that bridges the 4!-terpyridyl and porphyrin macrocycle meso-carbon positions. Pump–probe transient optical data demonstrate sub-picosecond excited singlet-to-triplet-state relaxation. The relaxed lowest-energy triplet (T 1 ) excited states of these chromophores feature absorption manifolds that span the 800–1200 nm spectral region, microsecond triplet-state lifetimes, and large absorptive extinction coefficients [ε(T 1 → T n ) > 4 × 10 4 M –1 cm –1 ]. Dynamic hyperpolarizability (β λ ) values were determined from hyper-Rayleigh light scattering (HRS) measurements carried out at several incident irradiation wavelengths over the 800–1500 nm spectral region. Relative to benchmark RuPZn and RuPZnRu chromophores which showed large βHRS values over the 1200–1600 nm range, RuPPd, RuPPt, RuPPdRu, and RuPPtRu displayed large βHRS values over the 850–1200 nm region. Generalized Thomas–Kuhn sum (TKS) rules and experimental hyperpolarizability values were utilized to determine excited state-to-excited state transition dipole terms from experimental electronic absorption data and thus assessed frequency-dependent β λ values, including two- and three-level contributions for both β zzz and β xzx tensor components to the RuPPd, RuPPt, RuPPdRu, and RuPPtRu hyperpolarizability spectra. Furthermore, these analyses qualitatively rationalize how the βzzz and βxzx tensor elements influence the observed irradiation wavelength-dependent hyperpolarizability magnitudes. The TKS analysis suggests that supermolecules related to RuPPd, RuPPt, RuPPdRu, and RuPPtRu will likely feature intricate dependences of experimentally determined βHRS values as a function of irradiation wavelength that derive from substantial singlet–triplet mixing, and complex interactions among multiple different β tensor components that modulate the long wavelength regime of the nonlinear optical response.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗