CyberSec Research

We present a pedagogical introduction to Floquet-Magnus theory through the classical example of Kapitza's pendulum \ -- a simple system exhibiting nontrivial dynamical stabilization under rapid periodic driving. By deriving the equations of motion and analyzing the system via Floquet theory and the Magnus expansion, we obtain analytical stability conditions and effective evolution equations. While grounded in classical mechanics, the techniques are directly applicable to periodically driven quantum systems as well. The approach is fully analytical, using only tools from theoretical mechanics and linear algebra, and is suitable for advanced undergraduate or graduate students.

The formation of planetesimals from cm-sized pebbles in protoplanetary disks faces significant barriers, including fragmentation and radial drift. We identify a previously unaccounted screening force, arising from mutual shielding of thermal gas particles between pebbles when their separation falls below the gas mean free path. This force facilitates pebble binding, overcoming key growth barriers under turbulent disk conditions. Unlike conventional mechanisms, screening forces operate independently of surface adhesion and complement streaming instability and pressure traps by enhancing aggregation in high-density regions. Our analysis predicts that screening interactions are most effective in the {middle disk regions ($ \sim 0.3$ to few AU),} consistent with ALMA observations (e.g., TW Hya) of enhanced dust concentrations. {Furthermore, we find that screening-induced pebble growth from centimeter to kilometer scales can occur on timescales significantly shorter than the disk lifetime ($\sim 10^5$ years). Importantly, this growth naturally terminates when particles smaller than the local gas mean free path are depleted, thereby avoiding runaway accretion.} Beyond planetary science, the screening forces have {potential} implications for high-energy astrophysics, dusty plasmas, confined particle suspensions and other relevant areas, suggesting a broader fundamental significance.

We derive an analytic solution for the electromagnetic vector potential in any gauge directly from Maxwell's equations for potentials for an arbitrary time-dependent charge-current distribution. No gauge condition is used in the derivation. Our solution for the vector potential has a gauge-invariant part and a gauge-dependent part. The gauge-dependent part is related to the scalar potential.

Braiding has attracted significant attention in physics because of its important role in describing the fundamental exchange of particles. Infusing the braiding with topological protection will make it robust against imperfections and perturbations, but such topological braiding is believed to be possible only in interacting quantum systems, e.g., topological superconductors. Here, we propose and demonstrate a new strategy of topological braiding that emerges from non-Abelian topological insulators, a class of recently discovered multi-band topological phase. We unveil a mathematical connection between braiding and non-Abelian quaternion invariants, by which Bloch eigenmodes under parallel transport produce braid sequences protected by the non-Abelian band topology. The braiding is also associated with geometric phases quantized over half the Brillouin zone. This new type of non-Abelian topological braiding is experimentally realized in acoustic systems with periodic synthetic dimensions. The results show that the principle discovered here is a new strategy towards topological braiding and can be extended for other types of classical waves and non-interacting quantum systems.

We introduce a mechanical control system for energy efficient and robust hoisting crane operations. The control system efficiently translates the harmonic motion of a spring loaded mediating system into the desired driving of the load, recyling most of the employed energy for subsequent operations. The control output is a shortcut-to-adiabaticity protocol borrowed from quantum mechanics. The control system reduces the single operation consuption in realistic working regimes, but it is in cyclical processes where the energetical advantage becomes substantial. The design of the control system and the control output is flexible enough to allow additional optimization of the robustness against perturbations.

The orbit of the S2 star around Sagittarius A* provides a unique opportunity to test general relativity and study dynamical processes near a supermassive black hole. Observations have shown that the orbit of S2 is consistent with a Schwarzschild orbit at a 10$\sigma$ confidence level, constraining the amount of extended mass within its orbit to less than 1200 M$_\odot$, under the assumption of a smooth, spherically symmetric mass distribution. In this work we investigate the effects on the S2 orbit of granularity in the mass distribution, assuming it consists of a cluster of equal-mass objects surrounding Sagittarius A*. Using a fast dynamical approach validated by full N-body simulations, we perform a large set of simulations of the motion of S2 with different realizations of the cluster objects distribution. We find that granularity can induce significant deviations from the orbit in case of a smooth potential, causing precession of the orbital plane and a variation of the in-plane precession. Interactions with the cluster objects also induce a sort of "Brownian motion" of Sagittarius A*. Mock data analysis reveals that these effects could produce observable deviations in the trajectory of S2 from a Schwarzschild orbit, especially near apocenter. During the next apocenter passage of S2 in 2026, astrometric residuals in Declination may exceed the astrometric accuracy threshold of GRAVITY of about 30 $\mu as$, as it happens in 35 to 60% of simulations for black holes of 20 to 100 M$_\odot$. This presents a unique opportunity to detect, for the first time, scattering effects on the orbit of S2 caused by stellar-mass black holes, thanks to the remarkable precision achievable with GRAVITY. We also demonstrate that any attempt to constrain the extended mass enclosed within the orbit of S2 must explicitly account for granularity in the stellar-mass black hole population.

