Papers

High harmonic generation (HHG) is a highly nonlinear optical process in which radiation from a strong driving field is up-converted into its high-order harmonics. In atomic systems, this nonlinearity manifests itself through the intensity scaling of the emitted harmonics with the driving field strength. Despite the highly nonlinear nature of HHG, when the driving field is prepared in a classical Gaussian state and atomic depletion remains negligible, the quantum statistical properties of the generated harmonics retains classical Gaussian quantum statistics. Driving HHG with bright squeezed vacuum (BSV) light challenges this paradigm, as its enhanced field fluctuations can modify the statistical properties of the generated harmonics. In this work, we investigate the conditions under which BSV-driven HHG gives rise to non-classical Gaussian states, and identify the regimes where this Gaussian description breaks down. For bichromatic driving by a strong coherent field at frequency $\omega$ and a perturbative BSV field at $2\omega$, the even-harmonic response is approximately linear in the BSV quadrature, leading to non-classical multimode Gaussian entanglement in the harmonic field. We show that this state can be described as a distributed collective squeezed mode over the even-harmonic manifold, and characterize its covariance matrix, entanglement structure, and quantum teleportation fidelity as an operational benchmark. Our results highlight the potential of non-classically driven HHG as a platform for engineering Gaussian and non-Gaussian states of light in the extreme ultraviolet regime.

The empirical bitstring distribution is the most accessible observable on Rydberg-atom arrays, but the bipartite von~Neumann entropy it constrains is far costlier to obtain. We present a six-term linear closed form for the entropy, built on bitstring-derivable physics scalars, and characterize its accuracy, portability, scaling behaviour, and calibration cost. The feature set is selected with guidance from a trained graph neural network: probing the network localizes its entropy prediction to the two-point correlators on the bipartition boundary, and an exhaustive ground-truth search restricted to those boundary correlators isolates the form. It reaches $0.024$~nats mean absolute error in distribution: $6.4$ times the network's error, but in a form a human can read and apply without retraining. Fit once and applied unchanged, it has lower error than the base network on five of six out-of-distribution pools and ties the sixth. An independent density-matrix renormalization-group study to one hundred atoms – five times the reach of exact diagonalization – settles the size-extrapolation question: coefficients frozen at small size fail at scale, but the failure is structured. Refit per size the form holds to $25$–$50$~mnat (cross-validated); two of its six slopes follow clean inverse-size laws, one a downward curving growth, and the others are trendless; the fitted laws deploy the form label-free at roughly $40$–$80$~mnat. The result fixes a label-budget rule: at large sizes, a few dozen labels recalibrate the closed form to match a fine-tuned in-distribution ensemble on the same features, while nonlinear ML models pull ahead only given large labelled datasets.

Higher-order topological insulators (HOTIs) have revolutionized our understanding of wave localization, extending the bulk-boundary correspondence to lower-dimensional hinges and corners. Thus far, the realization of mechanical HOTIs has relied exclusively on discretely engineered metamaterials or periodic phononic lattices. Here, we report a fundamental paradigm shift by demonstrating that continuous, homogeneous hyperelastic manifolds undergoing finite multiaxial deformations naturally harbor intrinsic higher-order topological phases. By extending the generalized Stroh-Lie impedance formalism into a fully coupled 3D finite-strain framework, we map the highly nonlinear orthotropic geometric frustration onto a four-band effective Dirac Hamiltonian spanned by Clifford $\Gamma$-matrices. We reveal that macroscopic orthogonal stretches act precisely as competing Dirac mass terms, driving the continuous spatial transitions of topological domain walls and triggering a breakdown of $C_{4v}$ spatial symmetry. Remarkably, we analytically prove that beyond classical 2D surface wrinkling (1st-order topology), concurrent multiaxial extreme compression unconditionally triggers the emergence of 1D hinge states (2nd-order) and completely localized 0D zero-energy corner states (3rd-order). We further extend this static bifurcation framework into the elastodynamic regime, proving the existence of mid-gap localized vibrational modes. The theoretically derived topological phase diagram, nested Wilson loops, and fractional corner charges are comprehensively verified. Finally, we propose a concrete experimental realization using electro-active dielectric elastomers, enabling the dynamic programming of 0D topological singularities.

