2025-04-03 | | Total: 9
We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn-Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method, while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method's accuracy by comparing with planewave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 seconds on 4608 processors.
The inverse design of metamaterial architectures presents a significant challenge, particularly for nonlinear mechanical properties involving large deformations, buckling, contact, and plasticity. Traditional methods, such as gradient-based optimization, and recent generative deep-learning approaches often rely on binary pixel-based representations, which introduce jagged edges that hinder finite element (FE) simulations and 3D printing. To overcome these challenges, we propose an inverse design framework that utilizes a signed distance function (SDF) representation combined with a conditional diffusion model. The SDF provides a smooth boundary representation, eliminating the need for post-processing and ensuring compatibility with FE simulations and manufacturing methods. A classifier-free guided diffusion model is trained to generate SDFs conditioned on target macroscopic stress-strain curves, enabling efficient one-shot design synthesis. To assess the mechanical response of the generated designs, we introduce a forward prediction model based on Neural Operator Transformers (NOT), which accurately predicts homogenized stress-strain curves and local solution fields for arbitrary geometries with irregular query meshes. This approach enables a closed-loop process for general metamaterial design, offering a pathway for the development of advanced functional materials.
Understanding and solving fluid dynamics equations efficiently remains a fundamental challenge in computational physics. Traditional numerical solvers and physics-informed neural networks struggle to capture the full range of frequency components in partial differential equation solutions, limiting their accuracy and efficiency. Here, we propose the Multi-stream Physics Hybrid Network, a novel neural architecture that integrates quantum and classical layers in parallel to improve the accuracy of solving fluid dynamics equations, namely Kovasznay flow problem. This approach decomposes the solution into separate frequency components, each predicted by independent Parallel Hybrid Networks, simplifying the training process and enhancing performance. We evaluated the proposed model against a comparable classical neural network, the Multi-stream Physics Classical Network, in both data-driven and physics-driven scenarios. Our results show that the Multi-stream Physics Hybrid Network achieves a reduction in root mean square error by 36% for velocity components and 41% for pressure prediction compared to the classical model, while using 24% fewer trainable parameters. These findings highlight the potential of hybrid quantum-classical architectures for advancing computational fluid dynamics.
We present a novel approach for simulating acoustic (pressure) wave propagation across different media separated by a diffuse interface through the use of a weak compressibility formulation. Our method builds on our previous work on an entropy-stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes/Cahn-Hilliard system \cite{manzanero2020entropyNSCH}, and incorporates a modified weak compressibility formulation that allows different sound speeds in each phase. We validate our method through numerical experiments, demonstrating spectral convergence for acoustic transmission and reflection coefficients in one dimension and for the angle defined by Snell's law in two dimensions. Special attention is given to quantifying the modeling errors introduced by the width of the diffuse interface. Our results show that the method successfully captures the behavior of acoustic waves across interfaces, allowing exponential convergence in transmitted waves. The transmitted angles in two dimensions are accurately captured for air-water conditions, up to the critical angle of 13∘. This work represents a step forward in modeling acoustic propagation in incompressible multiphase systems, with potential applications to marine aeroacoustics.
With the aim of generalizing histogram statistics to higher dimensional cases, density estimation via discrepancy based sequential partition (DSP) has been proposed [D. Li, K. Yang, W. Wong, Advances in Neural Information Processing Systems (2016) 1099-1107] to learn an adaptive piecewise constant approximation defined on a binary sequential partition of the underlying domain, where the star discrepancy is adopted to measure the uniformity of particle distribution. However, the calculation of the star discrepancy is NP-hard and it does not satisfy the reflection invariance and rotation invariance either. To this end, we use the mixture discrepancy and the comparison of moments as a replacement of the star discrepancy, leading to the density estimation via mixture discrepancy based sequential partition (DSP-mix) and density estimation via moments based sequential partition (MSP), respectively. Both DSP-mix and MSP are computationally tractable and exhibit the reflection and rotation invariance. Numerical experiments in reconstructing the d-D mixture of Gaussians and Betas with d=2,3,…,6 demonstrate that DSP-mix and MSP both run approximately ten times faster than DSP while maintaining the same accuracy.
Contemplating to register signature of a new vector boson unambiguously from the measured non-linear isotope shift (IS) effects in three recent experiments [Phys. Rev. X {\bf 12}, 021033 (2022); Phys. Rev. Lett. {\bf 128}, 163201 (2022) and Phys. Rev. Lett. {\bf 134}, 063002 (2025)], very precise values of quadrupole moments and IS constants for the 6s 2S1/2→5d 2D3/2 and 6s 2S1/2→5d 2D5/2 clock transitions of 171Yb+ are presented. This is accomplished by incorporating contributions from the computationally challenging triply excited configurations through the relativistic coupled-cluster (RCC) theory. Testament of quality atomic wave functions of states of the above transitions, obtained using the RCC theory, are gauged by comparing the calculated energies and magnetic dipole hyperfine structure constants with their measurements. The improved quadrupole moments from this work will be immensely useful to estimate quadrupole shifts of the clock transitions of Yb+. Complementary approaches are employed to ascertain accuracy and comprehend roles of orbital relaxation and correlation effects in evaluating the IS constants. Combining these constants with the IS measurements from Phys. Rev. Lett. {\bf 128}, 163201 (2022), differential nuclear charge radii of the Yb isotopes are inferred that deviate by 6-7\% from the literature data.
Convolutional operations are computationally intensive in artificial intelligence services, and their overhead in electronic hardware limits machine learning scaling. Here, we introduce a photonic joint transform correlator (pJTC) using a near-energy-free on-chip Fourier transformation to accelerate convolution operations. The pJTC reduces computational complexity for both convolution and cross-correlation from O(N4) to O(N2), where N2 is the input data size. Demonstrating functional Fourier transforms and convolution, this pJTC achieves 98.0% accuracy on an exemplary MNIST inference task. Furthermore, a wavelength-multiplexed pJTC architecture shows potential for high throughput and energy efficiency, reaching 305 TOPS/W and 40.2 TOPS/mm2, based on currently available foundry processes. An efficient, compact, and low-latency convolution accelerator promises to advance next-generation AI capabilities across edge demands, high-performance computing, and cloud services.
We present a novel approach to hard-constrain Neumann boundary conditions in physics-informed neural networks (PINNs) using Fourier feature embeddings. Neumann boundary conditions are used to described critical processes in various application, yet they are more challenging to hard-constrain in PINNs than Dirichlet conditions. Our method employs specific Fourier feature embeddings to directly incorporate Neumann boundary conditions into the neural network's architecture instead of learning them. The embedding can be naturally extended by high frequency modes to better capture high frequency phenomena. We demonstrate the efficacy of our approach through experiments on a diffusion problem, for which our method outperforms existing hard-constraining methods and classical PINNs, particularly in multiscale and high frequency scenarios.
Particle physics has an ambitious and broad global experimental programme for the coming decades. Large investments in building new facilities are already underway or under consideration. Scaling the present processing power and data storage needs by the foreseen increase in data rates in the next decade for HL-LHC is not sustainable within the current budgets. As a result, a more efficient usage of computing resources is required in order to realise the physics potential of future experiments. Software and computing are an integral part of experimental design, trigger and data acquisition, simulation, reconstruction, and analysis, as well as related theoretical predictions. A significant investment in computing and software is therefore critical. Advances in software and computing, including artificial intelligence (AI) and machine learning (ML), will be key for solving these challenges. Making better use of new processing hardware such as graphical processing units (GPUs) or ARM chips is a growing trend. This forms part of a computing solution that makes efficient use of facilities and contributes to the reduction of the environmental footprint of HEP computing. The HEP community already provided a roadmap for software and computing for the last EPPSU, and this paper updates that, with a focus on the most resource critical parts of our data processing chain.