Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bound state in the continuum”

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

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗

Quantum-Well Bound States in Graphene Heterostructure Interfaces

We present experimental evidence of electronic and optical interlayer resonances in graphene van der Waals heterostructure interfaces. Using the spectroscopic mode of a low-energy electron microscope (LEEM), we characterized these interlayer resonant states up to 10 eV above the vacuum level. Compared with nontwisted, AB-stacked bilayer graphene (AB BLG), an ≈ 0.2 Å increase was found in the interlayer spacing of 30° twisted bilayer graphene (30°-tBLG). In addition, we used Raman spectroscopy to probe the inelastic light-matter interactions. A unique type of Fano resonance was found around the D and G modes of the graphene lattice vibrations. This anomalous, robust Fano resonance is a direct result of quantum confinement and the interplay between discrete phonon states and the excitonic continuum.

36 MATERIALS SCIENCE↗

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Coronal Line Spectrum of the Iron Nova V3666 Ophiuchi

Optical and infrared spectroscopy are presented for the classical nova V3666 Ophiuchi (Nova Ophiuchi 2018 No.3) from 11 months after its outburst. At that epoch, V3666 Oph displayed an [Fe x] λ0.6375 emission line stronger than Hα and nearly a factor of 3 larger, relative to Hα, than in any nova observed previously. The [Fe x] feature, with an equivalent width exceeding 1300 Å, was enormously strong relative to the continuum emission as well. V3666 Oph also reached a state of high excitation for an Fe ii type nova, displaying coronal lines of [Fe xi], [S ix], and [Ar x], and permitted features that can be identified with the C vi (8–7), (9–8), and (10–9) transitions. In addition to [Fe x] and [Fe xi], V3666 Oph also displayed bright lines of [Fe vii], indicating that the prominent iron features were not confined to a single or narrow range of ionization states. Modeling the emission line spectrum of V3666 Oph with the photoionization code Cloudy indicated that the bulk of the emission line gas was matter-bounded, severely reducing the amount of low ionization gas and resulting in coronal lines that were very strong in comparison to lines such as H i and He i, and the familiar low excitation forbidden lines of nitrogen and oxygen. Such matter-bounded regions are observationally rare in novae since they are both faint and short-lived, offering a plausible explanation for why [Fe x] λ0.6375 had not be seen at this strength previously.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Photoionization from the Ground and Excited Vibrational States of H$^+_2$ and Its Deuterated Isotopologues

Photoionization cross sections and rate coefficients have been calculated for all bound vibrational levels of the 1sσ g state of H$^+_2$, HD+, and D$^+_2$. The Born–Oppenheimer approximation is employed in our calculation of vibrationally resolved photoionization cross sections. Vibrationally resolved and local thermal equilibrium photoionization rate coefficients are presented for photon temperatures less than 50,000 K and are found to be several orders of magnitude larger than previous results in the literature. Analytic fits for the vibrationally resolved and local thermal equilibrium photoionization rate coefficients are provided. Near-threshold oscillations in the vibrationally resolved photoionization are observed. A benchmark set of photoionization cross sections are presented. Fixed-nuclei photoionization cross sections are calculated using two-center true continuum wave functions and are verified by comparison with previous calculations and are found to be in excellent agreement in all cases. Data files for our set of benchmark cross sections, rate coefficients, and fitting parameters for H$^+_2$, HD + , and D$^+_2$ are available on Zenodo under an open-source Creative Commons Attribution license at doi:10.5281/zenodo.8304061.

74 ATOMIC AND MOLECULAR PHYSICS↗

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coma Morphology Due to an Extended Active Region and Implications for the Spin State of Comet Hale-Bopp

We show that the circular character of continuum structures observed in the coma of comet Hale-Bopp around the perihelion passage is most likely due to a dust jet from a large extended active region on the surface. Coma morphology due to a wide jet is different from that due to a narrow jet. The latter shows foreshortening effects due to observing geometry, wider jet produces more circular features. This circularization effect provides a self-consistent explanation for the evolution of near-perihelion coma morphology. No changes in the direction of the rotational angular momentum vector are required during this period in contrast to the models of Schleicher et al. This circularization effect also enables us to produce near-circular coma features in the S-E quadrant during 1997 late February and therefore questions the basic premise on which Sekanina bases his morphological arguments for a gravitationally bound satellite nucleus.