We construct an explicit Lewis-Ermakov-type dynamical invariant for time-dependent two-level systems by exploiting their algebraic correspondence with the time-dependent harmonic oscillator through the $su(2)$ and $su(1,1)$ algebras, which share the common complexification $sl(2,\mathbb{C})$. This invariant provides a closed-form evolution operator and an exact propagator for arbitrary time-dependent coupling $a(t)$ and detuning $b(t)$. We illustrate the method on representative scenarios, including Landau-Zener transitions, non-Hermitian dissipative processes, and adiabatic rapid passage, and we obtain inverse-engineered shortcuts to adiabaticity with fully explicit control fields when the auxiliary scaling function is taken to be real.

We propose a direct, model-independent search for irreversible coherence loss during a reversible measurement. A mesoscopic apparatus measures a microscopic two-state system by storing which-state information and is then returned to its pre-measurement state; unitary quantum mechanics predicts recovered coherence, while collapse leaves residual endpoint loss. In one representative device, a coherently controlled molecular force source displaces a charged nanoparticle that serves as the apparatus; during one $4.19\,{\rm ms}$ measurement-and-reversal cycle, its two states become almost fully distinguishable. Such a device would give a first direct bound in a measurement-and-reversal setting. Repeating this cycle $5\times10^5$ times would bound the corresponding irreversible coherence-loss rate beyond ordinary decoherence at $5.51\,{\rm s}^{-1}$. Continuous Spontaneous Localization is included as a secondary benchmark on the same apparatus history.

The Carruthers Geocorona Observatory is NASA's first mission dedicated to investigating the fundamental nature of Earth's exosphere and its dynamic response to space weather. Its primary payload, the GeoCoronal Imager, consists of two co-aligned broadband photometric imagers that support simultaneous, common-volume sensing of ultraviolet emission at 121.6 nm (Lyman-alpha) by exospheric hydrogen atoms. Accurate exospheric parameter retrieval from these images requires accurate knowledge of the instrument's optical responsivity, which enables conversion of measured exospheric signal rates into the scientifically relevant quantity of emission radiance. The Carruthers mission achieves absolute sensitivity calibration through the acquisition of photometric measurements of stars and subsequent inversion of the observed fluxes to retrieve the wavelength-dependent responsivity across the passband for each imaging configuration. The retrieval algorithm performance is enhanced by its incorporation of an objective, algorithm-driven ranking criterion to systematically select target calibration stars from a stellar spectral library. The end-to-end workflow, from the stellar selection criterion to the passband inversion, is validated using synthetically generated stellar measurements. These validation tests establish that our optical responsivity retrieval approach has high recovery fidelity, achieving error rates of <5% at the Lyman-alpha wavelength for all primary science imaging modes.

Time-varying electromagnetic media exhibit rich nonstationary wave phenomena, but the symmetry governing reversal of arbitrary temporal modulation sequences has remained unclear. We show that, in lossless, spatially homogeneous media with identical initial and final states, the scattering matrices of ordered and reversed sequences are related by inverse–conjugation, independent of the number of stages. This yields a classification of temporal reciprocity in bi-isotropic media: isotropic and chiral media are channel-preserving, whereas Tellegen media are channel-exchanging despite Lorentz nonreciprocity. Deterministic time rewinding follows directly. Our results provide a framework for predicting and designing temporal scattering responses in photonic media.

We present a numerical implementation of dissipative relativistic magnetohydrodynamics based on Bemfica-Disconzi-Noronha-Kovtun (BDNK) theory, in which first-order corrections render the equations causal without introducing additional dynamical variables. We show how these corrections can be incorporated in a finite-volume scheme describing the coupled dissipation of energy, momentum, and magnetic field, with the latter treated as a one-form charge. While a minimal set of BDNK terms can convert diffusive equations into telegrapher-type equations, we find that in the ultra-relativistic limit an additional correction is required for the system to behave in a stable and causal manner. With this set of equations, we develop an efficient method for primitive-variable recovery and validate the implementation through an analytical benchmark and various two-dimensional simulations.

We show that an oscillon can be understood as a localized discrete resonant (non-normalizable) mode. Specifically, oscillon in the vacuum arises from the threshold mode, which because of nonlinearity gets localized. Following this idea, we find wobblerons - nonlinear excitations of kinks, that is, oscillons-kink bound state. Now, the oscillon can also originate in an antibound mode, i.e., a discrete, positive energy but non-normalizable mode.

