2025-06-30 | | Total: 85
We propose a novel method, basis vector model material indexing (BVM-MI), for predicting atomic composition and mass density from two independent basis vector model weights derived from dual-energy CT (DECT) for Monte Carlo (MC) dose planning. BVM-MI employs multiple linear regression on BVM weights and their quotient to predict elemental composition and mass density for 70 representative tissues. Predicted values were imported into the TOPAS MC code to simulate proton dose deposition to a uniform cylinder phantom composed of each tissue type. The performance of BVM-MI was compared to the conventional Hounsfield Unit material indexing method (HU-MI), which estimates elemental composition and density based on CT numbers (HU). Evaluation metrics included absolute errors in predicted elemental compositions and relative percent errors in calculated mass density and mean excitation energy. Dose distributions were assessed by quantifying absolute error in the depth of 80% maximum scored dose (R80) and relative percent errors in stopping power (SP) between MC simulations using HU-MI, BVM-MI, and benchmark compositions. Lateral dose profiles were analyzed at R80 and Bragg Peak (RBP) depths for three tissues showing the largest discrepancies in R80 depth. BVM-MI outperformed HU-MI in elemental composition predictions, with mean RMSEs of 1.30% (soft tissue) and 0.1% (bony tissue), compared to 4.20% and 1.9% for HU-MI. R80 depth RMSEs were 0.2 mm (soft) and 0.1 mm (bony) for BVM-MI, vs. 1.8 mm and 0.7 mm for HU-MI. Lateral dose profile analysis showed overall smaller dose errors for BVM-MI across core, halo, and proximal aura regions. In conclusion, fully utilizing the two-parameter BVM space for material indexing led to significantly improved TOPAS MC dose calculations over the conventional HU-MI method, demonstrating BVM-MI's potential to enhance proton therapy planning.
In the realm of computational fluid dynamics, traditional numerical methods, which heavily rely on discretization, typically necessitate the formulation of partial differential equations (PDEs) in conservative form to accurately capture shocks and other discontinuities in compressible flows. Conversely, utilizing non-conservative forms often introduces significant errors near these discontinuities or results in smeared shocks. This dependency poses a considerable limitation, particularly as many PDEs encountered in complex physical phenomena, such as multi-phase flows, are inherently non-conservative. This inherent non-conservativity restricts the direct applicability of standard numerical solvers designed for conservative forms. This work aims to thoroughly investigate the sensitivity of Physics-Informed Neural Networks (PINNs) to the choice of PDE formulation (conservative vs. non-conservative) when solving problems involving shocks and discontinuities. We have conducted this investigation across a range of benchmark problems, specifically the Burgers equation and both steady and unsteady Euler equations, to provide a comprehensive understanding of PINNs capabilities in this critical area.
In Phase II of the Project 8 neutrino mass experiment, electrons from the decays of tritium or 83mKr are detected via their ≈26 GHz cyclotron radiation while contained within a circular waveguide. The signal from a given electron is characterized as a brief chirp, lasting ≲10 ms and changing in frequency by ≲1 MHz/ms. To detect these signals, the Project 8 collaboration developed a data acquisition (DAQ) system tailored to the signal properties. The DAQ is responsible for simultaneously selecting up to three 100 MHz-wide frequency windows to study, detect, and trigger on likely signals from different electron kinetic energies, and for writing the relevant data to disk. We describe the Phase II DAQ system in detail and address how the system was used for data-taking operations.
The institution of money can be seen as a foundational social mechanism, enabling communities to quantify collectively regulate economic processes. Money can be said, indeed, to constitute the micro-macro link in economics. This paper reviews influential views on the nature of money in economics and sociology, contrasting them to the relatively limited findings of recent agent-based models of "the emergence of money". Noting ample room for novel combinations of sociological and formal methods to drive insight into the many roles played by money in the economy, we conclude by indicating research directions in which we believe this combination can provide new answers to old questions in monetary theory
Correcting astigmatism and ellipticity in laser beams is critical for improving performance in many applications like microscopy, atomic physics, quantum information processing, and advanced manufacturing. Passive correction methods based on cylindrical lens telescopes require choosing lenses with precise focal lengths, effectively limiting the range of tunability when using standard catalog optics. Active solutions based on diffractive optical elements can achieve superior performance, but they are bulky, expensive, and suffer from finite diffraction efficiency and added complexity. Here, we introduce a simple method to convert astigmatic elliptical beams into circular Gaussian beams without astigmatism. Our method comprises three cylindrical lenses. The first lens focuses the beam along its major axis to create a plane where the intensity profile is radially symmetric. The second and third lenses are placed one behind the other in that plane at a relative angle, acting as a biaxial lens pair with tunable focal lengths. By adjusting the relative angle of the lenses, the two separate beam waists of the astigmatic beam can be overlapped, resulting in a circular Gaussian beam without astigmatism. We theoretically validate our method, numerically quantify its robustness to experimental imperfections, and experimentally demonstrate its ability to circularize the output beam of a commercial laser source. Our method corrects astigmatism and ellipticity in laser beams without requiring precise focal length matching, offering greater flexibility than other passive solutions and greater cost-effectiveness than active methods. Its simple and compact design makes it well suited for integration into both tabletop optical setups and industrial systems.
The metriplectic formulation of collisional guiding-center Vlasov-Maxwell-Landau theory is presented. The guiding-center Landau collision operator, which describes collisions involving test-particle and field-particle guiding-center orbits, is represented in terms of a symmetric dissipative bracket involving functional derivatives of the guiding-center Vlasov phase-space density Fgc and the electromagnetic fields (Dgc,B), where the guiding-center displacement vector Dgc≡E+4πPgc is expressed in terms of the electric field E and the guiding-center polarization Pgc. This dissipative Landau bracket conserves guiding-center energy-momentum and angular momentum, as well as satisfying a guiding-center H-theorem.
The mechanism by which odorants are recognized by olfactory receptors remains primarily unresolved. While charge transport is believed to play a significant role, its precise nature is still unclear. Here, we present a novel perspective by exploring the interplay between the intrinsic energy scales of odorant molecules and the gap states that facilitate intermolecular charge transport. We find that odorants act as weak tunneling conductors mainly because of the limited magnitude of electronic coupling. This behavior is further connected to electron-phonon interactions and reorganization energy, suggesting that physically meaningful values for these parameters emerge only in the deep off-resonant tunneling regime. These findings complement the swipe card model of olfaction, in which an odorant needs both the right shape to bind to a receptor and the correct vibrational frequency to trigger signal transduction. Moreover, they reveal that the underlying mechanisms are much more complex than previously assumed.
Three-body Förster resonances controlled by a dc electric field are of interest for the implementation of three-qubit quantum gates with single atoms captured in optical traps and laser-excited into strongly interacting Rydberg states. In Ref. [P. Cheinet et al., Quantum Electronics 50(3), 213 (2020)], we proposed and analyzed a new type of three-body Förster resonance 3×nP3/2→nS1/2+(n+1)S1/2+nP1/2 that can be realized with Rb Rydberg atoms for an arbitrary principal quantum number n. Its peculiarity is that the third atom goes into a state with a total angular moment J=1/2, which has no Stark structure, so two-body Förster resonances are completely absent. In the present work, an extended theoretical study of this three-body Förster resonance is performed for various spatial configurations of three interacting Rb Rydberg atoms and conditions for their experimental implementation are determined. It was found that one of the resonances has a weak dependence of the resonant electric field on the distance between atoms and is therefore most suitable for performing experiments to observe coherent oscillations of populations of collective three-body states and implement three-qubit quantum gates based on them.
Abnormal patterns of ventricular repolarisation contribute to lethal arrhythmias in various cardiac conditions, including inherited and acquired channelopathies, cardiomyopathies, and ischaemic heart disease. However, methods to detect these repolarisation abnormalities are limited. In this study, we introduce and assess a novel simulation-based method to infer ventricular activation and repolarisation times from the 12-lead electrocardiogram (ECG) and magnetic resonance-derived ventricular anatomical reconstruction, applicable for the first time to both healthy controls and cases with abnormal repolarisation. First, ventricular activation times were reconstructed through iterative refinement of early activation sites and conduction velocities, until the model and target QRS complexes matched. Then, ventricular repolarisation times were reconstructed through iterative refinement of ventricular action potential durations and an action potential shape parameter until the model and target T waves matched, including regularisation. Repolarisation inference was evaluated against 18 benchmark simulations with known repolarisation times, including both control and hypertrophic cardiomyopathy (HCM) cases with abnormal repolarisation. Inferred repolarisation times showed good agreement with the ground truth in control and HCM (Spearman r=0.63+/-0.11 and 0.65+/-0.19, respectively), with inferred model T waves closely matching the target T waves (r=0.81+/-0.05 and 0.78+/-0.08, respectively). The method further demonstrated flexibility in reconstructing the macroscopic patterns of delayed repolarisation across a range of abnormal ventricular repolarisation sequences, demonstrating applicability to a range of pathological cases. Simulation-based inference can accurately reconstruct repolarisation times from the 12-lead ECG in cases with both normal and abnormal repolarisation patterns.
The description of riblets and other drag-reducing devices has long used the concept of longitudinal and transverse protrusion heights, both as a means to predict the drag reduction itself and as equivalent boundary conditions to simplify numerical simulations by transferring the effect of riblets onto a flat virtual boundary. The limitation of this idea is that it stems from a first-order approximation in the riblet-size parameter s+, and as a consequence it cannot predict other than a linear dependence of drag reduction upon s+; in other words, the initial slope of the drag-reduction curve. Here the concept is extended to a full asymptotic expansion using matched asymptotics, which consistently provides higher-order protrusion heights and higher-order equivalent boundary conditions on a virtual flat surface. This procedure also allows us to explore nonlinearities of the Navier-Stokes equations and the way they enter the s+-expansion, with somewhat surprising results.
We present new and updated methods for reducing transport caused by electrostatic ion temperature gradient (ITG) driven turbulence in quasi-isodynamic (QI) configurations. We first show an updated model for the threshold (critical) gradient of localized toroidal ITG modes. It is then argued that it is desirable for ITG modes to "split" and localize in separate curvature drift wells, which is leveraged to produce a six-field-period QI configuration with a high critical gradient. We show that the destabilizing effect of kinetic electrons (Costello & Plunk 2025a) on localized ITG modes can be minimized in a magnetic field structure, which we dub "inverse mirror". Applying a general optimization target that improves ITG stability above the critical gradient yields an inverse mirror configuration, which produces heat fluxes below or equal to the W7-X high mirror configuration for a range of applied density gradients.
The relativistic second-order corrections at orders mα6 and mα6(m/M) in hydrogen molecular ions are calculated. Convergence of numerical results is studied, which allows estimating the relative numerical uncertainty to be less then 10−5. This accuracy is sufficient to enable future improvement of theoretical predictions beyond the 1 ppt (part-per-trillion) precision level for the ro-vibrational transition frequencies.
This paper introduces a novel multi-species Enskog-Vlasov solver. It is used to determine evaporation coefficients of fluids under high-pressure conditions, a critical factor for efficient fuel mixing in internal combustion engines. The solver handles collisions for different fluid species separately, thereby accurately capturing species-specific interactions essential for realistic evaporation modeling. A new pair correlation function based on the BMCSL equation of state is employed to enhance modeling accuracy, though compliance with Onsager relations remains to be explored. Validation through various numerical simulations demonstrates the solver's capability. Results from binary fluid relaxation simulations closely match molecular dynamics (MD) data, highlighting accurate collision frequency modeling. Comparative studies against state-of-the-art models verify the solver's precision in predicting liquid-vapor equilibria and evaporation coefficients across diverse pressure scenarios. Detailed simulations of argon-neon mixtures with realistic particle diameters and masses underscore the significant improvements achieved by employing this multi-species approach over traditional single-species methods. Comparisons with the SAFT-VRQ Mie and classical density functional theory confirm the solver's reliability in predicting liquid and vapor compositions and detailed density profiles at interfaces over a wide range of pressures and temperatures. Further analyses illustrate complex dependencies of evaporation and condensation coefficients on temperature and pressure, consistent with existing research findings. Overall, this work advances computational modeling of multi-species evaporation processes in high-pressure environments at different temperatures, providing essential data for improved combustion modeling and broader applications in multi-phase fluid dynamics.
We study double ionization in a one-dimensional model of a three-electron atom exposed to strong short laser fields. The corresponding ab initio simulations reveal a characteristic ``knee'' in the ionization yield at lower field intensities than in comparable two-electron systems. To understand the origin of this shift, a quantitative semi-analytic model for double ionization of a three-electron atom has been created. In this model, ionization is treated as a combination of different scenarios and channels. Each channel is simulated with either Quantitative Rescattering Theory (QRS) or two-electron ab initio model. Analysis of the semi-analytical results indicates that the shifted part of the ``knee'' is dominated by a recollision-excitation channel in which the rescattering electron initially escapes the parent atom from an inner orbital. This counterintuitive phenomenon is then discussed and explained.
This study examines the complex movement patterns and behavioral characteristics of wild boars using GPS telemetry data collected over a two-month period. Our methodological approach centers on the application of a Hidden Markov Model (HMM) to discern distinct behavioral states embedded within the trajectories. Furthermore, the study aimed to construct behavioral networks, derived from these segmented trajectories. The resultant network structures showed that the hidden behavioral patterns are mostly independent of geographical locations. While most locations have many behaviors occuring in them, our findings also suggest that Finally, the research incorporates a spatial trajectory analysis, complemented by raster data validation, to potentially delineate areas acting as repellents within the ecological context of Hainich National Park in Germany.
Comorbidity networks, which capture disease-disease co-occurrence usually based on electronic health records, reveal structured patterns in how diseases cluster and progress across individuals. However, how these networks evolve across different age groups and how this evolution relates to properties like disease prevalence and mortality remains understudied. To address these issues, we used publicly available comorbidity networks extracted from a comprehensive dataset of 45 million Austrian hospital stays from 1997 to 2014, covering 8.9 million patients. These networks grow and become denser with age. We identified groups of diseases that exhibit similar patterns of structural centrality throughout the lifespan, revealing three dominant age-related components with peaks in early childhood, midlife, and late life. To uncover the drivers of this structural change, we examined the relationship between prevalence and degree. This allowed us to identify conditions that were disproportionately connected to other diseases. Using betweenness centrality in combination with mortality data, we further identified high-mortality bridging diseases. Several diseases show high connectivity relative to their prevalence, such as iron deficiency anemia (D50) in children, nicotine dependence (F17), and lipoprotein metabolism disorders (E78) in adults. We also highlight structurally central diseases with high mortality that emerge at different life stages, including cancers (C group), liver cirrhosis (K74), subarachnoid hemorrhage (I60), and chronic kidney disease (N18). These findings underscore the importance of targeting age-specific, network-central conditions with high mortality for prevention and integrated care.
Optical beam steering has become indispensable in free-space optical communications, light detection and ranging (LiDAR), mapping, and projection. Optical phased array (OPA) leads this field, yet conventional versions still suffer from a narrow steering field of view (FOV), insufficient sidelobe suppression, and limited angular resolution. Thin-film lithium niobate (LN), with its strong Pockels electro-optic (EO) effect, offers a powerful integrated-photonics platform to overcome these limitations. Here we present a two-dimensional (2D) EO-steered OPA based on a non-uniformly spaced X-cut thin-film LN ridge-waveguide array. A superlattice ridge design suppresses optical crosstalk to -20 dB, enabling low-sidelobe far-field radiation. Using particle swarm optimization (PSO) method, we transform a uniformly spaced array into an optimized non-uniform design, largely improving angular resolution while maintaining sidelobe suppression. When combined with a single-radiating trapezoidal end-fire emitter incorporating an etched grating, the device produces a main-lobe beam width of 0.99 degree*0.63 degree from an aperture of only 140 um*250 um, achieving a wide 2D steering range of 47 degree*9.36 degree with a 20 dB sidelobe-suppression ratio. These results highlight thin-film LN OPA as a compelling route toward heterogeneous, compact, and high-performance EO beam-steering modules and ultra-miniaturized optical modulators.
Extreme weather events are projected to intensify global migration, increase resource competition, and amplify socio-spatial phenomena, including intergroup conflicts, socioeconomic inequalities, and unplanned displacements, among others. Addressing these challenges requires consolidating heterogeneous data to identify, estimate, and predict the dynamical process behind climate-induced movements. We propose a novel hybrid approach to reconstruct hazard-induced displacements by analysing the statistical properties of stochastic walks that explore the spatial network constructed from real movements. The likely trajectories produced by the walker inform the typical journey experienced by individuals, helping identify likely hazards encountered when fleeing high-risk areas. As a proof of concept, we apply this method to Somalia's detailed displacement tracking matrix, containing 20,220 movements dating from February 8th to June 18th, 2025. We reconstruct the likely routes that displaced persons could have taken when fleeing areas affected by conflict or climate hazards. We find that individuals using the most likely paths based on current flows would experience mainly droughts and conflicts, while the latter become less prominent at every next step. We also find that the probability of conflict and drought across all trajectories is widely dispersed, meaning that there is no typical exposure. This work provides an understanding of the mechanisms underlying displacement patterns and a framework for estimating future movements from areas expected to face increasing hazards. We also provide avenues to develop tools that can support policymakers implementing proactive responses and early interventions.
Elite football is believed to have evolved in recent years, but systematic evidence for the pace and form of that change is sparse. Drawing on event-level records for 13,067 matches in ten top-tier men's and women's leagues in England, Spain, Germany, Italy, and the United States (2020-2025), we quantify match dynamics with two views: conventional performance statistics and pitch-passing networks that track ball movement among a grid of pitch (field) regions. Between 2020 and 2025, average passing volume, pass accuracy, and the percent of passes made under pressure all rose. In general, the largest year-on-year changes occurred in women's competitions. Network measures offer alternative but complementary perspectives on the changing gameplay in recent years, normalized outreach in the pitch passing networks decreased, while the average shortest path lengths increased, indicating a wider ball circulation. Together, these indicators point to a sustained intensification of collective play across contemporary professional football.
We present a framework to calculate the scale-resolved turbulent Prandtl number Pr_t for the well-mixed and highly inertial bulk of a turbulent Rayleigh-Bénard mesoscale convection layer at a molecular Prandtl number Pr=10^{-3}. It builds on Kolmogorov's refined similarity hypothesis of homogeneous isotropic fluid and passive scalar turbulence, based on log-normally distributed amplitudes of kinetic energy and scalar dissipation rates that are coarse-grained over variable scales r in the inertial subrange. Our definition of turbulent (or eddy) viscosity and diffusivity does not rely on mean gradient-based Boussinesq closures of Reynolds stresses and convective heat fluxes. Such gradients are practically absent or indefinite in the bulk. The present study is based on direct numerical simulation of plane-layer convection at an aspect ratio of \Gamma=25 for Rayleigh numbers 10^5\leq Ra \leq 10^7. We find that the turbulent Prandtl number is effectively up to 4 orders of magnitude larger than the molecular one, \Pr_t\sim 10. This holds particularly for the upper end of the inertial subrange, where the eddy diffusivity exceeds the molecular value, \kappa_e(r)>\kappa. Highly inertial low-Prandtl-number convection behaves effectively as a high-Prandtl number flow, which also supports previous models for the prominent application case of solar convection.
In this study we investigate how hierarchical structures within the Roman Catholic Church shape the ideological orientation of its leadership. The full episcopal genealogy dataset comprises over 35,000 bishops, each typically consecrated by one principal consecrator and two co-consecrators, forming a dense and historically continuous directed network of episcopal lineage. Within this broader structure, we focus on a dataset of 245 living cardinals to examine whether genealogical proximity correlates with doctrinal alignment on a broad set of theological and sociopolitical issues. We identify motifs that capture recurring patterns of lineage, such as shared consecrators or co-consecrators. In parallel, we apply natural language processing techniques to extract each cardinal's publicly stated positions on ten salient topics, including LGBTQIA+ rights, women's roles in the Church, liturgy, bioethics, priestly celibacy, and migration. Our results show that cardinals linked by specific genealogical motifs, particularly those who share the same principal consecrator, are significantly more likely to exhibit ideological similarity. We find that the influence of pope John Paul II persists through the bishops he consecrated, who demonstrate systematically more conservative views than their peers. These findings underscore the role of hierarchical mentorship in shaping ideological coherence within large-scale religious institutions. Our contribution offers quantitative evidence that institutional lineages, beyond individual background factors, may have an impact on the transmission and consolidation of doctrinal positions over time.
We present a model and a rigorous method to calculate the transmission coefficient of silicon micro-rings with complex waveguide cross-section including non-linear effects and self-heating, with very short simulation times. The method is applied to the design of MRRs in the SISCAP platform with high Q and reduced non-linearity, namely due to two photon absorption and free carrier absorption. We demonstrate that the free carrier diffusion in rib waveguides and Shockely-Read-Hall recombination play a fundamental role in reducing the impact of non-linearities in the ring.
In this study, we systematically investigate how the accuracy of CIA spectra depends on the quality of the ab initio data used. We evaluate quantitatively the impact of different quantum chemical methods and basis sets on spectral features, finding that even the lowest-level calculations are accurate to approximately 10 % at room temperature, and better at higher temperatures. This study also reveals a previously unreported double-peak structure for the He-Ne complex, which cannot be described by simple but commonly used single-exponential models for the short-range dipole. Our analysis shows that the range of internuclear distances relevant for CIA spectra varies with temperature, with short-range interactions becoming increasingly important at high temperatures. The long-range van der Waals induced dipoles never contribute. These findings provide new insights into the temperature-dependent behavior of CIA spectra and emphasize the importance of accurate modeling of short-range interactions for reliable astronomical modeling.
Betatron radiation is produced in Laser Plasma Accelerators when the electrons are accelerated and simultaneously wiggle across the propagation axis. The mechanisms of electron acceleration and X-ray radiation production follow different scaling laws, and the brightest X-ray radiation is often produced for an electron beam with a lower quality in terms of energy and divergence. Here, we report a laser-driven Betatron X-ray source where the plasma density profile is tailored in order to separate the acceleration and wiggler stages, which allows for the independent optimizations of acceleration and X-ray production. We demonstrate this concept experimentally, and show that the Betatron photon energy can be controlled by adjusting the length of the plasma wiggler. This scheme offers a path to overcome the limitations of conventional Betatron sources, enabling the production of bright, stable, energetic, and collimated X-ray beams.
High-resolution processing of seismic signals is crucial for subsurface geological characterization and thin-layer reservoir identification. Traditional high-resolution algorithms can partially recover high-frequency information but often lack robustness, computational efficiency, and consideration of inter-trace structural relationships. Many deep learning methods use end-to-end architectures that do not incorporate prior knowledge or address data domain disparities, leading to limited generalization.To overcome these challenges, this paper presents the Domain-Adaptive Knowledge Distillation Network (DAKD-Net), which integrates a knowledge distillation strategy with a domain adaptation mechanism for high-resolution seismic data processing. Trained on datasets from forward modeling, DAKD-Net establishes physical relationships between low and high-resolution data, extracting high-frequency prior knowledge during a guided phase before detail restoration without prior conditions. Domain adaptation enhances the model's generalization to real seismic data, improving both generalization capability and structural expression accuracy.DAKD-Net employs a U-Net backbone to extract spatial structural information from multi-trace seismic profiles. The knowledge distillation mechanism enables prior knowledge transfer, allowing recovery of high-resolution data directly from low-resolution inputs. Domain-adaptive fine-tuning further enhances the network's performance in actual survey areas. Experimental results show that DAKD-Net outperforms traditional methods and classical deep networks in longitudinal resolution and complex structural detail restoration, demonstrating strong robustness and practicality.