Search

Lunar ejecta, produced by meteoroidal impacts, have been proposed for the origin of the near-Earth asteroid (469219) Kamo'oalewa, supported by its unusually Earth-like orbit and L-type reflectance spectrum (Sharkey et al., 2021). In a recent study (Castro-Cisneros et al. 2023), we found with N-body numerical simulations that the orbit of Kamo'oalewa is dynamically compatible with rare pathways of lunar ejecta captured into Earth's co-orbital region, persistently transitioning between horseshoe and quasi-satellite (HS-QS) states. Subsequently, Jiao et al. (2024) found with hydrodynamic and N-body simulations that the geologically young lunar crater Giordano Bruno generated up to 300 Kamo'oalewa-sized escaping fragments, and up to three of those could have become Earth co-orbitals. However, these results are based upon specific initial conditions of the major planets in the Solar System, close to the current epoch. In particular, over megayear time spans, Earth's eccentricity undergoes excursions up to five times its current value, potentially affecting the chaotic orbital evolution of lunar ejecta and their capture into Earth's co-orbital regions. In the present work, we carry out additional numerical simulations to compute the statistics of co-orbital outcomes across different launch epochs, representative of the full range of Earth's eccentricity values. Our main results are as follows: Kamo'oalewa-like co-orbital outcomes of lunar ejecta vary only slightly across the range of Earth's orbital eccentricity, suggesting no privileged ejecta launching epoch for such objects; the probability of co-orbital outcomes decreases rapidly with increasing launch speed, but long-lived HS-QS states are favored at higher launch speeds.

We give rigorous proofs and generalizations of partial fraction expansions for string amplitudes that were recently discovered by Saha and Sinha.

