2026-08-20 | | Total: 95
Seasonal precipitation anomalies are largely regulated by atmospheric circulation, which dynamical models predict with greater reliability than precipitation itself. Here, we employ a deep learning model that translates dynamical circulation predictions into precipitation estimates. Predictions initialized from March to May consistently indicate a dry anomaly over central China in summer 2026. Retrospective evaluations revealed higher predictive skill in the analogue years, which also tended to feature central equatorial Pacific warming persisting from the preceding winter into summer. This warming favors an anomalous cyclonic circulation over the western North Pacific-South China Sea-South China region, which induces northerly winds and moisture divergence that jointly suppress rainfall over central China. Supporting this mechanism, layer-wise relevance propagation (LRP) independently identifies these northerly winds as the dominant driver of the prediction among all model inputs. Perturbation tests supported this attribution: removing LRP-identified features effectively eliminates the dry anomaly. Our framework thus provides physically interpretable explanations for AI-derived regional climate projections, facilitating evidence-based assessment before observational data become available.
Silicon Photomultipliers (SiPMs) based on NUV-MT technology offer improved photo-detection efficiency and reduced correlated noise compared to earlier designs, making them increasingly attractive for particle and astroparticle physics applications. We present a systematic characterisation of commercially available Broadcom AFBR-S4N series NUV-MT SiPMs over a wide range of temperatures and overvoltages, from -85oC to 23oC and from 10-16V overvoltage (VoV). Key performance parameters are measured, including breakdown voltage, gain, signal-to-noise ratio, dark count rate, afterpulsing probability, and both internal and external optical crosstalk. At 12Vov and 23oC a dark noise of 107.2+/-2.1 kcps/mm^2 is measured, consistent with manufacturer specifications. At 12Vov and -30oC the measurements yielded an average dark noise of 1.69+/-0.12 kcps/mm^2, an average gain of (6.88+/-0.15)x10^6, an average signal to noise ratio of 13.6+/-0.9, an average afterpulsing probability of (0.59+/-0.06)%, an average direct crosstalk probability of (26.49+/-0.28)%, an average delayed crosstalk probability of (0.72+/-0.10)%, and an average external crosstalk probability of (0.10+/-0.04)%.
Horizontal earthquake ground motions exhibit substantial changes in amplitude with changes in orientation. While the effects of style of faulting on central tendency measures of ground motion intensity such as RotD50 are well established in ground motion models, less attention has been given to how faulting mechanism influences the complete directional characteristics of ground motions. This study addresses this gap by conducting a comprehensive probabilistic characterization of ground motion directionality for strike-slip and reverse-slip earthquakes using 3,862 records with moment magnitudes greater than five. Two directionality parameters were examined across 58 oscillator periods ranging from 0.01 to 10 s. The directionality parameters examined consist of pseudo-acceleration response spectral ordinates at all non-redundant orientations normalized by either RotD100 or RotD50 intensities (denoted as η and ν, respectively). Nonlinear regression was used to develop mechanism-specific models for the geometric mean, logarithmic standard deviations, and full probability distributions of these parameters. Results show that strike-slip earthquakes exhibit notably different directionality compared to reverse-slip earthquakes, particularly at longer periods. Strike-slip earthquakes are, on average, 1% to 44% more polarized than reverse-slip earthquakes at periods above 0.5 s, and have RotD100/RotD50 ratios that are 1% to 3.7% larger. The developed models enable engineers and engineering seismologists to obtain mechanism-specific estimates of pseudo-acceleration response spectral ordinates at specific orientations, for applications such as regional seismic risk assessments and orientation-specific probabilistic seismic hazard analyses.
Cosmic Explorer (CE) is a proposed next-generation gravitational-wave observatory that aims to extend our gravitational-wave vision to the edge of the observable universe. With a foundation of technology proven by the National Science Foundation's Laser Interferometer Gravitational-Wave Observatory (LIGO), CE will observe black holes and neutron stars across cosmic time, explore the nature of extreme matter with high fidelity, and probe the nature of gravity and fundamental physics. CE's reference design consists of two widely separated L-shaped detectors to be located in the conterminous United States, one with 20 km arms and one with 40 km arms. As of 2026, CE is in its design and site evaluation phase, with plans to begin observing in the early 2040s together with the Einstein Telescope in Europe. The size of CE observatories---up to an order of magnitude larger than the 4 km LIGO observatories---presents a significant challenge for identifying suitable candidate sites where CE will achieve its science goals, be built within cost boundaries, attract and retain a workforce, and align with community values. In this paper, we report on the design and use of a Python package, the Cosmic Explorer Location Search (CELS) package, to identify cost-favorable sites for CE. For a specified detector location and L-shaped geometry in the conterminous United States, CELS estimates site-preparation costs associated with excavation, land clearing, and land acquisition, while accounting for the scientific effects of detector tilt, arm length, and arm opening angle. After describing the package's methods, we present results for a national-level cost and positioning analysis that complements a recent national suitability analysis. We also discuss how future improvements to CELS will allow deeper, more local studies as the Cosmic Explorer team narrows its list of potential locations.
Antimicrobial resistance (AMR) complicates the treatment of diseases including lung and urinary tract infections (UTIs), which are among the most common bacterial infections worldwide. Motile pathogens can use rheotaxis to swim upstream against fluid flows, potentially promoting access to upper regions of anatomical tracts. However, it remains unclear how antibiotic exposure and resistance influence this transport process. Here, using single-cell tracking microscopy, we investigate how elongation induced by \b{eta}-lactam antibiotics affects the rheotactic migration of E. coli in confined microfluidic channels. Remarkably, we find that rheotaxis can be inhibited 100-fold by antibiotics, even if the susceptible elongated cells remain fully motile. However, resistant bacteria remain short and retain upstream migration under the same conditions. Using genetically engineered bacteria with tunable cell length, we show that the underlying mechanism that governs rheotaxis is the coupling between cell morphology and flow vorticity, where elongated cells are rapidly rotated downstream. Finally, we exploit this length-dependent transport difference to separate short and elongated cells under flow and enrich ampicillin-resistant cells from mixed populations. Together, these results establish bacterial elongation as a key control parameter for rheotactic transport, and provide a proof-of-concept strategy for enriching \b{eta}-lactam-resistant bacteria for potential use in rapid AMR detection.
Despite its renewed attention, the pathway to point-of-care portable low-field MRI systems remains challenging, limiting adoption across research groups. Incomplete documentation limits reproducibility, causing redesign and complicating cross-system comparison. Moreover, non-standardized testing and characterization complicates ethical approval for clinical studies. We present an open-source reference system for portable low-field MRI designed to support replication, reproducibility, benchmarking, and quantitative comparison. The system is fully open source, based on a ~50 mT permanent magnet, and integrated with a cloud-native acquisition platform. Pulseq-based calibration, characterization, and imaging sequences assessed noise level, eddy currents, image-based SNR, and geometric accuracy. Quantitative T1, T2, and B0 mapping sequences were developed and evaluated against reference values. Initial results from independent replications at two sites were compared. The system reached a noise level of 1.4 relative to the thermal noise floor and short eddy-current decay constants of 27-32 us across all gradient channels. Geometric deviations were below 2 mm over the field of view. Image-based SNR were consistent between the independent replications. Measured T1 values closely matched specified values, with an average absolute error of 3.1(1.8)%, while T2 values were overestimated by 10.4(5.8)%. Simulations showed only marginal errors for both quantities, suggesting experimental error sources for T2 mapping. The reference system combines openly documented hardware, software, calibration procedures, phantoms, quantitative MRI, and simulation tools in a reproducible ecosystem, aiming to support cross-site comparability, reproducible research, and collaborative development of future portable low-field MRI technologies.
Electric quadrupole absorption is driven by optical field gradients, not by field amplitudes alone. This distinction is crucial for structured light, where a paraxially small longitudinal field can generate a leading-order longitudinal gradient. For fully vectorial Laguerre-Gaussian modes, Maxwell transversality locks this longitudinal gradient to the transverse structured gradient of the leading paraxial field. This locking is therefore a general electromagnetic constraint, not a focusing-dependent correction that can be removed independently. Decomposing the symmetric-traceless optical field-gradient tensor into spherical components reveals a strongly channel-selective response: the $Δm=0$ transition requires the longitudinal gradient at the retained order, the $Δm=\pm1$ channels receive scalar and vectorial nonparaxial corrections on top of the carrier-scale longitudinal derivative of the transverse field, and the $Δm=\pm2$ channels remain transverse at the same order. These results show that longitudinal optical structure can enter structured-light quadrupole transitions as an indispensable leading-gradient contribution, even when the longitudinal field amplitude itself is paraxially small.
The Jiangmen Underground Neutrino Observatory (JUNO) experiment is a reactor antineutrino detector employing 20 kton of ultra-pure liquid scintillator to determine the neutrino mass ordering and to precisely measure oscillation parameters. The total singles background rate from radioactivity is required to be below 10 Hz in the energy range of 0.7-12 MeV within the fiducial volume for reactor neutrino analysis. The calibration system is designed to characterize the detector energy and position responses, while several of its components are located close to the target and may contribute to the background budget. Therefore, extensive material screening and selection are required to construct a low-background calibration system and to ensure that its contribution remains within the design requirements. In this work, a comprehensive study of the radioactive background induced by the calibration system is presented, including material radioactivity measurements using high-purity germanium detectors and neutron activation analysis techniques, detailed Monte Carlo simulations to evaluate the background, and comparisons with in-situ detector data to validate the predictions. In this data analysis, dedicated spatial selection methods are developed to isolate calibration-related contributions and to suppress the liquid scintillator background. The total radioactivity contribution from the calibration system is estimated to be less than 76 mHz, which satisfies the requirement of 200 mHz (2% of the total background budget). The results from in-situ data are found to be consistent with the expectations based on material assay and simulation within uncertainties. These results demonstrate that the calibration-induced background is well understood, in agreement between data and simulation, and negligible for reactor antineutrino measurements in JUNO.
No interatomic potential has offered universality across chemistry, near-first-principles accuracy and the speed of empirical potentials at once. Here we introduce DPA4C, an equivariant potential whose architecture and compressed CUDA operators are co-designed under deployment constraints to pursue accuracy and efficiency together. Five variants spanning a 49-fold parameter range form the high-throughput end of the measured accuracy--throughput frontier. The largest variant approaches the accuracy of the MACE-Omat models at about two orders of magnitude higher measured throughput. The most compact reduces the energy, force and stress errors of the fastest existing universal MLIP by 61.4%, 48.1% and 34.3% at 1.92 times its saturated throughput. All five variants complete multimillion-atom simulations on a single GPU and run molecular dynamics for 2.048 billion atoms on 1,024 16-GB NVIDIA V100 GPUs at 83.3--91.2% weak-scaling efficiency. Compared with the MEAM empirical potential, DPA4C-Nano reaches 1.8 and 2.5 times the saturated throughput in single-GPU scans on the same V100 hardware for diamond carbon and FCC copper, respectively. DPA4C therefore brings quantum-trained universal accuracy into a regime of speed and system size previously associated with empirical potentials.
Machine-learned exchange--correlation (XC) functionals offer a route to improve Kohn--Sham density-functional theory without incurring the cost of explicitly correlated electronic-structure methods. Their use in production simulation codes, however, requires a well-defined mapping between the learned model and the host-code density representation. We formulate and implement a Skala-1.1 interface in CP2K through the external GauXC library. CP2K supplies the geometry, Gaussian basis, spin-resolved atomic-orbital density matrix, and communicator, while GauXC evaluates the XC energy, atomic-orbital potential matrix, and available nuclear derivatives. The interface accepts both all-electron and valence-only density matrices. The latter may arise from separable dual-space pseudopotentials or molecular effective-core potentials. Implementation errors are isolated from functional differences by comparing the Perdew--Burke--Ernzerhof (PBE) functional evaluated through GauXC with native CP2K PBE. The resulting interface gives consistent energies, forces validated against finite-difference total-energy checks, and force-based molecular-virial diagnostics for representative molecular cases. The dietGMTKN55 benchmark suite is evaluated with an all-electron Gaussian augmented plane-wave treatment for elements up to bromine and def2 effective-core potentials for the heavier elements. The resulting aggregate mean absolute deviation of 1.255 kcal/mol is within 0.020 kcal/mol of the corresponding Skala reference value of 1.235 kcal/mol. This work establishes a validated molecular implementation of Skala in CP2K through GauXC.
Photonic microcavities offer high energy density and low mode volume and attract special attention due to wide range of applications. Subsequent developments in the area of plasmonic nanocavities have also opened up applications despite their high absorptive properties. In the last two decades, photonic-plasmonic hybrid microcavities have been investigated extensively. In this article, we have conducted a review of the literature on photonic and plasmonic systems individually as well as in hybrid systems. Along with the underlying physics, applications such as cavity quantum electrodynamics, lasing, sensing, nonlinear optics and Raman spectroscopy have been reviewed. The outputs of the systems for each applications have been discussed and compared for photonic and hybrid systems. This review will provide a thorough understanding of the effects of optical cavities and plasmonics and their applications.
High-refractive-index particles have recently attracted a growing interest in optical levitation experiments, offering the ability to further engineer optical forces through electromagnetic Mie resonances. Unlike standard silica particles, which are predominantly trapped in the dipole regime and exhibit trap frequencies mainly determined by material density, resonant meta-atoms formed by high-index particles enable qualitatively new trapping behaviors. In this work, we experimentally investigate the trapping of resonant silicon particles in an optical standing wave. A direct comparison of silicon and silica highlights the fundamental differences in their optical force scaling and trapping dynamics. Beyond conventional trapping at intensity-maxima, we demonstrate deterministic and stable three-dimensional trapping of silicon nanoparticles in optical intensity minima, a regime that remains inaccessible for silica particles. Drawing a mesoscopic analogy with blue-detuned atom trapping, our results establish meta-atoms as a versatile approach to further extend the optical manipulation tool box towards accessing novel trapping regimes e.g. in close proximity to a surface.
Scenario planning for rare, high-impact events often requires massive ensembles to stochastically sample relevant trajectories. Although autoregressive weather emulators can efficiently generate such ensembles, isolating trajectories of interest requires sifting through petabytes of data, a challenge that grows exponentially with lead time and rarity. In contrast, a non-autoregressive foundation model like Climate in a Bottle video (cBottle-video) can directly sample trajectories terminating in extremes, avoiding large-ensemble search. We use cBottle-video to generate 1000-member ensembles with start- and/or end-conditioning across three extreme events---the 2021 Pacific Northwest (PNW) heatwave, Superstorm Sandy, and Hurricane Ian. Antecedent 500 hPa geopotential height ($z_{500}$) spread at the free end of end-conditioned ensembles reaches 84--89\% of the final-state spread of start-conditioned ensembles, revealing substantial diversity consistent with each extreme event. For the 2021 PNW heatwave, end-conditioned ensemble members begin uniformly warmer than reanalysis and stay warm, replacing the observed rapid intensification with persistent antecedent heat. For Superstorm Sandy, the leading modes of $z_{500}$ at the antecedent end of the end-conditioned ensemble explain 44\% of the variance in track latitude, and roughly 10\% of ensemble members begin as stronger hurricanes than Sandy. For Hurricane Ian, variation in the first landfall location among end-conditioned trajectories underscores the importance of accounting for intermediate hazard exposure in risk planning.
High-frame rate pixel detectors can produce data volumes that exceed available off-chip bandwidth, yet in many applications only a sparse subset of each frame carries relevant information. Spatially localized events, including diffraction peaks in crystallography, particle hits in tracking detectors, fluorescence spots in biological imaging, and other applications, all require that clusters of above-threshold pixels be identified and extracted from an otherwise featureless background. Conventionally this is performed in software algorithms such as connected-component labeling on full frames, but at MHz frame rates the resulting throughput becomes prohibitive. We present a lightweight hardware architecture that performs peak localization and patch extraction directly in the sensor silicon, transmitting only small pixel patches rather than complete frames. FPGA-based testing on an AMD Alveo V80 validated the synthesizability and timing closure of the peak-finding module in real hardware. The design replaces global connected-component labeling with a cellular automaton that uses purely local, fixed-iteration neighborhood operations, eliminating the label storage and equivalence-resolution logic that make conventional approaches impractical on-chip. We validate the architecture on X-ray Bragg peak detection for far-field high-energy diffraction microscopy and show that every peak found by a software reference is recovered, with equivalent downstream reconstructions. The architecture sustains several-hundred-kHz frame rates in 130nm CMOS and exceeds 1MHz in 28nm.
Calibrating a stochastic car-following model is harder than its deterministic counterpart: the loss itself becomes a random variable, so a favorable random realization can be mistaken for a good parameter vector. This paper develops a structured framework for calibrating stochastic car-following models -- a completeness-controlled synthetic design, a corrected variance-based sensitivity analysis (VBSA), and the minimum-realization (MRMIN) calibration protocol -- across two structurally different stochastic mechanisms, QIDM and IDM2D. We test two claims from deterministic calibration -- that a small number of parameters, and the trajectory itself above all, dominates the sensitivity ranking, and that spacing calibration keeps dominating speed calibration once dynamics are stochastic -- and ask whether a model's noise term can be calibrated on its own. In a balanced synthetic experiment, driving-regime completeness has a mean total-effect index on par with the model's most influential parameter and roughly two orders of magnitude above pair identity, extending rather than reversing the deterministic finding on trajectory-identity dominance. Under MRMIN, calibrating only the noise parameter against a population-wide deterministic fit more than doubles median spacing error across 1644 NGSIM trajectories, but fitting the deterministic parameters per trajectory first and calibrating noise on top recovers it. Spacing calibration remains more cross-dimensionally robust than speed calibration on average, but the deterministic guarantee that this dominance can never reverse is violated in 19-26% of trajectories for both mechanisms. A multi-objective screen in relative-error space then favors joint spacing-speed goodness-of-fit functions over single-dimension spacing calibration. Deterministic calibration guarantees should therefore be re-tested, not assumed, once a model is stochastic.
Alternative divertor configurations (ADCs) must be evaluated under boundary plasma conditions approaching reactor-level values to be considered a reliable, physics-based solution for tokamak power exhaust. Most ADC experiments performed to date were at relatively low exhaust power. This work presents a high-power scenario on the TCV tokamak enabling the study of a wide variety of divertor magnetic shapes under an expanded SOL and power exhaust parameter space. The scenario is characterized by high power levels of electron cyclotron resonance heating ($2.5\,\text{MW}$ fully absorbed in a $\sim1\,\text{m}^{3}$ plasma) at high plasma current (edge safety factor $q_{95}\approx 2.5$), and low upstream separatrix densities ($n_{e,\text{u}}\approx1\times10^{19}\,\text{m}^{-3}$, Greenwald fraction $f_{\text{G}}\approx 0.1$). Stationary parallel heat fluxes up to $100\,\text{MW m}^{-2}$ are measured at the divertor target, an order of magnitude above previous TCV power exhaust studies. The obtained SOL collisionality and Lengyel detachment scaling metric lie within range of values expected in future reactors (SPARC, ITER, ARC).
Stefan Bergman worked in the field of pure and applied mathematics and became known primarily for his work on complex analysis of several variables. He studied the theory of solutions to linear partial differential equations via integral operators and introduced the concept of the Bergmann kernel function. He also worked on problems in electrical engineering, the theory of elasticity, potential theory, and hydrodynamics, particularly supersonic flows. His academic career in Germany was cut short by the Nazi purge, leading to a peripatetic exile through the USSR and France before finally establishing a distinguished career in the United States, notably at Harvard and Stanford. The paper sketches his biography based on archival sources and complements a bibliography of his writings compiled by M.Skwarczy'nski in 1981 with a number of references primarily from his early career.
The near-field shock structures of underexpanded sonic jets exiting into a subsonic coflow are investigated over a range of nozzle pressure ratio (NPR) and coflow-to-nozzle-exit velocity ratio ($U_c$), representative of a propulsive nozzle in subsonic flight. Time-averaged statistics from fully-resolved axisymmetric simulations and inviscid method-of-characteristics (MOC) analysis are used to understand how coflow alters the shock-cell structures, in particular the Mach-disk formation. It is well established that increasing NPR transitions the centerline reflection from regular (characterized by oblique shocks) to Mach reflection (characterized by a near-normal Mach disk). We find that coflow has the opposite influence: a strong coflow shrinks the Mach disk until it vanishes, reverting Mach reflection to regular reflection, so the NPR for this transition increases with $U_c$. This effect has previously been attributed to a reduction in the jet-boundary inclination at the nozzle lip, which confines the lip Prandtl-Meyer fan to a smaller angle and weakens the embedded shock. We show instead that this inclination is determined by the non-uniform pressure the coflow imposes along the jet boundary, which is the primary driver of the shock-structure transitions in coflowing jets. The non-uniform pressure weakens the boundary-reflected compression waves and orients them at shallower angles, so the embedded shock reflects regularly or fails to form. A simulation-informed MOC analysis with this non-uniform pressure boundary condition reproduces the transition behavior with increasing coflow. Coflow also lengthens the first shock cell linearly, which is accurately estimated by a simple correction to Prandtl classical shock-cell length scaling.
High harmonic generation (HHG) in solids offers a pathway to develop compact extreme ultraviolet (EUV) sources crucial for attosecond science and advanced spectroscopy. Here, we demonstrate theoretically that high pressure dramatically enhances HHG in superhard hexagonal tungsten nitride (h-WN6). Compared with solid-state systems at ambient pressure, the reshaped electronic environment under high pressure leads to a unique band-gap widening in h-WN6, which raises the material's damage threshold, allowing the use of stronger laser fields and enabling access to higher-energy bands. This pressure-induced band-gap widening offers a promising strategy to overcome the cutoff limitation of solid-state EUV light sources.
Nanometer-thick fluid films play a critical role in confined multiphase processes, yet the thermodynamics and stability of free gas nanofilms remain poorly understood compared with their liquid counterparts. Here, molecular dynamics (MD) simulations are employed to systematically investigate gas nanofilms confined between bulk liquid phases using Lennard-Jones argon as a model system. The results show that the surface tension decreases exponentially with decreasing film thickness, accompanied by an increasing magnitude of the negative disjoining pressure. Upon thinning, the planar gas film undergoes a distinct morphological transition from a stable planar state to a transient or persistent spherical bubble through the formation and growth of a liquid bridge. The film surface area strongly affects its thermodynamic properties, with larger areas producing stronger thickness dependence and larger deviations from classical density functional theory (cDFT) predictions. The closer agreement between MD and cDFT at smaller surface areas suggests that the discrepancy primarily arises from thermal capillary-wave fluctuations, which are included in MD but omitted in mean-field cDFT. Moreover, at small film thicknesses, the magnitude of the disjoining pressure increases with decreasing temperature, consistent with the enhanced sensitivity of the confined gas phase to thickness variations and contrasting with the trend generally reported for liquid nanofilms. These findings provide molecular insights into the thermodynamics and stability of gas nanofilms, with implications for confined multiphase transport and droplet coalescence.
In this work, particle-resolved direct numerical simulations are performed to investigate flow past the five Platonic solids, which represent a progression in particle sphericity with an increasing number of faces. The simulations cover particle Reynolds numbers in the range 0.1 <= Re_p <= 300 and multiple particle orientations relative to the incoming flow. Based on the numerical data, new correlations are developed for the drag, lift, and torque coefficients. The proposed drag correlation explicitly accounts for both Reynolds number and particle orientation, whereas the lift and torque coefficients are represented by orientation-dependent trigonometric and exponential basis functions whose coefficients vary with Reynolds number. The simulations are conducted using the immersed boundary method, and the resulting drag correlation accurately reproduces the numerical data. The lift and torque correlations capture the principal trends observed in the numerical simulations, including the strong dependence on particle orientation. The proposed correlations provide a computationally efficient framework for incorporating orientation-dependent hydrodynamic forces and torques into Euler--Lagrange and point-particle simulations, enabling a more realistic representation and predictions of non-spherical particle transport in multiphase flows.
We experimentally and theoretically study the hydrodynamic analog of the classical brachistochrone problem: the {\it time-} and {\it energy-minimizing} paths for a spherical particle rolling down an incline within a viscous fluid. We show that in the presence of viscous dissipation, the paths of minima diverge from the classical cycloid, into curves of opposing curvature for time and energy, and are characterized by an effective dimensionless parameter, $St_p$, representing the ratio of the particle's viscous response time scale to its gravitational time scale. Using a generalized variational framework, we show that the fastest path reduces to {nearly straight ramps}, however, beginning and terminating in localized cycloids of curvature, $κ_c \sim St_p^{-2}$. Remarkably, the path of fastest descent on a given energy budget requires navigating a non-monotonic path ({\it``S-shaped''}) with an interior point of inflection. Our findings reveal a unification of temporal and energetic optimality for transport through dissipative media, and expand the celebrated brachistochrone solutions to the hydrodynamic regime.
Magnetic resonance elastography (MRE) is established for noninvasive assessment of liver fibrosis. Conventional abdominal MRE is typically limited to 40-60 Hz. Lower frequencies remain largely unexplored, particularly with regard to hepatic inflammation. We developed ultra-wideband MRE covering 5-80 Hz to investigate frequency-resolved viscoelastic dispersion of the liver and spleen and to identify biomechanical markers of hepatic inflammation. Following phantom validation, nine healthy volunteers and nine patients with inflammatory liver disease were examined at 12 frequencies. Spatiotemporal phase unwrapping and frequency-adaptive wavefield preprocessing enabled reconstruction of shear wave speed (SWS), penetration rate (PR), and loss angle ($φ$). Six rheological models were evaluated. The largest inflammation-associated changes were observed at frequencies below 20 Hz: $φ$ increased by 63% (p<0.001), PR decreased by 37% (p=0.003), and SWS increased by 8% (p=0.008), indicating predominantly dissipative, rather than stiffness-related, changes and a shift toward fluid-like behavior with minor stiffness changes in the lower frequency regime. The rheological springpot model with serial dashpot provided the best fit and revealed distinct dispersion functions for liver and spleen. In patients, springpot elastic modulus increased (101%, p=0.001), while viscosity and springpot power-law exponent decreased (52%, p=0.002 and 58%, p<0.001) suggesting a shift from soft-fluid to stiff-solid liver properties. Ultra-wideband MRE revealed that inflammatory liver disease is associated with property shifts toward stronger dissipation and fluid-like behavior at low frequencies while displaying solid-like behavior at higher frequencies. Ultra-low frequency MRE may provide a diagnostic window into inflammation-associated liver viscoelasticity without full rheological modeling.
Terahertz (THz) radiation provides a powerful platform for ultrafast spectroscopy, imaging, and communication, yet deterministic control over its spatial and polarization structure remains challenging. Here we establish a unified framework for generating and synthesizing vectorial THz beams through coherent control of ultrafast photocurrents in semiconductors. By exploiting quantum interference between one- and two-photon excitation pathways driven by femtosecond vector beams, we demonstrate that the spatial phase and polarization structure of the optical fields can be directly mapped onto the magnitude and orientation of injected currents. This structured charge motion acts as a programmable THz antenna, enabling tailored far-field emission. Beyond forward modeling of THz generation from cylindrical vector beams and full Poincaré beams, we introduce an inverse-design methodology that reconstructs the required current distribution--and corresponding excitation beam profiles--from a desired THz field pattern. This invertible mapping transforms coherent photocurrent control into a systematic design strategy for THz beam shaping. Our results bridge structured light and THz photonics, providing a route toward compact, all-optical, and reconfigurable THz sources with engineered amplitude, phase, and polarization profiles.
The optical modes of mesoscale plasmonic waveguides (MPWs) are significantly affected by nonclassical quantum effects, which can be comprehensively described by the nonclassical electromagnetic boundary condition (NEBC) formulated with the surface-response Feibelman d-parameters. In this paper, a perturbation theory for the nonclassical waveguide modes (NWMs) supported by MPWs under the NEBC is proposed. In this theory, by adopting the classical waveguide modes (CWMs) under the classical electromagnetic boundary condition (CEBC) as the basis functions and treating the NEBC as a first-order perturbation, a general expression of the propagation constant of the NWM with an analytical dependence on the NEBC is derived. This theory transparently reveals the underlying general relation between the nonclassical effects and the propagation properties of the NWMs, thereby providing an effective tool for the understanding and design of MPW devices, as well as for the experimental measurement of the d-parameters.