Two-body Interaction Research Articles

Interaction-Induced Crystalline Topology of Excitons

We apply the topological theory of symmetry indicators to interaction-induced exciton band structures in centrosymmetric semiconductors. Crucially, we distinguish between the topological invariants inherited from the underlying electron and hole bands and those that are intrinsic to the exciton wave function itself. Focusing on the latter, we show that there exists a class of exciton bands for which the maximally localized exciton Wannier states are shifted with respect to the electronic Wannier states by a quantized amount; we call these excitons shift excitons. Our analysis explains how the exciton spectrum can be topologically nontrivial and sustain exciton edge states in open boundary conditions even when the underlying noninteracting bands have a trivial atomic limit. We demonstrate the presence of shift excitons as the lowest energy neutral excitations of the Su-Schrieffer-Heeger model in its trivial phase when supplemented by local two-body interactions. Published by the American Physical Society 2024

Landing Guidance Control Combining Powered Descending with Nonlinear Optimization and Vertical Descending with Model Prediction

In this study, we propose a landing guidance control method for the Moon using two control methods. The landing of a spacecraft on the Moon is divided into two phases: the powered descending phase and the vertical descending phase. The powered descending phase aims to optimize the entire trajectory to satisfy the termination constraint, instead of allowing for a relatively long control period. This phase performs feedback control using trajectory updates via nonlinear optimization. The vertical-descending phase aims to achieve accurate control over a short control period. This phase applies nonlinear model predictive control for fast optimization on finite time intervals. First, the landing of a spacecraft on the Moon is modeled as a two-body problem of the Moon and spacecraft, and equations of motion are derived. Subsequently, we formulated an optimization problem for each phase and developed the proposed landing guidance method by combining the two control methods. Furthermore, numerical simulations based on the derived equations of motion were performed to confirm the effectiveness of the proposed method and to compare its performance with another optimization method, the successive convexification (SCvx) algorithm.

Kepler’s problem of a two-body system perturbed by a third body

One of the most important problems in basic physics and astronomy is studying the motion of planets, satellites, and other celestial bodies. The solution to the two-body problem enables astronomers to predict the orbits of the Moon, satellites, and spaceships around the Earth. The general analytic solution for the three-body problem stands unsolved except in some special cases. This reduces the problem to a two-body problem. In this work, the authors present a closed-form approach to the three-body problem theoretically and numerically based on particle–particle vector analysis. The theoretical approach, which is based on the real Moon–Sun–Earth problem information, illustrates the perturbation of the Moon in the Sun–Earth problem and shows an expected orbital motion with a perturbation in the Sun–Earth orbit due to the revolution of the Moon. The numerical investigation uses the same information to study the same problem and calculate the angular momentums of each pair of objects. The two solutions show good agreement with the well-known Earth-Moon and Sun–Earth momentums. The Moon–Sun orbit is close to an elliptic shape with angular momentum of about 3.27 × 1038 J.s. This approach is the key to future studies for n-body problem solutions.

Spin Exchange-Enabled Quantum Simulator for Large-Scale Non-Abelian Gauge Theories

A central requirement for the faithful implementation of large-scale lattice gauge theories (LGTs) on quantum simulators is the protection of the underlying gauge symmetry. Recent advancements in the experimental realizations of large-scale LGTs have been impressive, albeit mostly restricted to Abelian gauge groups. Guided by this requirement for gauge protection, we propose an experimentally feasible approach to implement large-scale non-Abelian SU(N) and U(N) LGTs with dynamical matter in d+1D, enabled by two-body spin-exchange interactions realizing local emergent gauge-symmetry stabilizer terms. We present two concrete proposals for 2+1DSU(2) and U(2) LGTs, including dynamical bosonic matter and induced plaquette terms, that can be readily implemented in current ultracold-molecule and next-generation ultracold-atom platforms. We provide numerical benchmarks showcasing experimentally accessible dynamics, and demonstrate the stability of the underlying non-Abelian gauge invariance. We develop a method to obtain the effective gauge-invariant model featuring the relevant magnetic plaquette and minimal gauge-matter coupling terms. Our approach paves the way towards near-term realizations of large-scale non-Abelian quantum link models in analog quantum simulators. Published by the American Physical Society 2024

