Papers

We examine two major factors that influence the optimum projection angle: aerodynamic forces and the effect of run-up. With respect to aerodynamics, we consider not only the drag but also the lift generated by spin during flight. By linearizing the equations of motion that include these forces, we derive perturbation solutions with respect to drag and lift coefficients and clarify their qualitative effects. The results show that both drag and lift reduce the optimum projection angle, with the latter exerting a stronger influence. To investigate the effect of run-up, we use an extended projection model in which the initial speed depends on the initial angle. Analysis of this model reveals that a stronger run-up increases the relative projection angle but decreases the launch angle observed from the ground. These findings provide a mechanical explanation for the release angle in shot put and the takeoff angle in long jump. The present study establishes a simple theoretical framework for clarifying the respective roles of aerodynamic and run-up effects in determining the optimum projection angles in sports.

The paper presents a consistent thermomechanical framework for strongly nonlocal continua, a type of generalized media in which a material's response at a point depends on deformation and temperature gradients within its neighborhood. Such dependence is responsible for localization phenomena having an intrinsic size, and nonlocal modeling also serves as a practical regularization tool to prevent mesh dependence. Nonlocality is of the integral type, using a kernel function that gives the relative influence of neighboring points on a central point. The paper develops three sources of nonlocality: (1) nonlocal momentum balance equations; (2) nonlocal conservation laws (including the first and second laws of thermodynamics); and (3) a material's nonlocal constitutive structure. Balance equations are derived using the principle of virtual power, permitting non-smooth force and displacement fields and non-smooth boundary surfaces with indistinct normal directions. The balance equations differ from those of classical local continua. Conservation laws also differ from those of local continua, with the mechanical power and heating at a point being averaged over its neighborhood. These laws are developed for processes that are sufficiently slow to minimize additional meso-scale kinetic energy due to internal turbulence. Stress, entropy, and dissipative forces are obtained as averaged derivatives of the free energy function with respect to strain, temperature, and internal variables. An example is presented of a stretched elastoplastic bar with a small defect, and nonlocality is shown to impart a characteristic size to the deformation pattern.

Thermal radiation is discussed for a classical scalar field. Dilation and acceleration appear in conformal symmetry, and so are treated in connection with zero-point radiation and equilibrium thermal radiation. Classical zero-point radiation is completely invariant under dilation in a inertial frame. However, in a Rindler frame, zero-point radiation takes a different form which is related to Planck's thermal spectrum with temperature T related to acceleration. Traditionally, zero-point radiation is said to take a thermal (Unruh) form in a Rindler frame. However, the classical viewpoint indicates that zero-point radiation in a Rindler frame requires a divergent frequency spectrum and so is not thermal radiation. The Unruh form suggests an additional classical thermal spectrum satisfying the Stefan-Boltzmann relation.

Large plastic deformation can drive an initially uniform single crystal to spontaneously subdivide into misoriented grains separated by thin dislocation walls – a pattern-forming instability rooted in the loss of convexity of the crystal's elastic energy at large strain. We study this phenomenon for a ductile crystal in plane-strain simple shear within continuum dislocation theory, using a polyconvex (Ciarlet–Geymonat) elastic energy that guarantees existence of minimizers for the coupled deformation–slip problem. Minimizing over the plastic slip yields a condensed energy of double-well form whose non-quasiconvexity favours a lamellar microstructure; the gradient of the geometrically necessary dislocation density regularizes it, giving the grain boundaries a finite thickness and energy as functions of the misorientation angle. A block-coordinate finite element scheme – alternating a convex non-smooth solve for the slip with a Levenberg-regularized Newton solve for the deformation – resolves this microstructure numerically and detects its spontaneous onset, reproducing the lamellar grain structure in agreement with the closed-form analysis.

We solve in closed form a one-dimensional relativistic system: two masses interacting only through elastic collisions with a massless mediator bouncing between them. Momenta, times, and positions are hyperbolic functions of the collision index. The mediator energy, interpreted as the pair's effective potential, obeys an exact discrete Coulomb law, $V\propto 1/r$, with a Lorentz-invariant action as coupling. A massive Newtonian mediator instead transmits a $1/r^{3}$ force; one adiabatic invariant traces both laws to the mediator's dispersion relation. Continued to negative mediator energy, the closed forms turn trigonometric, binding a one-dimensional mechanical analog of the Coulomb atom.

