SEARCH · Engineering Papers

Results for “PHASE SPACE”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Parton distribution functions from scalar light-front parton gas model

Here, we propose an application of a microcanonical ensemble with light-front kinematics to model the phase-space distribution of relativistic constituents of a bound state. These constituents denoted by partons are treated as classical spin-zero particles confined inside the bound state with inter-parton collisions as their only interaction. The microcanonical molecular dynamics ensemble is applied to obtain the phase-space distribution of such a thermodynamic system. We sample this phase-space distribution using Monte Carlo algorithms to obtain the parton distribution functions (PDFs) in scenarios with 3, 4, and 5 identical partons. In addition PDFs when a selected number of massless partons are mixed with 3 massive partons are also presented.

Microcanonical ensemble↗

DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Neural refinement of sample weights

Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be negative, and several post hoc methods have been developed to eliminate or mitigate the negative values. All of these methods share the common strategy of approximating the average weight as a function of phase space. We introduce an alternative approach, which, instead of reweighting to the average, refines the initial weights with a scaling transformation, utilizing a phase space-dependent factor. Since this new refinement method does not need to model the full weight distribution, it can be more accurate. High-dimensional and unbinned phase space is processed using neural networks for the refinement method. In addition to the refinement method, we introduce a new resampling protocol, which can be used in conjunction with any weight transformation to not only preserve the average weight but also the statistical uncertainties of the initial distribution. Using both realistic and synthetic examples, we show that the new neural refinement method is able to match or exceed the accuracy of similar weight transformations and that the new resampling protocol is simpler in implementation than previous methods while exhibiting equivalent statistical properties.

Artificial neural networks↗

DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗

DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗

DOE OSTI2024

Pressure-Dependent Phase Behavior of Chemically Synthesized Ytterbium Hydride Nanoparticles and Confinement in CMK-3

In this work, we developed a chemical synthesis protocol to prepare both spherical and rod shaped YbH x nanoparticles and infiltrate a mesoporous CMK-3 structure with YbH x nanoparticles. The phase evolution of the YbH x nanoparticles was monitored as a function of the pressure in a diamond anvil cell. The cell volume of the Yb atom showed an unexpected positive deviation from the calculated equation of state (EoS) for the α-phase (space group Pnma) of YbH x . Additionally, the β-phase (space group Fm3̅m) YbH x nanoparticles appeared to become amorphous at approximately 27 GPa, although the same behavior was not observed in the Pnma (α-phase). The Pmna phase remained stable (measured up to 59.3 GPa) and showed a significantly higher than expected volume per Yb atom relative to the calculated EoS. This deviation could be attributed to the nanoscale effect on the thermodynamics of ytterbium hydride formation, which could provide a path forward to accessing higher order ytterbium hydrides (YbH x , x > 3) at pressures much lower than in the bulk metal.

diamond anvil cell↗

DOE OSTI2024

Quantum to classical parton evolution in the QGP

We study the time evolution of the density matrix of a high energy quark in the presence of a dense QCD background that is modeled as a stochastic Gaussian color field. At late times, we find that only the color singlet component of the quark’s reduced density matrix survives the in-medium evolution and that the density matrix becomes asymptotically diagonal in both transverse position and momentum spaces. In addition, we observe an accelerated entropy growth due to the larger phase space being explored by the quark and that the quantum and classical quark entropies converge at late times. We further observe that the quark state loses all memory of the initial condition. Combined with the fact that the reduced density matrix satisfies Boltzmann-diffusion transport, we conclude that the quark reduced density matrix can be interpreted as a classical phase space distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Analytic reconstruction with massive particles: one-loop amplitudes for $0\to \overline{q} qt\overline{t}H$

We present an analytic reconstruction of one-loop amplitudes for the process $0\to \overline{q} qt\overline{t}H$. Our calculation is a novel use of analytic reconstruction, retaining explicit covariance in the massive spin states through the massive spinor-helicity formalism. The analytic reconstruction relies on embedding the massive five-point kinematics in a fully massless eight-point phase space while still building a minimal ansatz directly in the five-point phase space. In order to obtain compact analytic expressions it is necessary to identify suitable partial fraction decompositions and extract common numerator factors, which we achieve through careful inspection of limits in which pairs of denominators vanish. We find that the resulting amplitudes are more numerically efficient than ones computed using automatic methods but that the gains are not as significant as in the massless case, at least at present. The method opens the door to applications at two-loop order, where numerical efficiency and improvements in the reconstruction methodology are more crucial, especially with regards to the number of free parameters in the ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

4D beam matrix reconstruction in particle accelerators