Near-Earth asteroid, Kamo'oalewa (469219), is one of a small number of known quasi-satellites of Earth; it transitions between quasi-satellite and horseshoe orbital states on centennial timescales, maintaining this dynamics over megayears. The similarity of its reflectance spectrum to lunar silicates and its Earth-like orbit both suggest that it originated from the lunar surface. Here we carry out numerical simulations of the dynamical evolution of particles launched from different locations on the lunar surface with a range of ejection velocities in order to assess the hypothesis that Kamo`oalewa originated as a debris-fragment from a meteoroidal impact with the lunar surface. As these ejecta escape the Earth-Moon environment, they face a dynamical barrier for entry into Earth's co-orbital space. However, a small fraction of launch conditions yields outcomes that are compatible with Kamo`oalewa's orbit. The most favored conditions are launch velocities slightly above the escape velocity from the trailing lunar hemisphere.

Zagier provided eleven conjectural rank two examples for Nahm's problem. All of them have been proved in the literature except for the fifth example, and there is no $q$-series proof for the tenth example. We prove that the fifth and the tenth examples are in fact equivalent. Then we give a $q$-series proof for the fifth example, which confirms a recent conjecture of Wang. This also serves as the first $q$-series proof for the tenth example, whose explicit form was conjectured by Vlasenko and Zwegers in 2011., Comment: 9 pages

We give new proofs of the twelve Rogers-Ramanujan-type identities due to Rogers and Slater that are traditionally associated with the moduli 7, 14 and 28., Comment: 11 pages; minor changes from previous version

Many potential mechanisms for delivering planetary debris to within a few Roche radii of white dwarfs rely on gravitational scattering events that feature perturbers which are giant planets or terrestrial planets. However, the population of these planets orbiting white dwarfs is still unknown, and for a substantial fraction of white dwarfs the largest planetary survivors of stellar evolution may be sub-terrestrial mass minor planets. Here, we attempt to identify the smallest mass perturbers that could pollute white dwarfs. Through computationally expensive numerical simulations of both unstable and stable configurations of minor planets, we find that this critical lower bound equals approximately one Luna mass. Further, we find that as this mass limit is approached from above, the typical cooling age at which white dwarf pollution occurs increases. Consequently, there is a two order-of-magnitude range of perturber masses between Earth and its moon that has remained largely unexplored in white dwarf pollution studies, despite the potential formation of thousands of such Luna-sized objects in these systems., Comment: Accepted for publication in MNRAS

We present a new proof of the Rogers-Ramanujan identities. Surprisingly, all its ingredients are available already in Rogers seminal paper from 1894, where he gave a considerably more complicated proof.

We study nearest-neighbour correlation functions for the ground state of the supersymmetric XYZ spin chain with odd length and periodic boundary conditions. Under a technical assumption related to the $Q$-operator of the corresponding eight-vertex model, we show that they can be expressed exactly in terms of the Painlev\'e VI tau functions $s_n$ and $\bar s_n$ introduced by Bazhanov and Mangazeev. Furthermore, we give an interpretation of the correlation functions in terms of the Painlev\'e VI Hamiltonian., Comment: 38 pages

We present three studies examining death in children's animated films. Study 1 is a content analysis of 49 films. We found that death is often portrayed in films, but many deaths occurred off-screen. Deaths were mostly portrayed in a biologically accurate manner, but some films portrayed biological misconceptions. Study 2 (n = 433) reports on parents' attitudes and parent-child conversations about death in films. Children's questions about death in animated films were similar to their questions about death more generally. Animated films may provide a context for parent-child conversations about death, as parents often watched these films with their children. However, it appeared that few parents took advantage of this opportunity to talk about death with their children. [This paper was published in "Cognitive Development" v59 Article 101063 2021.]

We present a detailed analysis of certain quantum spin systems with inhomogeneous (non-random) mean-field interactions. Examples include, but are not limited to, the interchange- and spin singlet projection interactions on complete bipartite graphs. Using two instances of the representation theoretic framework of Schur--Weyl duality, we can explicitly compute the free energy and other thermodynamic limits in the models we consider. This allows us to describe the phase-transition, the ground-state phase diagram, and the expected structure of extremal states., Comment: 55 pages, 9 figures

Children's understanding of death has been a topic of interest to researchers investigating the development of children's thinking and clinicians focusing on children's coping with the death of a loved one. Traditionally, researchers in cognitive development have mainly focused on death from a biological perspective. Current research suggests that exploring religious and spiritual conceptualizations might enrich our understanding of how children come to think about death. In particular, we review different methodological approaches that suggest that children form their understanding of death by engaging in conversations and question asking with family members, consuming cultural products, and actively participating in cultural rituals. We also provide some examples on how children combine different belief systems to form their understanding of death. Finally, we discuss recent research on how socialization with regards to death might be related to coping and bereavement after the death of a loved one. [This paper was published in "Child Development Perspectives" v14 n1 p55-60 2020.]

Using a mixed-methods approach, we examined how participants' memories of socialization regarding death might influence their self-reported coping with losses in childhood and adulthood. We recruited 318 adults to complete an online survey. Path analyses indicated that participants who remembered their parents shielding them less from issues related to death reported better coping as children and adults. Qualitative responses suggested participants wanted to receive more information about death from their parents as they went through the grieving process. We highlight the potential benefits of socializing children about death, and how it may aid in their coping with death-related events. [This paper was published in "Death Studies" v44 n2 p78-88 2020.]

A mixed-method approach was used to explore parent and child perspectives on death in Mexico. Parents' and children's death-related experiences and understanding of death were examined. While all children in this sample displayed a biological understanding of death, older children were less likely to endorse that all living things die. Children also displayed coexistence of beliefs related to death that can be attributed to both their biological and spiritual understanding of death. We also found that older children were more likely to report that a child should feel sad following the death of a loved one. These findings highlight how cultural practices shape the development of cognitive and affective processes related to death. [This paper was published in "Child Development" v91 n2 e491-e511 2020.]

People believe that treatments for illnesses are effective when they target the cause of the illness. Prior work suggests that biological essentialist explanations of mental illness lead people to prefer medications or other pharmacological treatments. However, prior work has not distinguished between biological and essentialist explanations. In three studies (total n = 517), we presented adults with vignettes about an individual with an artificial mental illness and manipulated the descriptions to emphasize or de-emphasize essentialist characteristics. Critically, none of the vignettes made reference to a biological basis for the disorder. Participants rated their willingness to interact with the person described in the vignettes and how effective they believed drug treatment and talk therapy would be on the mental illness. Across the three studies, describing mental illness with an essentialist framing led participants to think drug treatments would be more effective, but there was no effect for stigma or perceived effectiveness of talk therapy. This effect appears to be mediated by how much participants essentialized individuals with the disorder. The first framing that participants encountered seemed to shape their reasoning for the remainder of the study, even if they saw conflicting framing later on. The framing manipulation had similar effects for individuals with and without a mental illness. Results suggest that it is important to consider how mental illness is framed to the general public as it might impact people's treatment preferences. [This paper was published in "Memory & Cognition" v48 n8 p1317-1333 2020.]

People often have difficulty understanding processes of biological change, and they typically reject drastic life cycle changes such as metamorphosis, except for animals with which they are familiar. Even after a lesson about metamorphosis, people often do not generalize to animals not seen during the lesson. This might be partially due to the perceptual richness of the diagrams typically used during lessons on metamorphosis, which serves to emphasize the individual animal rather than a class of animals. In two studies, we examined whether the perceptual richness of a diagram influences adults' learning and transfer of knowledge about metamorphosis. One study was conducted in a laboratory setting and the other online. In both studies, adults who saw the bland diagram during the lesson accurately transferred more than adults who saw the rich diagram during the lesson. [This paper was published in "Applied Cognitive Psychology" v34 n5 p1101-1117 2020 (EJ1265961).]

Parents are typically in charge of purchasing the food that their children eat, but little is known about how parents decide if particular foods are healthy for their children and how their beliefs about nutrition influence their children's beliefs. In two studies, we investigated how parents of children ages 4 to 12 (N = 826) make decisions about the healthiness of foods, when presented with different representations of the same nutritional information. Providing parents with nutritional information did not influence their ratings of how healthy food items are, compared to when they are shown only pictures of the foods. Parents reported talking with their children about nutrition, believed they are the best source of information for children about nutrition, and believed their nutrition beliefs influence their child's beliefs. Our findings highlight the role of prior knowledge in food cognition and how beliefs about foods are transmitted from parents to children. [This paper was published in "Cognitive Development" v55 Article 100913 2020.]

We present four infinite families of mutually commuting difference operators which include the deformed elliptic Ruijsenaars operators. The trigonometric limit of this kind of operators was previously introduced by Feigin and Silantyev. They provide a quantum mechanical description of two kinds of relativistic quantum mechanical particles which can be identified with particles and anti-particles in an underlying quantum field theory. We give direct proofs of the commutativity of our operators and of some other fundamental properties such as kernel function identities. In particular, we give a rigorous proof of the quantum integrability of the deformed Ruijsenaars model., Comment: 29 pages

Kajihara obtained in 2004 a remarkable transformation formula connecting multiple basic hypergeometric series associated with $A$-type root systems of different ranks. By multiple principle specialisations of his formula, we deduce kernel identities for deformed Macdonald-Ruijsenaars (MR) and Noumi-Sano (NS) operators. The deformed MR operators were introduced by Sergeev and Veselov in the first order case and by Feigin and Silantyev in the higher order cases. As applications of our kernel identities, we prove that all of these operators pairwise commute and are simultaneously diagonalised by the super-Macdonald polynomials. We also provide an explicit description of the algebra generated by the deformed MR and/or NS operators by a Harish-Chandra type isomorphism and show that the deformed MR (NS) operators can be viewed as restrictions of inverse limits of ordinary MR (NS) operators., Comment: 29 pages

With the upcoming launch of large constellations of satellites in the low-Earth orbit (LEO) region it will become important to organize the physical space occupied by the different operating satellites in order to minimize critical conjunctions and avoid collisions. Here, we introduce the definition of space occupancy as the domain occupied by an individual satellite as it moves along its nominal orbit under the effects of environmental perturbations throughout a given interval of time. After showing that space occupancy for the zonal problem is intimately linked to the concept of frozen orbits and proper eccentricity, we provide frozen-orbit initial conditions in osculating element space and obtain the frozen-orbit polar equation to describe the space occupancy region in closed analytical form. We then analyze the problem of minimizing space occupancy in a realistic model including tesseral harmonics, third-body perturbations, solar radiation pressure, and drag. The corresponding initial conditions, leading to what we call minimum space occupancy (MiSO) orbits, are obtained numerically for a set of representative configurations in LEO. The implications for the use of MiSO orbits to optimize the design of mega-constellations are discussed., Comment: Submitted to The Journal of Guidance, Control, and Dynamics

Many random growth models have the property that the set of discovered sites, scaled properly, converges to some deterministic set as time grows. Such results are known as shape theorems. Typically, not much is known about the shapes. For first passage percolation on $\mathbb{Z}^d$ we only know that the shape is convex, compact, and inherits all the symmetries of $\mathbb{Z}^d$. Using simulated data we construct and fit a neural network able to adequately predict the shape of the set of discovered sites from the mean, standard deviation, and percentiles of the distribution of the passage times. The purpose of the note is two-fold. The main purpose is to give researchers a new tool for \textit{quickly} getting an impression of the shape from the distribution of the passage times -- instead of having to wait some time for the simulations to run, as is the only available way today. The second purpose of the note is simply to introduce modern machine learning methods into this area of discrete probability, and a hope that it stimulates further research.

Luna 3 (or Lunik 3 in Russian sources) was the first spacecraft to perform a flyby of the Moon. Launched in October 1959 on a translunar trajectory with large semi-major axis and eccentricity, it collided with the Earth in late March 1960. The short, 6-month dynamical lifetime has often been explained through an increase in eccentricity due to the Lidov-Kozai effect. However, the classical Lidov-Kozai solution is only valid in the limit of small semi-major axis ratio, a condition that is satisfied only for solar (but not for lunar) perturbations. We undertook a study of the dynamics of Luna 3 with the aim of assessing the principal mechanisms affecting its evolution. We analyze the Luna 3 trajectory by generating accurate osculating solutions, and by comparing them to integrations of singly- and doubly-averaged equations of motion in vectorial form. Lunar close encounters, which cannot be reproduced in an averaging approach, decisively affect the trajectory and break the doubly-averaged dynamics. Solar perturbations induce oscillations of intermediate period that affect the geometry of the close encounters and cause the singly-averaged and osculating inclinations to change quadrants (the orbital plane ``flips''). We find that the peculiar evolution of Luna 3 can only be explained by taking into account lunar close encounters and intermediate-period terms; such terms are averaged out in the Lidov-Kozai solution, which is not adequate to describe translunar or cislunar trajectories. Understanding the limits of the Lidov-Kozai solution is of particular significance for the motion of objects in the Earth-Moon environment and of exoplanetary systems., Comment: 17 pages, 7 figures. Accepted for publication in Celestial Mechanics and Dynamical Astronomy, topical collection "Toward the Moon and Beyond." Parts of this work were presented at the 2018 John L. Junkins Dynamical Systems Symposium and at the 2019 Meeting of the AAS Division on Dynamical Astronomy (DDA)

Multidimensional matrix inversions provide a powerful tool for studying multiple hypergeometric series. In order to extend this technique to elliptic hypergeometric series, we present three new multidimensional matrix inversions. As applications, we obtain a new $A_r$ elliptic Jackson summation, as well as several quadratic, cubic and quartic summation formulas.

We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first establish a baseline for evaluating various techniques for close-encounter prediction and collision-probability calculation (Hoots et al. 1984, Gronchi 2005, JeongAhn and Malhotra 2015) by carrying out brute-force numerical simulations and using a sequence of filters to greatly reduce the computational expense of the algorithm. Next, we estimate conjunction events in the orbital environment following the anticipated deployments of the OneWeb LEO and SpaceX Starlink mega-constellations. As a final step, we investigate Minimum Space Occupancy (MiSO) orbits (Bombardelli et al. 2018), a generalization of the well-known frozen orbits that account for the perturbed-Keplerian dynamics of the Earth-Moon-Sun-satellite system. We evaluate the ability of MiSO configurations of the proposed mega-constellations, as suggested by Bombardelli et al. 2018, to reduce the risk of endogenous (intra-constellation) collisions. The results indicate that the adoption of the MiSO orbital configuration can significantly reduce risk with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites., Comment: 32 pages, 19 figures; submitted for publication in Adv. Space Res.; comments welcome

The description of the long-term dynamics of highly elliptic orbits under third-body perturbations may require an expansion of the disturbing function in series of the semi-major axes ratio up to higher orders. To avoid dealing with long series in trigonometric functions, we refer the motion to the apsidal frame and efficiently remove the short-period effects of this expansion in vectorial form up to an arbitrary order. We then provide the variation equations of the two fundamental vectors of the Keplerian motion by analogous vectorial recurrences, which are free from singularities and take a compact form useful for the numerical propagation of the flow in mean elements.

Visualizations are commonly used in educational materials, however not all visualizations are equally effective at promoting learning. Prior research has supported the idea that both perceptually rich and bland visualizations are beneficial for learning and generalization. We investigated whether the perceptual richness of a life cycle diagram influenced children's learning of metamorphosis, a concept that prior work suggests is difficult for people to generalize. Using identical materials, Study 1 (n = 76) examined learning and generalization of metamorphosis in first and second grade students and Study 2 (n = 53) did so in fourth and fifth grade students. Bayesian regressions revealed that first and second grade students learned more from the lesson with the perceptually rich diagram, while fourth and fifth grade students generalized more broadly with the bland diagram. The results from the fourth and fifth grade students are similar to prior research with adults. This suggests that the effect of perceptual richness on learning and generalization changes over development. [This paper was published in "Journal of Experimental Child Psychology" v213 Article 10526 2022.]

We investigated children's (n = 120; 3- to 11-year-olds) and adults' (n = 18) reasoning about life-cycle changes in biological organisms by examining their endorsements of four different patterns of life- span changes. Participants were presented with two separate tasks: (a) judging possible adult versions of a juvenile animal and (b) judging possible juvenile versions of an adult animal. The stimuli enabled us to examine the endorsement of four different patterns of change: identical growth, natural growth, dramatic change, and speciation. The results suggest that endorsement of the different patterns is influenced by age and familiarity. Young children and individuals confronted with unfamiliar organisms often endorsed an identical growth that emphasizes the stability of features over the life span and between parents and offspring. The results are interpreted as supporting the idea that cognitive constraints influence individuals' reasoning about biological change and that the influence of these constraints is most notable when individuals are young or are presented with unfamiliar biological organisms. [This paper was published in "Journal of Experimental Child Psychology" v173 p205-221 2018.]

In this note, we consider the frog model on $\mathbb{Z}^d$ and a two-type version of it with two types of particles. For the one-type model, we show that the asymptotic shape does not depend on the initially activated set and the configuration there. For the two-type model, we show that the possibility for the types to coexist in that both of them activate infinitely many particles does not depend on the choice of the initially activated sets and the configurations there.

Kanade and Russell conjectured several Rogers-Ramanujan-type partition identities, some of which are related to level $2$ characters of the affine Lie algebra $A_9^{(2)}$. Many of these conjectures have been proved by Bringmann, Jennings-Shaffer and Mahlburg. We give new proofs of five conjectures first proved by those authors, as well as four others that have been open until now. Our proofs for the new cases use quadratic transformations for Askey-Wilson and Rogers polynomials. We also obtain some related results, including a new proof of a partition identity conjectured by Capparelli and first proved by Andrews., Comment: 21 pages

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term (120 years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes (400 km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded "ReDSHIFT" project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an "atlas" of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth's oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device., Comment: 25 pages, 10 figures

The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable DV cost. We simulated the dynamics over 120-200 years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90deg, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed ~40years, mainly around three specific inclination values: 46deg, 56deg, and 68deg; a result compatible with our understanding of MEO secular dynamics. For DV <= 300 m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with e<0.3), reentry times from GNSS altitudes exceed ~70 years, while low-cost (DV ~= 5-35 m/s) graveyard orbits, stable for at lest 200 years, are found for eccentricities up to e~0.018. This investigation was carried out in the framework of the EC-funded "ReDSHIFT" project., Comment: 39 pages, 23 figures

In this paper, we study the long-term dynamical evolution of highly-elliptical orbits (HEOs) in the medium-Earth orbit (MEO) region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realistic dynamical model, aimed primarily at recording the dynamical lifetime of such orbits (defined as the time needed for atmospheric reentry) and understanding its dependence on initial conditions and other parameters, such as the area-to-mass ratio (A/m). Our results are presented in the form of 2-D lifetime maps, for different values of inclination, A/m, and drag coefficient. We find that the majority of small debris (> 70%, depending on the inclination) can naturally reenter within 25-90 years, but these numbers are significantly less optimistic for large debris (e.g., upper stages), with the notable exception of those launched from high latitude (Baikonur). We estimate the reentry probability and mean dynamical lifetime for different classes of GTOs and we find that both quantities depend primarily and strongly on initial perigee altitude. Atmospheric drag and higher A/m values extend the reentry zones, especially at low inclinations. For high inclinations, this dependence is weakened, as the primary mechanisms leading to reentry are overlapping lunisolar resonances. This study forms part of the EC-funded (H2020) "ReDSHIFT" project., Comment: 16 pages, 10 figures

This paper is concerned with the comparison of semi-analytical and non-averaged propagation methods for Earth satellite orbits. We analyse the total integration error for semi-analytical methods and propose a novel decomposition into dynamical, model truncation, short-periodic, and numerical error components. The first three are attributable to distinct approximations required by the method of averaging, which fundamentally limit the attainable accuracy. In contrast, numerical error, the only component present in non-averaged methods, can be significantly mitigated by employing adaptive numerical algorithms and regularized formulations of the equations of motion. We present a collection of non-averaged methods based on the integration of existing regularized formulations of the equations of motion through an adaptive solver. We implemented the collection in the orbit propagation code THALASSA, which we make publicly available, and we compared the non-averaged methods to the semi-analytical method implemented in the orbit propagation tool STELA through numerical tests involving long-term propagations (on the order of decades) of LEO, GTO, and high-altitude HEO orbits. For the test cases considered, regularized non-averaged methods were found to be up to two times slower than semi-analytical for the LEO orbit, to have comparable speed for the GTO, and to be ten times as fast for the HEO (for the same accuracy). We show for the first time that efficient implementations of non-averaged regularized formulations of the equations of motion, and especially of non-singular element methods, are attractive candidates for the long-term study of high-altitude and highly elliptical Earth satellite orbits., Comment: 33 pages, 10 figures, 7 tables. Part of the CMDA Topical Collection on "50 years of Celestial Mechanics and Dynamical Astronomy". Comments and feedback are encouraged

The classical Selberg integral contains a power of the Vandermonde determinant. When that power is a square, it is easy to prove Selberg's identity by interpreting it as a determinant of one-variable integrals. We give similar proofs of summation and transformation formulas for continuous and discrete elliptic Selberg integrals. In the continuous case, the same proof was previously given by Noumi. Special cases of these identities have found applications in combinatorics., Comment: 10 pages; minor changes from previous version

The slow deformation of terrestrial orbits in the medium range, subject to lunisolar resonances, is well approximated by a family of Hamiltonian flow with $2.5$ degree-of-freedom. The action variables of the system may experience chaotic variations and large drift that we may quantify. Using variational chaos indicators, we compute high-resolution portraits of the action space. Such refined meshes allow to reveal the existence of tori and structures filling chaotic regions. Our elaborate computations allow us to isolate precise initial conditions near specific zones of interest and study their asymptotic behaviour in time. Borrowing classical techniques of phase- space visualisation, we highlight how the drift is mediated by the complement of the numerically detected KAM tori., Comment: 22 pages, 11 figures, 1 table, 52 references. Comments and feedbacks greatly appreciated. This article is part of the Research Topic `The Earth-Moon System as a Dynamical Laboratory', confer https://www.frontiersin.org/research-topics/5819/the-earth-moon-system-as-a-dynamical-laboratory