We apply a thin metal rod (copper rod) as a probe in the alternating magnetic field generated by a Helmholtz coil, and measure the variation of the induced electromotive force (EMF) on the metal rod at different radial positions of Helmholtz coil. Experimental results show that the induced EMF is zero when the center of the metal rod passes through the center of the cylindrical magnetic field inside the Helmholtz coil. When the metal rod is displaced from the center of field to different radial positions, the induced EMF gradually increases from zero, reaches a maximum at a certain position, and then decreases monotonically as the radial distance continues to increase. At the position where the induced EMF reaches its maximum, the metal rod intersects the radial cross-section of the internal magnetic field of the Helmholtz coil at two points, with a small central portion of the rod located inside the Helmholtz coil and the two end portions outside the coil. To explain the experimental phenomena, we construct four simplified models of the magnetic field distribution based on the actual field distribution of the Helmholtz coil to quantitatively investigate the radial variation of the induced EMF along the metal rod. Our theoretical results show that the four models yield similar results and can all qualitatively explain the experimental data curves, particularly reproducing well the variation trend of the induced EMF and the position of the extremum point. The result of piecewise function fitting model is quantitatively in good agreement with the experimental data. This work is also very much helpful and instructive for undergraduate-level students.

Incremental contact laws specify the normal contact force through a differential equation carrying an internal state, driven by the indentation and its rate. In some, the force is extinguished at a nonzero indentation, whether by plastic deformation or by an elastic aftereffect, so that a residual deformation remains at the separation. Such laws sit uneasily within rigid body dynamics, which admits no deformation. The tension is tolerable when the indentation is small relative to the bodies, so that it may be carried constitutively rather than geometrically. Even then, the contact law alone does not determine the interaction of the bodies. Because force and indentation no longer vanish together, conditions for the commencement and termination of contact must be supplied separately. So must the fate of the deformation and internal state at separation, neither of which the equations of motion contain. This article formulates the repeated direct collinear impact of two convex bodies under external forces as a hybrid dynamical system. The contact interface is modeled as a massless element carrying the contact law and its state, coupled to the bodies through relative velocity and an interaction force dictated by the contact state. Consequently, all switching and resets are confined to the interface, leaving the geometry and the equations of motion of the bodies unaltered. The principal analytical properties of the resulting formulations are established, among them passivity and completeness; the branching of solutions at the onset and termination of contact is also examined. The framework is demonstrated through numerical simulations.

The electromagnetic vacuum stress tensor in a system consisting of classical point-like atoms on a cubic lattice is derived directly from the electric field and, alternatively, using the known more phenomenological formula involving the (usually unknown) derivative of the dielectric tensor. The results agree. The stress tensor contains a nontrivial constant and therefore cannot be obtained from macroscopic electrodynamics alone. A system consisting of a mosaic of randomly oriented grains of a cubic crystal is used as a model of an amorphous solid. The electromagnetic stress tensor in this system is deduced by averaging the corresponding quantity in a cubic crystal over all orientations. Two components remain: the component in the direction of the electric field and the component perpendicular to it. The transverse component agrees with the analogous quantity for liquids derived by Peierls in an entirely different way.

We present a Hessian-based approach to numerically continue periodic orbits in dynamical systems. A loop (periodic orbit candidate) is parametrized as a Fourier series; a loss function is defined based on the deviation of the loop from the physical differential equations. Unlike previous work relying on hand-derived Jacobians, our method automates the process by leveraging automatic differentiation, a common machine learning technique. The continuation direction can be determined by the flat directions of the loss landscapes (directions with zero eigenvalues), making the search of periodic orbits efficient and guided. Our method is integrator-free, precisely initializes oscillations around unstable fixed points, and efficiently detects orbit family intersections and subharmonic bifurcations. As a demonstration, we present full continuations of periodic double pendulum oscillations from fixed points, showing bifurcations along orbit families and categorizing branches of periodic orbits. In particular, we find periodic orbits where both pendulum masses are never simultaneously at rest, which to our knowledge has been missing in the literature.

