Search NASA⌕ Search

SEARCH · Search NASA

Results for “Periodic orbits”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Finite-size effects in periodic coupled cluster calculations

Here, we provide the first rigorous study of the finite-size error in the simplest and representative coupled cluster theory, namely the coupled cluster doubles (CCD) theory, for gapped periodic systems. Given exact Hartree-Fock orbitals and their corresponding orbital energies, we demonstrate that the correlation energy obtained from the approximate CCD method, after a finite number of fixed-point iterations over the amplitude equation, exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Here $N_k$ is the number of discretization points in the Brillouin zone and characterizes the system size. Under additional assumptions ensuring the convergence of the fixed-point iterations, we demonstrate that the CCD correlation energy also exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Our analysis shows that the dominant error lies in the coupled cluster amplitude calculation, and the convergence of the finite-size error in energy calculation can be boosted to $\mathcal{O}(N^{-1}_k)$ with accurate amplitudes. This also provides the first proof of the scaling of the finite-size error in the third order Møller-Plesset perturbation theory (MP3) for periodic systems.

97 MATHEMATICS AND COMPUTING↗

TIC 168789840: A Sextuply Eclipsing Sextuple Star System

We report the discovery of a sextuply eclipsing sextuple star system from TESS data, TIC 168789840, also known as TYC 7037-89-1, the first known sextuple system consisting of three eclipsing binaries. The target was observed in Sectors 4 and 5 during Cycle 1, with lightcurves extracted from TESS Full Frame Image data. It was also previously observed by the WASP survey and ASAS-SN. The system consists of three gravitationally bound eclipsing binaries in a hierarchical structure of an inner quadruple system with an outer binary subsystem. Follow-up observations from several different observatories were conducted as a means of determining additional parameters. The system was resolved by speckle interferometry with a 0.″42 separation between the inner quadruple and outer binary, inferring an estimated outer period of ∼2 kyr. It was determined that the fainter of the two resolved components is an 8.217 day eclipsing binary, which orbits the inner quadruple that contains two eclipsing binaries with periods of 1.570 days and 1.306 days. Markov Chain Monte Carlo (MCMC) analysis of the stellar parameters has shown that the three binaries of TIC 168789840 are “triplets,” as each binary is quite similar to the others in terms of mass, radius, and T {sub eff}. As a consequence of its rare composition, structure, and orientation, this object can provide important new insight into the formation, dynamics, and evolution of multiple star systems. Future observations could reveal if the intermediate and outer orbital planes are all aligned with the planes of the three inner eclipsing binaries.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

TOI-503: The First Known Brown-dwarf Am-star Binary from the TESS Mission

We report the discovery of an intermediate-mass transiting brown dwarf (BD), TOI-503b, from the TESS mission. TOI-503b is the first BD discovered by TESS, and it has circular orbit around a metallic-line A-type star with a period of P = 3.6772 ± 0.0001 days. The light curve from TESS indicates that TOI-503b transits its host star in a grazing manner, which limits the precision with which we measure the BD’s radius (R{sub b}=1.34{sub −0.15}{sup +0.26}R{sub J}). We obtained high-resolution spectroscopic observations with the FIES, Ondřejov, PARAS, Tautenburg, and TRES spectrographs, and measured the mass of TOI-503b to be M {sub b} = 53.7 ± 1.2 M{sub J}. The host star has a mass of M {sub ⋆} = 1.80 ± 0.06 M {sub ⊙}, a radius of R {sub ⋆} = 1.70 ± 0.05R {sub ⊙}, an effective temperature of T {sub eff} = 7650 ± 160 K, and a relatively high metallicity of 0.61 ± 0.07 dex. We used stellar isochrones to derive the age of the system to be ∼180 Myr, which places its age between that of RIK 72b (a ∼10 Myr old BD in the Upper Scorpius stellar association) and AD 3116b (a ∼600 Myr old BD in the Praesepe cluster). Given the difficulty in measuring the tidal interactions between BDs and their host stars, we cannot precisely say whether this BD formed in situ or has had its orbit circularized by its host star over the relatively short age of the system. Instead, we offer an examination of plausible values for the tidal quality factor for the star and BD. TOI-503b joins a growing number of known short-period, intermediate-mass BDs orbiting main-sequence stars, and is the second such BD known to transit an A star, after HATS-70b. With the growth in the population in this regime, the driest region in the BD desert (35--55M{sub J}sin i) is reforesting.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectroscopic Orbits of Subsystems in Multiple Stars. VII