The size of the giant component in the configuration model, measured by the asymptotic fraction of vertices in the component, is given by a well-known expression involving the generating function of the degree distribution. In this note, we argue that the distribution over small degrees is more important for the size of the giant component than the precise distribution over very large degrees. In particular, the tail behavior of the degree distribution does not play the same crucial role for the size of the giant as it does for many other properties of the graph. Upper and lower bounds for the component size are derived for an arbitrary given distribution over small degrees $d\leq L$ and given expected degree, and numerical implementations show that these bounds are close already for small values of $L$. On the other hand, examples illustrate that, for a fixed degree tail, the component size can vary substantially depending on the distribution over small degrees.

We give a survey of elliptic hypergeometric functions associated with root systems, comprised of three main parts. The first two form in essence an annotated table of the main evaluation and transformation formulas for elliptic hypergeometric integeral and series on root systems. The third and final part gives an introduction to Rains' elliptic Macdonald-Koornwinder theory (in part also developed by Coskun and Gustafson)., Comment: 28 pages. Modulo minor edits this paper will appear as a chapter in the book "Multivariable Special Functions" (edited by Tom Koornwinder and Jasper Stokman) which is part of the Askey-Bateman project. This version: essential correction of equation (1.2.14) and some further minor corrections

A multi-type preferential attachment tree is introduced, and studied using general multi-type branching processes. For the $p$-type case we derive a framework for studying the tree where a type $i$ vertex generates new type $j$ vertices with rate $w_{ij}(n_1,n_2,\ldots, n_p)$ where $n_k$ is the number of type $k$ vertices previously generated by the type $i$ vertex, and $w_{ij}$ is a non-negative function from $\mathbb{N}^p$ to $\mathbb{R}$. The framework is then used to derive results for trees with more specific attachment rates. In the case with linear preferential attachment---where type $i$ vertices generate new type $j$ vertices with rate $w_{ij}(n_1,n_2,\ldots, n_p)=\gamma_{ij}(n_1+n_2+\dots +n_p)+\beta_{ij}$, where $\gamma_{ij}$ and $\beta_{ij}$ are positive constants---we show that under mild regularity conditions on the parameters $\{\gamma_{ij}\}, \{\beta_{ij}\}$ the asymptotic degree distribution of a vertex is a power law distribution. The asymptotic composition of the vertex population is also studied.