High-Throughput On-Demand Design Platform for Plasmonic Nanocavities: A Wavefunction Theory Approach.

Surface plasmon polaritons from plasmonic nanocavity have aroused great interest due to their applications in various fields, in which on-demand design is hindered by the lack of theoretical frameworks. Herein, based on its wave nature, we developed a wavefunction theory to explicitly describe individual surface plasmon polaritons and the resultant near-field and far-field behaviors, which serves as an efficient platform for high-throughput on-demand design of nanocavities. We found an applicative wavefunction form and proposed a two-body interaction function and a "shell" model for many-body interactions in surface plasmon polaritons' coupling. The wavefunction of individual surface plasmon polaritons and resultant near-field and far-field behaviors can be given explicitly and precisely. The theory provides a fundamental and quantitative understanding of surface plasmon polaritons and enables highly efficient on-demand design of plasmonic metamaterials and devices, leading to further methodological applications in numerous aspects.

Weighing Milky Way and Andromeda in an expanding ΛCDM Universe

The dynamics of the Local Group (LG), especially the contribution of the Milky Way (MW) and Andromeda (M 31) galaxies, is sensitive to the presence of dark energy. This work analyzes the evolution of the LG by considering it as a two-body problem in a homogeneous and isotropic expanding spacetime in a full Λcold dark matter (ΛCDM) background. Using the timing argument (TA), which links LG dynamics to LG mass, we find that the full ΛCDM background predicts a ∼10% lower mass for the LG; whereas Λ alone predicts a ∼10% higher mass. The TA mass is modified by (i) simulations and (ii) the effect of the Large Magellanic Cloud (LMC) to alleviate the poorly constrained internal mass distributions of M 31 and the MW, their time evolution, and the unknown distribution of dark matter between them. First, using IllustrisTNG simulations, we accounted for the effects of two extended halos and their environment (rather than point particles) and predicted their mass (3.89 ± 0.62)×1012 M⊙. Second, the LMC effectively changes the separation and velocities of M 31 towards the MW and reduces the predicted mass to (2.33 ± 0.72)×1012 M⊙. Despite the uncertainties around dark matter between these galaxies, the overall estimated mass is compatible with the mere sum of the MW and M 31 masses. The total mass of the TA is compatible with other estimates, such as the Hubble flow and the Virial Theorem with other dwarf galaxies. The combined result shows, for the first time, that a lower mass estimate can be obtained from the TA, with a consistent embedding and other systematic effects, and without an additional dark matter halo around the galaxies.

Classifying post-Minkowskian geometries for gravitational waves via loop-by-loop Baikov

We use the loop-by-loop Baikov representation to investigate the geometries in Feynman integrals contributing to the classical dynamics of a black-hole two-body system in the post-Minkowskian expansion of general relativity. These geometries determine the spaces of functions to which the corresponding Feynman diagrams evaluate. As a proof of principle, we provide a full classification of the geometries appearing up to three loops, i.e. fourth post-Minkowskian order, for all diagrams relevant to the conservative as well as the dissipative dynamics, finding full agreement with the literature. Moreover, we show that the non-planar top topology at four loops, which is the most complicated sector with respect to integration-by-parts identities, has an algebraic leading singularity and thus can only depend on non-trivial geometries through its subsectors.

Bertrand’s theorem and the double copy of relativistic field theories

Which relativistic field theories give rise to Kepler dynamics in the two-body problem? We consider a class of Hamiltonians that is the unique relativistic extension of the Kepler problem preserving its so(4) algebra, and have orbits related through time reparametrisation to orbits of the original Kepler problem. For three explicit examples, we give a natural interpretation in terms of spin-0,-1 and -2 interacting field theories in 5D. These are organically connected via the classical double copy, which therefore preserves maximal superintegrability.

