2025-04-07 | | Total: 12
Scale-invariance is a ubiquitous observation in the dynamics of large distributed complex systems. The computation of its scaling exponents, which provide clues on its origin, is often hampered by the limited available sampling data, making an appropriate mathematical description a challenge. This work investigates the behavior of correlation functions in fractal systems under conditions of severe subsampling. Analytical and numerical results reveal a striking robustness: the correlation functions continue to capture the expected scaling exponents despite substantial data reduction. This behavior is demonstrated numerically for the random 2-D Cantor set and the Sierpinski gasket, both consistent with exact analytical predictions. Similar robustness is observed in 1-D time series both synthetic and experimental, as well as in high-resolution images of a neuronal structure. Overall, these findings are broadly relevant for the structural characterization of biological systems under realistic sampling constraints.
Accurately predicting nonlinear transient thermal fields in two-dimensional domains is a significant challenge in various engineering fields, where conventional analytical and numerical methods struggle to balance physical fidelity with computational efficiency when dealing with strong material nonlinearities and evolving multiphysics boundary conditions. To address this challenge, we propose a novel cross-disciplinary approach integrating Green's function formulations with adaptive neural operators, enabling a new paradigm for multiphysics thermal analysis. Our methodology combines rigorous analytical derivations with a physics-informed neural architecture consisting of five adaptive hidden layers (64 neurons per layer) that incorporates solutions as physical constraints, optimizing learning rates to balance convergence stability and computational speed. Extensive validation demonstrates superior performance in handling rapid thermal transients and strongly coupled nonlinear responses, which significantly improves computational efficiency while maintaining high agreement with analytical benchmarks across a range of material configurations and boundary conditions.
Physical systems with complex unsteady dynamics, such as fluid flows, are often poorly represented by a single mean solution. For many practical applications, it is crucial to access the full distribution of possible states, from which relevant statistics (e.g., RMS and two-point correlations) can be derived. Here, we propose a graph-based latent diffusion (or alternatively, flow-matching) model that enables direct sampling of states from their equilibrium distribution, given a mesh discretization of the system and its physical parameters. This allows for the efficient computation of flow statistics without running long and expensive numerical simulations. The graph-based structure enables operations on unstructured meshes, which is critical for representing complex geometries with spatially localized high gradients, while latent-space diffusion modeling with a multi-scale GNN allows for efficient learning and inference of entire distributions of solutions. A key finding is that the proposed networks can accurately learn full distributions even when trained on incomplete data from relatively short simulations. We apply this method to a range of fluid dynamics tasks, such as predicting pressure distributions on 3D wing models in turbulent flow, demonstrating both accuracy and computational efficiency in challenging scenarios. The ability to directly sample accurate solutions, and capturing their diversity from short ground-truth simulations, is highly promising for complex scientific modeling tasks.
The supersolid phase of soft-core bosons in two dimensions is investigated using the self-consistent Hartree-Fock and quantum Monte Carlo methods. An approximate phase diagram at finite temperatures is initially constructed using the mean-field approach, which is subsequently validated through precise path-integral simulations, enabling a microscopic characterization of the various phases. Superfluid and melting/freezing transitions are analyzed through the superfluid density and the long-range behavior of correlation functions associated with positional and orientational order, in accordance with the general picture of Berezinskii-Kosterlitz-Thouless transitions. A broad region at low temperatures is identified where the supersolid phase exists, separating the uniform superfluid phase from the normal quasi-crystal phase. Additionally, a potential intermediate hexatic phase with quasi long-range orientational order is identified in a narrow region between the normal solid and fluid phases. These findings establish self-consistent Hartree-Fock theory beyond the local density approximation as an effective tool, complementary to computationally intensive quantum Monte Carlo simulations, for investigating the melting of the supersolid phase and the possible emergence of the hexatic superfluid phase in bosonic systems with various interaction potentials.
Linear-response time-dependent Density Functional Theory (LR-TDDFT) is a widely used method for accurately predicting the excited-state properties of physical systems. Previous works have attempted to accelerate LR-TDDFT using heterogeneous systems such as GPUs, FPGAs, and the Sunway architecture. However, a major drawback of these approaches is the constant data movement between host memory and the memory of the heterogeneous systems, which results in substantial \textit{data movement overhead}. Moreover, these works focus primarily on optimizing the compute-intensive portions of LR-TDDFT, despite the fact that the calculation steps are fundamentally \textit{memory-bound}. To address these challenges, we propose NDFT, a \underline{N}ear-\underline{D}ata Density \underline{F}unctional \underline{T}heory framework. Specifically, we design a novel task partitioning and scheduling mechanism to offload each part of LR-TDDFT to the most suitable computing units within a CPU-NDP system. Additionally, we implement a hardware/software co-optimization of a critical kernel in LR-TDDFT to further enhance performance on the CPU-NDP system. Our results show that NDFT achieves performance improvements of 5.2x and 2.5x over CPU and GPU baselines, respectively, on a large physical system.
A fully integrated waveguide-based, efficient surface plasmon coupler composed of a realistic non-tapered dielectric waveguide with graphene patches and sheet is designed and optimized for the infrared. The coupling efficiency can reach nearly 80\% for a coupler as short as 700 nm for an operating wavelength of 12 μm. This \gui{work} is carried out \gui{using rigorous numerical models} based on the finite element method \gui{taking into account 2D-materials as surface conductivities}. \gui{The} key numerical results are \gui{supported by} phdisplaycopyrightysical arguments based on modal approach or resonance condition.
Isotopically enriched silicon (Si) can greatly enhance qubit coherence times by minimizing naturally occurring 29Si which has a non-zero nuclear spin. Ultra-high fluence 28Si ion implantation of bulk natural Si substrates was recently demonstrated as an attractive technique to ultra-high 28Si isotopic purity. In this work, we apply this 28Si enrichment process to produce 28Si and 28Si-on-insulator (SOI) samples. Experimentally, we produced a 28Si sample on natural Si substrate with 29Si depleted to 7~ppm (limited by measurement noise floor), that is at least 100 nm thick. This is achieved with an ion energy that results in a sputter yield of less than one and a high ion fluence, as supported by simulations. Further, our simulations predict the 29Si and 30Si depletion in our sample to be less than 1~ppm. In the case of SOI, ion implantation conditions are found to be more stringent than those of bulk natural Si in terms of minimizing threading dislocations upon subsequent solid phase epitaxy annealing. Finally, we report the observation of nanoscopic voids in our 28SOI and 28Si samples located in the depth region between the surface and 70~nm. These voids appear to be stabilized by the presence of gold impurities and can be annealed out completely by a standard SPE annealing protocol at 620\deg C for 10 minutes in the absence of gold.
Converging flows are visible around bipolar magnetic regions (BMRs) on the solar surface, according to observations. Average flows are created by these inflows combined, and the strength of these flows depends on the amount of flux present during the solar cycle. In models of the solar cycle, this average flow can be depicted as perturbations to the meridional flow. In this article, we study the effects of introducing surface inflow to the surface flux transport models (SFT) as a possible nonlinear mechanism in the presence of latitude quenching for an inflow profile whose amplitude varies within a cycle depending on the magnetic activity. The results show that including surface inflows in the model in the presence of both LQ and tilt quenching (TQ) produced a polar field within a ±1σ of an average cycle polar field (σ is the standard deviation) and a correlation coefficient of 0.85. We confirm that including inflows produces a lower net contribution to the dipole moment (10\,--\,25\%). Furthermore, the relative importance of LQ vs. inflows is inversely correlated with the dynamo effectivity range (λR). With no decay term, introducing inflows into the model resulted in a less significant net contribution to the dipole moment. Including inflows in the SFT model shows a possible nonlinear relationship between the surface inflows and the solar dipole moment, suggesting a potential nonlinear mechanism contributing to the saturation of the global dynamo. For lower λR (⪅ 10 ^\circ), TQ always dominates LQ, and for higher \lambda_R LQ dominate. However, including inflows will make the domination a little bit earlier in case of having a decay term in the model.
Understanding how residue variations affect protein stability is crucial for designing functional proteins and deciphering the molecular mechanisms underlying disease-related mutations. Recent advances in protein language models (PLMs) have revolutionized computational protein analysis, enabling, among other things, more accurate predictions of mutational effects. In this work, we introduce JanusDDG, a deep learning framework that leverages PLM-derived embeddings and a bidirectional cross-attention transformer architecture to predict \Delta \Delta G of single and multiple-residue mutations while simultaneously being constrained to respect fundamental thermodynamic properties, such as antisymmetry and transitivity. Unlike conventional self-attention, JanusDDG computes queries (Q) and values (V) as the difference between wild-type and mutant embeddings, while keys (K) alternate between the two. This cross-interleaved attention mechanism enables the model to capture mutation-induced perturbations while preserving essential contextual information. Experimental results show that JanusDDG achieves state-of-the-art performance in predicting \Delta \Delta G from sequence alone, matching or exceeding the accuracy of structure-based methods for both single and multiple mutations.
We describe an experiment using superconducting transmon qubits that demonstrates wavefunction collapse consistent with Orchestrated Objective Reduction (Orch-OR) - the theory of consciousness proposed by Sir Roger Penrose and Stuart Hameroff. The experiment performs a partial measurement on a qubit. It uses the result of that measurement to move an estimated 10e-12 kg of mass in one of two locations separated by approximately 1 mm. In standard quantum mechanics, the partial measurement leaves the system in an improper mixture, a state that appears probabilistic but retains quantum coherence. This is mathematically indistinguishable from a proper mixture, where the state has genuinely collapsed. According to Penrose's theory, improper mixtures can lead to gravitationally induced collapse. Results of our experiment show a change in the circuit evolution that is consistent with wavefunction collapse driven by such an improper mixture. The experiment is performed on an IBM Eagle 127-qubit processor using IBM's programming framework Qiskit.
Simulation-to-Real (Sim2Real) transfer learning, the machine learning technique that efficiently solves a real-world task by leveraging knowledge from computational data, has received increasing attention in materials science as a promising solution to the scarcity of experimental data. We proposed an efficient transfer learning scheme from first-principles calculations to experiments based on the chemistry-informed domain transformation, that integrates the heterogeneous source and target domains by harnessing the underlying physics and chemistry. The proposed method maps the computational data from the simulation space (source domain) into the space of experimental data (target domain). During this process, these qualitatively different domains are efficiently bridged by prior knowledge of chemistry, the statistical ensemble and the relationship between source and target quantities. As a proof-of-concept, we predict the catalyst activity for the reverse water-gas shift reaction by using the abundant first-principles data in addition to the experimental data. Through the demonstration, we confirmed that the transfer learning model exhibits positive transfer in accuracy and data efficiency. In particular, a significantly high accuracy was achieved despite using a few (less than ten) target data in domain transformation, whose accuracy is one order of magnitude smaller than that of a full scratch model trained with over 100 target data. This result indicates that the proposed method leverages the high prediction performance with few target data, which helps to save the number of trials in real laboratories.
Large Language Models (LLMs) are playing an expanding role in physics research by enhancing reasoning, symbolic manipulation, and numerical computation. However, ensuring the reliability and interpretability of their outputs remains a significant challenge. In our framework, we conceptualize the collaboration between AI and human scientists as a dynamic interplay among three modules: the reasoning module, the interpretation module, and the AI-scientist interaction module. Recognizing that effective physics reasoning demands rigorous logical consistency, quantitative precision, and deep integration with established theoretical models, we introduce the interpretation module to improve the understanding of AI-generated outputs, which is not previously explored in the literature. This module comprises multiple specialized agents, including summarizers, model builders, UI builders, and testers, which collaboratively structure LLM outputs within a physically grounded framework, by constructing a more interpretable science model. A case study demonstrates that our approach enhances transparency, facilitates validation, and strengthens AI-augmented reasoning in scientific discovery.