We prove a multivariable elliptic extension of Jackson's summation formula conjectured by Spiridonov. The trigonometric limit case of this result is due to Gustafson and Rakha. As applications, we obtain two further multivariable elliptic Jackson summations and two multivariable elliptic Bailey transformations. The latter four results are all new even in the trigonometric case.

We study biorthogonal functions related to basic hypergeometric integrals with coupled continuous and discrete components. Such integrals appear as superconformal indices for three-dimensional quantum field theories and also in the context of solvable lattice models. We obtain explicit biorthogonal systems given by products of two of Rahman's biorthogonal rational ${}_{10}W_9$-functions or their degenerate cases. We also give new bilateral extensions of the Jackson and $q$-Saalsch\"utz summation formulas and new continuous and discrete biorthogonality measures for Rahman's functions., Comment: 22 pages

In these lecture notes I give an elementary introduction to elliptic hypergeometric functions. I focus on motivating the main ideas and constructions, rather than giving a comprehensive survey. The lectures include a brief explanation of the historical origin of elliptic hypergeometric functions in the context of solvable lattice models. In particular, I give a new proof of the fact that fused Boltzmann weights for the elliptic solid-on-solid model can be expressed as elliptic hypergeometric sums., Comment: Lecture notes for a course given at the summer school OPSF-S6, College Park, Maryland, 11-15 July 2016. Slight improvements compared to previous version