Quark core-gluon model for heavy hybrid baryons

Besides the ordinary hadrons, quantum chromodynamics (QCD) allows the existence of states in which excitations of the gluonic field can play the role of valence particles, either alone in a glueball, or coupled to quarks in a hybrid. So, hybrid baryons, made of three quarks and a gluon, can exist. Till now, there is no experimental evidence for such exotic hadrons but experimental efforts are being made to search for them at CEBAF Large Acceptance Spectrometer. In this work, a hybrid baryon is considered as a two-body system composed of a color octet three-quark core and a gluon, interacting via a QCD-inspired interaction. A semirelativistic potential model is built in which the dominant interaction is a potential simulating the flux tube confinement, and the Casimir scaling is assumed to link interactions between triplet and octet color sources. This picture is similar to the quark-diquark description for baryons. It is chosen in order to take properly into account the helicity of the gluon. Only cccg and bbbg states are considered because the strong mass asymmetry between the quark core and the gluon is expected to favor the formation of the core. As the results for heavy hybrid baryons seem relevant, we consider this paper as a proof of concept which can be extended for the study of light hybrid baryons. Published by the American Physical Society 2024

The Instability Mechanism of Compact Multiplanet Systems

To improve our understanding of orbital instabilities in compact planetary systems, we compare suites of N-body simulations against numerical integrations of simplified dynamical models. We show that, surprisingly, dynamical models that account for small sets of resonant interactions between the planets can accurately recover N-body instability times. This points toward a simple physical picture in which a handful of three-body resonances, generated by interactions between nearby two-body mean motion resonances, overlap and drive chaotic diffusion, leading to instability. Motivated by this, we show that instability times are well described by a power law relating instability time to planet separations, measured in units of fractional semimajor axis difference divided by the planet-to-star mass ratio to the 1/4 power, rather than the frequently adopted 1/3 power implied by measuring separations in units of mutual Hill radii. For idealized systems, the parameters of this power-law relationship depend only on the ratio of the planets’ orbital eccentricities to the orbit-crossing value, and we report an empirical fit to enable quick instability time predictions. This relationship predicts that observed systems comprised of three or more sub-Neptune-mass planets must be spaced with period ratios P≳1.35 and that tightly spaced systems ( P≲1.5 ) must possess very low eccentricities (e ≲ 0.05) to be stable for more than 109 orbits.

Post-Minkowskian theory meets the spinning effective-one-body approach for two-body scattering

Effective-one-body (EOB) waveforms employed by the LIGO-Virgo-KAGRA Collaboration have primarily been developed by resumming the post-Newtonian expansion of the relativistic two-body problem. Given the recent significant advancements in post-Minkowskian (PM) theory and gravitational self-force formalism, there is considerable interest in creating waveform models that integrate information from various perturbative methods in innovative ways. This becomes particularly crucial when tackling the accuracy challenge posed by upcoming ground-based detectors (such as the Einstein Telescope and Cosmic Explorer) and space-based detectors (such as LISA, TianQin, or Taiji) expected to operate in the next decade. In this context, we present the derivation of the first spinning EOB (SEOB) Hamiltonian that incorporates PM results up to three-loop order: the SEOB-PM model. The model accounts for the complete hyperbolic motion, encompassing nonlocal-in-time tails. To evaluate its accuracy, we compare its predictions for the conservative scattering angle, augmented with dissipative contributions, against numerical-relativity data of nonspinning and spinning equal-mass black holes. We observe very good agreement, comparable, and in some cases slightly better to the recently proposed wEOB-potential model, of which the SEOB-PM model is a resummation around the probe limit. Indeed, in the probe limit, the SEOB-PM Hamiltonian and scattering angles reduce to the one of a test mass in Kerr spacetime. Once complemented with nonlocal-in-time contributions for bound orbits, the SEOB-PM Hamiltonian can be utilized to generate waveform models for spinning black holes on quasicircular orbits. Published by the American Physical Society 2024