An alternative, streamlined methodology is presented for deriving the isotropic Cartesian tensor bases under the special orthogonal group $\text{SO}(3)$. By shifting the traditional interpretation of the isotropy condition to analyze it as an algebraic system of equations for the rotation matrices themselves rather than the tensor components, lower-rank bases can be explicitly established in just a few lines without requiring finite coordinate rotations or complicated contractions of infinitesimal generators. Furthermore, a transparent combinatorial interpretation is provided for the number of tensors in the general higher-rank spanning sets. Finally, the Gram-matrix method is advocated as an efficient computational sieve to resolve the subsequent linear dependence.

A symmetric distribution of charges or currents typically produces an electromagnetic field having the same symmetry group. However, under special circumstances we show that a discrete symmetry of the source distribution leads to a larger continuous symmetry for the resulting field in restricted regions. The conditions for this promotion from discrete to continuous symmetry are that the fields be affine linear and the symmetry of the source distribution act irreducibly on the ambient space. We illustrate this effect with many examples. This phenomenon is pedagogically interesting both as a source of inspiration in electro- and magnetostatics and as an accessible example of group representation theory.

In his seminal paper, Bell [Physics Physique Fizika 1, 195 (1964)] considers the correlations that result from space-like separated measurements on a pair of entangled particles. He uses relativity theory to motivate the Bell causal structure, then shows the existence of quantum correlations that cannot be explained classically within this causal structure. Classical explanations of such quantum correlations are possible in alternative causal structures, for instance, those that allow superluminal causal influences, but, as shown in [New Journal of Physics 17 033002 (2015)], all such alternative explanations require fine tuning (causation without correlation). Here we discuss the existence of spurious quantum correlations — correlations that look quantum in one causal structure, but have a natural classical explanation in another. More precisely, there are causal structures that admit non-classical quantum correlations, but for which the same correlations have a classical explanation in another causal structure without fine tuning. The realisation in the other causal structure can also be achieved without breaking any natural constraints on the causal structure that follow from relativity theory. However, similarly to non-classical quantum correlations in the Bell causal structure, we find other causal structures with non-classical quantum correlations that do not have a classical causal explanation in any alternative causal structure without fine tuning.

Classical macroscopic electrodynamics allows free and bound sources to cancel as distributions while the auxiliary equations retain nonzero free-source terms. A neutralized uniformly polarized sphere has $\mathbf{E}=\mathbf{0}$ and $\mathbf{D}=\mathbf{P}$, and a compensated uniformly magnetized sphere has $\mathbf{B}=\mathbf{0}$ and $\mathbf{H}=-\mathbf{M}$. These auxiliary fields are fixed by $\mathbf{P}$ and $\mathbf{M}$ rather than forming independent degrees of freedom. Time dependence separates charge cancellation from complete four-current cancellation. Canceling the charge alone leaves a divergence-free total current whose on-shell transverse transform controls radiation. For a neutralized polarized sphere of radius $R$ maintained by a tangential free-current sheet, the exterior field equals that of a point electric dipole whose effective moment is proportional to the spherical Bessel function $j_2(kR)$, where $k=\omega/c$. Every charge multipole and the ordinary magnetic dipole vanish, yet the source radiates at generic frequencies, and the complete exterior field vanishes at nonzero roots of $j_2$. Complete four-current cancellation removes all retarded source-generated fields. Among the compactly supported separable currents that maintain the same charge, it is the only choice that never radiates. In vacuum the auxiliary fields reduce to $\mathbf{D}=\varepsilon_0\mathbf{E}$ and $\mathbf{H}=\mathbf{B}/\mu_0$, so no radiation can be carried by $\mathbf{D}$ and $\mathbf{H}$ alone.

A charged test particle is confined to the round sphere in a uniform ambient magnetic field. We reduce its motion to quadrature by elementary means and solve explicitly the family of trajectories which reaches the poles. On that family the azimuth advances at a constant rate and the latitude obeys a pendulum equation. The winding gained per oscillation is a complete elliptic integral of the first kind. That integral is a strictly increasing bijection onto the positive numbers. Every positive rational winding is therefore carried by exactly one value of the half-cyclotron frequency. Then we ask which geometry has these trajectories. No affine connection has them among its geodesics. A Finsler metric of Randers type has all of them. The threshold which separates the two dynamical regimes is precisely the condition for that metric to exist. The closed trajectories are its closed geodesics. The example places a non-Riemannian geometry within reach of an upper-level undergraduate or beginning graduate course in classical mechanics. We supply the derivations, the numerical recipes and the exercises.