Samarasinha, Nalin H.↗

Insights into the emergence of mass from studies of pion and kaon structure

There are two mass generating mechanisms in the standard model of particle physics (SM). One is related to the Higgs boson and fairly well understood. The other is embedded in quantum chromodynamics (QCD), the SM’s strong interaction piece; and although responsible for emergence of the roughly 1 GeV mass scale that characterises the proton and hence all observable matter, the source and impacts of this emergent hadronic mass (EHM) remain puzzling. As bound states seeded by a valence-quark and -antiquark, pseudoscalar mesons present a simpler problem in quantum field theory than that associated with the nucleon. Consequently, there is a large array of robust predictions for pion and kaon properties whose empirical validation will provide a clear window onto many effects of both mass generating mechanisms and the constructive interference between them. This has now become significant because new-era experimental facilities, in operation, construction, or planning, are capable of conducting such tests and thereby contributing greatly to resolving the puzzles of EHM. These aspects of experiment, phenomenology, and theory, along with contemporary successes and challenges, are reviewed herein. In addition to providing an overview of the experimental status, we focus on recent progress made using continuum Schwinger function methods and lattice-regularised QCD. Further, advances made using other theoretical tools are also sketched. Our primary goal is to highlight the potential gains that can accrue from a coherent effort aimed at finally reaching an understanding of the character and structure of Nature’s Nambu–Goldstone modes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

On the origins of the continuum radiation of an underwater nanosecond pulsed discharge: an absolute-intensity optical emission spectroscopy study

The plasma properties of an underwater nanosecond pulsed discharge remain not fully understood despite being extensively studied for several decades. In this work, we focus on the continuum radiation generated in such discharges. The discharge is characterized by power measurements as well as by absolute emission spectroscopy. When observed, Stark broadenings of H α , H β and O (777 nm) are employed for electron number densities measurements. The discharge was generated by a 10 ns main voltage pulse followed by multiple secondary pulses, which last up to 4 μs after the primary pulse. It is shown that a peak power of 3.5 MW and energy of 35 mJ is coupled during the main voltage pulse. A quantitative estimation of the different possible continuum radiation sources is performed through analytical calculations. This includes emission (blackbody, free–bound and free–free bremsstrahlung radiations) and absorption (electron–ion and electron–neutral free–free inverse bremsstrahlung) mechanisms. Our results suggest that electron–neutral free–free bremsstrahlung is the principal mechanism responsible for the strong continuum radiation observed experimentally during the primary pulse. We also show that self-absorption through electron–neutral (and electron–ion) inverse bremsstrahlung plays an important role in the main discharge pulse. Further, our results indicate the non-negligible additional contribution of the H 2 continuum during the first reflected pulse which is likely ignited in bubbles generated by the first discharge pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗

Frameshifting Stimulatory Sequence Induces Large Structural Change of Ribosomal Proteins When Bound to E. coli Ribosomes

Biological macromolecular machines occupy a continuum of structural conformations to perform cellular tasks. Mapping this conformational space provides an insight into its functionality. While the cryo-electron microscopy resolution revolution has expanded our ability to characterize the conformational continuums, there are obstacles in structurally characterizing regions of high flexibility. These technical barriers have impeded characterization of flexible ribosomal proteins when the ribosome is interacting with mRNA stem-loop structures such as a frameshifting stimulatory sequence (FSS). Small-angle neutron/X-ray scattering and electron microscopy were used to study ribosomal samples and compared structural differences between a ribosome that is bound to an FSS stem-loop compared to a ribosome bound to linear mRNA. This comparison shows that a large protein stalk elongates by 22% when the 70S interacts with an mRNA stem-loop. Finally, our results suggest that ribosomal proteins have extensive flexibility and may influence important ribosomal mechanisms, such as those that involve FSS.

36 MATERIALS SCIENCE↗

Optimising the processing and storage of visibilities using lossy compression

