2026-09-09 | | Total: 301
In the Path Integral Monte Carlo (PIMC) simulation of fermions in a harmonic trap, with and without pairwise harmonic interactions, the partition functions for any discrete number of imaginary time slices (or beads) and for any choice of the short-time propagator can be analytically obtained from the contracted determinant form of the propagator. This work shows that the resulting recursion relation can be reformulated in the $λ$-ring language, yielding a closed-form finite-bead partition function in two dimensions in terms of permutation statistics. This closed-form partition function can be evaluated by a special algorithm in $O(n^2)$ time, providing an exact and numerically stable scheme for reproducing the energies of the original (undoable) fermion PIMC simulation for $n=10^4$ or more fermions. This result provides a concrete framework in which the numerical instability of the sign problem is completely bypassed, and serves as a counterexample to the prevailing view that all truly fermionic PIMC sign problems are NP-hard.
We present a quotient-space Bayesian formulation for DESI-era anomaly interpretation. Given a pattern-labeled catalog with map \(i\mapsto \Act(i)\), the induced posterior \(p(\Act\mid D)\), sector inclusion probabilities \(P_α\), co-activation probabilities \(P_{αβ}\), and grouped Bayes factors \(B_α(D)\) are exact summaries over predeclared physical activation events. Pairwise comparisons such as \(\lcdm\) versus \(\wacdm\) remain ordinary Bayes-factor tests between specified families; the quotient construction addresses the coarser question of which physical sector carries posterior support when different sectors are represented by unequal numbers of catalog elements. We derive a sector-resolved DESI-CMB-SN likelihood specification for late-time background, early-time ruler, supernova calibration, perturbation, and gravitational-wave propagation sectors. The construction includes an Alcock-Paczynski/isotropic-scale BAO decomposition, a pure-ruler projection, analytic marginalization of low-rank supernova calibration modes, Fisher-normalized sector priors, inactive-sector leakage tests, log-evidence uncertainty propagation, and prior/sector-partition diagnostics. The result is a quantitative procedure for reporting model-comparison support at the level of physically interpretable sectors.
In-situ measurements in space plasmas indicate that the velocity distributions of charged particles are not in thermal equilibrium, deviating from a standard Maxwellian mainly due to anisotropies and suprathermal populations which enhance high-energy tails. Although the Standard \k{appa}-Distribution (SKD) is well established in modeling these non-equilibrium distributions, its application is often viewed controversially due to certain unphysical implications, in particular divergent velocity moments. To address these issues, the Regularized \k{appa}-Distribution (RKD) was introduced. Such advanced, in general anisotropic RKDs are invoked here for the first time to investigate oblique firehose instabilities, including those induced by the temperature anisotropy of electrons and protons, the dominant species in space plasmas. In weakly collisional plasmas, both of these instabilities are expected to play significant roles in the self-regulation of the macroscopic properties of space plasmas (e.g. the expanding solar wind) reported by observations. Is not yet possible to resort to a general dispersion tensor related to RKD plasmas (whose derivation is still a challenge), instead the Arbitrary Linear Plasma Solver (ALPS) is exploited here. The unstable solutions obtained for the already established Maxwellian & SKDs successfully validate the capability of ALPS. For RKDs, the instabilities confirm the stimulating effect of suprathermal populations, with lower \k{appa} values generally enhancing the firehose growth rates. For the oblique electron firehose instability, RKDs substantially modify the competition between periodic & aperiodic branches only at intermediate angles & can sustain significantly increased growth rates of especially aperiodic modes at highly oblique propagation beyond both the Maxwellian unstable regime & the one that is accessible with SKDs.
Conventional semi-local density functional approximations for the exchange-correlation energy systematically underestimate fundamental band gaps. Advanced meta-generalized gradient approximations (meta-GGAs) can partially overcome this limitation through their orbital dependence via the kinetic energy density. In principle, the flexible meta-GGA form should allow for balanced performance across diverse properties. However, atomization energies, reaction barrier heights, and lattice constants, as well as band gaps, exhibit different sensitivities to the large parameter space of meta-GGA functionals. In this work, we explore the space of appropriate norms that guide the r$^2$SCAN-type functional construction through the interpolation function connecting two fundamental paradigms: the one-electron limit and the uniform electron gas limit. We introduce local modifications to a strongly-smoothed interpolation function and investigate their impact on functional performance. To further probe the limits of the meta-GGA framework, we include spin-unpolarized bonded systems among the set of norms. The resulting bonded-norm bn-r$^2$SCAN meta-GGA achieves balance: band-gap predictions approaching the accuracy of the ultra-nonlocal LAK meta-GGA and lattice constants of nearly r$^2$SCAN accuracy.
In aquatic environments, gradients in flow velocity and turbulence define a constantly shifting landscape that sets the physical constraints on propulsion and energy expenditure in fish. Successfully navigating these complex flows depends on the ability to sense and respond to subtle hydrodynamic cues in an energy-efficient manner. Revealing these mechanisms is central to understanding how fish identify and exploit favorable flow conditions. To this end, we designed a hydrodynamic maze with spatially varying mean flow and turbulence, creating a controlled heterogeneous flow landscape for understanding fish movement strategies. We found that fish escape energetically unfavorable regions by yawing their bodies relative to the flow. These yawed orientations appear to facilitate lateral migration through lift forces and may improve sensitivity to local flow variations. Notably, fish maintain body orientations near, but below, predicted stall conditions to exploit lift without a significant increase in drag. These findings provide insight into how fish navigate heterogeneous hydrodynamic environments, with potential relevance for both natural and engineered flow systems.
Stimulated Brillouin Scattering (SBS) is a nonlinear interaction between optical and acoustic waves in solids. First regarded as a parasitic process in optical fibers, it has gathered significant interest in microwave photonic applications such as optical signal processing and narrow linewidth lasers. While many Brillouin system demonstrations have been done on photonic chips, they struggle to provide simultaneous high Brillouin gain and narrow linewidth in a scalable platform, capable of integration with established photonic integrated circuits. Here, we present a novel integrated InGaP-on-SiO$_2$ platform, where a single crystalline InGaP waveguide layer offers superior material properties and strong nonlinearities to support SBS near 1550 nm wavelength. We demonstrate backward SBS with a measured high Brillouin gain coefficient of $588\, \mathrm{W}^{-1}\,\mathrm{m}^{-1}$ and a record narrow linewidth of 5.2 MHz at 9.346 GHz frequency shift. This work establishes a wafer-scale complementary metal-oxide semiconductor compatible fabrication process, paving a scalable path for high gain and narrow linewidth Brillouin photonics with applications such as precision signal processing.
We report simultaneous sub-Doppler cooling and optical dipole trapping of $^{39}$K-Cs and $^{41}$K-Cs mixtures. Both mixtures are cooled to temperatures $\sim10$ $\mathrmμ$K, achieving performance comparable to that obtained with each species individually. To the best of our knowledge, this constitutes the first realization of a laser-cooled and optically trapped $^{41}$K-Cs mixture. For the $^{39}$K-Cs mixture, we additionally implement parallel spin-resolved Feshbach spectroscopy enabled by Stern-Gerlach separation. For potassium isotopes, we implement sub-Doppler cooling with $D_1$-line gray molasses, while for cesium the entire cooling sequence is implemented using $D_2$ transitions. We spin-polarize the atoms and confine Bose-Bose mixtures in a 1064 nm optical dipole trap. Using the resulting $^{39}$K-Cs samples, we observe 14 heteronuclear Feshbach loss features. Five agree with resonances reported previously, while nine have, to the best of our knowledge, not been observed experimentally before, including $p$-wave features and resonances in additional spin channels. The shorter lifetime of the $^{41}$K-Cs mixture in the optical dipole trap currently hinders systematic Feshbach spectroscopy, which we therefore do not pursue in this work. To characterize this limitation quantitatively, we study the decay dynamics of $^{39}$K-Cs and $^{41}$K-Cs mixtures under comparable temperature and density conditions, revealing a substantially stronger nonexponential loss in the $^{41}$K-Cs mixture. The demonstrated preparation of ultracold $^{41}$K-Cs mixtures provides a starting point for Feshbach-resonance and photoassociation spectroscopy of this previously unexplored isotopologue. Such measurements are essential for identifying suitable pathways for magnetoassociation and coherent optical transfer, and ultimately for the production of ultracold ground-state $^{41}$KCs molecules.
Graph-based surrogate models offer a promising route to accelerate computational fluid dynamics (CFD) simulations on unstructured meshes. However, their development is limited by the scarcity of benchmark datasets spanning multiple flow regimes and standardized protocols for long-horizon autoregressive prediction. We introduce ONECYL (ONE CYLinder), a new benchmark for unsteady flow past a circular cylinder across laminar, transitional, and high-Reynolds-number regimes. The benchmark comprises 450 high-fidelity Variational Multiscale finite-element simulations (270,000 flow snapshots) with randomized cylinder geometries, providing time-resolved velocity and pressure fields together with mesh connectivity, geometric descriptors, Reynolds numbers, and integrated aerodynamic quantities. Beyond the dataset, ONECYL establishes a unified evaluation framework combining full-field rollout errors, virtual probes, and drag and lift predictions to assess numerical accuracy and physical fidelity. To accompany the benchmark, we develop a Graph Transformer as a reference baseline predicting velocity and pressure fields autoregressively on unstructured meshes. Using ONECYL, we investigate geometric representations and physics-based regularization across the three Reynolds-number regimes. The results show that explicitly encoding the cylinder geometry through a level-set representation consistently improves long-horizon prediction accuracy and generalization to unseen geometries, while divergence-based regularization becomes increasingly beneficial as flow complexity increases. The ONECYL benchmark and its Graph Transformer baseline provide a reproducible framework for evaluating graph-based surrogate models and establish a foundation for future research on long-horizon prediction of unsteady bluff-body flows.
Quasi-static sensing, i.e., measurement of a sensor signal change induced by analyte presence with respect to a baseline signal observed in a reference gas beforehand, is the state of the art for qualitative and quantitative analyte identification in chemiresistive sensing. However, this approach is prone to baseline drift and, especially for inexpensive miniaturized point-of-care sensors, repeated calibration with reference gas is not feasible. To access reliable baseline information, active switching between analyte and reference gas would be required, causing a limitation for miniaturized implementation. Dynamic excitation of sensors offers an alternative: Via external stimuli the sensor is reversibly driven out of equilibrium in a controlled way, and the resulting transient sensor response, which is dependent on analytes present, is recorded and interpreted. Correlation of excitation and response signals can reduce effects of baseline drift and furthermore, dynamic features in sensor responses may contain valuable information. Especially light activation is a powerful and potentially miniaturizable approach to agitate sorption processes and therefore induce dynamics. In this work we report photothermal heating of hybrid, chemiresistive gold nanoparticle composites via LED excitation matching the plasmon resonance, to dynamically shift the analyte sorption equilibrium and induce dynamic sensor responses in a highly controlled way. We demonstrate that these sensor responses enable the detection of volatile organic compounds (VOC) via interpretation of differential signal components and that this approach improves the chemiresistors' baseline stability.
Model Hamiltonians (e.g., the Hubbard or Heisenberg models) provide a simple yet physically meaningful description of complex electronic phenomena through a small number of parameters, such as electron hopping integrals or magnetic exchange couplings. Their construction typically exploits the local nature of electron correlation effects and therefore relies on localized orbitals. Although many localization schemes are available, the resulting localized orbitals depend not only on the chosen localization functional but also on the orbital space to which the localization procedure is applied. These choices are often treated as technical details and are rarely discussed explicitly. Here, we demonstrate that different, chemically reasonable choices of localized orbitals can lead to substantially different model Hamiltonian parameters, even when they produce nearly identical electronic energies. Using Density Matrix Downfolding, we systematically investigate how the localization functional and, more importantly, the choice of model orbital space affect effective Hamiltonians derived from ab initio calculations for representative $π$-conjugated systems and transition metal complexes.
We present a theoretical study of polarization dynamics in the elastic scattering of muons and antimuons by high-magnetic-moment atomic nuclei, using ${}_{83}^{209}$Bi as a benchmark. The analysis tracks the evolution of (anti)muon polarization and its transfer to the nucleus, emphasizing nuclear-structure corrections. These corrections significantly alter both the differential cross sections and the polarization observables. Our results demonstrate that measuring scattered muon polarization serves as a sensitive probe for nuclear charge radii and magnetic dipole moments.
Climate change poses a major scientific, environmental, and educational challenge, requiring approaches that connect physical concepts, empirical evidence, and socially relevant environmental issues. This study presents and analyzes a low-cost experimental investigation developed within the "Physics and the Environment" course of a Physics teacher training program at the State University of the Tocantina Region of Maranhão (UEMASUL), Brazil. The experiment involved constructing a simplified thermal simulator using accessible materials and comparing a confined system with an open control system, both exposed simultaneously to solar radiation. The results showed that the average temperature inside the simulator was $45.91^{\circ}\mathrm{C}$, compared to $38.58^{\circ}\mathrm{C}$ in the external control system, corresponding to a mean temperature difference of $7.33^{\circ}\mathrm{C}$ (among other analyses). These results demonstrate that the experimental setup was capable of producing and measuring a significant thermal difference between the two systems. However, the observed temperature rise should not be interpreted as a direct measure of the atmospheric greenhouse effect, given that the experimental apparatus involves both radiative and non-radiative heat transfer mechanisms. The study emphasizes the educational potential of low-cost experimentation as a strategy to integrate electromagnetic radiation, energy balance, temperature, heat transfer, and climate science into the initial training of physics teachers. It also highlights the importance of explicitly discussing the limitations of experimental models and distinguishing local experimental evidence from conclusions regarding global climate processes.
Significance: Most melanoma deaths arise from metastatic spread, yet current staging imperfectly identifies which primary tumors will progress. Melanin structure is altered during malignant transformation, and pump-probe microscopy (PPM) measures melanin excited-state dynamics in standard biopsy sections, offering molecular contrast that is complementary to morphology-based assessment. Aim: To determine whether PPM-derived melanin excited-state dynamics in primary cutaneous melanoma can predict metastatic outcome. Approach: We imaged unstained sections from 71 primary cutaneous melanomas and 18 melanoma metastases. Transient absorption curves were fit to a biexponential excited-state model, and the fit parameters were used as features for learning classifiers to predict metastatic outcome. Results: Melanin dynamics separated primary from metastatic tissue: both excited-state lifetimes were shorter in primaries, with common-language effect sizes up to 0.26. The same parameters did not separate primary tumors by outcome (effect sizes 0.41-0.61), and univariate logistic regression found no significant association for either the excited-state lifetimes or the ESA/GSB ratio. The best classifier reached 74% patient-level accuracy (AUC 0.73), a plateau reached by several unrelated architectures. Conclusions: PPM robustly distinguishes primary from metastatic melanoma tissue, consistent with progressive melanin disaggregation. While it is clear that pump-probe features of melanin thus correlate with tissue differences, melanin dynamics in single sections of the primary tumor, confounded by heterogeneity, do not give robust predictions of metastatic outcome.
In introductory physics laboratories, a central instructional goal is to help students construct, evaluate, and revise mathematical representations from experimental data rather than apply given formulas. We present a guided data-driven activity in which an RC circuit is treated as a black box defined through observable input-output relations. Experimental and graphical procedures are provided, while standard RC circuit theory is initially withheld. Students first examine discharging data obtained under different nominal values of R and C. They rescale time using T=t/RC and compare the normalized voltage v=V/E across circuit conditions, bringing the discharging data closer to a common trajectory. Graphical linearization through a plot of lnv versus T provides a basis for proposing an empirical representation and allows RC to be interpreted as a time scale rather than being assigned that meaning a priori. When the representation proposed for discharging is tested with charging data, logarithmic linearity is lost despite the continued organization of the data by T. Replacing v with the distance from the steady state, 1-v, restores linearity, and the resulting linear relation provides a basis for proposing a revised representation for charging. The activity thus uses a familiar physical system to make representation construction, breakdown, and reconstruction experimentally explicit. It also emphasizes that an empirical representation must be accompanied by a specified domain of validity and that standard circuit theory can be introduced a posteriori to interpret the physical meanings that emerge from representational constraints.
Molecular diffusion in fluctuating amorphous and macromolecular media governs key transport processes across soft-matter physics, energy storage, and biological membranes. Extracting localized trapping states from single-particle tracking trajectories remains a fundamental challenge; because thermal structural breathing continuously reconfigures pore boundaries, conventional geometric algorithms suffer from severe systematic biases, erroneously merging distinct localized states during cyclic molecular returns. Here, we address this deadlock by shifting the paradigm from local geometric recurrence to a structure-informed Bayesian regularization. Leveraging discrete Morse theory, we extract the time-invariant topological skeleton of the fluctuating host matrix to construct robust, gas-specific physical priors that account for individual molecular dimensions. Trajectory steps are sequentially partitioned via a two-stage probabilistic refinement that dynamically adapts to the transport landscape. Benchmarked against a rigorous environment where synthetic particles explore the actual interconnected matrix graph, our approach eliminates systemic biases, restricting macroscopic trapping parameter deviations to just a few percent under optimal linear $O(N)$ computational scaling. Applied to hydrogen and methane transport within a type-I kerogen matrix, serving as a prototype for highly tortuous, flexible macromolecular networks, the method successfully decodes the hidden microscopic mechanisms of confined diffusion. To ensure immediate broad impact, the documented open-source code and data are made publicly available, offering an accessible strategy readily adaptable to a broad spectrum of tracking phenomena, from ion transport in battery polymers to protein trafficking within cellular environments.
We present the cluster mean field(cMF) and its extensions into perturbation theory(cPT) and coupled cluster(cCC) methods for fermionic systems. We discuss how the commutation rules for fermionic systems modifies the method as well as other computational details. We then perform calculations on the Hubbard model and a few select conjugated molecules to demonstrate the efficacy of the method. We will also use these examples to compare cluster methods to CAS-SCF.
Silicon sensors are the foundational detection medium for X-rays and charged particles. While their bulk dopant distribution determines device performance, it is conventionally assumed homogeneous because traditional profiling is destructive, spatially restricted, and insensitive at the relevant concentrations. Here we introduce a non-destructive 3D doping imaging technique that turns the readout electronics of a charge-integrating hybrid pixel detector into a massively parallelized capacitance-voltage profiler. With a few tens of micrometres of 3D resolution over wafer-scale areas at concentrations on the order of $10^{11}$ cm$^{-3}$, we image the bulk doping concentration of operational sensors. Macroscopically, we resolve depth-evolving concentric doping rings; microscopically, we uncover scattered doping anomalies that distort local electric fields. The rings modulate the depletion voltage, while the anomalies disrupt local charge collection, a previously overlooked cause of pixel yield and performance degradation. By bridging manufacturing signatures with microscopic defects, this approach provides a non-destructive framework for sensor characterization and yield optimization.
Accurate simulation of glacier surface mass balance is essential for predicting sea level rise and freshwater resources, but it is constrained by uncertainties in meteorological forcing and model parameters. Here, we deploy glacier data assimilation strategies to assess the value of observations for improving surface mass balance simulation, focusing on observation quantity, quality, and timing. We perform synthetic twin experiments on Kongsvegen glacier, Svalbard, using a Particle Batch Smoother with 1000 ensemble members. Synthetic observations of albedo, snow depth, and surface temperature are assimilated at two quality levels, under two climatic scenarios, and over 12 years. Assimilation benefit is measured as the percentage improvement in the continuous ranked probability score of the posterior glacier surface mass balance relative to the prior. A single optimally timed high quality observation yields mean improvements of up to 80\%. Larger numbers of low quality observations partially compensate for lower improvement. In the accumulation zone, however, additional snow depth observations degrade performance through particle degeneracy. Optimal timing is governed by the seasonal transitions of the truth trajectory rather than by prior ensemble spread alone. The optimal windows shift by up to six weeks between early and late melting years. Joint assimilation adds value through temporal diversity rather than observational diversity, while independently timed observations outperform same day combinations. The asynchronously optimally timed combined assimilation of three variables sustains improvements of 85 to 97\% across all years in the ablation zone. These findings provide guidelines for adaptive observation scheduling in glacier monitoring and reanalysis.
Temporal modulation offers a fundamentally distinct degree of freedom for active wave control beyond static spatial structuring. Photonic time crystals (PTCs), based on periodic modulation of electromagnetic parameters in time, have expanded photonic band engineering from space to time by enabling controlled energy exchange between light and the modulation. Yet, the use of PTCs to synthesize rotational dynamics and thereby control chirality and circular dichroism (CD) remains largely unexplored. Here, we propose a spatiotemporal PTC whose central cylindrical element is driven by an azimuthally traveling-wave permittivity modulation. Although the structure is geometrically static, its dielectric profile evolves as an effectively rotating pattern in time. This synthetic rotation lifts a static modal degeneracy and produces two nondegenerate counter-rotating states with opposite orbital angular momenta. These chiral modes selectively couple to left- and right-circularly polarized light, giving rise to tunable CD. In addition, the spatiotemporal modulation induces orbital angular-momentum conversion between the Floquet replica bands. Our work reveals the microscopic origin of dynamic chiral response and establishes a strategy for reconfigurable chiral photonics without mechanical motion.
Aqueous redox-flow batteries based on TEMPO derivatives are promising for large-scale energy storage, but their practical use is limited by the chemical instability of the oxidized N -oxoammonium state. In this work, we investigate the degradation of five TEMPO derivatives using ab initio molecular dynamics combined with enhanced sampling. Two proposed degradation mechanisms, ring opening and Hofmann elimination, are examined and their corresponding activation free energies are compared. For all derivatives considered, ring opening exhibits a lower activation free energy than Hofmann elimination, identifying it as the kinetically preferred degradation pathway. The magnitude of the ring-opening barrier, however, varies significantly between molecules, showing that different functionalizations strongly influence its stability toward degradation. The predicted preference for ring opening is consistent with available experimental studies, which have identified or inferred ring-opening degradation for several TEMPO-based catholytes. These results provide an atomistic picture of degradation pathways that are difficult to resolve experimentally and highlight the importance of molecular structure in controlling the kinetic stability of TEMPO derivatives in aqueous electrolytes.
We consider Thouless pumping of two-dimensional quadratic solitons composed from coherently interacting fundamental frequency (FF) and second harmonic (SH) components propagating in a periodic $χ^{(2)}$ material. The pumping is induced by two mutually sliding two-dimensional lattices defined by shallow transverse and longitudinal periodic refractive index modulation. Focusing on solitons in the semi-infinite gap, we find three distinct pumping scenarios: the absence of transport for small-amplitude solitons, non-quantized transport in a transient regime at intermediate amplitudes, and stable quantized transport for solitons with relatively large amplitudes. Different directions of the sliding velocity were investigated resulting in different trajectories of soliton center. The transition to the regime of quantized transport is found to depend strongly on phase mismatch between FF and SH waves that also influences stability properties of two-dimensional $χ^{(2)}$ lattice solitons. We also show that pumping with longitudinal periods in the x and y directions, whose ratio approximates an irrational number, induces quantized transport whose direction rapidly converges, with increase of the accuracy of the approximation, to a limiting pumping direction determined by this irrational number and corresponding to truly incommensurate longitudinal periods.
Many models across the sciences require global distance-decay interactions on large sparse networks. Gravity models, accessibility measures, spatial economic models, and network influence processes all evaluate distance-weighted potential fields: each location accumulates contributions from every other location, weighted by a decaying function of the shortest-path travel cost between them. Computed directly, such aggregate fields require the dense matrix of all pairwise network costs, which scales quadratically in time and memory. Common approximations either discard long-range contributions or fail when interaction is governed by network distances rather than geometric proximity. We introduce HierX, a hierarchical sparse-plus-correction operator for distance-decay potential fields on networks. HierX constructs multi-scale representative layers with explicit correction terms ensuring each location pair contributes exactly once at the finest available resolution. Under bounded-growth assumptions common in spatially embedded networks, applying HierX scales as O(n log n). Systematic benchmarks confirm quasi-linear scaling to 100,000 nodes. In head-to-head comparison with distance cutoff truncation and Nystrom low-rank approximation on 25,000-zone networks, HierX achieves 5-9% RMSE at a fraction of the computational work; Nystrom degrades severely on steep decay kernels. Case studies compute population-weighted accessibility on the 2.58-million-node Great Britain driving network and the 1.77-million-node London pedestrian network: a one-time hierarchy construction (~1 hour) yields a compact reusable operator that then evaluates each national-scale accessibility field in under one second (~700 ms for Great Britain, ~150 ms for London). Open-source code and worked examples are provided.
Low Gain Avalanche Detectors (LGADs) achieve time resolutions below 50 ps through internal charge multiplication, yet the gain is not a fixed device property as the space-charge of the multiplied electron-hole density screens the gain layer field, which suppresses further impact ionization. While TCAD simulations reproduce the average gain suppression, their deterministic treatment of a continuous charge density cannot capture the event-by-event fluctuations that determine the timing performance. We present a Monte Carlo study of space-charge-induced gain suppression in LGADs using the Garfield++ framework, in which the dynamic field of the drifting carriers, computed from a coarse-grained system of charged rings, is superimposed on TCAD-computed static field maps of Centro Nacional de Microelectrónica (CNM) and Hamamatsu Photonics (HPK) produced sensors while preserving single carrier avalanche statistics. The simulated gain reproduces published two-photon absorption-transient current technique measurements across two orders of magnitude in injected charge, including the dependence on the injection depth. For minimum ionizing particles, the event-by-event coupling between deposited charge and gain compresses the Landau tail of the multiplied charge, narrowing the signal arrival time distribution and improving the intrinsic time resolution of the HPK-S1 sensor investigated in this work from 53.9 $\pm$ 1.0 ps to 44.5 $\pm$ 0.9 ps. Treating the primary ionization as a stochastic sequence of charge clusters, it is further shown analytically that any gain suppression that decreases with the cluster size narrows the centroid time distribution.
Pulsatile flow in corrugated tubes arises in hemodynamics and microfluidics, where oscillatory forcing and geometric constrictions jointly determine transport and wall loading. We extend Womersley's classical solution for oscillatory flow in a rigid circular tube to rigid axisymmetric tubes of slowly varying radius. Within the lubrication approximation, we derive closed-form expressions for the axial velocity profile, volumetric flow rate, phase lag, and wall shear stress at an arbitrary Womersley number. Three-dimensional lattice Boltzmann simulations are used to assess the regime of validity of the theory. The analysis shows that the local velocity profile varies strongly along the tube, ranging from plug-like in wide sections to more parabolic near bottlenecks. The cycle-maximum flow rate decreases with increasing corrugation, but this reduction weakens as pulsatility increases, reflecting a crossover from the quasi-steady scaling $\langle R^{-4}\rangle^{-1}$ to the high-frequency scaling $\langle R^{-2}\rangle^{-1}$. The wall shear stress is maximal at the bottleneck and decreases with Womersley number. For sinusoidal corrugations, the bottleneck wall shear stress depends non-monotonically on the corrugation because of the competition between local shear amplification and global hydraulic resistance. Closed-form expressions for the time-averaged wall shear stress and oscillatory shear index further connect the theory to standard hemodynamic metrics.
We address the challenge of scalable uncertainty quantification in large-scale scientific applications, where complex state-of-the-art machine learning methods are often computationally infeasible. Our primary contribution is a simple yet effective empirical Bayes method for automatically tuning the hyperparameters of a flexible, heteroscedastic Spectral-normalized Neural Gaussian Process. This approach retains the expressiveness and uncertainty-awareness of semi-Bayesian neural models while significantly reducing the computational burden by integrating hyperparameter learning directly into the training loop. We demonstrate the practical impact of our method on the task of estimating the dynamic aperture in circular particle accelerators, a fundamental problem in high-energy physics colliders and storage rings, using simulation data from the case of the Large Hadron Collider at CERN. Traditional approaches to DA estimation require extensive particle-tracking simulations, which are prohibitively time-consuming and resource-intensive. Our results show that the proposed method achieves competitive predictive performance and well-calibrated uncertainty estimates at much lower computational cost than state-of-the-art approaches. We stress that, beyond this application, the proposed empirical Bayes framework offers a general solution for training heteroscedastic neural models in situations where manual hyperparameter tuning is impractical. Accordingly, we anticipate that this framework can be applied to other domains that encounter comparable computational limitations.