Relative periodic orbits (RPOs) are solutions of the three-body problem that are periodic in a uniformly rotating reference frame and, in general, quasi-periodic in inertial coordinates. We present a numerical procedure for computing and continuing one-parameter families of RPOs of the planar Newtonian three-body problem. The method exploits consecutive syzygies, understood here as configurations in which the three bodies are aligned and their velocities satisfy the corresponding symmetry conditions. Matching the positions and momenta at two consecutive alignments reduces the computation of RPOs to a low-dimensional nonlinear problem. Its solutions are then numerically continued, and linear stability is determined from the nontrivial eigenvalues of the rotated monodromy matrix after removing the neutral directions associated with conserved quantities and continuous symmetries. The procedure is applied to several mass distributions and initial configurations, producing families of Poincaré, Hill, and binary-type solutions. These families exhibit transitions from nearly circular to highly eccentric motion, changes of stability near resonances and turning points, and absolute periodic solutions when the rotation angle is a rational multiple of 2{\pi}. In the Hill families, the continuation connects satellite configurations with circumstellar motion as the smallest body loses its gravitational binding to the intermediate body. Circumbinary and circumstellar configurations are also obtained in the binary regime. The results illustrate the dynamical diversity of RPOs and provide coherent three-body motions that can be used as prescribed trajectories in restricted four-body models.

The classical problem of a rigid body in an ideal fluid is a six-dimensional dynamical system with three conserved quantities. Despite its long history and its importance as a reduced model of many fluid-structure problems, basic questions remain, including the shapes and connectivities of the three-dimensional submanifolds of solutions, and the nature of chaotic motions on these. We explore the problem with a particular ellipsoidal body, beginning with integrable motions— steady linear translation and rotation and periodic planar tumbling and fluttering. The stability of these states depends on the relative magnitudes of linear and angular momentum or energy. Linear stability, including Floquet, analysis of lower-dimensional perturbed systems is consistent with direct integration of the full system. We document instabilities leading to a variety of regular and chaotic motions, including flipping and twirling, whose trajectories appear to follow paths near connections between the integrable states. Such indirect observation of these connections provides insight into the structure underlying the rich dynamics of this simple system. For the submanifold where linear and angular momentum are perpendicular, we further examine how solutions fit together to fill momentum space in different zones of stability of the integrable states, and thereby uncover additional bifurcations in connectivity and what appear to be new integrable solutions.

We investigate the statistical properties of the non-hierarchical planar three-body problem, relevant to astrophysical systems with nearly planar dynamics, including protoplanetary disks and AGN disks. Using the Flux-based statistical theory of the three-body problem, previously studied in three dimensions, we derive theoretical predictions for several statistical observables in the planar case. Focusing on ergodic disintegration outcomes, we compare these predictions directly with numerical simulations and contrast the results with the corresponding unconstrained three-dimensional case. Our analysis is based on one million simulations for each of eight distinct mass sets. We examine escape probabilities, properties of marginal escape events, eccentricity distributions, lifetime distributions, and relative prevalence of prograde and retrograde escapes. While several of the theoretical predictions are in good agreement with the simulations, others exhibit unexpected discrepancies. We discuss possible explanations for these deviations.

Most drivers slow down for speed bumps, sometimes to almost a stop, as if each speed bump were a stop or yield sign. The late Roy Schwitters and I shared the view (here called the Schwitters Doctrine) that one should take the bumps at a speed of at least 60 or 70 kph (35 to 45 mph), not because we were scofflaws, but because we found, based on informal experiments, that the traverse was less jarring at high speeds than at low. Here, the problem is worked in terms of the geometry of the wheels on the road and two coupled, damped harmonic oscillators representing the suspension of the car and the driver's cushioned seat. The Schwitters Doctrine is recovered, but not in the way that one expects.

