Computer Physics Communications最新文献

{"title":"Learning spatiotemporal dynamics from sparse data via a high-order physics-encoded network","authors":"Pu Ren ,&nbsp;Jialin Song ,&nbsp;Chengping Rao ,&nbsp;Qi Wang ,&nbsp;Yike Guo ,&nbsp;Hao Sun ,&nbsp;Yang Liu","doi":"10.1016/j.cpc.2025.109582","DOIUrl":"10.1016/j.cpc.2025.109582","url":null,"abstract":"<div><div>Learning unknown or partially known dynamics has gained significant attention in scientific machine learning (SciML). This research is mainly driven by the inherent sparsity and noise in scientific data, which poses challenges to accurately modeling spatiotemporal systems. While recent physics-informed learning strategies have attempted to address this problem by incorporating physics knowledge as soft constraints, they often encounter optimization and scalability issues. To this end, we present a novel physics-encoded learning framework for capturing the intricate dynamical patterns of spatiotemporal systems from limited sensor measurements. Our approach centers on a deep convolutional-recurrent network, termed Π<span>-block</span>, which hard-encodes known physical laws (e.g., PDE structure and boundary conditions) into the learning architecture. Moreover, the high-order time marching scheme (e.g., Runge-Kutta fourth-order) is introduced to model the temporal evolution. We conduct comprehensive numerical experiments on a variety of complex systems to evaluate our proposed approach against baseline algorithms across two tasks: reconstructing high-fidelity data and identifying unknown system coefficients. We also assess the performance of our method under various noisy levels and using different finite difference kernels. The comparative results demonstrate the superiority, robustness, and stability of our framework in addressing these critical challenges in SciML.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109582"},"PeriodicalIF":7.2,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679491","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Enhancing the Nektar++ spectral/hp element framework for parallel-in-time simulations","authors":"Jacques Y. Xing ,&nbsp;Chris D. Cantwell ,&nbsp;David Moxey","doi":"10.1016/j.cpc.2025.109584","DOIUrl":"10.1016/j.cpc.2025.109584","url":null,"abstract":"<div><div>We describe the efficient implementation of the Parareal algorithm in the <em>Nektar++</em> software, an open-source spectral/hp element framework for the solution of partial differential equations, which has been designed to achieve high-scalability on high-performance computing (HPC) clusters using distributed parallelism. Recently, time-parallel integration techniques are being recognized as a potential solution to further increase concurrency and computational speed-up beyond the limits of strong scaling obtained from a pure spatial domain decomposition. Amongst the various time-parallel approaches proposed in the literature, the Parareal algorithm is a non-intrusive and iterative approach, exploiting a fine and a coarse solvers to achieve time-parallelism, and can be applied to both linear and non-linear problems. We discuss the details of the implementation and discuss the specific techniques used to adapt the code to a time-parallel framework. We demonstrate the application of these methods to multiple linear and non-linear problems provided by the existing <em>Nektar++</em> solvers.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109584"},"PeriodicalIF":7.2,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679477","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Grid-free evaluation of phonon-limited electronic relaxation times and transport properties","authors":"Nenad Vukmirović","doi":"10.1016/j.cpc.2025.109583","DOIUrl":"10.1016/j.cpc.2025.109583","url":null,"abstract":"<div><div>Present calculations of electrical transport properties of materials require evaluations of electron-phonon coupling constants on dense predefined grids of electron and phonon momenta and performing the sums over these momenta. In this work, we present the methodology for calculation of carrier relaxation times and electrical transport properties without the use of a predefined grid. The relaxation times are evaluated by integrating out the delta function that ensures energy conservation and performing an average over the angular components of phonon momentum. The charge carrier mobility is then evaluated as a sum over appropriately sampled electronic momenta. We illustrate our methodology by applying to the Fröhlich model and to a real semiconducting material ZnTe. We find that rather accurate results can be obtained with a modest number of electron and phonon momenta, on the order of one hundred each, regardless of the carrier effective mass.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109583"},"PeriodicalIF":7.2,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143697066","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"GPU acceleration of overbridging boundary matching method without Green's functions based on real-space finite-difference method","authors":"Takanori Akamatsu ,&nbsp;Mitsuharu Uemoto ,&nbsp;Yoshiyuki Egami ,&nbsp;Tomoya Ono","doi":"10.1016/j.cpc.2025.109585","DOIUrl":"10.1016/j.cpc.2025.109585","url":null,"abstract":"<div><div>We present the graphics processing unit (GPU) acceleration of the overbridging boundary matching method for electron-transport property calculations, which is based on the density functional theory using the real-space finite-difference method. The execution of the implemented code using OpenACC and CUDA libraries on GPU is computationally more efficient and faster than a central processing unit. Furthermore, we achieve the ideal scalability in parallel execution from one to thirty-two nodes by adopting multiprocess parallelization schemes for two types of supercomputer with different configurations. To demonstrate the applicability of the accelerated code, the complex band structures of graphene and armchair carbon nanotubes with chiral indices (6, 6), (9, 9), and (12, 12) are calculated and compared with those obtained by the tight-binding method. We discuss the effects of the dispersion of evanescent waves on roll-up and axial strain in the carbon nanotubes.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109585"},"PeriodicalIF":7.2,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679584","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"SDGMPS: A spin-dependent Glauber model program for elastic proton-nucleus scattering","authors":"J.T. Zhang,&nbsp;X.L. Tu,&nbsp;Y. Huang,&nbsp;L.Y. Li,&nbsp;G.Q. Zhang,&nbsp;Z.H. Li","doi":"10.1016/j.cpc.2025.109587","DOIUrl":"10.1016/j.cpc.2025.109587","url":null,"abstract":"<div><div>SDGMPS is a Fortran program that calculates differential cross sections of elastic proton-nucleus scattering at intermediate energies based on the spin-dependent Glauber model. In the program, the Glauber model explicitly takes into account spin effects by using the spin-dependent nucleon-nucleon scattering amplitude, where the spin-orbit amplitude parameters are needed as input. It is particularly useful for analyses of the elastic proton scattering at both low and high momentum transfers and studies of the inner density distributions in nuclei. Such studies are an important part of the physics research program of the radiation beam facilities, such as the Heavy Ion Research Facility in Lanzhou (HIRFL).</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109587"},"PeriodicalIF":7.2,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679481","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"ALATDYN: A set of Anharmonic LATtice DYNamics codes to compute thermodynamic and thermal transport properties of crystalline solids","authors":"Keivan Esfarjani ,&nbsp;Harold Stokes ,&nbsp;Safoura Nayeb Sadeghi ,&nbsp;Yuan Liang ,&nbsp;Bikash Timalsina ,&nbsp;Han Meng ,&nbsp;Junichiro Shiomi ,&nbsp;Bolin Liao ,&nbsp;Ruoshi Sun","doi":"10.1016/j.cpc.2025.109575","DOIUrl":"10.1016/j.cpc.2025.109575","url":null,"abstract":"&lt;div&gt;&lt;div&gt;We introduce a lattice dynamics package which calculates elastic, thermodynamic and thermal transport properties of crystalline materials from data on their force and potential energy as a function of atomic positions. The data can come from density functional theory (DFT) calculations or classical molecular dynamics runs performed in a supercell. First, the model potential parameters, which are anharmonic force constants are extracted from the latter runs. Then, once the anharmonic model is defined, thermal conductivity and equilibrium properties at finite temperatures can be computed using lattice dynamics, Boltzmann transport theories, and a variational principle respectively. In addition, the software calculates the mechanical properties such as elastic tensor, Gruneisen parameters and the thermal expansion coefficient within the quasi-harmonic approximation (QHA). Phonons, elastic constants and thermodynamic properties results applied to the germanium crystal will be illustrated. Using the force constants as a force field, one may also perform molecular dynamics (MD) simulations in order to investigate the combined effects of anharmonicity and defect scattering beyond perturbation theory.&lt;/div&gt;&lt;/div&gt;&lt;div&gt;&lt;h3&gt;Program summary&lt;/h3&gt;&lt;div&gt;&lt;em&gt;Program Title:&lt;/em&gt; ALATDYN&lt;/div&gt;&lt;div&gt;&lt;em&gt;CPC Library link to program files:&lt;/em&gt; &lt;span&gt;&lt;span&gt;https://doi.org/10.17632/4jm4fh2nk2.1&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;&lt;em&gt;Developer's repository link:&lt;/em&gt; &lt;span&gt;&lt;span&gt;https://github.com/KeivanS/Anharmonic-lattice-dynamics/tree/main&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;&lt;em&gt;Licensing provisions:&lt;/em&gt; GPLv3&lt;/div&gt;&lt;div&gt;&lt;em&gt;Programming language:&lt;/em&gt; FORTRAN90&lt;/div&gt;&lt;div&gt;&lt;em&gt;Nature of problem:&lt;/em&gt; The ALATDYN suite of codes develops a lattice dynamical model of a crystalline solid. The FOCEX code extracts the model parameters from supercell calculations data on forces versus position and calculates the phonon spectrum, elastic constants and thermodynamic properties within the quasi-harmonic approximation. The SCOP8 code goes beyond QHA and implements the self-consistent phonon theory to minimize the free energy with respect to the strain tensor, atomic positions and harmonic force constants, and thus obtains the state of equilibrium at the given temperature along with the effective phonon bands structure. The THERMACOND code uses the cubic force constants and the crystal symmetries to solve the phonon Boltzmann equation (PhBE) efficiently and deduce the thermal conductivity. Finally, ANFOMOD uses the extracted force constants to perform a molecular dynamics simulation in a supercell.&lt;/div&gt;&lt;div&gt;&lt;em&gt;Solution method:&lt;/em&gt; Force constants are obtained from a singular value decomposition (or the ridge regression) method. PhBE is solved by first setting up the collision matrix and effectively inverting it using the conjugate-gradients method.&lt;/div&gt;&lt;div&gt;&lt;em&gt;Additional comments including restrictions and unusual features:&lt;/em&gt; This code has","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109575"},"PeriodicalIF":7.2,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679486","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"KARL - a Monte Carlo model for atomic and molecular processes in the tritium atmosphere of the KATRIN experiment","authors":"Christian Sendlinger ,&nbsp;Jonas Kellerer ,&nbsp;Felix Spanier","doi":"10.1016/j.cpc.2025.109580","DOIUrl":"10.1016/j.cpc.2025.109580","url":null,"abstract":"<div><div>A new parallelized simulation code is presented, which uses a Monte Carlo method to determine particle spectra in the KATRIN source. Reaction chains are generated from the decay of tritium within the source. The code includes all relevant processes: elastic scattering, ionization, excitation (electric, vibrational, rotational), recombination and various clustering processes. The main emphasis of the code is the calculation of particle spectra and particle densities and currents at specific points within the source. It features a new technique to determine these quantities. It also calculates target fields for the interaction of particles with each other as it is needed for recombination processes.</div><div>The code has been designed for the KATRIN experiment but is easily adaptable for other tritium based experiments like Project 8. Geometry and background tritium gas flow can be given as user input.</div><div>The code is parallelized using MPI and writes output using HDF5. Input to the simulation is read from a JSON description.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> KARL - <strong>KA</strong>trin WGTS elect<strong>R</strong>on and ion spectrum Monte Car<strong>L</strong>o</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/5bj3vwc6rg.1</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> GNU Public License v3</div><div><em>Programming language:</em> C++</div><div><em>External routines/libraries:</em> C++ compiler (tested with g++ 8.2 and 9.4.0), MPI 1.1 (tested with OpenMPI 3.1), HDF5 with support for parallel I/O (tested with version 1.10.0), Blitz++ (tested with version 1.0.2), Jansson (tested with version 2.12 and 2.13)</div><div><em>Nature of problem:</em> In the KATRIN experiment (and other experiments alike that feature large vessels filled with tritium) electrons are created from beta decay. These electrons interact with the ambient gas to produce secondary electrons through ionization. Subsequent processes include excitation, secondary ionization and collisions. The resulting electron and ion differential energy spectrum at various positions is relevant for further plasma analysis, and the current of charged particles to the ends of the experiments is an observable.</div><div><em>Solution method:</em> Semi-classical Monte Carlo.</div><div><em>Additional comments including restrictions and unusual features:</em> The geometry of the experiment is currently limited to the KATRIN experiment, but this may easily be changed. The configuration is stored in JSON files.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109580"},"PeriodicalIF":7.2,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644199","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"PolyMorph: Extension of PolyHoop for tissue morphogenesis coupled to chemical signaling","authors":"Nicolas Pascal Guido Müller ,&nbsp;Roman Vetter","doi":"10.1016/j.cpc.2025.109581","DOIUrl":"10.1016/j.cpc.2025.109581","url":null,"abstract":"&lt;div&gt;&lt;div&gt;We present PolyMorph, a lightweight standalone C++ program that extends its predecessor PolyHoop by a finite-difference solver for multi-component reaction-advection-diffusion equations. PolyMorph simulates two integral parts of tissue morphogenesis in two dimensions: 1) the mechanics of cellular deformation, growth and proliferation, and 2) transport and reaction of an arbitrary number of chemical species. Both of these components are bidirectionally coupled, allowing cells to base their behavior on local information on concentrations and flow, and allowing the chemical transport and reaction kinetics to depend on spatial information such as the local cell type. This bidirectional feedback makes PolyMorph a versatile tool to study a variety of cellular morphogenetic processes such as chemotaxis, cell sorting, tissue patterning with morphogen gradients, Turing patterning, and diffusion- or supply-limited growth with sub-cellular resolution.&lt;/div&gt;&lt;/div&gt;&lt;div&gt;&lt;h3&gt;Program summary&lt;/h3&gt;&lt;div&gt;&lt;em&gt;Program Title:&lt;/em&gt; PolyMorph&lt;/div&gt;&lt;div&gt;&lt;em&gt;CPC Library link to program files:&lt;/em&gt; &lt;span&gt;&lt;span&gt;https://doi.org/10.17632/4jscxhkd2s.2&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;&lt;em&gt;Licensing provisions:&lt;/em&gt; BSD 3-clause&lt;/div&gt;&lt;div&gt;&lt;em&gt;Programming language:&lt;/em&gt; C++11&lt;/div&gt;&lt;div&gt;&lt;em&gt;Supplementary material:&lt;/em&gt; Figure 1&lt;/div&gt;&lt;div&gt;&lt;em&gt;Journal reference of previous version:&lt;/em&gt; Comput. Phys. Commun. 299 (2024) 109128, &lt;span&gt;&lt;span&gt;https://doi.org/10.1016/j.cpc.2024.109128&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;&lt;em&gt;Does the new version supersede the previous version?:&lt;/em&gt; No&lt;/div&gt;&lt;div&gt;&lt;em&gt;Nature of problem:&lt;/em&gt; In tissue development and disease, morphogenesis and cell fate determination depends on mechanical processes as well as chemical signaling. PolyMorph couples the Newtonian mechanics of deformable cells (including growth and proliferation) in 2D with a customizable set of reaction-advection-diffusion equations to simulate problems that require an integrated approach with chemical-mechanical interactions. Typical use cases include the patterning of epithelial tissues with chemical signals (e.g., morphogen gradients or the Turing mechanism), chemotaxis and cell migration, wound healing, diffusion- or nutrition-limited growth, regulatory network dynamics in a spatial cellular environment, and other problems in tissue self-organization. PolyMorph enables the numerical solution of such problems with bidirectional feedback between mechanics and chemistry, in large monolayer tissues and with an arbitrary number of interacting species.&lt;/div&gt;&lt;div&gt;&lt;em&gt;Solution method:&lt;/em&gt; The off-lattice polygonal representation of cell boundaries in PolyHoop [1] is coupled to a lattice representation of diffusing chemical reactants. The reaction-advection-diffusion problem is solved with the finite difference method using the standard 5-point central difference stencil, and explicitly integrated in time. A scatter-gather approach inspired by the particle-in-ce","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109581"},"PeriodicalIF":7.2,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143642219","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"AFiD-Darcy: A finite difference solver for numerical simulations of convective porous media flows","authors":"Marco De Paoli ,&nbsp;Guru Sreevanshu Yerragolam ,&nbsp;Detlef Lohse ,&nbsp;Roberto Verzicco","doi":"10.1016/j.cpc.2025.109579","DOIUrl":"10.1016/j.cpc.2025.109579","url":null,"abstract":"<div><div>We present an efficient solver for massively-parallel simulations of convective, wall-bounded and incompressible porous media flows. The algorithm consists of a second-order finite-difference pressure-correction scheme, allowing the use of an efficient FFT-based solver in problems with different boundary conditions. The parallelization method is implemented in a two-dimensional pencil-like domain decomposition, which enables efficient parallel large-scale simulations. The original version of the code presented by van der Poel et al. (2015) <span><span>[35]</span></span> has been modified to solve the Darcy equation for the momentum transport, representative of porous media flows driven by buoyancy. Two schemes are implemented to treat the diffusive term of the advection-diffusion equation, namely a fully implicit and semi-implicit formulation. Despite exhibiting a higher computational cost per time step, the fully implicit scheme allows an efficient simulation of transient flows, leading to a smaller time-to-solution compared to the semi-implicit scheme. The implementation was verified against different canonical flows, and the computational performance was examined. To show the code's capabilities, the maximal driving strength explored has been doubled as compared to state-of-art simulations, corresponding to an increase of the associated computational effort of about 8 to 16 times. Excellent strong scaling performance is demonstrated for both schemes developed and for domains with more than 10¹⁰ spatial degrees of freedom.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> AFiD-Darcy</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/xhx3gzpj6n.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link:</em> <span><span>https://github.com/depaolimarco/AFiD-Darcy</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> CC BY 4.0</div><div><em>Programming language:</em> Fortran 90, MPI</div><div><em>External routines:</em> FFTW3, HDF5</div><div><em>Nature of problem:</em> Solving two- and three-dimensional Darcy equation coupled with a scalar field in a box bounded between two walls in one-direction and with periodic boundary conditions in the other two directions.</div><div><em>Solution method:</em> Second order finite difference method for spatial discretization, third order Runge–Kutta scheme in combination with Crank–Nicolson for the implicit terms for time advancement, two dimensional pencil distributed MPI parallelization. Implicit and semi-implicit formulations for the solution of the diffusive terms in the scalar transport equation.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109579"},"PeriodicalIF":7.2,"publicationDate":"2025-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679478","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Paths towards time evolution with larger neural-network quantum states","authors":"Wenxuan Zhang ,&nbsp;Bo Xing ,&nbsp;Xiansong Xu ,&nbsp;Dario Poletti","doi":"10.1016/j.cpc.2025.109577","DOIUrl":"10.1016/j.cpc.2025.109577","url":null,"abstract":"<div><div>In recent years, neural-network quantum states method in conjunction with the time-dependent variational Monte Carlo have been proposed to study the dynamics of many-body quantum systems. By interpreting the quantum dynamics problem as a ground state search of an effective Hamiltonian, we show that one can use stochastic reconfiguration (SR), a remarkable method that significantly boosts the efficiency and convergence of the variational training. Furthermore, since the vanilla SR method does not scale efficiently when the size of neural-network quantum states increases, we transfer to the study of time-dependent systems, or introduce altogether, three approaches that reduce the computational complexity of the SR method, and we compare their performance: Kronecker-factored approximate curvature (K-FAC), minimum-step stochastic reconfiguration (minSR), and sequential overlapping optimization (SOO). To demonstrate the generality of these approaches, we use both the restricted Boltzmann machine and the feed-forward neural network. We consider a titled Ising model and study the quantum quench from the paramagnetic to the anti-ferromagnetic phase. We show that the three approaches allow to use stochastic reconfigurations to describe the time evolution of a many-body quantum system using a neural network with more than 10000 parameters, which would be prohibitive otherwise. For systems up to 40 spins, we observe that minSR and SOO have similar performance and both provide better accuracy than K-FAC.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109577"},"PeriodicalIF":7.2,"publicationDate":"2025-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143620349","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}