Spectroscopic orbits of main-sequence stars HIP 3150A, 6873B, 11537A, 22531A, 22534B, 31089B, 49336A, 104833C, and 107731A belonging to eight multiple systems are determined from high-resolution spectra taken with CHIRON. Two of those are twins with mass ratios above 0.95. HIP 11537 is a young three-tier quadruple system with inner periods of 22.3 and 1146 days and the outer period of 3 kyr. HIP 22531 (ι Pic) is the brightest star in a hierarchical system with six components. It is a spectroscopic binary with periods of 1.56 days and 2.75 yr, as well as a γ Dor variable with a period of 0.67 day, possibly in a 7:3 resonance with the inner orbit. HIP 22534, also member of this system, is a double-lined binary with a period of 208 days. For HIP 31089, both the spectro-interferometric 32 yr outer orbit and the 213 day orbit of the subsystem are determined. HIP 107731 is a triple system with an inner period of 470 days and a fast spatial motion, likely metal-poor. New orbits contribute to the statistics of hierarchical multiplicity in the solar neighborhood.

79 ASTRONOMY AND ASTROPHYSICS↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗

Architecture of TOI-561 planetary system

ABSTRACT We present new observations from CHEOPS (CHaracterising ExOPlanet Satellite) and TESS (Transiting Exoplanet Survey Satellite) to clarify the architecture of the planetary system hosted by the old Galactic thick disc star TOI-561. Our global analysis, which also includes previously published photometric and radial velocity data, incontrovertibly proves that TOI-561 is hosting at least four transiting planets with periods of 0.44 d (TOI-561 b), 10.8 d (TOI-561 c), 25.7 d (TOI-561 d), and 77.1 d (TOI-561 e) and a fifth non-transiting candidate, TOI-561f with a period of 433 d. The precise characterization of TOI-561’s orbital architecture is interesting since old and metal-poor thick disc stars are less likely to host ultrashort-period super-Earths like TOI-561 b. The new period of planet -e is consistent with the value obtained using radial velocity alone and is now known to be $77.14399\pm 0.00025$ d, thanks to the new CHEOPS and TESS transits. The new data allowed us to improve its radius ($R_p = 2.517 \pm 0.045\,\mathrm{ R}_{\rm{\oplus }}$ from 5 per cent to 2 per cent precision) and mass ($M_p = 12.4 \pm 1.4\, \mathrm{ M}_{\rm{\oplus }}$) estimates, implying a density of $\rho _p = 0.778 \pm 0.097\, \rho _{\rm{\oplus }}$. Thanks to recent TESS observations and the focused CHEOPS visit of the transit of TOI-561 e, a good candidate for exomoon searches, the planet’s period is finally constrained, allowing us to predict transit times through 2030 with 20-min accuracy. We present an updated version of the internal structure of the four transiting planets. We finally performed a detailed stability analysis, which confirmed the long-term stability of the outer planet TOI-561 f.

Piotto, G. (ORCID:0000000299376387)↗

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Frustrated electron hopping from the orbital configuration in a two-dimensional lattice

Electron hopping on spatially periodic lattices gives rise to intriguing electronic behaviour. For example, hopping on the geometrically frustrated two-dimensional kagome, dice and Lieb lattices yields electronic band structures with both massless Dirac-like and perfectly dispersion-less, flat bands. As materials featuring the dice and Lieb lattice structures are scarce, an alternative approach proposes to leverage atomic orbitals to realize the characteristic electron hopping of geometrically frustrated lattices. This strategy promises to expand the list of candidate materials with frustrated electron hopping, but is yet to be shown in experiments. Here, in this study, we demonstrate frustrated hopping in the van der Waals intermetallic Pd 5 AlI 2 , emerging from the arrangement of atomic orbitals in a primitive square lattice. Using angle-resolved photoemission spectroscopy and quantum oscillation measurements, we reveal that the band structure of Pd 5 AlI 2 includes linear Dirac-like bands intersected at their crossing point by a locally flat band—an essential characteristic of frustrated hopping in Lieb and dice lattices. Moreover, this compound shows exceptional chemical stability, with its unusual bulk band structure and metallicity persisting in ambient conditions down to the monolayer limit. Hence, our results showcase a way to realize electronic structures characteristic of geometrically frustrated lattices in non-frustrated systems.

36 MATERIALS SCIENCE↗

TOI-2374 b and TOI-3071 b: two metal-rich sub-Saturns well within the Neptunian desert