Transverse beam parameters in particle accelerators are commonly described using the Twiss parameters, which are experimentally accessible yet inherently limited because they neglect correlations between different transverse coordinates. Such correlations frequently arise from uncompensated cathode magnetic fields or misaligned focusing quadrupoles, affecting beam quality and accelerator performance. To address this limitation, we propose and validate a novel diagnostic method for the complete four-dimensional (4D) transverse beam matrix. Our method involves passing the beam through a beamline comprising both conventional and skew quadrupole magnets, followed by downstream measurements of the resulting two-dimensional (2D) beam profiles. These measurements represent distinct 2D projections of the underlying 4D transverse phase–space distribution. By systematically varying quadrupole strengths, multiple independent projections of the beam phase space are obtained. We reconstruct the original 4D beam matrix from these measured projections using an optimization-based least-square fit, providing fast and robust reconstruction regardless of the specific beamline configuration. Through extensive numerical simulations and realistic particle-tracking studies, we demonstrate the diagnostic’s accuracy, robustness, and capability to achieve reconstruction uncertainties smaller than measurement errors, particularly when employing sufficient numbers of quadrupole scans. This method presents a powerful and flexible approach for comprehensive beam characterization and accelerator tuning.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Space charge dominated momentum spread and compensation strategies in the post-linac section of Proton Improvement Plan-II at Fermilab

The upcoming Proton Improvement Plan-II (PIP-II), designated for enhancements to the Fermilab accelerator complex, features a new 800 MeV superconducting linac and a Beam Transfer Line (BTL) to transport the beam to the existing Booster synchrotron. To mitigate the space charge tune shift associated with a high intensity accumulated beam, the low emittance linac beam is used to paint the ring phase space both transversely and longitudinally. To prevent losses caused by particles injected outside the rf separatrix while painting longitudinal phase space, the momentum spread of the incoming beam should not exceed 2.1 x 10^-4. Detailed simulations showed that due to space charge, the rms momentum spread increases to 4 x 10^-4 while it is transported in the BTL --about twice the allowable limit. In this paper, we outline a mitigation strategy involving a debuncher cavity. We discuss location, operating frequency, and gap voltage under both nominal and perturbed beam conditions, specifically accounting for momentum jitter. The impact of cavity misalignments is also assessed. The paper concludes by recommending an optimized configuration.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Calibration of a soft secondary vertex tagger using proton-proton collisions at s = 13 TeV with the ATLAS detector

Several processes studied by the ATLAS experiment at the Large Hadron Collider produce low-momentum b-flavored hadrons in the final state. This paper describes the calibration of a dedicated tagging algorithm that identifies b-flavored hadrons outside of hadronic jets by reconstructing the soft secondary vertices originating from their decays. The calibration is based on a proton-proton collision dataset at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 140 fb -1 . Scale factors used to correct the algorithm’s performance in simulated events are extracted for the b-tagging efficiency and the mistag rate of the algorithm using a data sample enriched in $t\overline{t}$ events. Several orthogonal measurement regions are defined, binned as a function of the multiplicities of soft secondary vertices and jets containing a b-flavored hadron in the event. The mistag rate scale factors are estimated separately for events with low and high average numbers of interactions per bunch crossing. The results, which are derived from events with low missing transverse momentum, are successfully validated in a phase space characterized by high missing transverse momentum and therefore are applicable to new physics searches carried out in either phase space regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Dosimetric and biological impact of activity extravasation of radiopharmaceuticals in PET imaging

The increasing use of nuclear medicine and PET imaging has intensified scrutiny of radiotracer extravasation. To our knowledge, this topic is understudied but holds great potential for enhancing our understanding of extravasation in clinical PET imaging. This work aims to (1) quantify the absorbed doses from radiotracer extravasation in PET imaging, both locally at the site of extravasation and with the extravasation location as a source of exposure to bodily organs and (2) assess the biological ramifications within the injection site at the cellular level. A radiation dosimetry simulation was performed using a whole-body 4D Extended Cardiac-Torso (XCAT) phantom embedded in the GATE Monte Carlo platform. A 10-mCi dose of 18 F-FDG was chosen to simulate a typical clinical PET scan scenario, with 10% of the activity extravasated in the antecubital fossa of the right arm of the phantom. The extravasation volume was modeled as a 5.5 mL rectangle in the hypodermal layer of skin. Absorbed dose contributions were calculated for the first two half-lives, assuming biological clearance thereafter. Dose calculations were performed as absorbed doses at the organ and skin levels. Energy deposition was simulated both at the local extravasation site and in multiple organs of interest and converted to absorbed doses based on their respective masses. Each simulation was repeated ten times to estimate Monte Carlo uncertainties. Biological impacts on cells within the extravasated volume were evaluated by randomizing cells and exposing them to a uniform radiation source of 18 F and 68 Ga. Particle types, their energies, and direction cosines were recorded in phase space files using a separate Geant4 simulation to characterize their entry into the nucleus of the cellular volume. Subsequently, the phase space files were imported into the TOPAS-nBio simulation to assess the extent of DNA damage, including double-strand breaks (DSBs) and single-strand breaks (SSBs). Organ-level dosimetric estimations are presented for 18 F and 68 Ga radionuclides in various organs of interest. With 10% extravasation, the hypodermal layer of the skin received the highest absorbed dose of 1.32 ± 0.01 Gy for 18 F and 0.99 ± 0.01 Gy for 68 Ga. The epidermal and dermal layers received absorbed doses of 0.07 ± 0.01 Gy and 0.13 ± 0.01 Gy for 18 F, and 0.14 ± 0.01 Gy and 0.29 ± 0.01 Gy for 68 Ga, respectively. In the extravasated volume, 18 F caused an average absorbed dose per nucleus of 0.17 ± 0.01 Gy, estimated to result in 10.58 ± 0.50 DSBs and 268.11 ± 12.43 SSBs per nucleus. For 68 Ga, the absorbed dose per nucleus was 0.11 ± 0.01 Gy, leading to an estimated 6.49 ± 0.34 DSBs and 161.24 ± 8.12 SSBs per nucleus. Absorbed doses in other organs were on the order of micro-gray (µGy). The likelihood of epidermal erythema resulting from extravasation during PET imaging is low, as the simulated absorbed doses to the epidermis remain below the thresholds that trigger such effects. Moreover, the organ-level absorbed doses were found to be clinically insignificant across various simulated organs. The minimal DNA damage at the extravasation site suggests that long-term harm, such as radiation-induced carcinogenesis, is highly unlikely.

