Search NASA⌕ Search

SEARCH · Search NASA

Results for “HARMONICS”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dimensionally Aligned Signal Projection Algorithms Library

Dimensionally aligned signal projection (DASP) algorithms are used to analyze fast Fourier transforms (FFTs) and generate visualizations that help focus on the harmonics for specific signals. At a high level, these algorithms extract the FFT segments around each harmonic frequency center, and then align them in equally sized arrays ordered by increasing distance from the base frequency. This allows for a focused view of the harmonic frequencies, which, among other use cases, can enable machine learning algorithms to more easily identify salient patterns. This work seeks to provide an effective open-source implementation of the DASP algorithms proposed by Vann et al. (2018) as well as functionality to help explore and test how these algorithms work with an interactive dashboard and signal-generation tool. The DASP library is implemented in Python and contains four types of algorithms for implementing these feature engineering techniques: fixed harmonically aligned signal projection (HASP), decimating HASP, interpolating HASP, and frequency aligned signal projection (FASP). Each algorithm returns a numerical array, which can be visualized as an image. The HASP algorithms are variations of the algorithms originally presented by Vann et al. (2018). For consistency, FASP, which is the terminology used for the short-time Fourier transform (STFT), has been implemented as part of the library to provide a similar interface to the STFT of the raw signal. Additionally, the library contains an algorithm to generate artificial signals with basic customizations such as the base frequency, sample rate, duration, number of harmonics, noise, and number of signals. Finally, the library provides multiple interactive visualizations, each of which is implemented using IPyWidgets and works in a Jupyter environment. A dashboard-style visualization is provided, which contains some common signal-processing visual components (signal, FFT, spectogram) updating in unison with the HASP functions (see Figure 1 below). Separate from the dashboard, an independent visualization is provided for each of the DASP algorithms as well as for the artifical signal generator. These visualizations are included in the library to aid in developing an intuitive understanding how the algorithms are affected by different input signals and parameter selections.

harmonics↗
DOE OSTI2025

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Microscopic examination of rf-cavity-quality niobium films through local nonlinear microwave response

The performance of superconducting radio-frequency (SRF) cavities is sometimes limited by local defects. To investigate the rf properties of these local defects, especially those that nucleate rf magnetic vortices, a near-field magnetic microwave microscope is employed. Local third-harmonic response ( P 3 f ) and its temperature dependence and rf power dependence are measured for one Nb / Cu film grown by direct current magnetron sputtering (DCMS) and six Nb / Cu films grown by high-power impulse magnetron sputtering (HiPIMS) with systematic variation of deposition conditions. Five out of the six HiPIMS Nb / Cu films show a strong third-harmonic response that is likely coming from rf vortex nucleation due to a low- T c surface defect with a transition temperature between 6.3 and 6.8 K, suggesting that this defect is a generic feature of air-exposed HiPIMS Nb / Cu films. A phenomenological model of surface-defect grain boundaries hosting a low- T c impurity phase is introduced and studied with time-dependent Ginzburg-Landau (TDGL) simulations of probe-sample interaction to better understand the measured third-harmonic response. The simulation results show that the third-harmonic response of rf vortex nucleation caused by surface defects exhibits the same general features as the data, including peaks in third-harmonic response with temperature, and their shift and broadening with higher microwave amplitude. We find that the parameters of the phenomenological model (the density of surface defects that nucleate rf vortices and the depth an rf vortex travels through these surface defects) vary systematically with film deposition conditions. From the point of view of these two properties, the Nb / Cu film that is most effective at reducing the nucleation of rf vortices associated with surface defects can be identified. Published by the American Physical Society 2024

Wang, Chung-Yang (ORCID:0000000184299633)↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2024

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗
DOE OSTI2025

Assessing methods in fusion and fitting for time series construction in remote sensing-based earth observations

This study evaluates the comparative performance of spatiotemporal fusion and time-series fitting methods for constructing high-spatiotemporal-resolution remote sensing time-series data. Due to in-class similarity of fusion methods and fitting methods, we employ the Fit-FC (Fitting, spatial Filtering, and residual Compensation) model as a representative fusion method and the linear harmonic fitting model as a representative fitting method. Both Fit-FC and the linear harmonic fitting are widely used for high-spatiotemporal-resolution time-series data construction, and we modify the original Fit-FC model to enable automatic time-series fusion. To ensure data representativeness, we use 3 years (2019–2021) of Harmonized Landsat and Sentinel-2 surface reflectance datasets and Terra MCD43A4 products. Eight experimental regions are selected worldwide to guarantee generalization of the comparative performance between fusion and fitting methods, covering diverse land-use types (cropland, developed land, forest, and grassland) and varying climatological conditions. Time-series of NDVI and surface reflectance are analyzed under both actual observations and simulated data-missing scenarios. The constructed time-series data reveals that (1) the modified Fit-FC and linear harmonic fitting model achieve excellent performance in constructing high-resolution time-series images; (2) the fusion method outperforms the fitting method in constructing time-series of NDVI and surface reflectance images in cropland-, forest-, and grassland-dominated regions; (3) both methods achieve comparable performance in developed-dominated regions; (4) the fusion method is more robust to missing data, and better captures abrupt phenological transitions under conditions of continuous missing data; (5) the fitting method is computationally more efficient, making it suitable for large-scale time-series image reconstruction. This study provides valuable insights for selecting optimal strategies to generate high-resolution time-series images across diverse application scenarios and lays a foundation for extensions to other vegetation indices or land surface variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Cybersecurity Standards, Certification, and Best Practices for DERs

