Related papers: Numerical solution of perturbed Kepler problem usi…
Quantum Hamiltonians containing nonseparable products of non-commuting operators, such as $\hat{\bf x}^m \hat{\bf p}^n$, are problematic for numerical studies using split-operator techniques since such products cannot be represented as a…
Multi-component polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures,…
The presence of A/F-type {\it Kepler} hybrid stars extending across the entire $\delta$ Sct-$\gamma$ Dor instability strips and beyond remains largely unexplained. In order to better understand these particular stars, we performed a…
The transmutation (transformation) operator associated with the perturbed Bessel equation is considered. It is shown that its integral kernel can be uniformly approximated by linear combinations of constructed here generalized wave…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
Using an operatorial formalism, we study the Kramers equation and its applications to numerical simulations. We obtain classes of algorithms which may be made precise at every desired order in the time step $\epsilon$ and with a set of free…
Over short time intervals planetary ephemerides have been traditionally represented in analytical form as finite sums of periodic terms or sums of Poisson terms that are periodic terms with polynomial amplitudes. Nevertheless, this…
Inspired by Quantum Mechanics, we reformulate Hilbert's tenth problem in the domain of integer arithmetics into either a problem involving a set of infinitely coupled differential equations or a problem involving a Shr\"odinger propagator…
The aim of this article is to present two different primal-dual methods for solving structured monotone inclusions involving parallel sums of compositions of maximally monotone operators with linear bounded operators. By employing some…
We present a new method for computing orbits in the perturbed two-body problem: the position and velocity vectors of the propagated object in Cartesian coordinates are replaced by eight orbital elements, i.e., constants of the unperturbed…
Determining the vibrational structure of a molecule is central to fundamental applications in several areas, from atmospheric science to catalysis, fuel combustion modeling, biochemical imaging, and astrochemistry. However, when significant…
The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear differential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the exact solutions are obtained. The rapid convergence…
Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of…
We propose a computationally efficient method to solve the dynamics of operators of bosonic quantum systems coupled to their environments. The method maps the operator under interest to a set of complex-valued functions, and its adjoint…
In a Keplerian system, a large number of bodies orbit a central mass. Accretion disks, protoplanetary disks, asteroid belts, and planetary rings are examples. Simulations of these systems require algorithms that are computationally…
Explicit and semi-explicit geometric integration schemes for dissipative perturbations of Hamiltonian systems are analyzed. The dissipation is characterized by a small parameter $\epsilon$, and the schemes under study preserve the…
We discuss, within the simplified context provided by the polymeric harmonic oscillator, a construction leading to a separable Hilbert space that preserves some of the most important features of the spectrum of the Hamiltonian operator.…
We propose a numerical method to solve the Wigner equation in quantum systems of spinless, non-relativistic particles. The method uses a spectral decomposition into $L^2(\mathbb{R}^d)$ basis functions in momentum-space to obtain a system of…
An extended Hamiltonian approach to conduct isothermal-isobaric molecular dynamics simulations with full cell flexibility is presented. The components of the metric tensor are used as the fictitious degrees of freedom for the cell, thus…
A hyperbolic integro-differential equation is considered, as a model problem, where the convolution kernel is assumed to be either smooth or no worse than weakly singular. Well-posedness of the problem is studied in the context of semigroup…