This paper presents a self-contained, vector-only formulation of gyroscope precession about a fixed pivot – gravity-driven, not the torque-free Euler-Poinsot case. The analysis uses only vector integrals, dot products, and cross products of intrinsic vectors: no tensor or dyadic formalism, no matrix representation, no Euler angles, and no component decomposition. Inertia enters through two mass integrals alone: the polar second moment of mass $J_0$ about the center of mass, and the vector operator $\vec F(\vec v)$, which carries directional inertia in place of an inertia tensor. Using material time derivatives referenced to the inertial observer, the torque about the pivot is assembled once from the complete acceleration of each mass element, giving an exact, unified vector equation for steady precession. The balance of the gravity torque against this inertia torque reproduces the classical steady-precession relation of the heavy symmetric top exactly at non-vertical tilts, on both the slow and fast branches, with the horizontal-shaft case exact rather than approximate. Closed forms are given for the disk, sphere, and circular and rectangular toruses. The constant tilt angle of dynamic equilibrium is explained as torque equalization: the two torques about the pivot balance exactly at that tilt and unequally on either side of it, treated throughout as a geometric function of the tilt under constant precession and self-spin rates, and never as a function of time. The formulation is informationally equivalent to the classical description; its contribution is methodological: economy of definition, a single algebra throughout, and geometric visualization preserved. The exact steady-precession relation predates this work; Brand obtained it in 1930. First derived by the author in 1975 and re-verified symbolically, it is intended for undergraduate readers, with every step retained.

Spinor structure is usually introduced as part of relativistic quantum theory. We show that its underlying algebra can instead emerge without quantisation from a covariant statistical description of a massive relativistic particle subject to an explicit finite-dimensional linear factorisation of the complete massive mass-shell constraint motivated by requiring the theory retain both sheets of the massive mass shell. A finite-dimensional factorisation of the complete mass-shell constraint, linear in all four momentum components, forces a Clifford algebra. In (3+1) dimensions its minimal complex representation is four-dimensional, while each on-shell factor selects a rank-two sector; in this representation statistical completeness requires an arbitrary $2\times2$ matrix distribution within that sector before quantisation. Projecting the Weyl-ordered matrix Liouvillian gives relativistic transport while preserving all internal populations and coherences. For universal deformations on flat canonical phase space satisfying the stated covariance, matrix-neutrality, homogeneity and associativity assumptions, the resulting product is the matrix Weyl–Moyal product. Its parameter has dimensions of action and, when identified with $\hbar$, supplies the scale of physical spin-$\tfrac12$. Requiring all-order constraint preservation gives a vanishing star commutator, while preservation of the selected star-spectral sector yields the two-sided Dirac–Wigner equations. Our central conclusion is therefore that relativistic mass-shell factorisation supplies spinor structure, while completion of the constraint-preservation hierarchy within the stated deformation class leads to the Weyl–Moyal phase-space formulation of quantum mechanics.

Magnetically coupled damped rotational degrees of freedom can give rise to tribological loss when excited by a relative sliding motion, linking friction to the dynamics of magnetic moments. We analytically investigate such magnetic friction by studying a recently introduced simplified model for a rigid magnetic rotor array sliding over a commensurate magnetic substrate. Here, the array rotors can rotate about an axis perpendicular to the sliding direction and are damped by microscopic shaft friction, while the substrate magnets have a fixed in-plane orientation. The simplified model reduces the collective rotor dynamics of this setup to a set of coupled nonlinear differential equations, which can be studied by the corresponding energy landscape in the quasi-static regime. We find that the parameter plane spanned by layer separation and time can be divided into distinct chambers, in which the structure of the energy landscape remains invariant. Using these chambers, we identify the layer separation intervals corresponding to different dynamical states and recover the emergence of a regime where the moment alignment alternates globally while exhibiting peak dissipation. Additionally, we derive the leading nontrivial orders of the collective magnetic friction in the regimes of small and large gaps between the array and the substrate. For the alternating regime, we find that friction can be computed by tracking energy jumps from unstable to stable critical points at the chamber boundaries. These results provide an analytical foundation for the observed dynamics and dissipation in sliding rotor arrays. Because the simplified model is scale-free, our conclusions transfer across a broad range of microscopic and macroscopic length scales.

