| Home | [11–20] << 21 22 23 24 25 26 27 28 29 30 >> [31–40] |
|
Ferrer-Sanchez, A., Martin-Guerrero, J., Ruiz de Austri, R., Torres-Forne, A., & Font, J. A. (2024). Gradient-annihilated PINNs for solving Riemann problems: Application to relativistic hydrodynamics. Comput. Meth. Appl. Mech. Eng., 424, 116906–18pp.
Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that forces the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss. The performance of our methodology is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with the GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative l(2) error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ''ground truth'', and with the corresponding error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by the GA-PINN model is comparable to that obtained with a shock-capturing scheme, achieving a performance superior to that of the baseline PINN algorithm in general. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved variables, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN methodology discussed in this paper could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.
|
|
Allanach, B. C., Martin, S. P., Robertson, D. G., & Ruiz de Austri, R. (2017). The inclusion of two-loop SUSYQCD corrections to gluino and squark pole masses in the minimal and next-to-minimal supersymmetric standard model: SOFTSUSY3.7. Comput. Phys. Commun., 219, 339–345.
Abstract: We describe an extension of the SOFTSUSY spectrum calculator to include two-loop supersymmetric QCD (SUSYQCD) corrections of order O(alpha(2)(s)) to gluino and squark pole masses, either in the minimal supersymmetric standard model (MSSM) or the next-to-minimal supersymmetric standard model (NMSSM). This document provides an overview of the program and acts as a manual for the new version of SOFTSUSY, which includes the increase in accuracy in squark and gluino pole mass predictions. Program summary Program title: SOFTSUSY Program Files doi: http://dx.doLorg/10.17632/sh77x9j7hs.1 Licensing provisions: GNU GPLv3 Programming language: C++, fortran, C Nature of problem: Calculating supersymmetric particle spectrum, mixing parameters and couplings in the MSSM or the NMSSM. The solution to the renormalization group equations must be consistent with theoretical boundary conditions on supersymmetry breaking parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters. Solution method: Nested fixed point iteration. Restrictions: SOFTSUSY will provide a solution only in the perturbative regime and it assumes that all couplings of the model are real (i.e. CP-conserving). If the parameter point under investigation is nonphysical for some reason (for example because the electroWeak potential does not have an acceptable minimum), SOFTSUSY returns an error message. The higher order corrections included are for the MSSM (R-parity conserving or violating) or the real R-parity conserving NMSSM only. Journal reference of previous version: Comput. Phys. Comm. 189 (2015) 192. Does the new version supersede the previous version?: Yes. Reasons for the new version: It is desirable to improve the accuracy of the squark and gluinos mass predictions, since they strongly affect supersymmetric particle production cross-sections at colliders. Summary of revisions: The calculation of the squark and gluino pole masses is extended to be of next-to next-to leading order in SUSYQCD, i.e. including terms up to O(g(s)(4)/(16 pi(2))(2)). Additional comments: Program obtainable from http://softsusy.hepforge.org/
|
|
Aebischer, J., Brivio, I., Celis, A., Evans, J. A., Jiang, Y., Kumar, J., et al. (2018). WCxf : An exchange format for Wilson coefficients beyond the Standard Model. Comput. Phys. Commun., 232, 71–83.
Abstract: We define a data exchange format for numerical values of Wilson coefficients of local operators parameterising low-energy effects of physics beyond the Standard Model. The format facilitates interfacing model-specific Wilson coefficient calculators, renormalisation group (RG) runners, and observable calculators. It is designed to be unambiguous (defining a non-redundant set of operators with fixed normalisation in each basis), extensible (allowing the addition of new EFTs or bases by the user), and robust (being based on industry standard file formats with parsers implemented in many programming languages). We have implemented the format for the Standard Model EFT (SMEFT) and for the weak effective theory (WET) below the electroweak scale and have added interfaces to a number of public codes dealing with SMEFT or WET. We also provide command-line utilities and a Python module for convenient manipulation of WCxf files, including translation between different bases and matching from SMEFT to WET. (C) 2018 Elsevier B.V. All rights reserved.
|
|
Consiglio, R., de Salas, P. F., Mangano, G., Miele, G., Pastor, S., & Pisanti, O. (2018). PArthENoPE reloaded. Comput. Phys. Commun., 233, 237–242.
Abstract: We describe the main features of a new and updated version of the program PArthENoPE, which computes the abundances of light elements produced during Big Bang Nucleosynthesis. As the previous first release in 2008, the new one, PArthENoPE2.0, is publicly available and distributed from the code site, http://parthenope.na.infn.it . Apart from minor changes, which will be also detailed, the main improvements are as follows. The powerful, but not freely accessible, NAG routines have been substituted by ODEPACK libraries, without any significant loss in precision. Moreover, we have developed a Graphical User Interface (GUI) which allows a friendly use of the code and a simpler implementation of running for grids of input parameters. New Version program summary Program Title: PArthENoPE2.0 Program Files doi : http://dx.doi.org/10.17632/wvgr7d8yt9.1 Licensing provisions: GPLv3 Programming language: Fortran 77 and Python Supplementary material: User Manual available on the web page http://parthenope.na.infn.it Journal reference of previous version: Comput. Phys. Commun. 178 (2008) 956 971 Does the new version supersede the previous version?: Yes Reasons for the new version: Make the code more versatile and user friendly Summary of revisions: (1) Publicly available libraries (2) GUI for configuration Nature of problem: Computation of yields of light elements synthesized in the primordial universe Solution method: Livermore Solver for Ordinary Differential Equations (LSODE) for stiff and nonstiff systems
Keywords: Primordial nucleosynthesis; Cosmology; Neutrino physics
|
|
KM3NeT Collaboration(Aiello, S. et al), Alves Garre, S., Calvo, D., Carretero, V., Colomer, M., Corredoira, I., et al. (2020). gSeaGen: The KM3NeT GENIE-based code for neutrino telescopes. Comput. Phys. Commun., 256, 107477–15pp.
Abstract: The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between tracktype and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project. Program summary Program Title: gSeaGen CPC Library link to program files: http://dx.doi.org/10.17632/ymgxvy2br4.1 Licensing provisions: GPLv3 Programming language: C++ External routines/libraries: GENIE [1] and its external dependencies. Linkable to MUSIC [2] and PROPOSAL [3]. Nature of problem: Development of a code to generate detectable events in neutrino telescopes, using modern and maintained neutrino interaction simulation libraries which include the state-of-the-art physics models. The default application is the simulation of neutrino interactions within KM3NeT [4]. Solution method: Neutrino interactions are simulated using GENIE, a modern framework for Monte Carlo event generators. The GENIE framework, used by nearly all modern neutrino experiments, is considered as a reference code within the neutrino community. Additional comments including restrictions and unusual features: The code was tested with GENIE version 2.12.10 and it is linkable with release series 3. Presently valid up to 5 TeV. This limitation is not intrinsic to the code but due to the present GENIE valid energy range. References: [1] C. Andreopoulos at al., Nucl. Instrum. Meth. A614 (2010) 87. [2] P. Antonioli et al., Astropart. Phys. 7 (1997) 357. [3] J. H. Koehne et al., Comput. Phys. Commun. 184 (2013) 2070. [4] S. Adrian-Martinez et al., J. Phys. G: Nucl. Part. Phys. 43 (2016) 084001.
|