It was recently found that the quasinormal modes (or Ruelle–Pollicott resonances) of the SYK model at infinite temperature form a Christmas tree shape, reminiscent of AdS black holes. We generalise this computation to finite temperature, allowing us to continuously connect the infinite temperature results to the low temperature regime dual to JT gravity. We contrast the movement of the quasinormal modes with a few examples: various AdS black holes, dynamical phase transitions, and the large $p$ SYK chain. We find that the relaxation rate increases monotonically with temperature only at strong coupling, corresponding to the gravitational regime. Byproducts of our investigations are new results on operator growth that may be of independent interest.

Skeletal editing enables precise structural modifications of molecules at late stages of a synthetic sequence, with applications in drug discovery and materials science. We recently demonstrated skeletal editing on the single-molecule scale. Voltage pulses applied by the tip of a scanning probe microscope to an oxygen-containing seven-membered heterocycle led to both oxygen deletion and ring-contraction rearrangement reactions. An open question is whether selective skeletal editing of a heterocyclic core can be achieved by an appropriate choice of the heteroatom. Here, we show that tip-induced reactions of an analogous sulfur-containing seven-membered ring results in sulfur deletion in virtually all cases. Our results demonstrate that the combination of tip-induced chemistry and heteroatom selection in the molecular design is a powerful strategy for single-molecule skeletal editing, with the potential to enable diverse structural transformations of heterocyclic frameworks.

We consider Bianchi IX dynamics of a Universe filled with a massless phantom field. Such an exotic matter source enables regimes impossible in vacuum or with a standard scalar field. In particular, two Kasner indices of BKL oscillations can be simultaneously negative, and the absolute value of a negative index can be large. We describe the consequences of these features and explain the nature of volume oscillations recently discovered in such a system by numerical methods.

We study the polar director structure in thin planar cells filled with bent-core liquid crystals in the ferroelectric smectic-A phase (SmAP$_\text{F}$). We analyze a continuum phenomenological model proposed in the physics literature and present rigorous proofs of the existence and uniqueness of the equilibrium solutions. We further investigate the qualitative properties of nontrivial solutions and examine the effects of a bias electric field, surface anchoring, and cell thickness on the polar director configuration. Our results are consistent with previous experimental and numerical simulations reported in the physics literature. In addition, our analysis reveals new parameter-dependent behaviors supported by our numerical simulations and extends results reported from previous literature.

Recent models of pulsar polarization predict that the position and ellipticity angles of the polarization vector can trace portions of a small or great circle on the Poincare sphere. A great circle can arise from a transition in dominance of orthogonal polarization modes, where the relative intensity of the modes changes with pulse longitude. A small circle may be caused by a rotation of the vector, where the phase difference between the modes changes with pulse longitude or wavelength. Observations of PSRs B1133+16 and B2016+28 are reanalyzed to search for these polarization features within their pulse profiles. The polarization angles observed in part of PSR B1133+16 are shown to follow a great circle on the Poincare sphere. The angles observed across the pulse of PSR B2016+28 follow an arc that resembles a portion of a great circle that has been altered by the pulsar's rotation. The observations are interpreted within the context of three different polarization models. All three models produce similar results for both pulsars and indicate that the observed geodesics are caused by mode transitions. The arc observed in PSR B2016+28 can also be interpreted as a vector rotation, provided the modes are elliptically polarized. The observations and accompanying analysis show that mode transitions are not restricted to the equatorial plane of the Poincare sphere and that arcs and partial circles may be more common than previously recognized.

We formulate four-dimensional finite-proper-time QED from a gauge-covariant source functional with a retained spectral lower boundary $s_0=\Lambda^{-2}>0$ and finite physical duration $\mathcal T$. Physically, the basic $s_0$ operation is a truncated Laplace transform: the arbitrarily small-proper-time part of the spectral integral, which carries the ultraviolet singularity, is excluded without imposing a cutoff on momentum space, while finite $\mathcal T$ expresses that a physical observable is prepared and detected during a finite experiment rather than by an infinite-time measurement. The retained $s_0$ belongs to the defining source construction and controls the correlated finite corrections after mass and charge matching. The charged interval and circle sectors are generated by the same covariant heat semigroup, while the photon covariance carries the same retained lower bound at the source level. At every fixed perturbative order, gauge covariance gives the Ward–Takahashi hierarchy and an assembled graph-level quadratic-form bound establishes ultraviolet finiteness. The physical electron, positron, and transverse-photon poles retain their canonical residues, and the fixed-order cutting relation uses the matched daughter domains on the stated regulated continuation domain. Mass and charge renormalisation are replaced by finite matching, while independent observables retain correlated $s_0$ dependence. At fixed $s_0>0$, the Euclidean source has a nonperturbative definition. Its physical standing-wave response is obtained order by order from the continued FPT amplitudes with the same finite-duration channel data. An exact all-coupling physical boundary realisation requires additional channel-norm control. The limits $s_0\to0^+$, $\mathcal T\to\infty$, and exact all-coupling free-space exhaustion are separate limiting questions and are not used to define the theory.