Many-Body Expansion-Based Quantum Mechanical Force Field for Cyclotrimethylene Trinitramine under High Pressure.

RDX undergoes pressures of approximately 30-50 GPa during detonation, leading to significant changes in intermolecular interactions. Accurately describing these interactions is crucial for understanding the energy transfer in the detonation process. To address this, this work introduces a many-body expansion-based quantum mechanical force field (MB-QMFF) to accurately describe RDX's intermolecular interactions under high pressures. Using MB-QMFF, we evaluated various density functionals and found that the M062X functional with GD3 dispersion correction provided the highest accuracy. Regarding intermolecular forces, two-body interactions were the most significant, with three-body interactions being negligible. Additionally, we investigated intermolecular energy variations at different densities (or pressures). The results clearly demonstrate an accurate description of intermolecular interactions by the MB-QMFF scheme. Therefore, we believe that the MB-QMFF scheme can serve as a foundation for the development of RDX-specific force fields and pave the way for future studies on the detonation process of RDX.

Variational approach for the two-body problem in a multiband extended-Hubbard model

Considering a spin-up and a spin-down fermion in a generic tight-binding lattice with a multi-site basis, we investigate the two-body problem using a multiband extended-Hubbard model with finite-ranged hopping and interaction parameters. We derive a linear eigenvalue problem for the entire two-body spectrum, alongside a nonlinear eigenvalue problem for the bound states in the form of a self-consistency equation. Our results, based on an exact variational approach, suggest potential applications across various lattice geometries. As an illustration, we apply them to the linear-chain model and show that the resultant spin singlet and triplet bound states align well with the existing literature.

Density matrix renormalization group study of the interacting Kitaev chain with quasi-periodic disorder

We document the ground state phase diagram of the one-dimensional Kitaev chain with quasi-periodic disorder in the presence of two-body interactions. Our data were obtained for systems of L=1000\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L=1000$$\\end{document} sites using large-scale density-matrix renormalization group numerics and is benchmarked against known results for the clean system. We demonstrate that moderate quasi-periodic disorder stabilizes the topological phase both for repulsive and attractive interactions. For larger disorder strengths, the system features re-entrance behavior and multiple phase transitions.Graphical abstractPhase diagram as a function of the chemical potential and the disorder strength for repulsive interactions