Distributed energy resources (DERs) are becoming increasingly important to the electric grid, including solar energy systems. However, DERs also introduce new cybersecurity risks, including those posed by cloud computing. Standards harmonization is essential for ensuring that DERs are secure and can be safely integrated into the grid. This panel will discuss cyber standards harmonization for solar security. The panel will feature experts from the S2G Program, National Labs and Industry who will discuss the following topics: the cybersecurity risks and future benefits posed by ubiquitous solar energy systems, the development and implementation of cloud-based security solutions for DERs, including solar energy systems, the challenges and opportunities for harmonizing DER cybersecurity standards, and Cyber Informed Engineering and the solar security implementations The panel will also discuss the following specific initiatives: the S2G Program's DER Cybersecurity Framework, UL's DER Cybersecurity Certification Program, and IEEE 1547 Updates. The panel will conclude with a discussion of the future of standards harmonization for DER cybersecurity.

14 SOLAR ENERGY↗
DOE OSTI2024

Role of gas composition in weakened nonlinear standing wave excitation and improved plasma radial uniformity in very-high-frequency asymmetric capacitive Ar/CF 4 discharges

Abstract The higher harmonics generated by nonlinear sheath motion would enhance the standing wave effect, and thus lead to center-peaked plasma density profile in very-high-frequency (VHF) capacitive discharges. In this work, a nonlinear transmission line (NTL) model introduced in Zhouet al(2021Plasma Sources Sci. Technol.30125017) has been extended, with radial transport of various particles and nonlinear sheath motion into account, to investigate the effects of CF 4 fraction α on the nonlinear standing wave excitation and plasma radial uniformity in VHF (60 MHz) capacitively coupled Ar/CF 4 plasmas at 3 Pa. The results indicate that for pure Ar discharges (i.e. α = 0%), the nonlinearly excited harmonics with short wavelength significantly enhance the electron power absorption at the radial center, resulting in a pronounced central-high plasma density profile. As α increases, the high-order harmonics are gradually damped due to the increase of resistance, as well as the longer wavelength caused by thicker sheath thickness. Thus, the radial profile of the electron absorbed power density shifts from center-peak to edge-high. Besides, at the radial center, the electron density and Ar + ion density decrease with α , CF 3 + ion density shows an increasing trend, while F − ion density initially rises and then decreases. Moreover, the density profiles of all the species become more uniform at higher CF 4 fraction, due to the suppressed nonlinear standing wave excitation and the longer wavelength of the nonlinear harmonics.

Physics↗
DOE OSTI2026

Inverse problem in the large momentum effective theory framework

One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed nonperturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using nonperturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of 𝑥 where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of 𝜆 = 𝑧⁢𝑃 𝑧 ∼ 5–15, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a “direct” computation of the 𝑥-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the 𝑥-dependence.