Noise can significantly influence thermoacoustic dynamics, yet the role of noise color in coupled thermoacoustic oscillators remains largely unexplored. Here, we examine the influence of colored noise on the dynamics of coupled thermoacoustic systems. The system consists of two coupled Rijke tube oscillators with time-delay and dissipative coupling. Stochastic forcing is modeled as an additive Ornstein-Uhlenbeck (OU) process, such that white and colored noise contain equal power within a band around the system's natural frequency. We find that noise influences the system most prominently near the transition between limit-cycle oscillations (LCO) and amplitude death (AD) states, where increasing noise amplitude smoothen the transition and reduces the extent of the AD regions. Our analysis reveals the emergence of coherence resonance near instability threshold under both white and colored noise. The peak coherence factor varies with the noise color, with the largest peak coherence observed for colored noise whose correlation time is much shorter than the acoustic time scale. White noise and the shortest-correlated OU noise exert the strongest influence on both the pressure amplitude response and the coherence resonance. Overall, our results show that, under both coupling mechanisms, colored noise induces qualitatively similar trends in the system response, governed by its amplitude and correlation time.

The precise realization of molecular electronic devices requires a comprehensive understanding of charge transport mechanisms and the specific interplay between electronic and nuclear degrees of freedom. While average current measurements (I-V characteristics) and conventional Inelastic Electron Tunneling Spectroscopy (IETS) offer valuable insights, they are fundamentally limited by temperature-dependent line-width broadening. This study presents a high-resolution spectroscopic methodology utilizing suspended-wire molecular junctions (SWMJs) based on self-assembled monolayers (SAMs) of 1-decanethiol (C10) and 1,1',4',1”-terphenyl-4-thiol (TPT). By systematically probing the voltage-dependent current noise (${\Delta}I$), we demonstrate that electronic noise spectroscopy circumvents thermal degradation by probing transition rates between vibrational manifolds rather than simple additions of conductance channels, which enables sub-thermal feature mapping. Leveraging a fast-convolution-based Landauer-Büttiker transport model fitted to experimental data, we map complex vibrational manifolds, including high-energy overtones. This allows for the direct extraction of crucial nanoscale molecular parameters, including mode energies, anharmonicities ($x_e$), dissociation energies ($D_e$), and local environment reorganization energies ($E_r$). These parameter-dense noise signatures act as a unique molecular fingerprint, establishing noise spectroscopy as a highly sensitive platform for chemical sensing and discrimination in advanced quantum devices.

The vorticity distribution associated to canonical Gibbs point vortex ensembles under mean field scaling has Gaussian fluctuations for small negative inverse temperature. The perturbative argument is based on a cluster expansion of the partition function.

Gradient-index (GRIN) media with freeform refractive-index distributions require ray tracing approaches capable of accurately modeling light propagation in complex three-dimensional environments. In this work, we develop a fully three-dimensional ray tracing method for arbitrary freeform GRIN media based on the local application of a set of vectorial relations of geometrical optics. At its core, the method constructs ray trajectories through successive local refractions using the vectorial form of Snell's law, while employing constant optical path length steps that naturally adapt the geometric step size within the medium. Unlike conventional continuous GRIN propagation approaches, the proposed method naturally handles refraction at the boundary of the medium as well as total internal reflection, enabling the treatment of both continuous media and discretized refractive-index distributions composed of isoindicial surfaces. The method is simple to implement and computationally efficient. Validation against an analytical solution demonstrates high accuracy, while applications to freeform GRIN configurations illustrate its capability not only for propagation analysis, but also for the direct design of complex tailor-made optical media.