We introduce several new identities combining basic hypergeometric sums and integrals. Such identities appear in the context of superconformal index computations for three-dimensional supersymmetric dual theories. We give both analytic proofs and physical interpretations of the presented identities., Comment: 25 pages, v2: minor corrections and comments

We consider Earth satellite orbits in the range of semi-major axes where the perturbing effects of Earth's oblateness and lunisolar gravity are of comparable order. This range covers the medium-Earth orbits (MEO) of the Global Navigation Satellite Systems and the geosynchronous orbits (GEO) of the communication satellites. We recall a secular and quadrupolar model, based on the Milankovitch vector formulation of perturbation theory, which governs the long-term orbital evolution subject to the predominant gravitational interactions. We study the global dynamics of this two-and-a-half degrees-of-freedom Hamiltonian system by means of the fast Lyapunov indicator (FLI), used in a statistical sense. Specifically, we characterize the degree of chaoticity of the action space using angle-averaged normalized FLI maps, thereby overcoming the angle dependencies of the conventional stability maps. Emphasis is placed upon the phase-space structures near secular resonances, which are of first importance to the space debris community. We confirm and quantify the transition from order to chaos in MEO, stemming from the critical inclinations, and find that highly inclined GEO orbits are particularly unstable. Despite their reputed normality, Earth satellite orbits can possess an extraordinarily rich spectrum of dynamical behaviors, and, from a mathematical perspective, have all the complications that make them very interesting candidates for testing the modern tools of chaos theory., Comment: 30 pages, 9 figures. Accepted for publication in the Astronomical Journal