Composite Pulses (CPs) are widely used in Nuclear Magnetic Resonance (NMR), optical spectroscopy, optimal control experiments and quantum computing to manipulate systems that are well-described by a two-level Hamiltonian. A careful design of these pulses can allow the refocusing of an ensemble at a desired state, even if the ensemble experiences imperfections in the magnitude of the external field or resonance offsets. Since the introduction of CPs, several theoretical justifications for their robustness have been suggested. In this work, we suggest another justification based on the classical mechanical concept of a stability matrix. The motion on the Bloch Sphere is mapped to a canonical system of coordinates and the focusing of an ensemble corresponds to caustics, or the vanishing of an appropriate stability matrix element in the canonical coordinates. Our approach highlights the directionality of the refocusing of the ensemble on the Bloch Sphere, revealing how different ensembles refocus along different directions. The approach also clarifies when CPs can induce a change in the width of the ensemble as opposed to simply a rotation of the axes. As a case study, we investigate the $90(x)180(y)90(x)$ CP introduced by Levitt, where the approach provides a new perspective into why this CP is effective.

This paper studies brachistochrone trajectories. Four rules are formulated as sufficient conditions. Two rules apply for a general conservative force. Two rules apply for a central force. A central force allows wire replacement. The wire is replaced by appropriate magnetic field. This enables solving motion equations directly. We replace Euler Lagrange with direct integration.

Gravitational systems in astrophysics often comprise a body -- the primary -- that far outweights the others, and which is taken as the centre of the reference frame. A fictitious acceleration, also known as the indirect term, must therefore be added to all other bodies in the system to compensate for the absence of motion of the primary. In this paper, we first stress that there is not one indirect term but as many indirect terms as there are bodies in the system that exert a gravitational pull on the primary. For instance, in the case of a protoplanetary disc with two planets, there are three indirect terms: one arising from the whole disc, and one per planet. We also highlight that the direct and indirect gravitational accelerations should be treated in a balanced way: the indirect term from one body should be applied to the other bodies in the system that feel its direct gravitational acceleration, and only to them. We point to situations where one of those terms is usually neglected however, which may lead to spurious results. These ideas are developed here for star-disc-planets interactions, for which we propose a recipe for the force to be applied onto a migrating planet, but they can easily be generalized to other astrophysical systems.

This work introduces a Hamiltonian approach to regularization and linearization of central force particle dynamics through a new canonical extension of the so-called "projective decomposition". The regularization scheme is formulated within the framework of classic analytical Hamiltonian dynamics as a redundant-dimensional canonical/symplectic coordinate transformation, combined with an evolution parameter transformation, on extended phase space. By considering a generalized version of the standard projective decomposition, we obtain a family of such canonical transformations which differ at the momentum level. From this family of transformations, a preferred canonical coordinate set is chosen that possesses a simple and intuitive connection to the particle's local reference frame. Using this transformation, closed-form solutions are readily obtained for inverse square and inverse cubic radial forces (or any superposition thereof) on any finite-dimensional Euclidean space. From these solutions, a new set of orbit elements for Kepler-Coulomb dynamics is derived, along with their variational equations for arbitrary perturbations (singularity-free in all cases besides rectilinear motion). Governing equations are numerically validated for the classic two-body problem, incorporating the J_2 gravitational perturbation.

The study guide (textbook) is part of a set of materials designed to support high-quality practical training in physics. It includes a collection of tasks for organizing both in-class and independent work. The guide serves as a foundation for further study in physics-related disciplines and aligns with current educational programs. This textbook presents a curated set of 120 physics problems with detailed solutions, structured according to the first-year bachelor's curriculum. Each section addresses common student questions and emphasizes conceptual understanding. Problem-solving is essential in physics education. It not only tests knowledge but also transforms theory into practical skills. Applying physical laws to real-world scenarios enhances comprehension and fosters analytical thinking. Through solving problems, students gain deeper insight into physical phenomena and develop effective strategies for analysis, and develop solutions to tasks-making the learning process truly comprehensive.

Symplectic numerical methods have become a widely-used choice for the accurate simulation of Hamiltonian systems in various fields, including celestial mechanics, molecular dynamics and robotics. Even though their characteristics are well-understood mathematically, relatively little attention has been paid in general to the practical aspect of how the choice of coordinates affects the accuracy of the numerical results, even though the consequences can be computationally significant. The present article aims to fill this gap by giving a systematic overview of how coordinate transformations can influence the results of simulations performed using symplectic methods. We give a derivation for the non-invariance of the modified Hamiltonian of symplectic methods under coordinate transformations, as well as a sufficient condition for the non-preservation of a first integral corresponding to a cyclic coordinate for the symplectic Euler method. We also consider the possibility of finding order-compensating coordinate transformations that improve the order of accuracy of a numerical method. Various numerical examples are presented throughout.

