Related papers: The computation of seismic normal modes with rotat…
There is recent interest in the inter and intra element interactions of metamaterial unit cells. To calculate the effects of these interactions which can be substantial, an ab initio general coupled mode equation, in the form of an…
Elliptic partial differential equations arise in many fields of science and engineering such as steady state distribution of heat, fluid dynamics, structural/mechanical engineering, aerospace engineering and seismology etc. In three…
We herein utilize the general three-body problem (GTBP) as a model, in order to simulate resonant systems consisting of a star and two planets, where at least one of them is highly eccentric. We study them in terms of their long-term…
We develop a general-purpose formulation, based on two-dimensional spectral integrals, for computing electromagnetic fields produced by arbitrarily-oriented dipoles in planar-stratified environments, where each layer may exhibit arbitrary…
Modern N-body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the timestep for…
We present a multiscale mixed finite element method for solving second order elliptic equations with general $L^{\infty}$-coefficients arising from flow in highly heterogeneous porous media. Our approach is based on a multiscale spectral…
The spectral transformation Lanczos method for the sparse symmetric definite generalized eigenvalue problem for matrices $A$ and $B$ is an iterative method that addresses the case of semidefinite or ill conditioned $B$ using a shifted and…
We propose an optimization approach to the solution of the partial quadratic eigenvalue assignment problem (PQEVAP) for active vibration control design with robustness (RPQEVAP). The proposed cost function is based on the concept of…
The eigenstates of a real or complex cubic anharmonic oscillator are investigated using original and alternative methods. The procedure consists of determining global solutions of the Schr\"odinger equation that comply with the pertinent…
One of the most suitable methods for modeling fully dynamic earthquake cycle simulations is the spectral boundary integral element method (sBIEM), which takes advantage of the fast Fourier transform (FFT) to make a complex numerical dynamic…
Aussel et al. (J Optim Theory Appl 170 818-837 2016) introduced the concept of projected solutions for the quasi-variational inequalities with a non-self constraint map, that is, the case where the constraint map may take values outside the…
Polynomial filtering can provide a highly effective means of computing all eigenvalues of a real symmetric (or complex Hermitian) matrix that are located in a given interval, anywhere in the spectrum. This paper describes a technique for…
We propose a novel numerical algorithm for computing the electronic structure related eigenvalue problem of incommensurate systems. Unlike the conventional practice that approximates the system by a large commensurate supercell, our…
Plasmonics has attracted much attention not only because it has useful properties such as strong field enhancement, but also because it reveals the quantum nature of matter. To handle quantum plasmonics effects, ab initio packages or…
We describe an algorithm to compute the extremal eigenvalues and corresponding eigenvectors of a symmetric matrix by solving a sequence of Quadratic Binary Optimization problems. This algorithm is robust across many different classes of…
Numerous methods for computing conformal mesh paramterizations has been developed due to the vast applications in the field of geometry processing. Spectral conformal parameterization (SCP) is one of these methods to computing a quality…
Determining low-energy eigenstates in electronic many-body quantum systems is a key challenge in computational chemistry and condensed-matter physics. Hybrid quantum-classical approaches, such as the Variational Quantum Eigensolver and…
Systems with an effectively non-Hermitian Hamiltonian display an enhanced sensitivity to parametric and dynamic perturbations, which arises from the nonorthogonality of their eigenstates. This enhanced sensitivity can be quantified by the…
The spatial mode is an essential component of an electromagnetic field description, yet it is challenging to characterize it for optical fields with low average photon number, such as in a squeezed vacuum. We present a method for…
Modeling the internal structure of self-gravitating solid and liquid bodies presents a challenge, as existing approaches are often limited to either overly simplistic constant-density approximations or more complex numerical equations of…