Susanna Terracini
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Susanna Terracini has a total of 159 co-authors. Here are the top co-authors by number of shared publications:
- Gianmaria Verzini
Dipartimento di Matematica, Politecnico di Milano, Milano, Italy
27 shared publications - Veronica Felli
Università degli Studi di Milano-Bicocca, Milano, Italy
22 shared publications - Nicola Soave
Universita degli Studi di Torino, Italy
17 shared publications - Hugo Tavares
Departamento de Matemática do Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
15 shared publications - Stefano Vita
Dipartimento di Matematica “F. Casorati”, Università di Pavia, Pavia, Italy
10 shared publications - Vivina Barutello
Dipartimento di Matematica, Università degli Studi di Torino, Turin, Italy
10 shared publications
Latest Publications
- 16 September 2026
- Published by European Mathematical Society - EMS - Publishing House GmbH in Annales de l'Institut Henri Poincaré C, Analyse non linéaire
Chaotic Boltzmann billiard systems at positive energy
This paper deals with the Boltzmann billiard , that is, a billiard subjected to a central force of the type V(r)=-\alpha/r-\beta/r^{2} , \alpha and \beta being real constants, and with a straight reflection table. In the particular case of \alpha and \beta positive, we prove the presence of a symbolic dynamics, and hence of positive topological entropy, at positive energy for \beta sufficiently small.
- 16 September 2026
- Published by European Mathematical Society - EMS - Publishing House GmbH in Annales de l'Institut Henri Poincaré C, Analyse non linéaire
Chaotic Boltzmann billiard systems at positive energy
This paper deals with the Boltzmann billiard , that is, a billiard subjected to a central force of the type V(r)=-\alpha/r-\beta/r^{2} , \alpha and \beta being real constants, and with a straight reflection table. In the particular case of \alpha and \beta positive, we prove the presence of a symbolic dynamics, and hence of positive topological entropy, at positive energy for \beta sufficiently small.
- Open Access
- 4 August 2026
- Published by Springer Nature in Archive for Rational Mechanics and Analysis
- Vol. 250 (5) , 75
On the Regularity of Solutions to the Hamilton–Jacobi Equations for the N-Body Problem
We prove that certain renormalized value functions associated with the d -dimensional ( $$d\geqq 2$$ ) N -body problem corresponding to different limiting shapes of expanding solutions, under the assumption that the center of mass is at the origin, are viscosity solutions of the associated Hamilton–Jacobi equation. We analyze their singularities, defined as the initial configurations for which the minimizer of the associated variational problem is not unique. Moreover, we estimate the size of the closure of the singular set by proving its $$\mathcal {H}^{d(N-1)-1}$$ -rectifiability, and we provide an upper bound on the Hausdorff dimension of the set of regular conjugate points.- 4 August 2026
- Published by Springer Nature in Archive for Rational Mechanics and Analysis
- Vol. 250 (5) , 75
Open AccessOn the Regularity of Solutions to the Hamilton–Jacobi Equations for the N-Body Problem
We prove that certain renormalized value functions associated with the d -dimensional ( $$d\geqq 2$$ ) N -body problem corresponding to different limiting shapes of expanding solutions, under the assumption that the center of mass is at the origin, are viscosity solutions of the associated Hamilton–Jacobi equation. We analyze their singularities, defined as the initial configurations for which the minimizer of the associated variational problem is not unique. Moreover, we estimate the size of the closure of the singular set by proving its $$\mathcal {H}^{d(N-1)-1}$$ -rectifiability, and we provide an upper bound on the Hausdorff dimension of the set of regular conjugate points. - Open Access
- 27 March 2026
- Published by Springer Nature in La Matematica
- Vol. 5 (2) , 25
Bifurcation for a Three-Component Competition System
- 27 March 2026
- Published by Springer Nature in La Matematica
- Vol. 5 (2) , 25
Open AccessBifurcation for a Three-Component Competition System
- Open Access
- 30 October 2025
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 406 (12) , 1-40
Frozen Planet Orbits for the n-Electron Atom
We investigate periodic trajectories in a classical system consisting of multiple mutually repelling electrons constrained to move along a half-line, with an attractive nucleus fixed at the origin. Adopting a variational framework, we seek critical points of the associated Lagrangian action functional using a modified Lusternik-Schnirelmann theory for manifolds with boundary. Furthermore, in the limit where the electron charges vanish, we demonstrate that frozen planet orbits converge to segments of a brake orbit in a Kepler-type problem, thereby drawing a strong analogy with Schubart orbits in the gravitational N -body problem.
- 30 October 2025
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 406 (12) , 1-40
Open AccessFrozen Planet Orbits for the n-Electron Atom
We investigate periodic trajectories in a classical system consisting of multiple mutually repelling electrons constrained to move along a half-line, with an attractive nucleus fixed at the origin. Adopting a variational framework, we seek critical points of the associated Lagrangian action functional using a modified Lusternik-Schnirelmann theory for manifolds with boundary. Furthermore, in the limit where the electron charges vanish, we demonstrate that frozen planet orbits converge to segments of a brake orbit in a Kepler-type problem, thereby drawing a strong analogy with Schubart orbits in the gravitational N -body problem.
- Open Access
- 7 August 2025
- Published by Elsevier in Communications in Nonlinear Science and Numerical Simulation
- Vol. 152
Equivariant optimisation for the gravitational n-body problem: A computational factory of symmetric orbits
- 7 August 2025
- Published by Elsevier in Communications in Nonlinear Science and Numerical Simulation
- Vol. 152
Open AccessEquivariant optimisation for the gravitational n-body problem: A computational factory of symmetric orbits
- Open Access
On the regularity of solutions to the Hamilton-Jacobi equations for the N-body problem
We prove that certain suitably renormalized value functions associated with the $d$-dimensional ($d\geq2$) $N$-body problem corresponding to different limiting shapes of expanding solutions, under the assumption that the center of mass is at the origin, are viscosity solutions of the associated Hamilton-Jacobi equation. We analyze their singularities, defined as the initial configurations for which the minimizer of the associated variational problem is not unique. Moreover, we estimate the size of the closure of the singular set by proving its $\mathcal{H}^{d(N-1)-1}$-rectifiability, and we provide an upper bound on the Hausdorff dimension of the set of regular conjugate points.
Open AccessOn the regularity of solutions to the Hamilton-Jacobi equations for the N-body problem
We prove that certain suitably renormalized value functions associated with the $d$-dimensional ($d\geq2$) $N$-body problem corresponding to different limiting shapes of expanding solutions, under the assumption that the center of mass is at the origin, are viscosity solutions of the associated Hamilton-Jacobi equation. We analyze their singularities, defined as the initial configurations for which the minimizer of the associated variational problem is not unique. Moreover, we estimate the size of the closure of the singular set by proving its $\mathcal{H}^{d(N-1)-1}$-rectifiability, and we provide an upper bound on the Hausdorff dimension of the set of regular conjugate points.
- Open Access
A priori Hölder estimates for equations degenerating on nodal sets
We prove a priori Hölder bounds for continuous solutions to degenerate equations with variable coefficients of type $$ \mathrm{div}\left(u^2 A\nabla w\right)=0\quad\mathrm{in \ }Ω\subset\mathbb R^n,\qquad \mbox{with}\qquad \mathrm{div}\left(A\nabla u\right)=0, $$ where $A$ is a Lipschitz continuous, uniformly elliptic matrix (possibly $u$ has non-trivial singular nodal set). Such estimates are uniform with respect to $u$ in a class of normalized solutions that have a bounded Almgren frequency. As a consequence, a boundary Harnack principle holds for the quotient of two solutions vanishing on a common set. This analysis relies on a detailed study of the associated weighted Sobolev spaces, including integrability of the weight, capacitary properties of the nodal set, and uniform Sobolev inequalities yielding local boundedness of solutions.Open AccessA priori Hölder estimates for equations degenerating on nodal sets
We prove a priori Hölder bounds for continuous solutions to degenerate equations with variable coefficients of type $$ \mathrm{div}\left(u^2 A\nabla w\right)=0\quad\mathrm{in \ }Ω\subset\mathbb R^n,\qquad \mbox{with}\qquad \mathrm{div}\left(A\nabla u\right)=0, $$ where $A$ is a Lipschitz continuous, uniformly elliptic matrix (possibly $u$ has non-trivial singular nodal set). Such estimates are uniform with respect to $u$ in a class of normalized solutions that have a bounded Almgren frequency. As a consequence, a boundary Harnack principle holds for the quotient of two solutions vanishing on a common set. This analysis relies on a detailed study of the associated weighted Sobolev spaces, including integrability of the weight, capacitary properties of the nodal set, and uniform Sobolev inequalities yielding local boundedness of solutions. - 22 July 2025
- Published by European Mathematical Society - EMS - Publishing House GmbH
- p. 219-247
Bifurcation for the Lotka–Volterra competition model
We analyze the bifurcation phenomenon for the following two-component competition system: $$ \begin{cases} -\Delta u\_1=\mu u\_1(1-u\_1)-\beta \alpha u\_1u\_2,& \text{in } B\_1\subset \mathbb{R}^N,\\ -\Delta u\_2=\sigma u\_2(1-u\_2)-\beta \gamma u\_1u\_2,& \text{in } B\_1\subset \mathbb{R}^N,\\ \frac{\partial u\_1}{\partial n}= \frac{\partial u\_2}{\partial n} =0,&\text{on } \partial B\_1, \end{cases} $$ where $N\ge 2$, $\alpha>\gamma>0$, $\sigma\ge\mu>0$, and $\beta>\frac{\sigma}{\gamma}$. More precisely, treating $\beta$ as the bifurcation parameter, we initially perform a local bifurcation analysis around the positive constant solutions, obtaining precise bifurcation locations and determining the direction of bifurcation. As a byproduct, the instability of the constant solution is provided. Furthermore, we extend our exploration to the global bifurcation analysis. Lastly, under the condition $\sigma=\mu$, we demonstrate the limiting configuration on each bifurcation branch as the competition rate $\beta\rightarrow+\infty$.- 22 July 2025
- Published by European Mathematical Society - EMS - Publishing House GmbH
- p. 219-247
Bifurcation for the Lotka–Volterra competition model
We analyze the bifurcation phenomenon for the following two-component competition system: $$ \begin{cases} -\Delta u\_1=\mu u\_1(1-u\_1)-\beta \alpha u\_1u\_2,& \text{in } B\_1\subset \mathbb{R}^N,\\ -\Delta u\_2=\sigma u\_2(1-u\_2)-\beta \gamma u\_1u\_2,& \text{in } B\_1\subset \mathbb{R}^N,\\ \frac{\partial u\_1}{\partial n}= \frac{\partial u\_2}{\partial n} =0,&\text{on } \partial B\_1, \end{cases} $$ where $N\ge 2$, $\alpha>\gamma>0$, $\sigma\ge\mu>0$, and $\beta>\frac{\sigma}{\gamma}$. More precisely, treating $\beta$ as the bifurcation parameter, we initially perform a local bifurcation analysis around the positive constant solutions, obtaining precise bifurcation locations and determining the direction of bifurcation. As a byproduct, the instability of the constant solution is provided. Furthermore, we extend our exploration to the global bifurcation analysis. Lastly, under the condition $\sigma=\mu$, we demonstrate the limiting configuration on each bifurcation branch as the competition rate $\beta\rightarrow+\infty$.- Open Access
New solutions for the symmetrical n-body problem through variational approach and optimisation techniques
Advances in the variational approach to the $n$-body problem have led to significant progress in celestial mechanics, uncovering new types of possible orbits. In this paper, critical points of the Lagrangian action associated with the $n$-body problem are analysed using evolutionary algorithms to identify periodic and symmetrical solutions of the discretised system. A key objective is to locate minimum points of the action functional, as these correspond to feasible periodic solutions that satisfy the system's differential equations. By employing both stochastic and deterministic algorithms, we explore the solution space and obtain numerical representations of these orbits. Next, we examine the stability of these orbits by treating them as critical points. One approach is to compute their discrete Morse index to distinguish between minimum points and saddle points. Another is to classify them based on their action levels. Finally, analysing the boundaries of their attraction basins allows us to identify non-minimal critical points via the Ambrosetti-Rabinowitz Mountain Pass Theorem. This leads to an updated version of the algorithm that provides a constructive proof of the theorem, yielding new orbits in specific cases. This paper builds upon and extends the results presented in \cite{nostro}, providing a more detailed theoretical framework and deeper insights into the formulation. Additionally, we present new numerical results and an extended analysis of the critical points found, further enhancing the findings of the previous study.
Open AccessNew solutions for the symmetrical n-body problem through variational approach and optimisation techniques
Advances in the variational approach to the $n$-body problem have led to significant progress in celestial mechanics, uncovering new types of possible orbits. In this paper, critical points of the Lagrangian action associated with the $n$-body problem are analysed using evolutionary algorithms to identify periodic and symmetrical solutions of the discretised system. A key objective is to locate minimum points of the action functional, as these correspond to feasible periodic solutions that satisfy the system's differential equations. By employing both stochastic and deterministic algorithms, we explore the solution space and obtain numerical representations of these orbits. Next, we examine the stability of these orbits by treating them as critical points. One approach is to compute their discrete Morse index to distinguish between minimum points and saddle points. Another is to classify them based on their action levels. Finally, analysing the boundaries of their attraction basins allows us to identify non-minimal critical points via the Ambrosetti-Rabinowitz Mountain Pass Theorem. This leads to an updated version of the algorithm that provides a constructive proof of the theorem, yielding new orbits in specific cases. This paper builds upon and extends the results presented in \cite{nostro}, providing a more detailed theoretical framework and deeper insights into the formulation. Additionally, we present new numerical results and an extended analysis of the critical points found, further enhancing the findings of the previous study.
- Open Access
- 5 May 2025
- Published by Springer Nature in Aerotecnica Missili & Spazio
- Vol. 105 (2) , 153-170
New Solutions for the Symmetrical n-Body Problem Through Variational Approach and Optimisation Techniques
Advances in the variational approach to the n -body problem have led to significant progress in celestial mechanics, uncovering new types of possible orbits. In this paper, critical points of the Lagrangian action associated with the n -body problem are analysed using evolutionary algorithms to identify periodic and symmetrical solutions of the discretised system. A key objective is to locate minimum points of the action functional, as these correspond to feasible periodic solutions that satisfy the system’s differential equations. By employing both stochastic and deterministic algorithms, we explore the solution space and obtain numerical representations of these orbits. Next, we examine the stability of these orbits by treating them as critical points. One approach is to compute their discrete Morse index to distinguish between minimum points and saddle points. Another is to classify them based on their action levels. Finally, analysing the boundaries of their attraction basins allows us to identify non-minimal critical points via the Ambrosetti–Rabinowitz Mountain Pass Theorem. This leads to an updated version of the algorithm that provides a constructive proof of the theorem, yielding new orbits in specific cases. This paper builds upon and extends the results presented in Introna et al. (75th International Astronautical Conference, Milan, Italy), providing a more detailed theoretical framework and deeper insights into the formulation. Additionally, we present new numerical results and an extended analysis of the critical points found, further enhancing the findings of the previous study.
- 5 May 2025
- Published by Springer Nature in Aerotecnica Missili & Spazio
- Vol. 105 (2) , 153-170
Open AccessNew Solutions for the Symmetrical n-Body Problem Through Variational Approach and Optimisation Techniques
Advances in the variational approach to the n -body problem have led to significant progress in celestial mechanics, uncovering new types of possible orbits. In this paper, critical points of the Lagrangian action associated with the n -body problem are analysed using evolutionary algorithms to identify periodic and symmetrical solutions of the discretised system. A key objective is to locate minimum points of the action functional, as these correspond to feasible periodic solutions that satisfy the system’s differential equations. By employing both stochastic and deterministic algorithms, we explore the solution space and obtain numerical representations of these orbits. Next, we examine the stability of these orbits by treating them as critical points. One approach is to compute their discrete Morse index to distinguish between minimum points and saddle points. Another is to classify them based on their action levels. Finally, analysing the boundaries of their attraction basins allows us to identify non-minimal critical points via the Ambrosetti–Rabinowitz Mountain Pass Theorem. This leads to an updated version of the algorithm that provides a constructive proof of the theorem, yielding new orbits in specific cases. This paper builds upon and extends the results presented in Introna et al. (75th International Astronautical Conference, Milan, Italy), providing a more detailed theoretical framework and deeper insights into the formulation. Additionally, we present new numerical results and an extended analysis of the critical points found, further enhancing the findings of the previous study.