Search results for " Computational Physics"
showing 10 items of 116 documents
On the p-length of some finite p-soluble groups
2014
The main aim of this paper is to give structural information of a finite group of minimal order belonging to a subgroup-closed class of finite groups and whose $p$-length is greater than $1$, $p$ a prime number. Alternative proofs and improvements of recent results about the influence of minimal $p$-subgroups on the $p$-nilpotence and $p$-length of a finite group arise as consequences of our study
A New Three-Dimensional Track Fit with Multiple Scattering
2017
Modern semiconductor detectors allow for charged particle tracking with ever increasing position resolution. Due to the reduction of the spatial hit uncertainties, multiple Coulomb scattering in the detector layers becomes the dominant source for tracking uncertainties. In this case long distance effects can be ignored for the momentum measurement, and the track fit can consequently be formulated as a sum of independent fits to hit triplets. In this paper we present an analytical solution for a three-dimensional triplet(s) fit in a homogeneous magnetic field based on a multiple scattering model. Track fitting of hit triplets is performed using a linearization ansatz. The momentum resolution…
PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE
2021
Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kind of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, stand-alone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming lan…
LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories
2021
We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples.
The Ferroelectric Photo-Groundstate of SrTiO$_3$: Cavity Materials Engineering
2021
Significance Controlling collective phenomena in quantum materials is a promising route toward engineering material properties on demand. Strong THz lasers have been successful at inducing ferroelectricity in S r T i O 3 . Here we demonstrate, from atomistic calculations, that cavity quantum vacuum fluctuations induce a change in the collective phase of S r T i O 3 in the strong light–matter coupling regime. Under these conditions, the ferroelectric phase is stabilized as the ground state, instead of the quantum paraelectric one. We conceptualize this light–matter hybrid state as a material photo ground state: Fundamental properties such as crystal structure, phonon frequencies, and the col…
On numerical broadening of particle size spectra: a condensational growth study using PyMPDATA
2020
This work discusses the numerical aspects of representing the diffusional (condensational) growth in particulate systems such as atmospheric clouds. It focuses on the Eulerian modeling approach, in which the evolution of the particle size spectrum is carried out using a fixed-bin discretization associated with inherent numerical diffusion. Focus is on the applications of MPDATA numerical schemes (variants explored include: infinite-gauge, non-oscillatory, third-order-terms and recursive antidiffusive correction). Methodology for handling coordinate transformations associated with both particle size distribution variable choice and numerical grid layout are expounded. Analysis of the perform…
On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method
2017
The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations ($Au=b$) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of $P$ relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of $M$ consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method. …
Scheduled Relaxation Jacobi method: improvements and applications
2016
Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficien…
Efficient numerical integration of neutrino oscillations in matter
2016
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.
Self-consistent field theory based molecular dynamics with linear system-size scaling
2012
We present an improved field-theoretic approach to the grand-canonical potential suitable for linear scaling molecular dynamics simulations using forces from self-consistent electronic structure calculations. It is based on an exact decomposition of the grand canonical potential for independent fermions and does neither rely on the ability to localize the orbitals nor that the Hamilton operator is well-conditioned. Hence, this scheme enables highly accurate all-electron linear scaling calculations even for metallic systems. The inherent energy drift of Born-Oppenheimer molecular dynamics simulations, arising from an incomplete convergence of the self-consistent field cycle, is circumvented …