2026-03-23 | | Total: 67
Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse datasets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000 times faster than DFT. While previous MLIP training datasets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) dataset provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pre-trained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation-solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion-correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations.
In this letter, we numerically investigate a long-haul coherent data transmission system with a cascade of semiconductor optical amplifiers (SOAs). We exploit low-complexity neural networks that can be implemented in real time to compensate for the accumulated distortions induced by a cascade of SOAs. This equalization provides an order-of-magnitude reduction in bit error rate at low dispersion (in the O-band), whereas higher dispersion degrades performance.
Quasi-axisymmetric stellarators are the stellarator analogue of the axisymmetric tokamak, retaining many of its favorable confinement properties, its compacity and its relative coil simplicity, while avoiding its principal limitation, the need for an inductively driven plasma current. Despite these attractive physics properties, the development of quasi-axisymmetric configurations has been severely constrained by the absence of an experimentally validated divertor concept compatible with their large bootstrap current. In this Letter, approximately quasi-axisymmetric fields, complemented with piecewise omnigenous perturbations, are proposed as the basis for a new strategy towards a stellarator reactor that simultaneously achieves simple coil geometries, tokamak-like confinement properties and, through tailoring of the bootstrap current, compatibility with an island divertor. Implications for attaining a high bootstrap current fraction in tokamak devices are also discussed.
This work presents a variational physics-informed deep learning framework for phase-field modelling of brittle crack propagation in anisotropic media. Previous Deep Ritz Method (DRM) approaches have focused on second-order, isotropic phase-field fracture formulations. In contrast, the present work introduces, for the first time within a variational deep learning setting, a family of higher-order anisotropic phase-field models through a generalised crack density functional. The resulting fracture problem is solved by minimising the total energy using the DRM. The trial space is enriched with higher-order B-spline basis functions to represent higher-order gradients accurately and stably, thereby eliminating the need for conventional automatic differentiation. The methodology is assessed for isotropic, cubic, and orthotropic fracture surface energy densities. Numerical examples demonstrate direction-dependent crack growth in anisotropic cases, highlighting the capability of the method to accurately capture this behaviour.
The Microstrip Silicon Detector (MSD) is one of the subsystems of the FragmentatiOn Of Target (FOOT) experiment whose goal is to measure double differential nuclear fragmentation cross sections for applications in particle therapy and radioprotection in space. The MSD, composed of six 150 $μ$m-thick silicon sensors arranged in three X-Y measuring planes, is part of the FOOT experiment tracking region. In this work, we propose a new method to set the MSD thresholds for the clustering algorithm independently of the other detectors of the experiment. The obtained values can be used as reference values from which to perform a threshold scan in order to evaluate the single-ion detection efficiency of each MSD sensor and to orient the clustering analysis using tracking information.}
A new library VecAmpFit for multidimensional amplitude analyses in high-energy physics has been developed for an ongoing amplitude analysis at Belle II experiment. It includes a fitter performing likelihood calculation and explicitly-vectorized subprograms for amplitude implementation. The fitter supports explicit gradient calculation and simultaneous fitting of multiple data sets.
Accurately calculating time delays between signals is pivotal in many modern physics applications. One approach to estimating these delays is computing the cross-spectrum in the time-frequency domain. Linear time-frequency representations, such as the continuous wavelet transform (CWT), are widely used to construct these cross-spectra. However, it is well known that the frequency resolution is inherently limited by the localized nature of the convolving wavelet. Moreover, the functional form of the CWT cross-spectrum is not a proper correlation measure and typically requires post-processing smoothing. Conversely, quadratic representations achieve joint time-frequency resolution approaching the Gabor-Heisenberg limit while also providing an adequate measure of similarity between the signals. Motivated by these advantages, we propose a time-delay estimation method based on the Wigner-Ville Distribution (WVD). Considering nonstationary signals arising from two typical wave-physics scenarios, we show that the WVD yields more accurate time-delay estimates with lower uncertainty, particularly in the most energetic frequency bands.
The $Δ$SCF DFT approach defines the system energy as a function of orbital occupancy. Inspired by Landau Fermi liquid theory, we develop an occupancy extrapolation (OE) method that captures excited-state energies via a Taylor expansion of the energy with respect to occupation fluctuation from a reference state. OE retains the physics of $Δ$SCF while offering a physical interpretation of excitation energies as sums of quasiparticle energies and their generalized screened interactions. It yields accurate valence, Rydberg, and charge-transfer excitation energies at $O(N^3)$ cost, avoids separate SCF calculations for each excited state, and enables efficient large-scale excited-state simulations from ground-state calculations.
Ensembles of regional-global climate model combinations show substantial spread in projected wind and solar resources. Using 31 RCM-GCM pairs, we quantify the sources of this spread with a spatially and seasonally resolved variance decomposition, separating contributions from RCMs and GCMs. For both wind speed and solar radiation, RCMs dominate the variability in the absolute historical fields. In contrast, projected changes in wind speed are largely controlled by the driving GCMs, except in mountainous regions where RCM-induced variance becomes larger than that induced by GCMs. For solar radiation, contributions are strongly season-dependent, with RCMs dominating in summer and GCMs in winter. Our findings support that GCM and RCM variability together define the uncertainty of wind and solar climate projections. This provides guidance for designing climate model ensembles that better support uncertainty-aware energy system decisions under climate change.
Water-based Liquid Scintillator (WbLS) was proposed over a decade ago as a novel detector medium that might allow the separation and tuning of the relative ratio of the Cherenkov and Scintillation signals. A detector deploying this technology could combine GeV-scale and MeV-scale neutrino detection at scale. Furthermore, the metal-loading capability of such a material enables neutron tagging and allows the effective particle containment to be tuned. WbLS is attractive both for the potential to use it in large detectors and the ability to modify the configuration in situ. At Brookhaven National Laboratory (BNL), two prototypes have been built for understanding WbLS properties and stability, with masses of 1-ton and 30-ton, respectively. We present here the 30-ton prototype detector design, installation, and operation. Results from the analysis of data collected in the two detectors will follow in future publications.
Machine-learning electronic Hamiltonians achieve orders-of-magnitude speedups over density-functional theory, yet current models omit long-range Coulomb interactions that govern physics in polar crystals and heterostructures. We derive closed-form long-range Hamiltonian matrix elements in a nonorthogonal atomic-orbital basis through variational decomposition of the electrostatic energy, deriving a variationally consistent mapping from the electron density matrix to effective atomic charges. We implement this framework in HamGNN-LR, a dual-channel architecture combining E(3)-equivariant message passing with reciprocal-space Ewald summation. Benchmarks demonstrate that physics-based long-range corrections are essential: purely data-driven attention mechanisms fail to capture macroscopic electrostatic potentials. Benchmarks on polar ZnO slabs, CdSe/ZnS heterostructures, and GaN/AlN superlattices show two- to threefold error reductions and robust transferability to systems far beyond training sizes, eliminating the characteristic staircase artifacts that plague short-range models in the presence of built-in electric fields.
Localized surface plasmon resonances are self-sustained, collective oscillations of free electrons in metallic nanostructures. They have a wide range of applications. The most common plasmonic metals are noble metals, such as gold and silver. However, there are applications, such as surface-enhanced Raman spectroscopy, in which using non-noble metals is advantageous. This review summarizes the investigation of localized surface plasmons in non-noble metal nanoparticles, providing an overview of the plasmonic properties of non-noble metals. We cover the following metals: aluminium (Al), antimony (Sb), bismuth (Bi), chromium (Cr), copper (Cu), gallium (Ga), indium (In), lead (Pb), magnesium (Mg), molybdenum (Mo), nickel (Ni), potassium (K), selenium (Se), sodium (Na), tellurium (Te), tin (Sn), titanium (Ti), tungsten (W), and zinc (Zn). Our summary therefore compares the plasmonic properties of non-noble metals and briefly introduces their potential to the readers.
Localized cavitation in liquids and soft tissues, typically initiated by the rarefaction phase of high-amplitude ultrasound waves, is leveraged in several biomedical applications such as ablation techniques and drug delivery with vaporizing agents. However, safety considerations aimed at avoiding unwanted bubble activity outside the targeted region pose a limit to the maximum allowed peak rarefaction pressure, which on the other hand can hinder the therapeutic efficacy of these techniques. This study shows that a purely compressive shock wave can generate localized, negative pressure and initiate cavitation inside a sub-millimetric perfluorohexane droplet, without requiring any externally applied rarefaction wave. The Gouy phase shift is identified as the physical mechanism responsible for the conversion of positive pressure into tension during shock focusing, and its occurrence is demonstrated through numerical simulations and direct experimental measurements. Comparison of the regions affected by cavitation, visualized \emph{in-situ} by means of high-speed x-ray phase-contrast imaging, with prediction from Classical Nucleation Theory suggests homogeneous nucleation as the underlying mechanism behind bubble formation. The presented findings offer valuable insights into the physics of shock wave propagation which can inspire the development of novel acoustic driving strategies for cavitation generation, facilitating the reduction of negative pressures outside the target region and improving the safety and precision of biomedical treatments.
The SWUF-3D drone fleet is used in the atmospheric boundary layer (ABL) for in situ turbulence measurements of complex flows, such as in mountainous terrain or wind turbine wakes. Previous calibrations for measuring vertical wind speed $w$ using the drones' avionics data were performed on field data, limiting applicability to low winds ($\leq 8~\mathrm{m}\,\mathrm{s}^{-1}$) and being prone to high uncertainties. To overcome these limitations, we calibrate $w$ measurement in a wind tunnel and validate it in field measurements. Calibration is performed in a wind tunnel with an active grid used to deflect horizontal flow into the vertical. This wind is measured with a multi hole probe, while wind forces acting on the drone are determined from the avionics data, allowing an empirical relationship between these quantities. For validation, we conduct comparative fleet measurements with up to 10 drones simultaneously around an array of meteorological masts equipped with sonic anemometers. The results show high accuracy for turbulence statistics: the variance determination for $w$ has a root mean square error (RMSE) of 0.12~$\mathrm{m^2\,s^{-2}}$ and a normalized RMSE (nRMSE) of 17.0~\%, for the horizontal wind components the RMSEs are $\sim$0.3~$\mathrm{m^2\,s^{-2}}$ and nRMSEs $\sim$25~\%. The RMSEs for the covariances of the components are $<\,$0.3~$\mathrm{m^2\,s^{-2}}$. The variance spectra of $w$ measured with drones and reference sensors agree in all frequency ranges, the RMSE for covariances between different measurement points is $\sim$0.1~$\mathrm{m^2\,s^{-2}}$. Accurate $w$ retrieval at all wind speeds sustainable by the drone enables studies of strongly three-dimensional flows, supports eddy-covariance flux estimation, enables resolving diurnal turbulence evolution in the ABL, and improves spatial turbulence characterization.
We present the positron coupled cluster singles and doubles (POS-CCSD) method to calculate positron binding energies in molecules. This framework treats electrons and positrons on an equal footing and includes up to simultaneous double-electron-single-positron excitations. We benchmark the approach by computing binding energies for atomic anions and several polar and non-polar polyatomic systems, comparing the results with independent theoretical studies and, where available, experimental data. The fully converged results for H$^{-}$ are in excellent agreement with quantum Monte Carlo and multi-reference configuration interaction results. Quantitative agreement with experiments is not reached in the present study due to the slow convergence of the binding energy with respect to the size of the orbital bases for the electrons and the positron. However, the POS-CCSD results underscore the critical role of electron correlation in the description of electron-positron systems required for a balanced description of these complex systems. In addition, we examine nuclear relaxation effects following positron attachment in LiH.
Quantitative high-resolution transmission electron microscopy (HRTEM) provides an indispensable means to understand the structure-property relationships of a material in atomic dimensions. Successful quantification requires reliable retrieval of essential atomic structural information despite artifacts arising from unwanted but practically unavoidable imaging imperfections. Experimental observation carried out in tandem with model-based iterative image simulation shows vast applications in quantitative structural and chemical determination of objects spanning zero to three dimensions [Prog. Mater. Sci. 133, 101037, 2023]. However, the large number of parameters involved in the simulations make the current multi-step, user-guided iterative approach highly time consuming, thereby restricting its application primarily to small sample areas and to experienced users. In this work, we implement and apply a physics-informed Bayesian optimization (BO) framework to advance HRTEM quantification towards full automation and large-field-of-view analysis. Unlike conventional optimization approaches, our method adopts a stepwise strategy that fully leverages the strength of BO in handling high-dimensional parameters, while its probabilistic engine rigorously and efficiently refines the parameter space to enable rapid quantification. Using a BaTiO3 single crystal that contains heavy, medium and light elements as a model system, we demonstrate that the three-dimensional crystal structure can be determined from a single HRTEM image with a three to four order-of-magnitude improvement in time efficiency. This approach thus opens new avenues for fast and automated image quantification over larger sample volumes and, potentially, in the time domain.
Determination of orientation in the imaged sample/scene has a large application potential when the anisotropy of properties is analysed, usually, under a linearly polarised illumination. This study combined several improvements of microscopy imaging: use of an incoherent white illumination source (a lamp) with a spectral filter to define a spectral window, a plastic circular polariser to image with circularly polarised light (instead of linear), and a 4-pol. camera with integrated polarisers for simultaneous acquisition of images at four pi/4-azimuth shifts. In the transmission mode, a high-fidelity readout of form birefringent optical elements, q-plates, was achieved using a fitting procedure based on the analytical expression of S0 (intensity) Stokes parameters at the pixel level. In the reflection mode, the S0 fit was used to determine the azimuth orientation of the q-plates, as well as a generic Amp x cos(2Q1 - 2Q0) + Offset fit (at pixel level) applied to images taken at four Q1 azimuths. The 4-pol. analysis in reflection under circularly polarised illumination is discussed.
Time domain airborne electromagnetic (AEM) surveying is a mature geophysical tool for imaging the Earth's shallow subsurface. It produces images of the electromagnetic conductivity structure of the earth, down to depths of a few hundred metres. The AEM method is fast, with aircraft acquiring data at speeds of 100-300 km/hr, making it an ideal near-surface reconnaissance tool. The physics of the AEM method are sensitive primarily to the subsurface conductivity, which is influenced by a range of geological factors such as mineral content, porosity, and water content and chemistry. In addition, the inferred subsurface conductivity depends on the accurate measurement and modelling of airborne transmitter and receiver geometries. In this work, we present inferences of the subsurface conductivity over Lake Menindee, New South Wales, Australia, using data from various AEM systems over the period 2014-2024. The lake storage has varied dramatically over this time and while this difference in storage volume undoubtedly influences the near surface conductivity, a remarkably consistent interpretation of the regional geology emerges. While the upper ten metres of the modelled depth sections exhibit the greatest time-variability in inferred electromagnetic conductivity, a correlation of lakebed near-surface conductivity with the lake water volume cannot robustly be established. We also provide information theoretic calculations for each inversion result to aid in their quantitative comparison. The implications of our study are that subtle, shallow, hydrogeological changes are difficult to image with repeat overflights. Conversely, we establish that different AEM systems robustly image the regional geo-electric structure of the near surface, validated by known stratigraphy and borehole conductivity logs.
In this contribution paper, we construct a two-dimensional non-Hermitian (NH) photonic crystal (PhC) to prototype its NH skin effect for experimental proposal. Based on the tight-binding model for Lieb lattice with NH coupling, a nontrivial spectral winding number is pinpointed for certain eigenstates, which translates to geometry-dependent skin modes in a PhC slab with tilt boundaries. For ease of implementation, complex refractive indices are employed for the Lieb unit cell of PhC to emulate the NH coupling. Validated by full wave simulation, our work under scores the boundary dependence of the geometric skin effect, and provides a concrete prototype design of NH skin effect easily implemented in classical wave systems by state-of-the-art of topological metamaterial platforms.
Objective: This research focuses on the design and the evaluation of an ultrathin cylindrical metasurface for improving the transmit efficiency of travelling-wave magnetic resonance imaging (MRI) of the human brain. Methods: The cylindrical metasurface was designed as a compact and lightweight replacement of a high-permittivity dielectric waveguide previously proposed for the same purpose. Based on the optimized unit cell geometry, full numerical model of the cylindrical metasurface in the presence of a voxel human body model was constructed. We compared the proposed metasurface to the dielectric waveguide in the travelling-wave setup experimentally, including in-vivo measurements performed with a healthy volunteer. Results: The proposed metasurface showed an improved B1 + homogeneity (by 17.3%), transmit efficiency (by 27.4%), and SAR-efficiency (by 23%) compared to the dielectric waveguide. Conclusion: The proposed cylindrical metasurface, optimized for field enhancement in the human brain at 7 T in the travelling-wave excitation regime, can further improve the transmit efficiency and homogeneity in the region of interest compared to state-of-the-art structures for travelling-wave MRI, at the same time, granting the advantages of light weight and compactness. Significance: The proposed metasurface can be used for increasing the transmit efficiency of single-channel travelling-wave excitation setups suitable for prospective research and clinical MRI tasks.
Structured-core thulium-doped fibers were developed to reduce heat load, enable shorter-wavelength operation, and achieve a pedestal-free design. In a proof-of-principle experiment, laser slope efficiencies of 52% at 1907 nm and 54% at 1940 nm were achieved with respect to absorbed power.
Transition-metal complexes (TMCs) are promising photosensitizers for Type~I photodynamic therapy (PDT), where electron-transfer processes can generate reactive oxygen species under hypoxic conditions. Yet identifying candidates with the required ground- and excited-state redox energetics remains challenging across the vast chemical space of TMCs. Here, we develop a data-efficient active learning (AL) framework for the discovery of Type~I active TMC photosensitizers by combining a chemically structured design space of over 2.1 million Ru(II), Os(II), and Ir(III) complexes with targeted DFT calculations and pretrained atomistic representations. With only 300 quantum-chemical evaluations, the approach efficiently enriches candidates within a mechanistically defined optimal redox region. Analysis of the viable complexes reveals chemical design principles linking metal identity, ligand framework, substituent pattern, and physicochemical properties to Type~I photoreactivity, including a pronounced preference for Os(II)-based complexes and electronically asymmetric ligand environments along with combination of electronic donating and accepting substituents. More broadly, the strategy presented herein provides a scalable, mechanism-guided route for the rational design of transition-metal photocatalysts for applications spanning biomedicine, solar energy conversion, and photoredox chemistry.
Electromagnetic (EM) coupling between molecular exciton and plasmon has been studied using in Rayleigh scattering or extinction spectroscopy. However, evaluating EM coupling involving subradiant plasmon is challenging because this resonance does not manifest clearly in far-field spectra. In this study, we developed a method to evaluate such coupling using EM enhancement factors (FR) derived from ultrafast surface-enhanced fluorescence (ultrafast SEF). This SEF, which appears as a broad background in surface-enhanced resonant Raman scattering (SERRS) spectra, were measured using silver nanoparticle dimers containing dye molecules within their nanogaps. Our results show that the spectral peaks of FR for subradiant resonances appear near the dips in Rayleigh scattering spectra. Furthermore, these FR peaks exhibit blue-shifts during the quenching processes of both ultrafast SEF and SERRS. We examined these static and temporal spectral properties using a coupled oscillator model composed of radiant plasmons, subradiant plasmons, and molecular excitons. The static properties were reproduced by increasing the linewidths of the radiant plasmon resonance, while the temporal properties were captured by decreasing the EM coupling energies between the exciton and both plasmon oscillators. These findings indicate that this methodology is a powerful tool for evaluating EM coupling between subradiant plasmons and molecular excitons.
Understanding the complex anisotropic acoustic propagation in crystals is crucial for optimizing the performance of surface and bulk acoustic wave devices. Here, we investigate the anisotropy and coupling of GHz acoustic modes in (001)-cut gallium arsenide through theory and experiment. We first numerically calculate the angle-dependent phase velocities for surface and bulk modes, and we provide a code which can easily be adapted to different material systems. We validate our theoretical model experimentally by exciting surface modes with an interdigital transducer, and achieve omnidirectional acoustic propagation through random scattering of the acoustic waves. We measure the complex acoustic field with a scanning optical interferometer, and extract the angle-dependent velocities of surface and bulk modes using Fourier domain analysis. Our method could be used for the optimization of GHz-range classical and quantum acoustic devices, by studying losses of surface and bulk modes.
Understanding the formation of social ties requires disentangling the roles of individual traits and local network structure. We analyse signed social relationships among 3,395 students using an interpretable machine learning model -- the Explainable Boosting Machine (EBM) -- to predict link polarity from individual attributes (prosociality, cognitive reflection, and gender) and a structural metric, triadic influence. Our results show that triadic influence overwhelmingly dominates link prediction, confirming that local network structure is the primary driver of social relationships. Nevertheless, a small subset of links (0.24\%) is primarily explained by individual-level traits. A detailed characterisation of this subset reveals that these links do not arise from distinct structural conditions, but rather correspond to weaker and less structurally embedded relationships. In particular, they are more likely to be negative ties and exhibit lower levels of structural balance, whereas triadic-dominant links are strongly associated with positive relationships and highly balanced configurations. Furthermore, we find that links without indirect structural paths are not explained by individual traits, but by the absence of structural reinforcement itself. These findings support a layered view of social tie formation, in which structural mechanisms dominate globally, while individual-level effects emerge in specific, less constrained contexts. More broadly, our work highlights the value of explainable machine learning for uncovering the mechanisms underlying social network formation.