The navigation satellite constellations in medium-Earth orbit exist in a background of third-body secular resonances stemming from the perturbing gravitational effects of the Moon and the Sun. The resulting chaotic motions, emanating from the overlapping of neighboring resonant harmonics, induce especially strong perturbations on the orbital eccentricity, which can be transported to large values, thereby increasing the collision risk to the constellations and possibly leading to a proliferation of space debris. We show here that this transport is of a diffusive nature and we present representative diffusion maps that are useful in obtaining a global comprehension of the dynamical structure of the navigation satellite orbits., Comment: 8 pages, 5 figures, conference Chaos, complexity and transport (Marseille, France)

We introduce and study twelve multivariable theta functions defined by pfaffians with elliptic function entries. We show that, when the crossing parameter is a cubic root of unity, the domain wall partition function for the eight-vertex-solid-on-solid model can be written as a sum of two of these pfaffians. As a limit case, we express the domain wall partition function for the three-colour model as a sum of two Hankel determinants. We also show that certain solutions of the TQ-equation for the supersymmetric eight-vertex model can be expressed in terms of elliptic pfaffians., Comment: 34 pages

In this article we introduce a dynamic Erd\H{o}s-R\'enyi graph model, in which, independently for each vertex pair, edges appear and disappear according to a Markov on-off process. In studying the dynamic graph we present two main results. The first being on how long it takes for the graph to reach stationarity. We give an explicit expression for this time, as well as proving that this is the fastest time to reach stationarity among all strong stationary times. The second result concerns the time it takes for the dynamic graph to reach a certain number of edges. We give an explicit expression for the expected value of such a time, as well as study its asymptotic behavior. This time is related to the first time the dynamic Erd\H{o}s-R\'enyi graph contains a cluster exceeding a certain size.