ABSTRACT We report the discovery of two transiting planets detected by the Transiting Exoplanet Survey Satellite (TESS), TOI-2374 b and TOI-3071 b, orbiting a K5V and an F8V star, respectively, with periods of 4.31 and 1.27 d, respectively. We confirm and characterize these two planets with a variety of ground-based and follow-up observations, including photometry, precise radial velocity monitoring and high-resolution imaging. The planetary and orbital parameters were derived from a joint analysis of the radial velocities and photometric data. We found that the two planets have masses of (57 ± 4) M⊕ or (0.18 ± 0.01) MJ, and (68 ± 4) M⊕ or (0.21 ± 0.01) MJ, respectively, and they have radii of (6.8 ± 0.3) R⊕ or (0.61 ± 0.03) RJ and (7.2 ± 0.5) R⊕ or (0.64 ± 0.05) RJ, respectively. These parameters correspond to sub-Saturns within the Neptunian desert, both planets being hot and highly irradiated, with Teq ≈ 745 K and Teq ≈ 1812 K, respectively, assuming a Bond albedo of 0.5. TOI-3071 b has the hottest equilibrium temperature of all known planets with masses between 10 and 300 M⊕ and radii less than 1.5 RJ. By applying gas giant evolution models we found that both planets, especially TOI-3071 b, are very metal-rich. This challenges standard formation models which generally predict lower heavy-element masses for planets with similar characteristics. We studied the evolution of the planets’ atmospheres under photoevaporation and concluded that both are stable against evaporation due to their large masses and likely high metallicities in their gaseous envelopes.

Hacker, Alejandro (ORCID:0009000994520181)↗

Tidal Circularization of Binaries by Resonance Locking. I. The Importance of the Pre-main Sequence

Although tidal dissipation in binary stars has been studied for over a century, theoretical predictions have yet to match the observed properties of binary populations. This work quantitatively examines the recent proposal of tidal circularization by resonance locking, where tidal dissipation arises from resonances between the star’s natural oscillation frequencies and harmonics of the orbital frequency, and where resonances are “locked” for an extended period of time due to concurrent stellar evolution. We focus on tidal resonances with axi-symmetric gravity modes, and examine binaries with primary masses from one to two solar masses. We find that orbital evolution via resonance locking occurs primarily during the star’s pre-main-sequence phase, with the main-sequence phase contributing negligibly. Resonance locking, ignoring nonlinearity, can circularize binaries with peri-center distances out to ∼10 stellar radii, corresponding to circular periods of ∼4–6 days. However, we find resonantly excited gravity modes will become nonlinear in stellar cores, which prevents them from reaching their full, linear amplitudes. We estimate that such a “saturated resonance lock” reduces the circularization period by about a third, but resonance locking remains much more effective than the cumulative actions of equilibrium tides. In a companion paper, we examine recent binary data to compare against theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

WIYN Open Cluster Study. LXXXII. Radial-velocity Measurements and Spectroscopic Binary Orbits in the Open Cluster NGC 7789

We present an extensive time-series radial-velocity (RV) survey of stars in the rich open cluster NGC 7789 (1.6 Gyr, [Fe/H] = +0.02). The stellar sample lies within an 18′ circular radius from the cluster center (10 pc in projection, or about 2 core radii), and includes giants, red clump stars, blue stragglers, red stragglers, sub-subgiants, and main-sequence stars down to 1 mag below the turnoff. Our survey began in 2005 and comprises more than 9000 RV measurements from the Hydra Multi-Object Spectrograph on the WIYN 3.5 m telescope. We identify 624 likely cluster members and present the orbital solutions for 81 cluster binary stars with periods between 1.45 and 4200 days. From the main-sequence binary solutions we fit a circularization period of 7.2{sub −1.1}{sup +0.6} days. We calculate an incompleteness-corrected main-sequence binary frequency of 31% ± 4% for binaries with periods less than 10{sup 4} days, similar to other WIYN Open Cluster Survey (WOCS) open clusters of all ages. We detect a blue straggler binary frequency of 33% ± 17%, consistent with the similarly aged open cluster NGC 6819. We also find one secure, rapidly rotating sub-subgiant and one red straggler candidate in our sample.

79 ASTRONOMY AND ASTROPHYSICS↗

MEGASIM: Lifetimes and Resonances of Earth Trojan Asteroids—The Death of Primordial ETAs?