Nanosized hafnia-zirconia (HfxZr1-xO2) in the form of thin films, multilayers, and nanoparticles is one of the most promising CMOS-compatible ferroelectric materials for advanced electronic memories and logic devices. Using the Landau-Ginzburg-Devonshire free energy functional with trilinear and biquadratic couplings of polar, nonpolar, and antipolar order parameters, we analyze the pyroelectric and electrocaloric properties in an ensemble of spherical core-shell HfxZr1-xO2 nanoparticles. To test the theoretical model, we experimentally measured the temperature dependence of the electric charge accumulated in pressed powders consisting of oxygen-deficient core-shell Hf0.5Zr0.5O2 nanoparticles with an average size of 7 nm. The observed temperature-dependent behavior of the accumulated charge and its derivative are in qualitative agreement with the calculated polarization and pyroelectric coefficient for the ensemble of densely packed spherical core-shell HfxZr1-xO2 nanoparticles; this suggests that the theoretical model captures the physical mechanisms responsible for the experimentally observed charge accumulation. The combined theoretical and experimental results provide a physical foundation for the future development of CMOS-compatible HfxZr1-xO2 nanoparticles for pyroelectric and electrocaloric applications.

We give the complete solution of the one-parameter deformation of the six-vertex model on dynamical lattice introduced in [1] and dubbed gravitational seven-vertex model. The statistical model in question is mapped to a gas of self- and mutually avoiding loops on dynamical triangulations, with a temperature coupling controlling the volume not occupied by loops. The phase diagram is characterised by massive, dilute and dense critical phases, similarly to the gravitational O(n) loop model. There is however an important difference – in our model the weights of the loops are not topological but depend on the form of the loop and on the curvature defects of the lattice via lattice spin connection. The seven-vertex model on dynamical lattice is nevertheless exactly solvable after being reformulated as a large-N matrix model, which we will refer to as 7vMM, and the solution in the scaling limit was found in [1]. Here we derive the full solution in terms of Jacobi theta functions and present the (non-algebraic) spectral curve of 7vMM in a parametric form. We obtain the phase diagram in the space of the two coupling constants – the cosmological constant and the temperature – and identify the critical phases along the boundary of the physical domain. We derive the scaling solution of [1] as the asymptotic of the full solution in the vicinity of the tricritical point separating the phases of dense and massive loops.

Optical imaging in complex samples such as biological tissues is fundamentally challenging due to random light scattering that degrades resolution and contrast. When imaging realistic targets that contain natural dynamics such as flowing blood, the temporal variability introduces an additional obstacle, as the leading computational scattering-compensation methods require the target to remain stationary during a multi-frame acquisition process. Here we show that instead of struggling to perform rapid acquisitions, the target dynamics themselves can serve as an intrinsic information source for scattering compensation, replacing multiple controlled illuminations. The mathematical equivalence between target dynamics and conventional varying illumination patterns allows to demonstrate this approach in coherent holographic imaging and incoherent fluorescence microscopy using established matrix and model-based scattering-compensation techniques. Our general framework enables reconstruction of dynamic scenes with a number of acquisitions equal to the number of reconstructed frames, without the use of any spatial light modulators or illumination control.

Understanding dynamical properties of quantum many-body systems remains a central challenge because excitations generally require information beyond a ground-state wave function. Here we introduce a neural polaron ansatz that represents quasiparticle excitations by neural many-body operators acting on a correlated ground state. Instead of learning an independent excited-state wave function, the method parameterizes a local dressing operator through a compact neural head defined on the feature map of a pretrained ground-state network. This operator-based construction builds in translation symmetry, momentum resolution, and quasiparticle locality, while separating ground-state correlations from excitation-specific dressing. We benchmark the method on the square-lattice $J_1$-$J_2$ Heisenberg model, where it accurately reproduces magnon dispersions and spectral weights over a broad range of frustration. In particular, it captures nontrivial many-body features including the $(\pi,0)$ anomaly and its progressive softening with increasing $J_2/J_1$. These results establish neural operators as a physically transparent route for extending neural quantum states from ground-state properties to dynamical response.

