Related papers: A new method for computing self-gravity in an isol…
In this paper, a new model is proposed for the inverse random source scattering problem of the Helmholtz equation with attenuation. The source is assumed to be driven by a fractional Gaussian field whose covariance is represented by a…
A new numerical method for an inverse problem for an elliptic equation with unknown potential is proposed. In this problem the point source is running along a straight line and the source-dependent Dirichlet boundary condition is measured…
A new approximation scheme for nonperturbative renormalisation group equations for quantum gravity is introduced. Correlation functions of arbitrarily high order can be studied by resolving the full dependence of the renormalisation group…
We extend the work of Yen et al. (2012) and develop 2nd order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a…
Multiple solutions are common in various non-convex problems arising from industrial and scientific computing. Nonetheless, understanding the nontrivial solutions' qualitative properties seems limited, partially due to the lack of efficient…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
Prompt isolated leptons are essential in many analyses in high-energy particle physics but are subject to fake-lepton background, i.e. objects that mimic the lepton signature. The fake-lepton background is difficult to estimate from…
We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. in [6] for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic…
We present a novel deep learning approach to approximate the solution of large, sparse, symmetric, positive-definite linear systems of equations. These systems arise from many problems in applied science, e.g., in numerical methods for…
This thesis introduces a set of methods for testing models of modified gravity using galaxy clusters. In particular, a technique for constraining models with a chameleon screening is introduced. In addition, the outlined technique is…
The envelope theory is a method to compute approximate eigensolutions of quantum $N$-body Hamiltonians with a quite general structure in $D$ dimensions. The advantages of the method are that it is easy to implement and that $N$ is treated…
The entropy computation of Gaussian mixture distributions with a large number of components has a prohibitive computational complexity. In this paper, we propose a novel approach exploiting the sphere decoding concept to bound and…
The outcomes of projective measurements on a quantum many-body system in a chosen basis are inherently probabilistic. The Shannon entropy of this probability distribution (the "diagonal entropy") often reveals universal features, such as…
The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system,…
A simple reweighting scheme is proposed for Monte Carlo simulations of interacting particle systems, permitting one to study various parameter values in a single study, and improving efficiency by an order of magnitude. Unlike earlier…
We develop and analyse finite volume methods for the Poisson problem with boundary conditions involving oblique derivatives. We design a generic framework, for finite volume discretisations of such models, in which internal fluxes are not…
This paper presents an iterative method suitable for inverting semilinear problems which are important kernels in many numerical applications. The primary idea is to employ a parametrization that is able to reduce semilinear problems into…
We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical…
A viable weak-field and slow-motion approximation method is constructed in $F(R,R_{\mu\nu}R^{\mu\nu}, R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma})$ gravity, a general class of fourth-order theories of gravity. By applying this method, the…
We present a new rational approximation algorithm based on the empirical interpolation method for interpolating a family of parametrized functions to rational polynomials with invariant poles, leading to efficient numerical algorithms for…