The Looping pendulum phenomenon was first introduced in 2019 at the 32nd edition of the IYPT, wherein a lighter bob sweeps around a cylindrical rod to support the weight of a heavier bob. In this paper, the phenomenon was divided based on rotating and non-rotating forces, and differential equations were derived for each. To verify the theoretical derivation, an experimental analysis was done, varying the mass ratio with the vertical distance travelled by the heavier bob. (Tracked using tracker) Experimental findings fit a logarithmic curve fit -- falling succinctly with a similar trend with the simulation run with MATLAB solving the derived differential equations. Furthermore, to verify the simulation, the trajectory of both the lighter and heavier mass was also compared for the simulation and experimental findings. The experimental findings fit very closely to the simulation findings, accrediting the validity and accuracy of the derived theory.

A linear elastic circular disc is analyzed under a self-equilibrated system of loads applied along its boundary. A distinctive feature of the investigation, conducted using complex variable analysis, is the assumption that the material is incompressible (in its linearized approximation), rendering the governing equations formally identical to those of Stokes flow in viscous fluids. After deriving a general solution to the problem, an isoperimetric constraint is introduced at the boundary to enforce inextensibility. This effect can be physically realized, for example, by attaching an inextensible elastic rod with negligible bending stiffness to the perimeter. Although the combined imposition of material incompressibility and boundary inextensibility theoretically prevents any deformation of the disc, it is shown that the problem still admits non-trivial solutions. This apparent paradox is resolved by recognizing the approximations inherent in the linearized theory, as confirmed by a geometrically nonlinear numerical analysis. Nonetheless, the linear solution retains significance: it may represent a valid stress distribution within a rigid system and can identify critical conditions of interest for design applications.

We present a detailed analysis of all possible regular precessions of a heavy asymmetric body with a fixed point not coinciding with the center of mass. The calculations are done in terms of the rotation matrix, by writing the Euler-Poisson equations with all involved vectors parameterized in the Laboratory frame. It is shown that a regular precession is possible if the suspension point is chosen on the straight lines (lying in the principal plane) which are frontiers of the regions where, as the distance from the center of mass increases, the permutation of intermediate and largest moments of inertia occurs. Like the spin of an electron in quantum mechanics, the frequency of regular precession in classical mechanics turns out to be rigidly fixed by two values, i.e., quantized.

This paper is concerned with non-Hermitian degeneracy and exceptional points associated with resonances in an acoustic scattering problem with sound-hard obstacles. The aim is to find non-Hermitian degenerate (defective) resonances using numerical methods. To this end, we characterize resonances of the scattering problem as eigenvalues of a holomorphic integral operator-valued function. This allows us to define defective resonances and associated exceptional points based on the geometric and algebraic multiplicities. Based on the theory on holomorphic Fredholm operator-valued functions, we show fractional-order sensitivity of defective resonances with respect to operator perturbation. This property is particularly important in physics and associated with intriguing phenomena, e.g., enhanced sensing and dissipation. A defective resonance is sought based on the perturbation analysis and Nystr\"om discretization of the boundary integral equation. Numerical evidence of the existence of a defective resonance is provided. The numerical results combined with theoretical analysis provide a new insight into novel concepts in non-Hermitian physics.

We provide a statistical and correlational analysis of the spatial and energetic properties of equilibrium configurations of a few-body system of two to eight equally charged classical particles that are confined on a one-dimensional helical manifold. The two-body system has been demonstrated to yield an oscillatory effective potential, thus providing stable equilibrium configurations despite the repulsive Coulomb interactions. As the system size grows, the number of equilibria increases, approximately following a power-law. This can be attributed to the increasing complexity in the highly non-linear oscillatory behavior of the potential energy surface. This property is reflected in a crossover from a spatially regular distribution of equilibria for the two-body system to a heightened degree of disorder upon the addition of particles. However, in accordance with the repulsion within a helical winding, the observed interparticle distances in equilibrium configurations cluster around values of odd multiples of half a helical winding, thus maintaining an underlying regularity. Furthermore, an energetic hierarchy exists based on the spatial location of the local equilibria, which is subject to increasing fluctuations as the system size grows.

In this work, we extend a fractional-dimensional space model for anisotropic solids by incorporating a q-deformed derivative operator, inspired by Tsallis' nonadditive entropy framework. This generalization provides an analytical framework for exploring anisotropic thermal properties, within a unified and flexible mathematical formalism. We derive modified expressions for the phonon density of states and specific heat capacity, highlighting the impact of the deformation parameters on thermodynamic behavior. We apply the model to various solid-state materials, achieving excellent agreement with experimental data, across a wide temperature range and demonstrating its effectiveness in capturing anisotropic and subextensive effects in real systems.