Dutrieux, Hervé [Aix-Marseille Université, Marseil↗
DOE OSTI2025

Nonlinear van der Waals Metasurfaces with Resonantly Enhanced Light Generation

Efficient nonlinear wave mixing is of paramount importance for a wide range of applications. However, weak optical nonlinearities pose significant challenges for accessing nonlinear light–matter interaction in compact systems. Here, we experimentally study second harmonic generation in deeply subwavelength 3R-MoS 2 metasurfaces (<λ/13 thick). Our measurements, supported by theoretical analysis, reveal a complex interplay and coupling between geometric resonances, optical extinction, and exciton-driven strong nonlinear susceptibility dispersion. We further demonstrate >150-fold enhancement in second harmonic signal at 740 nm mediated by the A exciton resonance. Additionally, our theoretical studies predict an enhancement of more than 10 6 in second harmonic generation in <100 nm thick structures exhibiting bound states in the continuum resonance. These findings provide insight into accessing and harnessing the unprecedented 3R-MoS 2 nonlinearities at a subwavelength scale, paving the way to ultracompact nonlinear photonic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Digital Real-Time Simulation and Power Quality Analysis of a Hydrogen-Generating Nuclear-Renewable Integrated Energy System

This paper investigates the challenges and solutions associated with integrating a hydrogen-generating nuclear-renewable integrated energy system (NR-IES) under a transactive energy framework. The proposed system directs excess nuclear power to hydrogen production during periods of low grid demand while utilizing renewables to maintain grid stability. Using digital real-time simulation (DRTS) in the Typhoon HIL 404 model, the dynamic interactions between nuclear power plants, electrolyzers, and power grids are analyzed to mitigate issues such as harmonic distortion, power quality degradation, and low power factor caused by large non-linear loads. A three-phase power conversion system is modeled using the Typhoon HIL 404 model and includes a generator, a variable load, an electrolyzer, and power filters. Active harmonic filters (AHFs) and hybrid active power filters (HAPFs) are implemented to address harmonic mitigation and reactive power compensation. The results reveal that the HAPF topology effectively balances cost efficiency and performance and significantly reduces active filter current requirements compared to AHF-only systems. During maximum electrolyzer operation at 4 MW, the grid frequency dropped below 59.3 Hz without filtering; however, the implementation of power filters successfully restored the frequency to 59.9 Hz, demonstrating its effectiveness in maintaining grid stability. Future work will focus on integrating a deep reinforcement learning (DRL) framework with real-time simulation and optimizing real-time power dispatch, thus enabling a scalable, efficient NR-IES for sustainable energy markets.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Can Restoring Tidal Wetlands Reduce Estuarine Nuisance Flooding of Coasts Under Future Sea‐Level Rise?

Wetland restoration is an increasingly popular nature‐based method for flood risk mitigation in coastal communities. In this study, we present a novel method using hydrodynamic modeling and harmonic analysis to quantify wetlands' ability to reduce future nuisance flooding. The method leverages a hydrodynamic model calibrated to present day data and was run for a range of future sea‐level rise (SLR) and wetland restoration scenarios to quantify changes to tidal harmonic amplitudes and phases. The harmonic constituents are used to generate water surface elevations over a time period of interest (e.g., one year) and compared to critical exceedance thresholds such as levee elevations. Then, changes to nuisance flooding are calculated by counting the number of hours critical thresholds are exceeded under different SLR and wetland restoration scenarios. We applied the method to Coos Bay, Oregon, USA as a test case. We found restoration reduces the number of hours nuisance flooding occurs in downtown Coos Bay from 15 hr (present day conditions) to 0 hr (fully restored condition) under median SLR (82 cm by 2100). Restoration had spatially variable impacts on reducing peak flood elevations with minimal impacts near the estuary mouth and greatest impact 32 km inland. The effectiveness of restoration was heavily dependent on future SLR. Restoration was maximally effective in 2050 under all SLR scenarios, less effective in 2100 under median SLR, and not effective under high SLR. Modeling results suggest increased tidal prism and accommodation space are driving restoration‐associated reductions in tidal amplitudes.

Brand, Matthew W. [Louisiana State Univ., Baton Ro↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: I. Stability

Ion cyclotron emission (ICE) driven by fast ions has the potential to be a reactor relevant diagnostic of confined and escaping fast ions; moreover, controlled experiments in magnetic fusion devices can shed light on similar magnetospheric instabilities below and above the ion cyclotron frequency f ci . Dedicated DIII-D experiments studied instabilities below f ci and at cyclotron harmonics l f ci by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H, D, and 3He at five different values of toroidal field. High-frequency magnetic loops diagnose mode properties. Each beam type is individually injected for an approximate slowing-down time (~100 ms), resulting in a database of 2529 conditions that are analyzed in frequency bands corresponding to ICE harmonics and < f ci global Alfvén eigenmodes. The sub-cyclotron modes are more unstable at low field and in plasmas with large H concentration. For D injection, first harmonic ICE at f cD is more unstable in hydrogen plasmas but, for H injection, the converse is true (ICE at f cH is more unstable in D), suggesting that thermal same-species ions damp fundamental ICE. ICE stability also depends upon the pitch-angle anisotropy gradient and upon the ‘bump-on-tail’ gradient ∂f/∂v. Although the frequency of unstable ICE in these plasmas always occurs near the central cyclotron frequency, the precise value also depends on the Alfvén speed and beam direction, consistent with the idea that instability occurs for modes that simultaneously satisfy ω ≃ 2πlf ci , the magnetosonic dispersion relation, and the Doppler-shifted resonance condition.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2025

Analytical expression of a finite, long, conical canted-cosine-theta coil for particle collider interaction regions

Magnets in the accelerator interaction region (IR) present significant challenges because of high field requirements and limited available space. Conical-shaped magnets offer advantages in these environments by allowing closer placement to the interaction point while maintaining clearance from synchrotron radiation. Interestingly, numerical studies have shown that conical canted-cosine-theta (CCT) designs produce a constant field distribution along the axial direction in the IR quadrupoles for the Electron-Ion Collider (EIC) at Brookhaven National Laboratory. However, the field harmonics generated by conical CCT windings are not yet fully understood. This paper presents an analytical approach to describe the magnetic field produced by a conical surface current and proposes a method for designing conical CCT magnets for accelerator applications. First, we begin with a surface current sheet having a general cosine-theta distribution in spherical coordinates and solve the vector potential using the Green’s function. The magnetic fields generated by the conical current sheet are expressed using associated Legendre polynomials. These results are then related to circular field harmonics and integral field harmonics for designing a coil that produces a pure multipole field. Next, a single layer of the conical CCT winding path is produced based on the cosine-theta current distribution. Finally, the magnetic field quality of dipole and quadrupole conical CCT coils with multiple layers is verified using the Biot-Savart law.

Yang, Ye↗
DOE OSTI2024

Adiabatic approximation and Aharonov-Casher bands in twisted homobilayer transition metal dichalcogenides

Topological flat moiré bands with nearly ideal quantum geometry have been identified in homobilayer transition metal dichalcogenide moiré superlattices, and are thought to be crucial for understanding the fractional Chern insulating states recently observed therein. Previous study proposed viewing the system using an adiabatic approximation that replaces the position dependence of the layer spinor with a nonuniform periodic effective magnetic field. When the local zero-point kinetic energy of this magnetic field cancels identically against that of an effective Zeeman energy, a Bloch-band version of Aharonov-Casher zero-energy modes, which we refer to as Aharonov-Casher band, emerges leading to ideal quantum geometry. Here, we critically examine the validity of the adiabatic approximation and identify the parameter regimes under which Aharonov-Casher bands emerge. We show that the adiabatic approximation is accurate for a wide range of parameters including those realized in experiments. Furthermore, we show that while the cancellation leading to the emergence of Aharonov-Casher bands is generally not possible beyond the leading Fourier harmonic, the leading harmonic is the dominant term in the Fourier expansions of the zero-point kinetic energy and Zeeman energy. Finally, the leading harmonic expansion accurately captures the trend of the bandwidth and quantum geometry, although it may fail to quantitatively reproduce more detailed information about the bands such as the Berry curvature distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microwave microscope studies of trapped vortex dynamics in superconductors

Trapped vortices in superconductors introduce residual resistance in superconducting radio-frequency (SRF) cavities and disrupt the operation of superconducting quantum and digital electronic circuits. Understanding the detailed dynamics of trapped vortices under oscillating magnetic fields is essential for advancing these technologies. We have developed a near-field magnetic microwave microscope to study the dynamics of a limited number of trapped vortices under the probe when stimulated by a localized rf magnetic field. By measuring the local second-harmonic response (𝑃 2⁢f ) at subfemto-Watt levels, we isolate signals exclusively arising from trapped vortices, excluding contributions from surface defects and Meissner screening currents. Toy models of niobium superconductor hosting vortex pinning sites are introduced and studied with time-dependent Ginzburg-Landau (TDGL) simulations of probe/sample interaction to better understand the measured second-harmonic response. The simulation results demonstrate that the second-harmonic response of trapped vortex motion under a localized rf magnetic field shares key features with the experimental data. Here, this measurement technique provides access to vortex dynamics at the micrometer scale, such as depinning events and spatially resolved pinning properties, as demonstrated in measurements on a niobium film with an antidot flux pinning array.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Subharmonic Shapiro steps in depinning dynamics of a two-dimensional solid dusty plasma modulated by one-dimensional nonlinear deformed periodic substrates

Langevin dynamical simulations are performed to investigate the depinning dynamics of a two-dimensional (2D) solid dusty plasma, which is modulated by one-dimensional (1D) nonlinear deformed periodic substrates, and also driven by the combination of the DC and AC forces. As the DC driving force increases gradually, pronounced subharmonic and harmonic Shapiro steps are discovered under the combined modulation of the deformed substrate and the external AC drive. Here, these observed subharmonic and harmonic Shapiro steps are attributed to the dynamic mode locking. The data analysis results indicate that the nonlinear deformed substrate strongly influences these subharmonic Shapiro steps, which can be accurately diagnosed using the fraction of sixfold coordinated particles. Furthermore, the diagnostic of the kinetic temperature clearly indicates the difference between harmonic and subharmonic Shapiro steps, i.e., the particle motion is slightly less synchronized at subharmonic Shapiro steps, caused by the unstable locations of the deformed substrate.

36 MATERIALS SCIENCE↗