Magnetic levitation is one of the most captivating phenomena in physics, combining a striking visual appeal with the rich underlying physics of forces and interactions in equilibrium. Bringing it into undergraduate laboratories in a quantitative way is nevertheless often hindered by the intrinsic instability of static magnetic systems and by the high cost of active control systems. Here we present a simple, safe, and very low-cost experiment that uses ring-shaped ferrite magnets together with a mechanical stabilization system fabricated by 3D printing in a male–female configuration. The apparatus lets undergraduate students explore directly and quantitatively the relationship between the magnitude of the magnetic force and the separation distance. We assess the robustness of the device through experimental tests, whose results show excellent quantitative agreement with the equivalent current-loop model. Furthermore, we validate this model through a complementary, independent first-principles-based methodology based on axial magnetic field profiling. To support reproducibility, the 3D design files are made available in a public, open-access repository.

Nanoelectromechanical actuators operating at sub-100-nm gaps collapse through a pull-in instability set by competing electrostatic and Casimir forces. The quasi-static fold that bounds their safe operating range has been known in closed form for three decades, together with the Casimir ceiling above which no static equilibrium survives, which fixes the smallest gap a given stiffness and area can hold open against the quantum vacuum. That fold does not give the threshold reached from rest, its dependence on damping, or the design sensitivities of either. We train a physics-informed neural network in a rapidity coordinate that maps the movable pull-in pole to infinity, which keeps the residual bounded across the collapse threshold where fixed-step Runge-Kutta integration steps into unphysical states. Differentiating the trained surrogate returns pull-in-voltage sensitivities that match the closed-form fold to a relative error of $3\times10^{-6}$ and inverts a device specification to a gap of 97.036 nm at a target actuation voltage. Applied to the from-rest boundary, which carries no closed form once the damping is finite, it supplies the same sensitivities where no analytic root exists. We prove that this boundary is bracketed by two closed-form curves, that it is nondecreasing in the damping ratio, that it merges with the fold once the damping ratio exceeds $2^{-1/4}$, and that the gap closes as $(\tau_*-\tau)^{2/5}$. Numerically the merger already occurs at $0.396$, and the growth of the collapse time changes there from logarithmic to inverse square root. The classical-limit bound on the thermal Lifshitz derating is at the percent level, and the physical shift at these gaps lies orders of magnitude below it.

Within the Clifford algebra of spacetime, a covariant electrodynamics is formulated in which the derivative of the four-potential contains a scalar field (S) in addition to the electromagnetic bivector. A generalized system of Maxwell equations is derived that reduces to the classical Maxwell system when S=0. The energy-momentum tensor, the corresponding conservation laws, and a generalized Lorentz force density are obtained. It is shown that, in free space, the field S admits wave solutions, termed scalar S-waves. The structure and energy characteristics of a free S-wave, as well as its possible sources and the conditions for its experimental detection, are determined.

We clarify some mathematical features of the density of a dilute particle flow on a body surface in a two-fluid model, with a view toward its application to icing simulations. Inviscid solutions in the vicinity of a front stagnation point are obtained in exact and perturbed forms for small and large Stokes numbers, respectively, and the behavior of the density of the dispersed phase is clarified over the entire range of Stokes numbers. The viscous effect of the laminar boundary layer in the carrier phase on the particle-free thin layer appearing in the low-Stokes-number regime is also investigated based on the incompressible Navier-Stokes equations. An order-of-magnitude estimation of the velocity field of the dispersed phase reveals the cause of this layer and the criterion for its appearance. In particular, the latter result provides an improved form of Michael's criterion.

We demonstrate a physical mechanism for causal information filtering in a physical reservoir computing (PRC) by exploiting asymmetric coupling anisotropy. Using a network of coupled Duffing oscillators, we show that the directionality of internal coupling induces a spatial gradient in the effective potential, establishing a deterministic upstream-to-downstream information flow. This anisotropy allows for the selective amplification of semantic drifts, triggering a macroscopic saddle-node bifurcation as a physical interlock before global computational failure. Through spatiotemporal analysis of a 50-node system under traveling wave inputs, we confirm that local phase transitions effectively purge anomalous information while preserving the computational integrity of the remaining nodes. The results suggest that the intrinsic causality of the reservoir's topology provides a robust framework for autonomous reliability and fault-tolerant physical intelligence.

