Related papers: Accuracy Controlled Schemes for the Eigenvalue Pro…
Reliable and efficient computation of the pseudospectral abscissa in the large-scale setting is still not settled. Unlike the small-scale setting where there are globally convergent criss-cross algorithms, all algorithms in the large-scale…
Linear kinetic transport equations play a critical role in optical tomography, radiative transfer and neutron transport. The fundamental difficulty hampering their efficient and accurate numerical resolution lies in the high dimensionality…
Interatomic potentials are essential to go beyond ab initio size limitations, but simulation results depend sensitively on potential parameters. Forward propagation of parameter variation is key for uncertainty quantification, whilst…
With increasing demand for accurate calculation of isotope shifts of atomic systems for fundamental and nuclear structure research, an analytic energy derivative approach is presented in the relativistic coupled-cluster theory framework to…
In medical SPECT imaging, we seek to simultaneously obtain the internal radioactive sources and the attenuation map using not only ballistic measurements but also first order scattering measurements. The problem is modeled using the…
The current capacity of computers makes it possible to perform simulations of small systems with portable, explicit-solvent potentials achieving high degree of accuracy. However, simplified models must be employed to exploit the behaviour…
In this paper we present a multilevel projection-based iterative scheme for solving thermal radiative transfer problems that performs iteration cycles on the high-order Boltzmann transport equation (BTE) and low-order moment equations.…
We propose a control-oriented optimal experimental design (cOED) approach for linear PDE-constrained Bayesian inverse problems. In particular, we consider optimal control problems with uncertain parameters that need to be estimated by…
We develop a "hybrid approximative scheme" in the framework of the post-Newtonian approximation for computing general-relativistic polytropic models simulating neutron stars in critical rigid rotation. We treat the differential equations…
The paper considers a linear system of Boltzmann transport equations modelling the evolution of three species of particles, photons, electrons and positrons. The system is coupled because of the collision term (an integral operator). The…
A substantial fraction of systematic uncertainties in neutrino oscillation experiments stem from the lack of precision in modeling the nuclear target in neutrino-nucleus interactions. Whilst this has driven significant progress in the…
Interpolation models are critical for a wide range of applications, from numerical optimization to artificial intelligence. The reliability of the provided interpolated value is of utmost importance, and it is crucial to avoid the…
We introduce a notion of nuclear numerical range defined as the set of expectation values of a given operator $A$ among normalized pure states, which belong to the nucleus of an auxiliary operator $Z$. This notion proves to be applicable to…
We investigate how unitary control can improve parameter estimation by designing the effective spectrum of the imprinting Hamiltonian. We show that, for commuting Hamiltonians, the general problem of spectral manipulation via unitary…
Uncertainties in the proton range in tissue during proton therapy limit the precision in treatment delivery. These uncertainties result in expanded treatment margins, thereby increasing radiation dose to healthy tissue. Real-time range…
Random walks are frequently used as a model for very diverse physical phenomena. The Monte Carlo method is a versatile tool for the study of the properties of systems modelled as random walks. Often, each walker is associated with a…
This work proposes a general strategy for solving possibly nonlinear problems arising from implicit time discretizations as a sequence of explicit solutions. The resulting sequence may exhibit instabilities similar to those of the base…
We study in this work an integral formulation for the radiative transfer equation (RTE) in anisotropic media with truncated approximation to the scattering phase function. The integral formulation consists of a coupled system of integral…
From characterizing the speed of a thermal system's response to computing natural modes of vibration, eigenvalue analysis is ubiquitous in engineering. In spite of this, eigenvalue problems have received relatively little treatment compared…
Many problems in physics, chemistry and other fields are perturbative in nature, i.e. differ only slightly from related problems with known solutions. Prominent among these is the eigenvalue perturbation problem, wherein one seeks the…