Recent studies have shown that the medium-Earth orbit (MEO) region of the Global Navigation Satellite Systems is permeated by a devious network of lunisolar secular resonances, which can interact to produce chaotic and diffusive motions. The precarious state of the four navigation constellations, perched on the threshold of instability, makes it understandable why all past efforts to define stable graveyard orbits, especially in the case of Galileo, were bound to fail; the region is far too complex to allow of an adoption of the simple geosynchronous disposal strategy. We retrace one such recent attempt, funded by ESA's General Studies Programme in the frame of the GreenOPS initiative, that uses a systematic parametric approach and the straightforward maximum-eccentricity method to identify long-term stable regions, suitable for graveyards, as well as large-scale excursions in eccentricity, which can be used for post-mission deorbiting of constellation satellites. We then apply our new results on the stunningly rich dynamical structure of the MEO region toward the analysis of these disposal strategies for Galileo, and discuss the practical implications of resonances and chaos in this regime. We outline how the identification of the hyperbolic and elliptic fixed points of the resonances near Galileo can lead to explicit criteria for defining optimal disposal strategies., Comment: 19 pages, 12 figures. Submitted to MNRAS. Comments are greatly appreciated

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula's expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane's expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long--term motion of the highly elliptical and critically--inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics., Comment: 30 pages, 4 figures