Skip to content
Dr. Denys Dutykh

Dr. Denys Dutykh

Associate Professor · Associate Dean of Graduate Studies
Applied Mathematics

Applied mathematician working on water waves, tsunami modeling and spectral numerical methods. Since 2024 these tools also serve general relativity: quasinormal modes of black holes and wormholes.

∇²u = ∂²u/∂t²
∫∫∫ ∇ × F · dS
λ = 2πc/ω

Recent Publications

View all
CComput. Appl. Math.
Comput. Appl. Math.

Error estimation for numerical approximations of ODEs via composition techniques. Part I: one-step methods

In this study, we introduce a refined method for estimating errors in numerical simulations of dynamical systems through an innovative application of composition techniques. Our approach involves a dual application: a basic one-step numerical method of order p in this part, and a class of Backward Difference Formulas (BDF) schemes in Part II (Deeb et al. 2026). This dual application uses complex coefficients, resulting in outputs in the complex plane. The method’s innovation lies in demonstrating that the real parts of these outputs correspond to approximations of the solutions with an enhanced order of p + 1 , while the imaginary parts serve as error estimates of the same order, a novel proof presented herein. The linear stability of the resulting scheme is improved over that of the basic one. The performance of the composition in computing the approximation is also compared. The results show that the proposed technique attains higher accuracy with reduced computational time relative to the basic integrators; compared with established methods of the same order it remains competitive in cost while additionally supplying a built-in error estimate, and it is most advantageous for integrators that lack a native error estimator. This dual composition technique has been rigorously applied to a variety of dynamical problems, demonstrating its efficacy in adapting the time step, particularly in situations where numerical schemes lack theoretical error estimates. Consequently, the technique has the potential to advance adaptive time-stepping strategies in numerical simulations.

A. Deeb, D. Dutykh2027
View
AAnn. Phys.
Ann. Phys.

An exact Kerr-like rotating Morris–Thorne wormhole: Throat, ergoregion and equatorial shadow slice

We construct an exact Kerr-like rotating extension of the zero-redshift Morris–Thorne wormhole within a Teo-type stationary and axisymmetric ansatz. Imposing the Einstein equations for an anisotropic source comoving with zero angular momentum observers fixes the circumferential metric function to K(r) = r² + a² and reduces the frame-dragging function to a single radial quadrature. For b(r) = r₀²/r, this quadrature is evaluated in terms of incomplete elliptic integrals and normalised by the ADM angular momentum. The spacetime is asymptotically flat but has vanishing ADM mass, so the Kerr relation a = J_ADM/M_ADM is not applicable. The field equations leave the Kerr-like oblateness length a and the physical angular momentum J_ADM independent. The exact family is therefore labelled by (r₀, a, J_ADM). The relation J_ADM = a r₀ is used only to select a one-parameter slice for the numerical ergoregion and equatorial-capture illustrations, and not as a field-equation constraint. We obtain r_th = √(r₀² − a²), the reality bound |a| ⩽ r₀, and the regular range 0 ⩽ |a| < r₀. Within the canonical stationary foliation, the throat is characterised quasi-locally as the unique member S₀ of the closed two-surface family S_ℓ with vanishing mean curvature and positive area second variation. The surfaces r = r_th and ℓ = 0 are only coordinate representations of this surface. A signed throat-adapted coordinate displays two asymptotically flat ends, excludes closed timelike curves, and shows that the throat is timelike rather than a horizon. These geometrical and causal conclusions do not require J_ADM = a r₀. The supporting anisotropic stress tensor is reconstructed from the Einstein tensor and is interpreted here as an effective phenomenological source. No microscopic matter Lagrangian or realistic equation of state is claimed. We also determine the ergoregion onset. Because the full Hamilton–Jacobi equation is not separable, the optical result is only the one-dimensional equatorial capture interval, not the complete two-dimensional shadow.

M. E. Sukaiti, D. Batic, D. Dutykh2026
View
PPhys. Rev. D
Phys. Rev. D

Quasinormal Modes of Gauss–Bonnet Black Holes via the Spectral Method: Scalar, Vector, and Tensor Perturbations

We present a unified study of scalar, vector, and tensor quasinormal modes (QNMs) of Schwarzschild black holes corrected by a Gauss—Bonnet (GB) term in higher dimensions. Using a high-precision Chebyshev spectral method, we map the QNM spectra across $D\in\{5,6,7,8,10,11,12,26\}$ well beyond the regime where sixth-order WKB and characteristic-integration techniques remain reliable. Across the three spin sectors, we find several robust signatures of higher-curvature dynamics: the appearance of overdamped purely imaginary modes, non-monotonic behaviour in the real parts of higher overtones, and a strong amplification of the dimensionless QNM frequencies in string-motivated dimensions. In the scalar and vector sectors, we uncover an exact isospectrality between the scalar monopole ($\ell=0$) and vector dipole ($\ell=1$) at vanishing GB coupling, and we provide an analytic proof based on a Darboux factorisation of the corresponding Hamiltonians. In the tensor sector, we obtain the first numerical confirmation of the long-predicted instability in six dimensions; its onset is sharply captured by the Cohn—Calogero bound and leads to the mass threshold $GM\leq 158.1\,α^{3/2}$. No analogous instability is found for $D\geq 7$, and no tensor isospectrality occurs. Converting the dimensionless frequencies to physical units suggests that the amplified modes in higher dimensions may enter the sensitivity window of future space-based detectors such as DECIGO. The merged analysis provides a comprehensive benchmark for QNMs in Einstein—Gauss—Bonnet gravity and highlights the limitations of standard approximation schemes in the strong-coupling and high-overtone regimes.

D. Batic, D. Dutykh2026
View

Latest Blog Posts

View all
dds@spy
Cloning into 'academic-web-genesis'...
✔ Receiving objects: 100% (842/842), done.
Switched to a new branch 'feature/terminal-component'
[feature/terminal-component 7fa91b3] feat: add interactive terminal component3 files changed, 125 insertions(+)
✔ Branch 'feature/terminal-component' set up to track remote branch.
ℹ Create a pull request for 'feature/terminal-component' on GitHub:https://github.com/dutykh/academic-web-genesis/pull/new
✔ Packages installed successfullyDependencies: 57, devDependencies: 28
▲ Next.js 15.3.4 (Turbopack)- Local: http://localhost:3000✓ Ready in 875ms
Dr. Denys Dutykh
Dr. Denys Dutykh
Dr. Denys Dutykh
Dr. Denys Dutykh
Dr. Denys Dutykh
Dr. Denys Dutykh
Dr. Denys Dutykh

Dr. Denys Dutykh

Associate Professor of Applied Mathematics

Always learning, always building. 

Stay Updated

Subscribe to receive notifications about new blog posts and research updates