In this work, we study a $\mathbb{Z}_3$-symmetric extension of the Standard Model with three hyperchargeless $SU(2)_L$ fermion triplets responsible for neutrino mass generation $\textit{via}$ the Type-III seesaw mechanism together with a complex scalar singlet $\chi$ whose vacuum expectation value spontaneously breaks the $\mathbb{Z}_3$ symmetry. Radiative corrections induced by the Yukawa interactions between the $SU(2)_L$ fermion triplets and the complex scalar singlet $\chi$ generate a Coleman-Weinberg vacuum bias that lifts the degeneracy among the $\mathbb{Z}_3$ vacua, leading to the annihilation of unstable domain-walls. Consequently, the degeneracy among the $\mathbb{Z}_3$ vacua is lifted radiatively through the Coleman-Weinberg effective potential, generating a dynamical bias term that triggers the annihilation of unstable domain walls. We perform a numerical analysis consistent with current neutrino oscillation data and identify viable regions of parameter space accommodating the observed neutrino masses and leptonic mixing parameters. The observed baryon asymmetry of the Universe is generated through thermal leptogenesis $\textit{via}$ the out-of-equilibrium decay of the lightest fermion triplet for masses around $\mathcal{O}(10^{9})\,\mathrm{GeV}$, consistent with the Type-III seesaw framework. Depending on the choice of model parameters, the predicted gravitational-wave spectrum can fall within the sensitivity reach of future space-based and ground-based gravitational-wave detectors. Our framework therefore establishes a correlation between neutrino mass generation, leptogenesis, radiative domain-wall instability, and gravitational-wave phenomenology.

We extend the nullifier-based certification framework of catability from parity-defined coherent cat states to parity-less Kerr cat states. Since Kerr cats cannot be distinguished by parity alone, we replace the parity component of the original nullifier with a displaced-parity operator that captures their characteristic interference structure. The resulting generalized catability remains directly observable and can be evaluated from a finite number of photon-number measurements without full state tomography. Numerical benchmarks show that the method faithfully certifies Kerr-cat features, produces accurate approximations of ideal Kerr cat states, and is more resilient to optical loss than fidelity-based certification. We also identify an exact anti-linear nullifier based on complex conjugation, providing an ideal algebraic description of Kerr cat states. Our results broaden the scope of catability and provide an experimentally practical approach to the characterization of Kerr cat states.

A phase-field-based reaction-diffusion corrosion model is developed to predict microbially influenced corrosion (MIC) in metal alloys, with a focus on anaerobic conditions and sulfate-reducing bacteria (SRB). The formulation couples microbial sulfate reduction, sulfate transport, electrochemical kinetics, material dissolution, and mechanical effects. Microbial activity is modelled using a Monod-type expression for sulfate consumption, whereas the mechano-chemical coupling is incorporated through an enhanced mobility relationship that captures the influence of mechanical fields on corrosion kinetics. The model is calibrated against experiments and shows strong agreement in predicting pitting kinetics under SRB activity. Sensitivity analyses quantify the competing roles of microbial kinetics, transport, and thermodynamic driving forces in governing corrosion behaviour. The capability of the formulation to capture both MIC-induced pitting and stress-assisted corrosion across multiple length scales is demonstrated through case studies that include microstructure-sensitive simulations and structural-scale coupling with a cathodic protection (CP) model. Results show that finer grain sizes reduce pitting severity but promote faster defect propagation under mechanical loading. At the structural scale, coupling with the CP model enables predictions of defect growth under varying electrochemical conditions and over engineering-relevant length scales, as exemplified with the analysis of an offshore wind turbine monopile. CP delays pitting and suppresses crack propagation, although its effectiveness diminishes as sacrificial anodes degrade. The framework provides a predictive and computationally efficient tool for assessing MIC-induced damage over extended times, with potential applications in the integrity and life assessment of metallic structures operating in aggressive microbial environments.

A methodology for passive mitigation of vortex-induced vibrations (VIVs) in subsea dynamic power cables is developed using optimized, strongly nonlinear bi-stable mass-spring-damper attachments, termed bi-stable nonlinear energy sinks (BNESs), within the open-source MoorDyn library. A fully three-dimensional time-domain framework captures cable dynamics, Morison-type hydrodynamic forcing, nonlinear vibration-mitigation mechanisms, and spatially and temporally varying currents. The BNESs are consistently integrated into the cable model, allowing treatment of highly non-stationary VIVs. A data-driven optimization study samples the BNES design space across multiple current profiles and shows that properly tuned configurations substantially reduce peak-to-peak cable vibration amplitudes. The BNESs also produce significant and robust reductions in cumulative VIV-induced energy intake from the surrounding flow and in drag energy. To the authors' knowledge, passive nonlinear attachments are shown for the first time to reduce both energy intake and drag amplification in a subsea cable, with potential benefits for short- and long-term fatigue. Time-frequency wavelet analysis reveals targeted nonlinear energy transfers and scattering from dominant high-amplitude, low-frequency cable modes to lower-amplitude, higher-frequency modes. This modal redistribution promotes rapid dissipation through hydrodynamic damping and internal structural losses, explaining the simultaneous reductions in vibration amplitude and cumulative energy intake. The results demonstrate that BNESs can provide effective and robust VIV mitigation for subsea power cables under realistic unsteady operating conditions and motivate future studies involving combined current-wave loading, platform-induced motion, and fatigue-life assessment.

