2025-06-19 | | Total: 83
We study gravitational waves from a stellar-mass binary orbiting a spinning supermassive black hole, a system referred to as a binary extreme mass ratio inspiral (b-EMRI). We use Dixon's formalism to describe the stellar-mass binary as a particle with internal structure, and keep terms up to quadrupole order to capture the generation of gravitational waves by the inner motion of the stellar-mass binary. The problem of emission and propagation of waves is treated from first principles using black hole perturbation theory. In the gravitational waveform at future null infinity, we identify for the first time Doppler shifts and beaming due to the motion of the center of mass, as well as helicity breaking gravitational lensing, and resonances with ringdown modes of the supermassive black hole. We establish that previously proposed phenomenological models inadequately capture these effects.
Neural network potentials (NNPs) offer a powerful alternative to traditional force fields for molecular dynamics (MD) simulations. Accurate and stable MD simulations, crucial for evaluating material properties, require training data encompassing both low-energy stable structures and high-energy structures. Conventional knowledge distillation (KD) methods fine-tune a pre-trained NNP as a teacher model to generate training data for a student model. However, in material-specific models, this fine-tuning process increases energy barriers, making it difficult to create training data containing high-energy structures. To address this, we propose a novel KD framework that leverages a non-fine-tuned, off-the-shelf pre-trained NNP as a teacher. Its gentler energy landscape facilitates the exploration of a wider range of structures, including the high-energy structures crucial for stable MD simulations. Our framework employs a two-stage training process: first, the student NNP is trained with a dataset generated by the off-the-shelf teacher; then, it is fine-tuned with a smaller, high-accuracy density functional theory (DFT) dataset. We demonstrate the effectiveness of our framework by applying it to both organic (polyethylene glycol) and inorganic (L$_{10}$GeP$_{2}$S$_{12}$) materials, achieving comparable or superior accuracy in reproducing physical properties compared to existing methods. Importantly, our method reduces the number of expensive DFT calculations by 10x compared to existing NNP generation methods, without sacrificing accuracy.
Determining onflow parameters is crucial from the perspectives of wind tunnel testing and regular flight and wind turbine operations. These parameters have traditionally been predicted via direct measurements which might lead to challenges in case of sensor faults. Alternatively, a data-driven prediction model based on surface pressure data can be used to determine these parameters. It is essential that such predictors achieve close to real-time learning as dictated by practical applications such as monitoring wind tunnel operations or learning the variations in aerodynamic performance of aerospace and wind energy systems. To overcome the challenges caused by changes in the data distribution as well as in adapting to a new prediction task, we propose a transfer learning methodology to predict the onflow parameters, specifically angle of attack and onflow speed. It requires first training a convolutional neural network (ConvNet) model offline for the core prediction task, then freezing the weights of this model except the selected layers preceding the output node, and finally executing transfer learning by retraining these layers. A demonstration of this approach is provided using steady CFD analysis data for an airfoil for i) domain adaptation where transfer learning is performed with data from a target domain having different data distribution than the source domain and ii) task adaptation where the prediction task is changed. Further exploration on the influence of noisy data, performance on an extended domain, and trade studies varying sampling sizes and architectures are provided. Results successfully demonstrate the potential of the approach for adaptation to changing data distribution, domain extension, and task update while the application for noisy data is concluded to be not as effective.
How can we detect when global events fundamentally reshape public discourse? This study introduces a topological framework for identifying structural change in media narratives using persistent homology. Drawing on international news articles surrounding major events - including the Russian invasion of Ukraine (Feb 2022), the murder of George Floyd (May 2020), the U.S. Capitol insurrection (Jan 2021), and the Hamas-led invasion of Israel (Oct 2023) - we construct daily co-occurrence graphs of noun phrases to trace evolving discourse. Each graph is embedded and transformed into a persistence diagram via a Vietoris-Rips filtration. We then compute Wasserstein distances and persistence entropies across homological dimensions to capture semantic disruption and narrative volatility over time. Our results show that major geopolitical and social events align with sharp spikes in both H0 (connected components) and H1 (loops), indicating sudden reorganization in narrative structure and coherence. Cross-correlation analyses reveal a typical lag pattern in which changes to component-level structure (H0) precede higher-order motif shifts (H1), suggesting a bottom-up cascade of semantic change. An exception occurs during the Russian invasion of Ukraine, where H1 entropy leads H0, possibly reflecting top-down narrative framing before local discourse adjusts. Persistence entropy further distinguishes tightly focused from diffuse narrative regimes. These findings demonstrate that persistent homology offers a mathematically principled, unsupervised method for detecting inflection points and directional shifts in public attention - without requiring prior knowledge of specific events. This topological approach advances computational social science by enabling real-time detection of semantic restructuring during crises, protests, and information shocks.
We have successfully fabricated the nitrogen-doped cuprous oxide thin films on the amorphous standard soda-lime glass by reactive high-power impulse magnetron sputtering. The energy of film-forming particles was controlled by the value of pulse-averaged target power density, which has a significant impact on the elemental composition, structure and optoelectrical properties of the films. We have shown that the high-energy regime is more suitable for preserving Cu2O structure and leads to continuous substitution of oxygen by nitrogen compared with the low-energy regime. Moreover, in the high-energy regime, it is possible, to some extent, to independently control the electrical resistivity and optical properties. The electrical resistivity decreases down to 5 x10-2 ohm.cm at the optical band gap 2.0-2.3 eV. Special attention is paid to the formation of nitrogen molecules and their ability to form shallow acceptor states. Experimental results supported by our DFT calculations indicate that N2 replacing Cu in the Cu2O lattice is one possible (but not the only possible) acceptor. We have also found that the formation of nitrogen molecules is preferred in a high-energy regime.
In classical probability theory, the convergence of empirical frequencies to theoretical probabilities: as captured by the Law of Large Numbers (LLN): is treated as axiomatic and emergent from statistical assumptions such as independence and identical distribution. In this work, a novel dynamical framework is constructed in which convergence arises as a consequence of structured evolution in outcome space, rather than a statistical postulate. Through this formalism, statistical convergence is derived dynamically, revealing an internal structure to randomness and exposing entanglement between successive trials. The system recovers classical LLN behaviour in the large-number limit, while predicting deviations, transient fluctuations, and geometric asymmetries in the early regime. This work inaugurates a new paradigm: dynamical probability mechanics: in which randomness is modelled not as a sequence of disconnected stochastic events, but as a physically structured, feedback-driven process. The theory provides a novel explanatory layer beneath statistical laws and opens pathways toward a mechanistic foundation of probability itself.
The computational discovery and design of new crystalline materials, particularly metal-organic frameworks (MOFs), heavily relies on high-quality, computation-ready structural data. However, recent studies have revealed significant error rates within existing MOF databases, posing a critical data problem that hinders efficient high-throughput computational screening. While rule-based algorithms like MOSAEC, MOFChecker, and the Chen and Manz method (Chen-Manz) have been developed to address this, they often suffer from inherent limitations and misclassification of structures. To overcome this challenge, we introduce MOFClassifier, a novel machine learning approach built upon a positive-unlabeled crystal graph convolutional neural network (PU-CGCNN) model. MOFClassifier learns intricate patterns from perfect crys-tal structures to predict a crystal-likeness score (CLscore), effectively classifying MOFs as computation-ready. Our model achieves a ROC value of 0.979 (previous best 0.912) and, importantly, can identify subtle structural and chemical errors that are fundamentally undetectable by current rule-based methods. By accurately recovering previously misclassified false-negative structures, MOFClassifier reduces the risk of overlooking promising material candidates in large-scale computational screening efforts. This user-friendly tool is freely available and has been integrated into the preparation workflow for the updated CoRE MOF DB 2025 v1, contributing to accelerated computational discovery of MOF materials.
The shallow water equations often assume a constant velocity profile along the vertical axis. However, this assumption does not hold in many practical applications. To better approximate the vertical velocity distribution, models such as the shallow water moment expansion models have been proposed. Nevertheless, under non-slip bottom boundary conditions, both the standard shallow water equation and its moment-enhanced models struggle to accurately capture the vertical velocity profile due to the stiff source terms. In this work, we propose modified shallow water equations and corresponding moment-enhanced models that perform well under both non-slip and slip boundary conditions. The primary difference between the modified and original models lies in the treatment of the source term, which allows our modified moment expansion models to be readily generalized, while maintaining compatibility with our previous analysis on the hyperbolicity of the model. To assess the performance of both the standard and modified moment expansion models, we conduct a comprehensive numerical comparison with the incompressible Navier--Stokes equations -- a comparison that is absent from existing literature.
We present a general and automated approach for computing model gradients for PDE solvers built on sparse spectral methods, and implement this capability in the widely used open-source Dedalus framework. We apply reverse-mode automatic differentiation to symbolic graph representations of PDEs, efficiently constructing adjoint solvers that retain the speed and memory efficiency of this important class of modern numerical methods. This approach enables users to compute gradients and perform optimization for a wide range of time-dependent and nonlinear systems without writing additional code. The framework supports a broad class of equations, geometries, and boundary conditions, and runs efficiently in parallel using MPI. We demonstrate the flexibility and capabilities of this system using canonical problems from the literature, showing both strong performance and practical utility for a wide variety of inverse problems. By integrating automatic adjoints into a flexible high-level solver, our approach enables researchers to perform gradient-based optimization and sensitivity analyses in spectral simulations with ease and efficiency.
Optical observations of solar corona provide key information on its magnetic geometry. The large-scale open field of the corona plays an important role in shaping the ambient solar wind and constraining the propagation dynamics of the embedded structures, such as interplanetary coronal mass ejections. Rigorous analysis of the open-flux coronal regions based on coronagraph images can be quite challenging because of the depleted plasma density resulting in low signal-to-noise ratios. In this paper, we present an in-depth description of a new image segmentation methodology, the Quasi-Radial Field-line Tracing (QRaFT), enabling a detection of field-aligned optical coronal features approximating the orientation of the steady-state open magnetic field. The methodology is tested using synthetic coronagraph images generated by a three-dimensional magnetohydrodynamic model. The results of the numerical tests indicate that the extracted optical features are aligned within $\sim 4-7$ degrees with the local magnetic field in the underlying numerical solution. We also demonstrate the performance of the method on real-life coronal images obtained from a space-borne coronagraph and a ground-based camera. We argue that QRaFT outputs contain valuable empirical information about the global steady-state morphology of the corona which could help improving the accuracy of coronal and solar wind models and space weather forecasts.
We apply machine learning methods to demonstrate range superresolution in remote sensing radar detection. Specifically, we implement a denoising autoencoder to estimate the distance between two equal intensity scatterers in the subwavelength regime. The machine learning models are trained on waveforms subject to a bandlimit constraint such that ranges much smaller than the inverse bandlimit are optimized in their precision. The autoencoder achieves effective dimensionality reduction, with the bottleneck layer exhibiting a strong and consistent correlation with the true scatterer separation. We confirm reproducibility across different training sessions and network initializations by analyzing the scaled encoder outputs and their robustness to noise. We investigate the behavior of the bottleneck layer for the following types of pulses: a traditional sinc pulse, a bandlimited triangle-type pulse, and a theoretically near-optimal pulse created from a spherical Bessel function basis. The Bessel signal performs best, followed by the triangle wave, with the sinc signal performing worst, highlighting the crucial role of signal design in the success of machine-learning-based range resolution.
There exist elegant methods of aligning point clouds in $\mathbb R^3$. Unfortunately, these methods rely on the positive definite property of the Euclidean metric, and do not easily extend to the indefinite Minkowski metric. In this paper, we propose two solutions to the following problem: given inertial reference frames $A$ and $B$, and given (possibly noisy) measurements of a set of 4-vectors $\{v_i\}$ made in those reference frames with components $\{v_{A,i}\}$ and $\{v_{B,i}\}$, find the optimal Lorentz transformation $\Lambda$ such that $\Lambda v_{A,i}=v_{B,i}$. The method we outline is conceptually simple and easily extends to alignment problems in other matrix Lie groups.
At room temperature, low frequency vibrations at far-infrared frequencies are thermally excited ($k_B T > h \nu$) and not restricted to harmonic fluctuations around a single potential energy minimum. For folded proteins, these intrinsically anharmonic vibrations can contain information on slow conformational transitions. Recently, we have developed FREquency-SElective ANharmonic (FRESEAN) mode analysis, a method based on time correlation functions that isolates low-frequency vibrational motions from molecular dynamics simulation trajectories without relying on harmonic approximations. We recently showed that low-frequency vibrations obtained from FRESEAN mode analysis are effective collective variables in enhanced sampling simulations of conformational ensembles. However, FRESEAN mode analysis is based on velocity time correlations between all degrees of freedom, which creates computational challenges for large biomolecules. To facilitate future applications, we demonstrate here how coarse-graining of all-atom simulation trajectories can be combined with FRESEAN mode analysis to extract information on low-frequency vibrations at minimal computational cost.
Nearly a decade ago it was discovered that the spherical cell body of the alga $Chlamydomonas~reinhardtii$ can act as a lens to concentrate incoming light onto the cell's membrane-bound photoreceptor and thereby affect phototaxis. Since many nearly transparent cells in marine environments have complex, often non-axisymmetric shapes, this observation raises fundamental, yet little-explored questions in biological optics about light refraction by the bodies of microorganisms. There are two distinct contexts for such questions: the $absorption$ problem for $incoming$ light, typified by photosynthetic activity taking place in the chloroplasts of green algae, and the $emission$ problem for $outgoing$ light, where the paradigm is bioluminescence emitted from scintillons within dinoflagellates. Here we examine both of these aspects of ``algal optics" in the special case where the absorption or emission is localized in structures that are small relative to the overall organism size, taking into account both refraction and reflections at the cell-water boundary. Analytical and numerical results are developed for the distribution of light intensities inside and outside the body, and we establish certain duality relationships that connect the incoming and outgoing problems. For strongly non-spherical shapes we find lensing effects that may have implications for photosynthetic activity and for the angular distribution of light emitted during bioluminescent flashes.
The issue of gender bias in scientific publications has been the subject of ongoing debate. One aspect of this debate concerns whether women receive equal credit for their contributions compared to men. Conventional wisdom suggests that women are more likely to be acknowledged than listed as co-authors, a role that carries greater prestige. Here, we analyze data from hundreds of thousands of scientists across nine disciplines and a broad range of publications. Our results confirm persistent gender disparities: women are more frequently acknowledged than credited as co-authors, especially in roles involving investigation and analysis. To account for status and disciplinary effects, we examined collaboration pairs composed of highly cited (high-status) and less cited (low-status) scientists. In such collaborations, the highly cited scientist is more likely to be listed as a co-author, regardless of gender. Notably, highly cited women in these pairs are even more likely to be listed as co-authors than their male counterparts. These findings suggest that power dynamics and perceived success heavily influence how credit is distributed in scientific publishing. The results underscore the role of status in shaping authorship and call for a more nuanced understanding of how gender, power, and recognition interact. This research offers valuable insights for scientists, editors, and funding agencies committed to advancing equity in science.
While cryogenic-temperature subthreshold swing (SS) in crystalline semiconductors has been widely studied, a careful study on the temperature-dependent SS in amorphous oxide semiconductors remains lacking. In this paper, a comprehensive analysis of the SS in thin-film transistors with an amorphous indium gallium zinc oxide (IGZO) channel at temperatures from 300 K down to 4 K is presented. Main observations include: 1) At room temperature (300 K), the devices exhibit a SS of 61 mV/dec, and a low interface trap density (<1011 cm-2), among the best reported values for IGZO devices. 2) A SS saturation around 40 mV/dec is observed between 200 K and 100 K. It is well explained by the electron transport via band tail states with exponential decay (Wt) of 13 meV. 3) At deep-cryogenic temperature, SS increase significantly exceeding 200 mV/dec at 4 K. Such high SS values are actually limited by the measurement current range, confirmed by Id-Vg simulations based on the variable range hopping (VRH) model. This work not only elucidates the SS behavior in amorphous IGZO devices but also provides a deep understanding of the physical mechanisms of electron transport in amorphous semiconductors.
Quantum neural networks (QNNs) are widely employed as ansätze for solving variational problems, where their expressivity directly impacts performance. Yet, accurately characterizing QNN expressivity remains an open challenge, impeding the optimal design of quantum circuits. In this work, we introduce the effective rank, denoted as $\kappa$, as a novel quantitative measure of expressivity. Specifically, $\kappa$ captures the number of effectively independent parameters among all the variational parameters in a parameterized quantum circuit, thus reflecting the true degrees of freedom contributing to expressivity. Through a systematic analysis considering circuit architecture, input data distributions, and measurement protocols, we demonstrate that $\kappa$ can saturate its theoretical upper bound, $d_n=4^n-1$, for an $n$-qubit system when each of the three factors is optimally expressive. This result provides a rigorous framework for assessing QNN expressivity and quantifying their functional capacity. Building on these theoretical insights, and motivated by the vast and highly structured nature of the circuit design space, we employ $\kappa$ as a guiding metric for the automated design of highly expressive quantum circuit configurations. To this end, we develop a reinforcement learning framework featuring a self-attention transformer agent that autonomously explores and optimizes circuit architectures. By integrating theoretical characterization with practical optimization, our work establishes $\kappa$ as a robust tool for quantifying QNN expressivity and demonstrates the effectiveness of reinforcement learning in designing high-performance quantum circuits. This study paves the way for building more expressive QNN architectures, ultimately enhancing the capabilities of quantum machine learning.
The generation of the ground state of a spin vortex in a $^{87}$Rb Bose-Einstein condensate with the assistance of an optical plug has been studied. However, gravity is everywhere, and this potential linear dependence on the spatial position will destroy the axisymmetric structure of the system with the optical plug. In this case, the question of whether the spin vortex ground state still exists remains unresolved. The present study aims to explore the impact of microgravity on the formation of the spin vortex state with the assistance of an optical plug. To this end, a simple model has been employed to provide a comprehensive understanding of the phenomenon. The Gross-Pitaevskii equations are solved by setting the optical plug intensity, adjusting the optical plug width, and adjusting the microgravity strength. This process results in the phase diagram for the single-mode state and spin vortex state. Under microgravity situations, we observe an off-axis structure of the spin vortex state. Our calculations offer a reliable approach to generating spin vortex states in a microgravity environment.
Physics-informed neural networks (PINNs) are extensively used to represent various physical systems across multiple scientific domains. The same can be said for cardiac electrophysiology, wherein fully-connected neural networks (FCNNs) have been employed to predict the evolution of an action potential in a 2D space following the two-parameter phenomenological Aliev-Panfilov (AP) model. In this paper, the training behaviour of PINNs is investigated to determine optimal hyperparameters to predict the electrophysiological activity of the myocardium in 3D according to the AP model, with the inclusion of boundary and material parameters. An FCNN architecture is employed with the governing partial differential equations in their strong form, which are scaled consistently with normalization of network inputs. The finite element (FE) method is used to generate training data for the network. Numerical examples with varying spatial dimensions and parameterizations are generated using the trained models. The network predicted fields for both the action potential and the recovery variable are compared with the respective FE simulations. Network losses are weighed with individual scalar values. Their effect on training and prediction is studied to arrive at a method of controlling losses during training.
Shear-thickening fluids (STFs) become more viscous under shear stress, which makes them useful for many engineering and scientific applications. However, their behavior under normal forces, especially when these forces are applied harmonically, is less understood. Here, we examine the dynamic response of STFs under harmonic excitation. Our experimental setup features an unbalanced rotor attached to a vibrating plate submerged in an STF-filled container. We monitored the rotor's speed, the system's displacement, and the STF force. Data from experiments without STF were used to identify system parameters, and measurements of the STF force revealed the dynamic nature of the STF force. By comparing responses with and without STF, we identified three regions of behavior: in the pre-resonance region, the STF force is negligible; at resonance, it acts as a significant damping force; and in the post-resonance region, it behaves like an on-off force. Overall, the STF effectively reduces resonance amplitudes. These results can inform the design of complex structures incorporating STFs.
We combine classical heuristics with partial shadow tomography to enable efficient protocols for extracting information from correlated ab initio electronic systems encoded on quantum devices. By proposing the use of a correlation energy functional and sampling of a polynomial set of excitation amplitudes of the quantum state, we can demonstrate an almost two order of magnitude reduction in required number of shots for a given statistical error in the energy estimate, as well as observing a linear scaling to accessible system sizes. Furthermore, we find a high-degree of noise resilience of these estimators on real quantum devices, with up to an order of magnitude increase in the tolerated noise compared to traditional techniques. While these approaches are expected to break down asymptotically, we find strong evidence that these large system arguments do not prevent algorithmic advantage from these simple protocols in many systems of interest. We further extend this to consider the extraction of beyond-energetic properties by mapping to a coupled cluster surrogate model, as well as a natural combination within a quantum embedding framework. This embedding framework avoids the unstable self-consistent requirements of previous approaches, enabling application of quantum solvers to realistic correlated materials science, where we demonstrate the volume-dependence of the spin gap of Nickel Oxide.
Quantum key distribution (QKD) with deterministic single photon sources has been demonstrated over intercity fiber and free-space channels. The previous implementations relied mainly on polarization encoding schemes, which are susceptible to birefringence, polarization-mode dispersion and polarization-dependent loss in practical fiber networks. In contrast, time-bin encoding offers inherent robustness and has been widely adopted in mature QKD systems using weak coherent laser pulses. However, its feasibility in conjunction with a deterministic single-photon source has not yet been experimentally demonstrated. In this work, we construct a time-bin encoded QKD system employing a high-brightness quantum dot (QD) single-photon source operating at telecom wavelength. Our proof-of-concept experiment successfully demonstrates the possibility of secure key distribution over fiber link of 120 km, while maintaining extraordinary long-term stability over 6 hours of continuous operation. This work provides the first experimental validation of integrating a quantum dot single-photon source with time-bin encoding in a telecom-band QKD system. In addition, it demonstrates the highest secure key rate among the time-bin QKDs based on single-photon sources. This development signifies a substantial advancement in the establishment of a robust and scalable QKD network based on solid-state single-photon technology
Equilibrium and nonequilibrium states of matter can exhibit fundamentally different behavior. A key example is the Kardar-Parisi-Zhang universality class in two spatial dimensions (2D KPZ), where microscopic deviations from equilibrium give rise to macroscopic scaling laws without equilibrium counterparts. While extensively studied theoretically, direct experimental evidence of 2D KPZ scaling has remained limited to interface growth so far. Here, we report the observation of universal scaling consistent with the KPZ universality class in 2D exciton-polariton condensates -- quantum fluids of light that are inherently driven and dissipative, thus breaking equilibrium conditions. Using momentum-resolved photoluminescence spectroscopy as well as space- and time-resolved interferometry, we probe the phase correlations across microscopically different systems, varying drive conditions in two distinct lattice geometries. Our analysis reveals correlation dynamics and scaling exponents in excellent agreement with 2D KPZ predictions. These results establish exciton-polariton condensates as a robust experimental platform for exploring 2D nonequilibrium universality quantitatively, and open new avenues for investigating the emergence of coherence in interacting quantum systems far from equilibrium.
Neutral atom quantum processors are a promising platform for scalable quantum computing. An obstacle to implementing deep quantum circuits is managing atom loss, which constitutes a significant fraction of all errors. Current approaches are either not capable of replacing lost atoms in the middle of a circuit -- and therefore restricted to fixed, short circuit depths -- or require more than an order of magnitude longer time than gate and measurement operations to do so. In this work, we demonstrate fast, continuous atom replacement leveraging the metastable $^{171}$Yb qubit. A continuously loaded reservoir near the computation zone enables on-demand atom extraction with tweezers up to 500 times per second. New qubit arrays can be prepared 30 times per second when including single-atom preparation, non-destructive imaging and initialization. Importantly, existing qubits are completely undisturbed by the reloading process, owing to the extreme isolation of the metastable qubit from cooling and imaging light. This work establishes a complete foundation for implementing fast quantum circuits with unlimited depth, removing a final roadblock for fault-tolerant quantum computing with neutral atoms.
In recent decades, it has been emphasized that the evolving structure of networks may be shaped by interaction principles that yield sparse graphs with a vertex degree distribution exhibiting an algebraic tail, and other structural traits that are not featured in traditional random graphs. In this respect, through a mean-field approach, this review tackles the statistical physics of graph models based on the interaction principle of duplication-divergence. Additional sophistications extending the duplication-divergence model are also reviewed as well as generalizations of other known models. Possible research gaps and related prior results are then discussed.