Skip to main content

Computer Physics Communications

ISSN: 0010-4655

Visit Journal website

Datasets associated with articles published in Computer Physics Communications

Filter Results
1970 2024
5954 results
  • CUDA-based focused Gaussian beams second-harmonic generation efficiency calculator
    We present an object-oriented programming (OOP) CUDA-based package for fast and accurate simulation of second-harmonic generation (SHG) efficiency using focused Gaussian beams. The model includes linear as well as two-photon absorption that can ultimately lead to thermal lensing due to self-heating effects. Our approach speeds up calculations by nearly 40x (11x) without (with) temperature profiles with respect to an equivalent implementation using CPU. The package offers a valuable tool for experimental design and study of 3D field propagation in nonlinear three-wave interactions. It is useful for optimization of SHG-based experiments and mitigates undesired thermal effects, enabling improved oven designs and advanced device architectures, leading to stable, efficient high-power SHG.
    • Dataset
  • FlexibleSUSY extended to automatically compute physical quantities in any beyond the standard model theory: Charged lepton flavor violation processes, Higgs decays, and user-defined observables
    FlexibleSUSY is a framework for the automated computation of physical quantities (observables) in models beyond the Standard Model (BSM). This paper describes an extension of FlexibleSUSY which allows to define and add new observables that can be enabled and computed in applicable user-defined BSM models. The extension has already been used to include Charged Lepton Flavor Violation (CLFV) observables, but further observables can now be added straightforwardly. The paper is split into two parts. The first part is non-technical and describes from the user's perspective how to enable the calculation of predefined observables, in particular CLFV observables. The second part of the paper explains how to define new observables such that their automatic computation in any applicable BSM model becomes possible. A key ingredient is the new NPointFunctions extension which allows to use tree-level and loop calculations in the model-independent setup of observables. Three examples of increasing complexity are fully worked out. This illustrates the features and provides code snippets that may be used as a starting point for implementation of further observables.
    • Dataset
  • The TDHF code Sky3D version 1.2
    The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.
    • Dataset
  • Multem 3: An updated and revised version of the program for transmission and band calculations of photonic crystals
    We present here Multem 3, an updated and revised version of Multem 2, which syntax has been upgraded to Fortran 2018, with the source code being divided into modules. Multem 3 is equipped with LAPACK, the state-of-the art Faddeeva complex error function routine, and the Bessel function package AMOS. The amendments significantly improve both the speed, convergence, and precision of Multem 2. Increased stability allows to freely increase the cut-off value LMAX on the number of spherical vector wave functions and the cut-off value RMAX controlling the maximal length of reciprocal vectors taken into consideration. An immediate bonus is that Multem 3 can be reliably used to describe bound states in the continuum (BICs). To ensure convergence of the layer coupling scheme, it appears that appreciably larger values of convergence paramaters LMAX and RMAX are required than those reported in numerous published work in the past using Multem 2. We hope that Multem 3 will become a reliable and fast alternative to generic commercial software, such as COMSOL Multiphysics, CST Microwave Studio, or Ansys HFSS, and that it will become the code of choice for various optimization tasks for a large number of research groups. The improvements concern the core part of Multem 2, which is common to the extensions of Multem 2 for acoustic and elastic multiple scattering and to the original layer-Kohn-Korringa-Rostocker (LKKR) code. Therefore, the enhancements presented here can be readily applied to the above codes as well. The previous version of this program may be found at
    • Dataset
  • HYPIC: A fast hybrid EM PIC-MCC code for ion cyclotron resonance energization in cylindrical coordinate system
    Ion cyclotron resonance energization (ICRE) such as ion cyclotron resonance heating (ICRH) is widely applied to magnetic confinement fusion and high-power electric propulsion. Since ICRE involves cyclotron resonance processes, a kinetic model is required. Both conventional particle-in-cell (PIC) simulations and solving the Boltzmann equation require enormous computation and memory. The hybrid simulation incorporating of adiabatic electrons and PIC ions provides a viable solution for both a substantial reduction in computation and the inclusion of cyclotron resonance effects. Under the adiabatic electron approximation, we have developed a two-dimensional (r, z) hybrid electromagnetic (EM) PIC-MCC (Monte-Carlo collision) simulation program, named HYPIC. The advantages of HYPIC are the inclusion of ion kinetic effects, electrostatic (ES) and EM effects, and collisional effects of ions and electrons, with a small computation. The HYPIC program is able to fast simulate the antenna-plasma interactions and the ion cyclotron resonance energization and/or ion cyclotron resonance heating processes in linear devices, such as high-power electric propulsion, magnetic mirror, and field-reversed-configuration (FRC), etc.
    • Dataset
  • FeAmGen.jl: A Julia program for Feynman Amplitude Generation
    FeAmGen.jl is a Julia package designed to generate Feynman diagrams and their corresponding amplitudes for various processes in particle physics. Utilizing the models in the Universal Feynman Output (UFO) format and Qgraf for diagram generation, it also employs SymEngine.jl and Form for amplitude generation. Additionally, the package offers functions for constructing Feynman integral topologies. In conclusion, FeAmGen.jl provides usability and versatility for the high precision calculations of the perturbative quantum field theory in the Standard Model or in the extensions beyond it. The corresponding codes are available at
    • Dataset
  • New version of Hyper-Fractal Analysis application for estimating the fuzzy fractal dimension of hyperspectral satellite ocean color images
    In this work we present a new version of the Hyper-Fractal Analysis C# .Net application for estimating the fuzzy fractal dimension of hyperspectral images. As an example, we also present some initial results on 116D algal images obtained from the HYPSO-1 satellite.
    • Dataset
  • H-COUP Version 3: A program for one-loop corrected decays of any Higgs bosons in non-minimal Higgs models
    The H-COUP program is provided as a package of Fortran codes, which can compute observables related to Higgs bosons including radiative corrections in various extended Higgs sectors. We give a manual for the latest version of H-COUP (H-COUP_3.0), in which decay rates and branching ratios of all the Higgs bosons can be calculated at one-loop level in EW and Higgs interactions with QCD corrections in the Higgs singlet model, four types of the two Higgs doublet model with a softly-broken Z_2 symmetry, and the inert doublet model. The previous version (H-COUP_2.0) can evaluate those only for the standard model like Higgs boson with the mass of 125 GeV (h). In H-COUP_3.0, renormalized quantities are computed based on the gauge-independent on-shell renormalization scheme. The source code of H-COUP_3.0 can be downloaded via the following link: By using H-COUP_3.0, we can compare the precise measurements of the properties of h and direct searches for additional Higgs bosons with their predictions at one-loop level, by which we can reconstruct the structure of the Higgs sector.
    • Dataset
  • ModelFLOWs-app: Data-driven post-processing and reduced order modelling tools
    This article presents an innovative open-source software named ModelFLOWs-app,1 written in Python, which has been created and tested to generate precise and robust hybrid reduced order models (ROMs) fully data-driven. By integrating modal decomposition and deep learning in diverse ways, the software uncovers the fundamental patterns in dynamic systems. This acquired knowledge is then employed to enrich the comprehension of the underlying physics, reconstruct databases from limited measurements, and forecast the progression of system dynamics. The hybrid ROMs produced by ModelFLOWs-app combine experimental and numerical databases, serving as highly accurate alternatives to numerical simulations. As a result, computational expenses are significantly reduced, and the models become powerful tools for optimization and control in various applications. The exceptional capability of ModelFLOWs-app in developing reliable data-driven hybrid ROMs has been demonstrated across a wide range of applications, making it a valuable resource for understanding complex nonlinear dynamical systems and providing insights in diverse domains. This article presents the mathematical background, as well as a review of some examples of applications.
    • Dataset
  • GeoTaichi: A Taichi-powered high-performance numerical simulator for multiscale geophysical problems
    This study introduces GeoTaichi, an open-source high-performance numerical simulator designed for addressing multiscale geophysical problems. By leveraging the power of the Taichi parallel language, GeoTaichi maximizes the utilization of modern computer resources on multicore CPU and GPU architectures. It offers robust and reliable modules for the discrete element method (DEM), material point method (MPM), and coupled material point-discrete element method (MPDEM). These modules enable efficient solving of large-scale problems while being implemented in pure Python. The design philosophy of GeoTaichi focuses on creating a framework that is readable, extensible, and user-friendly. This paper highlights the coupling procedure of MPDEM, the code structures, and the most important features of GeoTaichi. Rigorous benchmark tests have been conducted to verify the validity and robustness of GeoTaichi. Additionally, the performance of GeoTaichi is compared with similar software tools in the field, underscoring a notable improvement in both computational efficiency and memory savings when compared to existing alternatives.
    • Dataset