2026-08-14 | | Total: 104
The orbital angular momentum (OAM) of light is a discrete, unbounded degree of freedom that underpins mode-multiplexed communications and high-dimensional quantum photonics. Yet, dynamic OAM control remains dependent on bulky free-space optics or cascaded architectures that separate switching from wavefront shaping, hindering nanoscale integration. Here, we engineer artificial van der Waals crystals from rhombohedrally stacked (3R) MoS$_2$, in which spatial control of the local crystal orientation imprints a nonlinear geometric phase onto the second-harmonic (SH) field, enabling background-free generation of SH vortex beams in an ultrathin (46 nm) van der Waals platform. Leveraging the C$_{3v}$ symmetry of 3R-MoS$_2$, we demonstrate monolithic, all-optical switching with sub-optical-cycle precision between Hermite-Gauss-like and Laguerre-Gaussian vortex SH beams with opposite topological charges ($l=\pm1$). Our results establish artificial 3R-MoS$_2$ crystals as a monolithic platform for the generation and all-optical reconfiguration of nonlinear structured light at the nanoscale, advancing active nanophotonic sources for integrated classical and quantum photonic technologies.
Complex multipath environments are usually avoided in wave-based information processing because repeated scattering creates many interfering propagation paths, obscuring controllability and generating extreme sensitivity to perturbations. The addition of nonlinear mechanisms fundamentally alters the wave-control landscape by breaking the superposition principle that underpins most wave-management strategies. Here, we show that these two apparent impediments -- multipath complexity and nonlinearity -- can instead be harnessed as key resources for physical optimization. We demonstrate an in-situ adjoint optimization protocol in a wave-chaotic platform incorporating a single localized nonlinear defect, in which the system itself performs both the forward and the adjoint propagations required for gradient evaluation. Recurrent multipath returns repeatedly expose the wave to the defect, producing from a minimal hardware a rich nonlinear input-output map with many pathway-mediated degrees of freedom. At the same time, a suitable adjoint excitation enables direct extraction of the sensitivities from measurements alone, without a digital twin or conventional numerical backpropagation. We experimentally validate the protocol on a minimal nonlinear multipath platform composed of incommensurate coaxial cables connected via T-junctions, one of which hosts a diode-loaded cavity. Our approach opens a route to adaptive wireless communications, imaging and analog intelligence in complex, partially unknown environments where conventional modeling is impractical.
AI weather emulators benefit from clear objectives and metrics, which have led to the rapid development of models that outperform traditional benchmarks. In contrast, long-term climate emulators must reliably reproduce forced responses over months to centuries, while relying on training objectives that span a small number of model time steps. We assess autoregressive, full-depth ocean emulators using data from the midHolocene experiment of a numerical climate model to examine their skill in responding to surface forcings from an in-distribution, out-of-sample climate. We demonstrate that these emulators generalize to new orbital forcings, reproducing the spatial structure of the large-scale response as well as changes in seasonal patterns and in the spatial structure of ocean variability, while underestimating their amplitude. Baselines that infer the ocean state directly from the boundary forcings also recover much of the large-scale pattern, but only near the surface, and capture neither the seasonal nor the variability changes, indicating that these require some representation of dynamics. Despite these successes, the emulators fail to reproduce the slow, internally driven evolution of the ocean interior. We then show that the emulators' total forced response is well reconstructed by linearly composing their independent responses to each forcing component. Tracking response across training epochs, we find that convergence on mean state metrics in the training climate does not guarantee that the emulators capture the dynamics necessary for a skillful response. Together, these experiments establish the midHolocene as a controlled, ground-truthed setting for diagnosing forced-response failures before emulators are pushed to out-of-distribution climates.
Deep learning surrogates for RANS flow prediction around airfoils face two persistent bottlenecks. A single neural architecture cannot simultaneously resolve sharp near-wall boundary layers and smooth far-field potential flow. Additionally, force prediction is undermined by the numerical instability of computing wall-normal velocity gradients from continuous-field approximations. Both of these points are addressed with a DD-RNO (domain-decomposed routed neural operator), combining a spectral geometry encoder with two physics-guided innovations: (a) a differentiable domain routing mechanism that partitions the flow field into inviscid, boundary-layer, and wake regimes---dispatching query points to specialized regional decoders, and (b) learned canonical quadrature (LCQ), which replaces unstable pressure integration with flow-conditioned, learned integration weights that predict lift and drag directly from surface pressure. On the AirfRANS benchmark, DD-RNO cuts velocity field mean-square error (MSE) by 17$x$ ($u_x$) and 12$x$ ($u_y$) over the strongest baseline, widening to 23$x$ under out-of-distribution Reynolds extrapolation---evidence that the routing mechanism generalizes with the physics it encodes rather than merely fitting the training distribution. LCQ reduces drag MSE by 7.5$x$ relative to conventional pressure integration and raises drag rank correlation from $ρ= 0.250$ to $ρ= 0.997$. Ablations confirm that both components are indispensable to performance: removing domain routing increases velocity error by 8.2$x$, and removing LCQ increases relative drag error more than 40-fold. At ~144 ms per sample---a 10,000$x$ speedup over conventional RANS solvers---DD-RNO offers a surrogate accurate and fast enough for real-time aerodynamic design and optimization loops.
The macroscopic ideal-MHD stability of an axisymmetric mirror device with sonic levels of plasma rotation is analyzed by approximating the plasma equilibrium as a rotating theta pinch possessing an artificial gravity. An eigenmode equation is derived that governs the stability of the equilibrium to small perturbations in the case of an arbitrary plasma angular velocity profile. The stability of the m=1 and m=2 modes is investigated. The plasma is found to be stable to these two modes provided that it is sufficiently short in the axial direction. The critical axial length of the device below which the modes are stabilized first decreases with increasing plasma rotation, attains a minimum value when the rotation is roughly sonic, and then increases with increasing plasma rotation. The value of the plasma rotation off the magnetic axis is found to have a significantly stronger effect on the stability of the modes than the value on the magnetic axis.
We report the formation of an original dome-shaped air-water cavity at the impact site of closely-packed plunging water jets. The injector assembly consists of concentric rings of circular jets, similar to a shower-head configuration. We infer from high-speed imaging, LASER-induced fluorescence, and optical phase detection probes that the cavity consists of multiple stems and sheets of air that stretch, retract and pinch off to produce bubbles. We also illustrate various bubble production mechanisms in the cavity, along with coalescence and bubble breakup scenarios in the subsequent bubble cloud. Thereby, we describe the wide range of bubble sizes generated by such jets, from a few micrometers up to about $5$~mm. Measurements of the cavity size are carried out for varying impact velocity, number of jets, and jet spacing. We propose that this distinctive air-water cavity is generated by liquid entrainment when successive rings of plunging jets progressively mix with the pool water.
A lattice Boltzmann scheme for the three-dimensional compressible Navier--Stokes--Fourier equations, derived automatically from the declared system by a symbolic compiler, is validated against exact solutions and published reference data. The declared system carries the viscous stress and the heat flux as transported state, and is discretized on a D3Q7 lattice in single precision. Against the exact Sod and Becker solutions the captured shock thickness converges at first order. On the supersonic Taylor-Green vortex at $M_0 = 1.25$ the scheme at $512^3$ matches the reference dilatational dissipation more closely than seven compared solvers, by thirty percent over the next best. Every operator in this solver, the generated collision and constitutive closure and the added shock sensor alike, reads only the cell it acts on, and data reaches a neighbor only by streaming along the lattice characteristics.
Metallic nanostructures confine electromagnetic fields at subwavelength scales, making them attractive as plasmonic nanoantennas. At these scales, the response of metals becomes nonlocal, and the hydrodynamic model is widely used to capture this response. However, existing solvers provide only the response to a prescribed excitation and do not directly reveal the intrinsic resonances of the structure. This work extends the characteristic mode analysis to plasmonic nanostructures to enable excitation-independent modal analysis of their resonant behavior. For simple metals, the coupled hydrodynamic and volume integral equations are reduced to a single hydrodynamic volume integral equation in terms of the induced current. The equation is discretized and cast as a generalized eigenvalue problem within the characteristic mode analysis framework, whose solution yields the characteristic mode currents and modal significance curves of the structure. The proposed framework is validated through three metallic nanostructures: a nanosphere, a nanorod, and a nanodimer. The results show that the method identifies the intrinsic resonances of each structure, including resonances not excited by a given source and additional resonances arising from the nonlocal response, which are absent in local models. The proposed framework provides physical insight into the modal mechanisms of plasmonic nanostructures and serves as a practical tool for their analysis and design.
Bond-dependent magnetic interactions, exemplified by the Kitaev model, are known to arise from the interplay between spin-orbit coupling(SOC) and specific coordination geometries, yet their existence has so far been predominantly associated with edge-sharing octahedral systems. Whether analogous interactions survive in edge- sharing tetrahedral environments - relevant to iron-based superconductor parent compounds - remains an open question. Here, we construct a Kitaev-like model for monolayer FeTe and FeSe and demonstrate the presence of a previously unrecognized bond-dependent Ising-type interaction, induced jointly by chalcogen-mediated SOC and the tetrahedral crystal-field geometry. By mapping the first-principles calculations derived magnetic anisotropy energy mapping across representative linear magnetic orders, we disentangle the bond-dependent contributions from single-ion anisotropy. We reveal that the Kitaev-like interaction dominates the magnetic anisotropy in FeTe, whereas in FeSe, it fiercely competes with a single-ion anisotropy of opposite sign. The resulting noncollinear local anisotropy axes generate intrinsic single-site spin frustration, providing a microscopic mechanism for magnetic disorder beyond isotropic exchange models. Our results establish edge-sharing tetrahedral magnets as a new setting for bond-dependent interactions and extend the scope of Kitaev physics beyond octahedral coordination.
Aerodynamic unsteadiness is inherent to the operation of many engineering applications, especially those that involve large-scale rotating blades, such as modern wind turbines. Over the past few decades, wind turbine rotors have grown rapidly and are now exceeding 200 meters in diameter, causing them to operate in conditions where limited empirical data are available, and where models have not been validated. Here, we use a highly pressurized wind tunnel to probe these conditions and evaluate the quasi-steady approximation, often used to simplify the modeling of the unsteady aerodynamic response of an airfoil to slow changes in inflow conditions. We find that the quasi-steady approximation is not valid for a large range of frequencies where it is normally applied. This finding suggests that the loading on wind turbine blades will be significantly underestimated when using conventional models, especially near stalling conditions, which modern wind turbines typically encounter.
Intermittent departures from nominal Kelvin--Helmholtz shedding can produce rare and pronounced drag excursions in transitional airfoil flow. We examine these events using two-dimensional direct numerical simulations of flow over a NACA0012 airfoil at an angle of attack of $5^\circ$, a freestream Mach number of $0.4$, and chord-based Reynolds numbers of $5\times10^4$ and $5\times10^5$. At the lower Reynolds number, event-resolved analysis shows that individual primary vortices are released from the separated shear layer through the eruption of wall-generated, opposite-signed secondary vorticity. Each eruption interrupts the connection between a developing primary vortex and its feeding shear layer, releasing the vortex downstream. During nominal shedding, the vortex reaching the trailing-edge region is associated with a single such release and remains sufficiently isolated to pass the trailing edge without strong collective interaction. Extreme events instead arise through clustered vortex release, in which several secondary-vorticity eruptions occur within a short interval and produce a compact group of primary vortices with small initial streamwise spacing. Differential convection further reduces their spacing and promotes strong near-trailing-edge interactions, where the combined pressure footprint of these vortices produces a localized suction peak and a sharp increase in drag. These interactions range from prolonged deformation and filamentation to rapid core coalescence. Similar compact vortex organization and near-trailing-edge interactions are recovered at $Re=5\times10^5$, indicating that the downstream event pathway persists despite the smaller vortical scales.These findings suggest that controlling vortex-release timing through secondary-vorticity dynamics may provide a route to disrupt clustered release and mitigate extreme aerodynamic loading.
Time-domain inverse design in photonics is known to be suitable for maximizing the efficiency of optical devices over broad frequency ranges. Objectives commonly used in this context include time-integrated field quantities derived from Poynting's theorem, such as energy, flux, or dissipated power, which can be directly linked to integrated frequency-domain responses via Parseval's theorem. While computationally efficient, these objectives measure only the total response over the targeted bandwidth and, as we demonstrate, are insufficient to capture undesired in-band ripple, narrow spectral features, or sidelobes. We overcome this limitation by introducing a time-domain metric quantifying such spectral variations based on the weighted long-lag autocorrelation energy of the optical response. We incorporate this metric into an FDTD-based topology-optimization framework and demonstrate its beneficial effect on the example of inverse designing one-dimensional dielectric Bragg mirrors via the time-domain adjoint method.
Alchemical relative binding free energy (RBFE) calculations are limited by the fixed-charge approximation of classical force fields. Hybrid machine learning interatomic potential/molecular mechanics (MLIP/MM) schemes correct ligand strain, but under mechanical embedding still describe ligand--environment electrostatics with static point charges. Electrostatic embedding schemes coupling machine-learned charges to the MM environment have been proposed and validated against QM/MM for simple systems, but not tested in a production alchemical workflow. We take the electrostatic embedding scheme of Semelak et al.\ and evaluate it on protein--ligand RBFE. We trained a TensorNet2 model, \texttt{AceFF-2-RESP-1}, on $10^{6}$ conformations from the AceFF dataset, jointly predicting energies, forces and Restrained Electrostatic Potential (RESP) charges. We chose RESP over MBIS for commensurability with the AMBER-family force field it couples to. The predicted charges enter the short-range direct-space part of the particle mesh Ewald sum, with Thole damping to prevent polarization catastrophes during alchemical transformations. We tested the scheme across five targets from the Wang et al.\ benchmark set, fixed in advance by a prior study, with three replicates per edge and matched protocols. Electrostatic embedding improved every accuracy and correlation metric for TYK2 ($ΔΔG$ RMSE $0.86 \rightarrow 0.45$~kcal/mol against GAFF2), but performed comparably to the classical and mechanical-embedding baselines for CDK2, thrombin, p38 and JNK1. Standard single-molecule energy and charge benchmarks were not good predictors of this target-dependent outcome. TYK2 combined good $ΔΔG$ accuracy with the lowest force error on the Schrödinger benchmark, but this pattern did not hold for the other targets.
Phenomenological models of non-Fourier heat conduction often lack a strict microstructural foundation, leading to ambiguities when modeling complex heterogeneous materials. In this study, we derive a continuum heat equation beyond Fourier's law using spatial volume averaging for a two-component system. We analytically prove that the experimentally observed static and dynamic thermal diffusivity arise directly from the distinct material properties, concluding that heterogeneous media are inherently over-diffusive. The resulting heat equation is thermodynamically compatible, and the microstructural origin allows the calculation of non-Fourier transport coefficients. Furthermore, we demonstrate that finite-sample boundaries introduce higher-order spatial non-localities, thereby explaining the size dependence of over-diffusion. We validate the model against experimental data across metal and carbon foams, rocks, and metal-organic frameworks.
Geometric parametric instability (GPI) is the resonant growth of discrete spectral sidebands enabled by longitudinally periodic multimode evolution and has been studied primarily in graded-index fibers. Here we show theoretically and numerically that GPI can occur in gas-filled nonlinear multipass cells (MPCs). By mapping a mode-matched MPC onto an equivalent waveguide, we derive a Floquet quasi-phase-matching condition governed by the single-pass Gouy-phase imbalance of the signal--idler pair relative to the pump pair. The theory predicts the small-signal gain and bandwidth. A pump-depleted coupled-mode model (CMM) further relates the maximum converted fraction to the residual phase mismatch. The CMM predicts multiple geometrically tunable sideband pairs associated with different radial indices and Floquet orders. For argon at $5$~bar, varying the cavity geometry shifts the sideband detuning from approximately $96$ to $46$~THz when the $p_{\mathrm{s}}=1$, $h=0$ branch is considered. A truncated multimode generalized nonlinear Schrödinger equation (MMGNLSE) model is used for numerical simulations with a semiclassical stochastic seed corresponding to one photon per spectral mode. The MMGNLSE simulations reproduce the predicted sideband frequencies and reveal pump depletion and competition among the retained radial channels. GPI in MPCs may therefore limit spatial beam quality in nonlinear pulse compression while providing a tunable mechanism for broadband multicolor generation.
SiN photonic crystal (PhC) cavities are promising platforms for nonlinear and quantum photonics because of their wide transparency window, CMOS compatibility, and negligible two-photon absorption. However, realizing high-Q/V cavities remains challenging because of the relatively low refractive index of SiN. Here, we employ inverse design to optimize a two-dimensional SiN PhC cavity and experimentally demonstrate a quality factor of approximately 80,000, the highest reported for a near-stoichiometric SiN 2D PhC cavity. Furthermore, both second- and third-harmonic generation are observed from the same cavity, providing experimental evidence of strong optical confinement and large Q/V. Our results establish inverse-designed SiN PhC cavities as a promising platform for nonlinear photonics and future heterogeneous integrated photonic devices.
High-precision phase measurement of microwave fields underpins a wide range of applications, including wireless communications, distributed radar, plasma diagnostics, and antenna metrology. Existing Rydberg-atom-based approaches, however, often face trade-offs among phase resolution, measurement range, and system complexity. Here we demonstrate a Rydberg-atom-based microwave Mach-Zehnder-type interferometer using a dual-local-oscillator configuration. The two local oscillators establish two coherent interferometric pathways in the Rydberg medium. Their coherent mixing with the signal field produces an interferometric intermediate-frequency output governed by a phase-to-intensity transfer characteristic that enables critical-point enhancement. This scheme supports direct phase retrieval with a resolution exceeding $0.1^\circ$ and unambiguous full $360^\circ$ phase coverage with the reconfigurable dual-LO architecture. Moreover, near the critical interference point, the system exhibits a sharply enhanced phase-to-amplitude transduction, where weak amplitude variations are converted into pronounced phase responses, yielding a sensitivity enhancement exceeding 25 dB. Besides, the same interferometric transfer mechanism enables microwave propagation-distance and polarization metrology, achieving a propagation-distance precision below 20 $μ$m at 5.7 GHz together with a polarization-angle resolution exceeding $0.1^\circ$. This approach eliminates the need for complex optical configurations and lock-in detection, providing a simple, scalable, and reconfigurable Mach-Zehnder-type quantum microwave interferometry framework for multifunctional high-precision microwave metrology.
Earth's feedback parameter measures how the Earth system responds to forcing and is inversely proportional to climate sensitivity. Sea-surface temperature (SST) patterns can modulate the value of the feedback parameter. Differences between observed and simulated SSTs have raised the question how the observed SSTs evolution impacts the global feedback. The standard method for estimating this effect uses observed SSTs prescribed to an atmospheric model with fixed pre-industrial atmospheric forcing. This method makes two assumptions: first, that the observed SSTs capture all relevant effects from the forcing, so that prescribing a time-varying forcing is unnecessary; second, that the temporal variations in the feedback parameter are driven by the evolving SST pattern and can be estimated via moving-window regressions. We test these assumptions by running controlled experiments in which SSTs from fully-coupled historical simulations are prescribed to an atmospheric model. We find that the prescribed-SST experiments fail to capture the coupled feedback evolution. This is explained by two effects: First, the absence of prescribed atmospheric forcing, and second, statistical noise arising from the computation of moving-window regressions. We find no evidence of any significant relationship between evolving SST patterns and changes in the feedback time series in a 4000-year pre-industrial control simulation. Any trends in the feedback parameter detected on timescales shorter than ~100 years are indistinguishable from statistical noise, making their attribution to the evolving SST pattern extremely difficult. Our results imply that prescribed SST simulations offer limited potential for inferring temporal changes in Earth's parameter over the observational period.
In this work, a three-dimensional computational fluid dynamics model is developed for a vortex-stabilized microwave CO$_2$ plasma reactor operating over the pressure range of 100-400 mbar. The model combines experimentally constrained, emission-based plasma sizes and volumetric heat-source distributions with thermally dominated finite-rate heavy-particle chemistry for a multi-component mixture. Turbulent flow and transport are described using the SST k-omega model. The model reproduces the measured radial gas-temperature profiles in the plasma core and the non-monotonic pressure dependence of CO$_2$ conversion, including a maximum at 150 mbar and a pronounced decrease at 400 mbar. A vortex-driven recirculation region redistributes gas upstream. Turbulent mixing and cooling are strongest near the upper reactor boundary, but their contribution decreases as pressure increases. The pressure dependence of conversion is determined by the competition between CO$_2$ dissociation and CO recombination. CO$_2$ direct dissociation reaction dominates in the high-temperature plasma core, whereas O-assisted conversion reaction contributes near the plasma edges and in the surrounding hot region. At 150 mbar, enhanced CO$_2$ dissociation is accompanied by limited CO loss, resulting in the highest conversion. With pressure increasing to 400 mbar, slower cooling and more frequent three-body collisions promote CO recombination in the afterglow, causing more than 60% of the CO formed near the plasma to be lost downstream. Moreover, additional CO loss occurs in the upper region of the reactor at higher pressures because of the reduced cooling rate.
Extreme manufacturing with ultrashort-pulse (USP) lasers at the physical limit of precision and efficiency requires understanding ablation dynamics on the picosecond-to-nanosecond timescale. Pump-probe reflectometry (PPR) provides direct access to photomechanical spallation through Newton ring (NR) interference, but this signature vanishes when the spallation layer becomes optically opaque or the ablated material strongly attenuates the probe. Here, we combine PPR with phase-sensitive interferometric pump-probe (PPI) measurements to track the spallation layer in bulk steel, aluminum, and gold. PPI resolves the propagating layer even when the reflected probe signal is suppressed by >95%. Joint PPR/PPI analysis with transfer-matrix modelling (TMM) yields spallation layer thickness, vapor layer absorption, and the layer disintegration times. These quantities are key determinants of the energy coupling of subsequent pulses in GHz burst processing.
Many traditional models of electoral dynamics emphasize opinion change, social influence and persuasion mechanisms. In contrast, contemporary polarized societies often exhibit relatively stable partisan blocs coexisting with an electorally relevant mobile population. We propose a minimal dynamical model in which electoral competition is governed not by direct persuasion between opposing blocs, but by the dynamics and allocation of a mobile electoral interface. The electorate is partitioned into two stable partisan blocs and a mobile fraction, while direct transitions between the polarized blocs are strongly suppressed. Mobile voters are allocated between the competing blocs through a Fermi-like probabilistic rule governed by rejection asymmetries. Analytical calculations show that the electoral susceptibility to political shocks is proportional to the stationary size of the mobile electoral interface, identifying electoral mobility as the key quantity controlling the macroscopic response of polarized systems to external perturbations. The model naturally predicts a continuum of mobility regimes, ranging from frozen polarization to highly responsive electoral states characterized by a broad mobile interface. These results suggest that, in strongly polarized elections, aggregate electoral changes may be governed primarily by fluctuations at the mobile electoral interface rather than by large-scale ideological conversion.
Focusing of synchrotron light across extreme ultraviolet (XUV) and soft X-ray regimes is increasingly desired for photoemission-based techniques where reduced beam width gives access to smaller samples such as microscopic single crystals and functioning two-dimensional (2D) heterostructures and devices. Many existing focusing methods, however, are not able to take full advantage of the synchrotron beam due to limited photon energy range or low transmission. Modern capillary optics have enabled achromatic, high transmission focusing of XUV and X-ray light. Here, we present a detailed characterisation of such an achromatic capillary optic installed at the AU-SGM4 beamline for spatial- and angle-resolved photoemission spectroscopy (ARPES) experiments at the ASTRID2 light source. The transmission of the capillary as a function of photon energy is given, and the dependence of the beam width, position, and transmission are measured against the source size. Analysis of the far-field image of the beam allows for slope errors on the inner surface of the capillary to be overcome by selectively aperturing the beam, resulting in a minimum beam width of 900 nm measured in a photoemission geometry.
Inverse problems are ubiquitous in physics, chemistry, and engineering, arising when reconstructing hidden quantities from indirect measurements. A key example is the analytic continuation of imaginary-time correlation functions (iTCFs) to the real-frequency domain. This process requires an inverse Laplace transform, which is inherently ill-posed and highly sensitive to small input variations. Recent neural network (NN)-based methods have shown promising results by learning mappings from imaginary-time to real-frequency spectra, often outperforming traditional techniques such as maximum entropy. However, because the problem is ill-posed, many spectra fit the same iTCF. Regression-based approaches output a single solution, which approximates an average over the true solution space, and therefore fail to capture the full distribution of plausible power spectra. To address this issue, we introduce a diffusion-based framework for analytic continuation that learns the distribution of spectra consistent with a given iTCF. It offers two key advantages. First, it quantifies uncertainty directly from the learned distribution. Second, by analyzing the spread and structure of this distribution, we can quantitatively assess the intrinsic hardness of each inversion problem. We measure this hardness with a new metric, the uncertainty pseudo-volume. Applying the framework to an iTCF from a path-integral molecular dynamics simulation of liquid parahydrogen, we obtain the self-diffusion coefficient with an error bar and flag a secondary high-frequency peak as a possible spurious artifact. In contrast to previous attempts at uncertainty quantification, our generative approach rests on a concrete probabilistic basis, providing a more theoretically grounded measure of confidence in the reconstructed power spectra.
Fast pixelated detectors are driving 4D-STEM toward microsecond dwell times, a regime in which the finite response of the scan deflection coils becomes comparable to the dwell time itself. Using direct probe imaging and sub-frame diffraction analysis, we document a significant intra-dwell scan-coil delay that systematically smears the recorded signal anisotropically along the fast scan direction, with a settling timescale of several tens of microseconds. We present a phase-correlation-based sub-frame alignment procedure that measures and corrects this smearing, and we assess its impact on focused and defocused 4D-STEM reconstructions over a range of scan step sizes. The correction restores signal across a broad range of spatial frequencies, with the largest gains at the large step sizes required for low-dose biological imaging. Because it operates on existing data with no modification to the microscope, the method offers a practical route to recovering signal that would otherwise be lost to scan-coil delay.
Dual-comb spectroscopy, constructed from two frequency combs with slightly different repetition frequencies, enables real-time and high-precision measurements without mechanical scanning. In the terahertz spectral range, dual-comb techniques provide a powerful tool for high-resolution and rapid spectroscopy. Quantum cascade lasers (QCLs), owing to their compact footprint and favorable size, weight, and power, have emerged as promising sources for THz dual-comb systems. In QCL frequency combs, the repetition frequency generated through intrinsic four-wave mixing is typically equal to the cavity round-trip frequency, corresponding to the fundamental comb. Although this repetition frequency can be tuned via current and temperature control, the accessible tuning range remains limited. Recently, harmonic frequency combs, whose repetition frequencies are integer multiples of the cavity round-trip frequency, have attracted increasing attention, offering enhanced single-line signal-to-noise ratios and providing direct insight into the strong nonlinearity of QCLs. Here, we demonstrate a reconfigurable multi-harmonic dual-comb system realized on a single self-detected THz QCL platform. By precisely controlling the driving current and thermal conditions, we achieve and switch between multiple dual-comb configurations, including fundamental-fundamental, fundamental-second-harmonic, second-harmonic-second-harmonic, and second-harmonic-third-harmonic dual-combs. These results establish harmonic order as an additional degree of freedom for dual-comb operation within a single laser system. The demonstrated platform enables high-resolution and high-sensitivity measurements while significantly simplifying the system architecture. This reconfigurable multi-band dual-comb approach based on a single QCL opens new opportunities for compact THz spectroscopy and frequency metrology.