DNA strand breaks↗

DOE OSTI2024

Transverse emittance reduction in muon beams by ionization cooling

Accelerated muon beams have been considered for the next-generation studies of high-energy lepton–antilepton collisions and neutrino oscillations. However, high-brightness muon beams have not yet been produced. The main challenge for muon acceleration and storage stems from the large phase-space volume occupied by the beam, derived from the production mechanism of muons through the decay of pions. The phase-space volume of the muon beam can be decreased through ionization cooling. Here we show that ionization cooling leads to a reduction in the transverse emittance of muon beams that traverse lithium hydride or liquid hydrogen absorbers in the Muon Ionization Cooling Experiment. Our results represent a substantial advance towards the realization of muon-based facilities that could operate at the energy and intensity frontiers.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Quantifying carrier density in monolayer MoS 2 by optical spectroscopy

The successful design and device integration of nanoscale heterointerfaces hinges upon precise manipulation of both ground- and excited-state charge carrier (electron and hole) densities. However, it is particularly challenging to quantify these charge carrier densities in nanoscale materials, leading to uncertainties in the mechanisms of many carrier density-dependent properties and processes. Here, we demonstrate a method that utilizes steady-state and transient absorption spectroscopies to correlate monolayer MoS 2 electron density with the easily measured metric of excitonic optical absorption quenching in a variety of mixed-dimensionality s-SWCNT/MoS 2 heterostructures. By employing a 2D phase-space filling model, the resulting correlation elucidates the relationship between charge density, local dielectric environment, and concomitant excitonic properties. The phase-space filling model is also able to describe existing trends from the literature on transistor-based measurements on MoS 2 , WS 2 , and MoSe 2 monolayers that were not previously compared to a physical model, providing additional support for our method and results. The findings provide a pathway to the community for estimating both ground- and excited-state carrier densities in a wide range of TMDC-based systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗

DOE OSTI2024

Shearless effective barriers to chaotic transport induced by even twin islands in nontwist systems

For several decades now it has been known that systems with shearless invariant tori, nontwist Hamiltonian systems, possess barriers to chaotic transport. These barriers are resilient to breakage under perturbation and therefore regions where they occur are natural places to look for barriers to transport. Here we describe a kind of effective barrier that persists after the shearless torus is broken. Because phenomena are generic, for convenience we study the standard nontwist map (SNM), an area-preserving map that violates the twist condition locally in the phase space. The barrier occurs in nontwist systems when twin even period islands are present, which happens for a broad range of parameter values in the SNM. With a phase space composed of regular and irregular orbits, the movement of chaotic trajectories is hampered by the existence of shearless curves, total barriers, and a network of partial barriers formed by the stable and unstable manifolds of the hyperbolic points. Being a degenerate system, the SNM has twin islands and, consequently, twin hyperbolic points. We show that the structures formed by the manifolds intrinsically depend on period parity of the twin islands. For this even scenario the structure that we call a torus free barrier occurs because the manifolds of different hyperbolic points form an intricate chain atop a dipole configuration and the transport of chaotic trajectories through the chain becomes a rare event. This structure impacts the emergence of transport, the escape basin for chaotic trajectories, the transport mechanism, and the chaotic saddle. The case of odd periodic orbits is different: we find for this case the emergence of transport immediately after the breakup of the last invariant curve, and this leads to a scenario of higher transport, with intricate escape basin boundary and a chaotic saddle with nonuniformly distributed points.

classical mechanics↗