We report the fabrication and characterization of an integrated quantum photonic device consisting of an electrically driven whispering-gallery-mode micropillar laser evanescently coupled to a ridge waveguide, both incorporating InGaAs quantum dots (QDs). The lasing characteristics of microlasers are systematically investigated as a function of the pillar-waveguide gap distance. Coherent emission from the whispering-gallery-mode microlaser coupled into the waveguide enables on-chip optical excitation of QDs embedded in an electrically contacted micropillar at the end of the waveguide. Under continuous-wave on-chip excitation, we observe single-photon emission with $g^{(2)}(0) = (3.49 \pm 0.01) \%$ for a QD integrated in the outcoupling micropillar which can be spectrally tuned-by the quantum confined Stark effect. These results constitute an important step toward low-footprint, deterministic, and scalable single-photon sources for QD-based integrated quantum photonic circuits.

Building on resistive communication, this paper presents a physics-based design of an on-chip neural network with multi-level memristive synapses supporting a dense spectrum of conductance states. Derived from ionic transport physics, we develop a state-variable model and quantify storable sub-levels under thermal noise, drift, and quantized conductance. We assemble these devices into a 1T1R crossbar fabric, derive the linear algebra of analog vector-matrix multiplication (VMM) under wire resistance, and design a differential synapse for signed weights. A multilayer pipeline executes inference, backpropagation, and weight updates physically in the analog domain. We derive the in-situ outer-product learning rule, its discretization onto the conductance lattice, and the resulting quantization noise. We provide energy, area, capacity, and inter-tile models, showing this substrate is ideally suited for large language models (LLMs). Our design eliminates weight movement, surpassing binary ReRAM and traditional CMOS. We detail the material stack (HfO_2-based), the FEOL/BEOL CMOS foundry-integration flow, a self-contained SPICE model, the complete memristive-FPGA neuromorphic system, and an in-memory self-attention engine with current-mode translinear softmax. Finally, a ternary BitNet datapath shows projected per-token efficiency orders of magnitude better than advanced CPUs or GPUs. The result is a self-contained hardware-native blueprint for a high-density, analog, in-memory neural processor.

Luminous red novae are astrophysical transients associated with unstable mass transfer in interacting binaries and are commonly interpreted as outcomes of common-envelope evolution, possibly ending in merger. These interactions can liberate large amounts of gas, part of which may later condense into dust. We study the luminous red novae AT2021biy and AT2021blu to constrain their binary progenitors, estimate the mass ejected during the outbursts, and infer the dust mass from the infrared evolution of their remnants. We computed two grids of binary stellar-evolution tracks with the MESA binary module, constrained by pre-outburst photometry. We applied mass-transfer instability criteria to select progenitors able to merge on timescales compatible with archival observations. From these models, we estimated the gas mass lost during the mass-transfer phase and lower and upper bounds on the envelope mass that could be ejected during common-envelope evolution using the available orbital energy. We compared these ejecta-mass estimates with values inferred from light-curve models. Finally, we modeled mid-infrared NEOWISE data to derive dust masses up to ~3 years post-outburst, providing an additional lower limit on the total ejecta mass. We constrained the donor masses to Md = 18-23 Msun for AT2021biy and Md = 14 +/- 0.5 Msun for AT2021blu. Lower limits on the ejected envelope mass are 0.03-2.98 Msun for AT2021biy and 0.02-0.1 Msun for AT2021blu. Comparison with light-curve models favors intermediate mass ratios, q = 3-10 for AT2021biy and q = 5-15 for AT2021blu. The inferred dust masses are 1-5 orders of magnitude below the estimated ejected envelope masses, implying that only a small fraction of the gas condenses into dust. Their evolution is consistent with shock interaction and suggests pre-existing circumstellar material, in line with the pre-outburst mass loss predicted by our MESA models.