We consider the motion of a charged test particle confined to the round unit sphere in a uniform ambient magnetic field. We use an extension of Noether's theorem for systems with magnetic forces to reduce the problem to a quadrature. The rotation number of a trajectory, the azimuth it gains over one oscillation in latitude, is a complete elliptic integral of the third kind in Legendre form. A trajectory closes if and only if that advance is a rational multiple of a full turn. The criterion holds on every level set of the two first integrals. Then, we differentiate the rotation number with respect to the half-cyclotron frequency and sign that derivative on each of the two branches into which a level set divides. The rotation number sees the charge, the mass, the field and the speed through a single dimensionless ratio, that of the half-cyclotron frequency to the speed. A closed trajectory therefore fixes a value of that ratio. We say how many values carry a given rational winding. The poles are attainable on a single level set, where the motion reduces to a pendulum. The field has no flux through the sphere, so it is globally exact. We compute the Mañé strict critical value of the resulting exact magnetic flow. Above that value the closed trajectories are closed Reeb orbits of an explicit contact form.

Measurements of Newton's gravitational constant remain limited by environmental and apparatus-dependent systematics whose treatment is often less explicit than that of instrumental noise. Among these, atmospheric density fluctuations generate unshieldable gravity gradients that couple directly to torsion-balance observables as atmospheric Newtonian noise. Here, we develop a GUM-consistent framework for propagating this contribution through the torque estimator into the uncertainty budget of torsion-balance measurements. We derive closed-form spatial transfer functions for two benchmark geometries: a two-mass dumbbell, which provides an upper-coupling reference, and a perfectly symmetric cross, which serves as an idealized rejection limit and exposes the trade-off between suppressing low-order environmental coupling and preserving signal response. We also separate stationary correlated inputs from non-stationary baseline drift, using the Ornstein–Uhlenbeck process only as a benchmark for the former. Atmospheric pressure benchmarks indicate that the resulting background contribution is below present reference uncertainty levels, but can become relevant as systematic floors approach the part-per-million regime. This framework provides a practical route for incorporating site-specific environmental gravity gradients into future torsion-balance uncertainty budgets.

Can a source radiate when its total charge density vanishes identically? Consider a uniformly polarized sphere coated with free surface charge that cancels the bound surface charge at every point and time. Then $\rho_{\mathrm{tot}}=0$, so every electric charge multipole vanishes. Continuity determines only $\nabla\cdot\mathbf{J}$, which allows the same charge history to be supported by different conserved currents. A compensating interior current gives $\mathbf{J}_{\mathrm{tot}}=\mathbf{0}$ and produces no $\mathbf{E}$ or $\mathbf{B}$ at any frequency. The minimum-norm tangential sheet current instead leaves $\mathbf{J}_{\mathrm{tot}}$ nonzero and divergence-free. Its radiation-zone field is exact in $kR$ and proportional to $j_2(kR)$. At long wavelength the radiated power is suppressed by $(kR)^4/100$, and it vanishes exactly at the positive roots of $j_2$. The same calculation gives the interior field and a closed-form energy balance. The average work supplied by driving equals the radiated power and falls to zero at those roots even though interior fields remain. For comparison, the bare sphere has the factor $3j_1(kR)/(kR)$ and is silent at the roots of $j_1$.

The thermodynamics of cavity radiation is usually introduced through the quantum theory of a photon gas. In contrast, we investigate which features of thermal radiation follow from classical electrodynamics alone. We show, using two complementary approaches, that the proportionalities between energy density, radiation pressure, and energy flux arise without invoking photons or quantum statistics. The first approach is based on explicit solutions of Maxwell's equations and spatio-temporal averaging of electromagnetic waves. The second approach develops the classical statistical mechanics of the electromagnetic field from first principles by constructing the canonical Hamiltonian, implementing the gauge constraints, formulating the partition function as a functional integral, and evaluating the corresponding field correlation functions. By combining wave analysis, statistical mechanics, constrained Hamiltonian dynamics, and functional methods within a single classical framework, this work provides a coherent classical formulation of electromagnetic radiation thermodynamics.