Combining Gaia and GRAVITY: Characterising five new directly detected substellar companions

Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which – detected at an angular separation of (34.01 ± 0.15) mas from the host – is the closest substellar companion ever imaged. In combination with the system’s distance and the orbital elements, this translates to a semi-major axis of (0.938 ± 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 ± 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia’s upcoming DR4 and the advent of GRAVITY+.

Two mass-imbalanced atoms in a hard-wall trap: Deep learning integrability of many-body systems.

The study of integrable systems has led to significant advancements in our understanding of many-body physics. We design a series of numerical experiments to analyze the integrability of a mass-imbalanced two-body system through energy-level statistics and deep learning of wave functions. The level spacing distributions are fitted by a Brody distribution and the fitting parameter ω is found to separate the integrable and nonintegrable mass ratios by a critical line ω=0. The convolutional neural network built from the probability density images could identify the transition points between integrable and nonintegrable systems with high accuracy, yet in a much shorter computation time. A brilliant example of the network's ability is to identify a new integrable mass ratio 1/3 by learning from the known integrable case of equal mass, with a remarkable network confidence of 98.22%. The robustness of our neural networks is further enhanced by adversarial learning, where samples are generated by standard and quantum perturbations mixed in the probability density images and the wave functions, respectively.

Limits on an improved action for contact effective field theory in two-body systems

We consider a possible resummation of subleading effects in two-body systems with a large scattering length as described by a short-range effective field theory (EFT). In particular, we investigate the consequences of a resummation of part of the range corrections. Explicit calculations of the two-body phase shifts and charge form factor indicate that, except for extreme choices, resummations do not alter the convergence of the EFT expansion and are often beneficial at the lowest orders. We have considered the expansion when the regulator cutoff is removed as well as when it is finite, and find that the cutoff is not an important factor for resummations. Our results connect with other works where the partial resummation is induced by potentials with finite cutoffs or interaction ranges.

One-Dimensional Relativistic Self-Gravitating Systems.

One of the oldest problems in physics is that of calculating the motion of N particles under a specified mutual force: the N-body problem. Much is known about this problem if the specified force is non-relativistic gravity, and considerable progress has been made by considering the problem in one spatial dimension. Here, I review what is known about the relativistic gravitational N-body problem. Reduction to one spatial dimension has the feature of the absence of gravitational radiation, thereby allowing for a clear comparison between the physics of one-dimensional relativistic and non-relativistic self-gravitating systems. After describing how to obtain a relativistic theory of gravity coupled to N point particles, I discuss in turn the two-body, three-body, four-body, and N-body problems. Quite general exact solutions can be obtained for the two-body problem, unlike the situation in general relativity in three spatial dimensions for which only highly specified solutions exist. The three-body problem exhibits mild forms of chaos, and provides one of the first theoretical settings in which relativistic chaos can be studied. For N≥4, other interesting features emerge. Relativistic self-gravitating systems have a number of interesting problems awaiting further investigation, providing us with a new frontier for exploring relativistic many-body systems.

Extended Lambert-Based Approach for Ground-Track Adjustment with Specified Overflight Altitude

This paper presents an extended Lambert-based approach for ground-track adjustment under J2 perturbation, aiming to achieve a satellite flyby over a designated ground location at a predetermined altitude. Extended equations are developed in the orbital plane based on Lagrange’s form for the unperturbed Lambert’s problem to address the J2-perturbed two-body motion. By combining them with the out-of-plane flight time equation constraint regarding the ground-track adjustment, a nonlinear system of equations with two unknowns and two equations is established and split into different modes. To solve these modes, a modified Newton–Raphson method is used, and a necessary condition is established to determine the maximum allowable number of revolutions. To enhance computational efficiency, a radius restriction method is proposed to narrow down the allowable range of the number of revolutions. The proposed algorithm is applicable to various scenarios, including single- or multiple-day flights, ascending or descending overflights, coplanar or non-coplanar transfers, and initial or reverse-initial directional transfers. The simulation results demonstrate that the proposed algorithm is effective and accurate for the perturbed adjustment problem.

Viscous Current-Induced Forces.

We study the motion (translational, vibrational, and rotational) of a diatomic impurity immersed in an electron liquid and exposed to electronic current. An approach based on the linear response time-dependent density functional theory combined with the Ehrenfest dynamics leads to a system of linear algebraic equations, which account for the competing and counteracting effects of the current-induced force (electron wind) and the electronic friction. We find and emphasize the coupling between the center of mass motion and that of the nuclei relative to each other, the feature due to the mediation of the two-body interaction by the environment. The current-induced forces, by means of the dynamic exchange-correlation (xc) kernel f_{xc}(r,r^{'},ω), include the electronic viscosity contribution. Starting from the ground state at the equilibrium internuclear distance and applying a current pulse, we observe three phases of the motion: (i)acceleration due to the prevalence of the current-induced force, (ii)stabilization upon balancing of the two forces, and (iii)deceleration due to the friction after the end of the pulse. At lower, but still metallic, electron densities, the dynamic xc contribution to the force significantly affects the acceleration (deceleration) at the first (third) phase of the process. For the Cs density (r_{s}≈6 a.u.), this correction amounts up to 40% in the rotation regime.