The next-generation radio astronomy instruments are providing a massive increase in sensitivity and coverage, largely through increasing the number of stations in the array and the frequency span sampled. The two primary problems encountered when processing the resultant avalanche of data are the need for abundant storage and the constraints imposed by I/O, as I/O bandwidths drop significantly on cold storage. An example of this is the data deluge expected from the SKA Telescopes of more than 60 PB per day, all to be stored on the buffer filesystem. While compressing the data is an obvious solution, the impacts on the final data products are hard to predict. In this paper, we chose an error-controlled compressor – MGARD – and applied it to simulated SKA-Mid and real pathfinder visibility data, in noise-free and noise-dominated regimes. As the data have an implicit error level in the system temperature, using an error bound in compression provides a natural metric for compression. MGARD ensures the compression incurred errors adhere to the user-prescribed tolerance. To measure the degradation of images reconstructed using the lossy compressed data, we proposed a list of diagnostic measures, exploring the trade-off between these error bounds and the corresponding compression ratios, as well as the impact on science quality derived from the lossy compressed data products through a series of experiments. We studied the global and local impacts on the output images for continuum and spectral line examples. We found relative error bounds of as much as 10%, which provide compression ratios of about 20, have a limited impact on the continuum imaging as the increased noise is less than the image RMS, whereas a 1% error bound (compression ratio of 8) introduces an increase in noise of about an order of magnitude less than the image RMS. For extremely sensitive observations and for very precious data, we would recommend a 0.1% error bound with compression ratios of about 4. These have noise impacts two orders of magnitude less than the image RMS levels. At these levels, the limits are due to instabilities in the deconvolution methods. We compared the results to the alternative compression tool DYSCO, in both the impacts on the images and in the relative flexibility. MGARD provides better compression for similar error bounds and has a host of potentially powerful additional features.

Techniques: interferometric↗

Core spectroscopy of oxazole

We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second order and time-dependent density functional levels of theory. We demonstrate that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. Furthermore, the resonant AES, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.

74 ATOMIC AND MOLECULAR PHYSICS↗

Plasma screening in mid-charged ions observed by K-shell line emission

Dense plasma environment affects the electronic structure of ions via variations of the microscopic electrical fields, also known as plasma screening. This effect can be either estimated by simplified analytical models, or by computationally expensive and to date unverified numerical calculations. We have experimentally quantified plasma screening from the energy shifts of the bound-bound transitions in matter driven by the x-ray free electron laser (XFEL). This was enabled by identification of detailed electronic configurations of the observed Kα, Kβ and Kγ lines. This work paves the way for improving plasma screening models including connected effects like ionization potential depression and continuum lowering, which will advance the understanding of atomic physics in the Warm Dense Matter regime.

Physics - Plasma physics↗

Binding energy of the 𝑇 𝑏⁢𝑏 tetraquark from lattice QCD with relativistic and nonrelativistic heavy-quark actions

We present a new determination of the $b\bar{b}$𝑢⁢𝑑 (𝐽 𝑃 = 1 + , 𝐼 = 0) tetraquark binding energy using lattice quantum chromodynamics (QCD) with domain-wall light quarks and a nonperturbatively tuned three-parameter anisotropic-clover “relativistic” action for the 𝑏 quarks. We also perform a direct comparison with a reanalysis of data generated in prior work using a lattice-nonrelativistic QCD (NRQCD) action for the 𝑏 quarks and otherwise identical parameters. Using the new data with relativistic 𝑏 quarks from seven different ensembles with multiple lattice spacings and pion masses, we perform combined chiral and continuum extrapolations and obtain (𝑚 𝑇 𝑏⁢𝑏 −𝑚 𝐵 −𝑚 𝐵* ) RHQ =(−76 ±23) MeV. For the NRQCD data from five ensembles, we perform chiral-only extrapolations and obtain (𝑚 𝑇 𝑏⁢𝑏 −𝑚 𝐵 −𝑚 𝐵* ) NRQCD = (−74 ±17 ±10) MeV. The lower magnitude of the results obtained here, compared to the original analysis in [Phys. Rev. D 100, 014503 (2019)], is due to the use of the symmetric parts of the correlation matrices with local four-quark operators only.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Asymptotic consistency of the WSINDy algorithm in the limit of continuum data

In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.

asymptotic consistency↗