Abstract We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set. Trojan asteroids co-orbit the Sun with a planet, but remain bound to the Lagrange points, L4 (60° leading the planet) or L5 (60° trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and the central mass. For the inner planets, the range of stability becomes increasingly small, so perturbations from the planets have made primordial Trojans rare. To date, there have been just two ETAs (2010 TK 7 and 2020 XL 5 ), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are shorter than a million years, suggesting they are interlopers rather than members of a stable and long-lasting population. With the largest ETA n -body simulation to date, we are able to track their survival across a wide initialized parameter space. We find that the remaining fraction of ETAs over time is well fit with a stretched exponential function that, when extrapolated beyond our simulation run time, predicts zero ETAs by 2.33 Gyr. We also show correlations between ETA ejections and the periods of the Milankovitch cycles. Though Earth’s orbital dynamics dominate the instabilities of ETAs, we provide evidence that ETA ejections are linked to resonances found in the variation of the orbital elements of many if not all of the planets.

79 ASTRONOMY AND ASTROPHYSICS↗

TESS giants transiting giants V – two hot Jupiters orbiting red giant hosts

ABSTRACT In this work, we present the discovery and confirmation of two hot Jupiters orbiting red giant stars, TOI-4377 b and TOI-4551 b, observed by Transiting Exoplanet Survey Satellite in the Southern ecliptic hemisphere and later followed-up with radial-velocity (RV) observations. For TOI-4377 b, we report a mass of $0.957^{+0.089}_{-0.087} \ M_\mathrm{J}$ and a inflated radius of 1.348 ± 0.081 RJ orbiting an evolved intermediate-mass star (1.36 M⊙ and 3.52 R⊙; TIC 394918211) on a period of of 4.378 d. For TOI-4551 b, we report a mass of 1.49 ± 0.13 MJ and a radius that is not obviously inflated of $1.058^{+0.110}_{-0.062} \ R_\mathrm{J}$, also orbiting an evolved intermediate-mass star (1.31 M⊙ and 3.55 R⊙; TIC 204650483) on a period of 9.956 d. We place both planets in context of known systems with hot Jupiters orbiting evolved hosts, and note that both planets follow the observed trend of the known stellar incident flux-planetary radius relation observed for these short-period giants. Additionally, we produce planetary interior models to estimate the heating efficiency with which stellar incident flux is deposited in the planet’s interior, estimating values of $1.91 \pm 0.48~{{\ \rm per\ cent}}$ and $2.19 \pm 0.45~{{\ \rm per\ cent}}$ for TOI-4377 b and TOI-4551 b, respectively. These values are in line with the known population of hot Jupiters, including hot Jupiters orbiting main-sequence hosts, which suggests that the radii of our planets have re-inflated in step with their parent star’s brightening as they evolved into the post-main sequence. Finally, we evaluate the potential to observe orbital decay in both systems.

Pereira, Filipe (ORCID:0000000221577146)↗

Magnetic and singlet phases in the three-dimensional periodic Anderson Model

Heavy-fermion materials are compounds in which localized f orbitals hybridize with delocalized d ones, leading to quasiparticles with large renormalized masses. The presence of strongly correlated f electrons at the Fermi level may also lead to long-range order, such as magnetism, or unconventional superconductivity. From a theoretical point of view, the “standard model” for heavy-fermion compounds is the periodic Anderson model (PAM). Despite being extensively scrutinized, its thermodynamic properties in three-dimensional (3D) lattices have not been carefully addressed by unbiased methodologies. Here we investigate the 3D PAM employing state-of-the-art finite-temperature auxiliary-field quantum Monte Carlo simulations. Furthermore, we present the behavior of the kinetic energy, the entropy, the specific heat, and the double occupancy as functions of the temperature and the hybridization strength. From these quantities, and by the analysis of the spin-spin correlation functions, we investigate the occurrence of magnetic phase transitions at finite temperatures, and we determine the phase diagram of the model, including the behavior of the Néel temperature as a function of the external parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Five percent measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids’ periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler’s third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be G LMC /G SS = $0.93{^{+0.05}_{-0.04}}$. We discuss the implications of this 5% measurement of Newton’s constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers G LMC < G SS at 2.6σ. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α = 0.90{^{+0.36}_{-0.26}}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

79 ASTRONOMY AND ASTROPHYSICS↗

Proton and electron RLA optics design

We describe optics designs of the key components of proton and electron Recirculating Linear Accelerators (RLAs). They are presented in the context of a high-power hadron accelerator being considered at ORNL and a CEBAF electron energy doubling study, FFA@CEBAF, being developed at Jefferson Lab. Both concepts rely on the Fixed-Field Alternating gradient (FFA) arc optics designs where multiple beam passes are transported by a single beam line. Some of the main challenges of the optics design are in synchronization and orbital and linear optics match of the individual passes with the linac straight. To minimize the complexity of the matching sections and make them more practical, it is preferable to complete as much of the orbital and optical match as possible within a single beam line using the FFA approach. The adiabatic matching technique proposed earlier is one such solution, but it requires a large number of cells. We present a much shorter matching section design based on another approach. A compact dog-leg-type section of several bending dipoles and quadrupoles is used to suppress the orbital offset and dispersion coming out of a periodic cell. Additional quadrupoles are used in the section to keep the beta functions under control and simplify their subsequent match to the linac. Initial beam tracking results are presented.

Benesch, J.↗