SEARCH · Engineering Papers

Results for “monotonicity”

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 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Stellar migration and chemical enrichment in the milky way disc: a hybrid model

ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

SDSS-IV MaNGA: global properties of kinematically misaligned galaxies

ABSTRACT We select 456 gas–star kinematically misaligned galaxies from the internal Product Launch-10 of MaNGA survey, including 74 star-forming (SF), 136 green-valley (GV), and 206 quiescent (QS) galaxies. We find that the distributions of difference between gas and star position angles for galaxies have three local peaks at ∼0°, 90°, and 180°. The fraction of misaligned galaxies peaks at log (M*/M⊙) ∼ 10.5 and declines to both low- and high-mass end. This fraction decreases monotonically with increasing star formation rate and specific star formation rate. We compare the global parameters including gas kinematic asymmetry Vasym, H i detection rate and mass fraction of molecular gas, effective radius Re, Sérsic index n as well as spin parameter $\lambda _{R_e}$ between misaligned galaxies and their control samples. We find that the misaligned galaxies have lower H i detection rate and molecular gas mass fraction, smaller size, higher Sérsic index, and lower spin parameters than their control samples. The SF and GV misaligned galaxies are more asymmetric in gas velocity fields than their controls. These pieces of observational evidence point to the gas accretion scenario followed by angular momentum redistribution from gas–gas collision, leading to gas inflow and central star formation for the SF and GV misaligned galaxies. We propose three possible origins of the misaligned QS galaxies: (1) external gas accretion, (2) merger, and (3) GV misaligned galaxies evolve into QS galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Constructing high-fidelity halo merger trees in abacussummit

ABSTRACT Tracking the formation and evolution of dark matter haloes is a critical aspect of any analysis of cosmological N-body simulations. In particular, the mass assembly of a halo and its progenitors, encapsulated in the form of its merger tree, serves as a fundamental input for constructing semi-analytic models of galaxy formation and, more generally, for building mock catalogues that emulate galaxy surveys. We present an algorithm for constructing halo merger trees from abacussummit, the largest suite of cosmological N-body simulations performed to date consisting of nearly 60 trillion particles, and which has been designed to meet the Cosmological Simulation Requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. Our method tracks the cores of haloes to determine associations between objects across multiple time slices, yielding lists of halo progenitors and descendants for the several tens of billions of haloes identified across the entire suite. We present an application of these merger trees as a means to enhance the fidelity of abacussummit halo catalogues by flagging and ‘merging’ haloes deemed to exhibit non-monotonic past merger histories. We show that this cleaning technique identifies portions of the halo population that have been deblended due to choices made by the halo finder, but which could have feasibly been part of larger aggregate systems. We demonstrate that by cleaning halo catalogues in this post-processing step, we remove potentially unphysical features in the default halo catalogues, leaving behind a more robust halo population that can be used to create highly accurate mock galaxy realizations from abacussummit.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Accuracy of power spectra in dissipationless cosmological simulations

Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, of scale-free models with a range of spectral indices n, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At 1 per cent level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, ${\sim}2{-}3 k_\Lambda$ at the very latest times we evolve to. At the $5{{\ \rm per\ cent}}$ level, accuracy extends up to wavenumbers of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple n-dependence, allowing a extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale-free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small-scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_\Lambda$ may need revision in light of our convergence analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The frequency and sizes of inner bars and nuclear rings in barred galaxies and their dependence on galaxy properties

ABSTRACT Using a volume- and mass-limited (D < 30 Mpc, $\log \, (M_{\star }/M_{\odot })\ge 9.75$) sample of 155 barred S0–Sd galaxies, I determine the fraction with secondary structures within their bars. Some 20 ± 3 per cent have a separate inner bar, making them double-barred; an identical fraction have nuclear rings, with $11^{+3}_{-2}$ per cent hosting both. The inner-bar frequency is a strong, monotonic function of stellar mass: only $4^{+3}_{-2}$ per cent of barred galaxies with $\log \, (M_{\star }/M_{\odot })= 9.75$–10.25 are double-barred, while 47 ± 8 per cent of those with $\log \, (M_{\star }/M_{\odot })\gt 10.5$ are. The nuclear-ring frequency is a strong function of absolute bar size: only $1^{+2}_{-1}$ per cent of bars with semimajor axes <2 kpc have nuclear rings, while $39^{+6}_{-5}$ per cent of larger bars do. Both inner bars and nuclear rings are absent in very late-type (Scd–Sd) galaxies. Inner bar size correlates with galaxy stellar mass, but is clearly offset to smaller sizes from the main population of bars. This makes it possible to define ‘nuclear bars’ in a consistent fashion, based on stellar mass. There are eight single-barred galaxies where the bars are nuclear-bar-sized; some of these may be systems where an outer bar failed to form, or previously double-barred galaxies where the outer bar has dissolved. Inner bar size is even more tightly correlated with host bar size, which is likely the primary driver. In contrast, nuclear ring size is only weakly correlated with galaxy mass or bar size, with more scatter in size than is true of inner bars.

Astronomy & Astrophysics↗

DOE OSTI2024

The structure of the stellar halo of the Andromeda galaxy explored with the NB515 for Subaru/HSC – I. New insights on the stellar halo up to 120 kpc

ABSTRACT We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and Canada France Hawaii Telescope/MegaCam g and i bands. We succeed in separating M31’s halo stars from foreground contamination with $\sim$90 per cent accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with the Giant Southern Stream (GSS) based on their photometric metallicity estimates. We also produce the distance and photometric metallicity estimates for the known substructures. While these quantities for the GSS are reproduced in our study, we find that the north-western stream shows a steeper distance gradient than found in an earlier study, suggesting that it is likely to have formed in an orbit closer to the Milky Way. For two streams in the eastern halo (Stream C and D), we identify distance gradients that had not been resolved. Finally, we investigate the global halo photometric metallicity distribution and surface brightness profile using the NB515-selected halo stars. We find that the surface brightness of the metal-poor and metal-rich halo populations, and the all population can be fitted to a power-law profile with an index of $\alpha =-1.65\pm 0.02$, $-2.82\pm 0.01$, and $-2.44\pm 0.01$, respectively. In contrast to the relative smoothness of the halo profile, its photometric metallicity distribution appears to be spatially non-uniform with non-monotonic trends with radius, suggesting that the halo population had insufficient time to dynamically homogenize the accreted populations.

Ogami, Itsuki (ORCID:0000000182394549)↗

DOE OSTI2024

Fundamental physics with ESPRESSO: a new determination of the D/H ratio towards PKS1937-101

ABSTRACT Primordial abundances of light elements are sensitive to the physics of the early Universe and can directly constrain cosmological quantities, such as the baryon-to-photon ratio $\eta _{10}$, the baryon density, and the number of neutrino families. Deuterium is especially suited for these studies: its primordial abundance is sensitive and monotonically dependent on $\eta _{10}$, allowing an independent measurement of the cosmic baryon density that can be compared, for instance, against the Planck satellite data. The primordial deuterium abundance can be measured in high H i column density absorption systems towards distant quasars. We report here a new measurement, based on high-resolution ESPRESSO data, of the primordial D i abundance of a system at redshift $z \sim 3.572$, towards PKS1937-101. Using only ESPRESSO data, we find a D /H ratio of $2.638\pm 10^{-5}$, while including the available UVES data improves the precision, leading to a ratio of $2.608 \pm 10^{-5}$. The results of this analysis agree with those of the most precise existing measurements. We find that the relatively low column density of this system ($\log {N_{\rm H_I}/ {\rm cm}^{-2}}\sim 18$) introduces modelling uncertainties, which become the main contributor to the error budget.

Guarneri, Francesco (ORCID:0000000347409762)↗

DOE OSTI2025

The role of electric dominance for particle injection in relativistic reconnection

ABSTRACT Magnetic reconnection in relativistic plasmas – where the magnetization $\sigma \gg 1$ – is regarded as an efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. We employ two-dimensional (2D) particle-in-cell simulations of relativistic pair-plasma reconnection with vanishing guide field and outflow boundaries to quantify the impact of the energy gain occurring in regions of electric dominance ($E\gt B$) for the early stages of particle acceleration (i.e. the ‘injection’ stage). Given an injection threshold energy $\epsilon ^\ast =\sigma /4$ for the particles that eventually attain energy $\epsilon _{\rm T}$ by time T, we calculate the mean fractional contribution $\zeta (\epsilon ^\ast ,\epsilon _{\rm T})$ by $E\gt B$ fields to particle energization at the time when the threshold $\epsilon ^\ast$ is reached. We find that $\zeta$ monotonically increases with $\sigma$ and $\epsilon _{\rm T}$; for $\sigma \gtrsim 50$ and $\epsilon _{\rm T}/\sigma \gtrsim 8$, we find that $\gtrsim 80~{{\ \rm per\ cent}}$ of the energy gain obtained before reaching $\epsilon ^\ast =\sigma /4$ occurs in $E\gt B$ regions. We find that $\zeta$ is independent of simulation box size $L_x$, as long as $\epsilon _{\rm T}$ is normalized to the maximum particle energy, which scales as $\epsilon _{\rm max}\propto L_{\rm x}^{1/2}$ in 2D. The distribution of energy gains $\epsilon _{\chi }$ acquired in $E\gt B$ regions can be modelled as $\mathrm{ d}N/\mathrm{ d}\epsilon _{\chi }\propto \epsilon _{\chi }^{-0.35}\exp [-(\epsilon _{\chi }/0.06\, \sigma)^{0.5}]$. Our results help assess the role of electric dominance in relativistic reconnection with vanishing guide fields, which is realized in the magnetospheres of black holes and neutron stars.

Gupta, Sanya (ORCID:0000000151944384)↗

DOE OSTI2020

A 16 deg2 survey of emission-line galaxies at z < 1.6 from HSC-SSP PDR2 and CHORUS

Abstract We have conducted a comprehensive survey of emission-line galaxies at z ≲ 1.6 based on narrow-band (NB) imaging data taken with Hyper Suprime-Cam (HSC) on the Subaru telescope. In this paper, we update the catalogs of Hα, [O iii], and [O ii] emission-line galaxies using the data from the second Public Data Release (PDR2) of the Subaru Strategic Program (SSP) of the HSC and Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey along with the spectroscopic redshifts for 2019 emission-line galaxies selected with the PDR1 data. The wider effective coverage of NB816 and NB921, 16.3 deg2 and 16.9 deg2, respectively, are available in the Deep and UltraDeep layers of HSC-SSP from the PDR2. The CHORUS survey provides us with data with additional three NBs (NB527, NB718, and NB973) in the COSMOS field in the UltraDeep layer (1.37 deg2). The five NB datasets allow us to investigate the star-forming galaxies presenting emission-lines at 14 specific redshifts ranging from z ∼ 1.6 down to z ∼ 0.05. We revisit the distribution of large-scale structures and luminosity functions (LFs) for the emission-line galaxies with the large samples of 75377 emission-line galaxies selected. The redshift revolution of LFs shows that the star formation rate densities (SFRDs) decreases monotonically from z ∼ 1.6, which is consistent with the cosmic SFRD known to-date. Our samples of emission-line galaxies covering a sufficiently large survey volume are useful to investigate the evolution of star-forming galaxies since the cosmic noon in a wide range of environments including galaxy clusters, filaments, and voids.

Hayashi, Masao (ORCID:0000000293217406)↗

DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗

DOE OSTI2026

Rapidity-dependent spin decomposition of the nucleon

We revisit the two-dimensional Fourier transform of generalized parton distributions (GPDs) at nonzero skewness. At 𝜂 = 0 it reduces to the standard impact-parameter density, while at 𝜂 ≠ 0 it is an off-forward amplitude that we interpret as a genuine parton–nucleon correlation. Its overall strength (the transverse-plane integral of the density) is fixed by the GPD at the kinematic point 𝑡 =−𝑐 𝜂 =−4⁢𝜂 2 ⁢𝑚$^{2}_{𝑁}$/(1 − 𝜂 2 ) and decreases monotonically with the rapidity gap Δ⁢𝑦 = ln⁡[(1+𝜂)/(1−𝜂)] = 2 artanh⁡(𝜂). This rapidity dependence implies rapidity-modified Ji identities that connect helicity, orbital, and total angular momenta of the correlation in closed form. To quantify these effects, we construct leading-twist quark and gluon GPDs in a string-based conformal framework: conformal moments are parametrized by linear open- and closed-string Regge trajectories with slopes constrained by parton distribution functions (PDFs), hadron/glueball spectroscopy, and form-factor data, and GPDs are reconstructed over the full (𝑥,𝜂,𝑡) domain by Mellin-Barnes inversion with next-to-leading order evolution. We find qualitative agreement (and fair quantitative agreement within quoted uncertainties) for several moments and selected nonsinglet 𝑥-space channels at 𝜇 = 2 GeV when compared with lattice QCD, while we also identify channels with visible tension and discuss likely sources (PDF priors and 𝑡-slope systematics).

Gauge-gravity dualities↗

DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Mass radius and D-term of atomic nuclei in relativistic mean field theory

Based on relativistic mean field theory for atomic nuclei, we compute the mass radius and other radii associated with the energy momentum tensor for dozens of spin-0 nuclei across the nuclear chart. We also compute the D-term of these nuclei, the forward limit of the gravitational form factor 𝐷⁡(𝑡=0)=𝐷. The dependence on the neutron number 𝑁 is systematically studied for calcium (Ca), nickel (Ni), zirconium (Zr), tin (Sn), and lead (Pb) isotopes. Remarkably, |𝐷| does not monotonically increase with 𝑁. Instead, it exhibits local maxima and minima when 𝑁 equals a magic number and even a submagic number. This results in characteristic kinks in the mass, scalar, tensor, and shear radii of these isotopes. Our work for the first time elucidates the strong sensitivity of the various mechanical properties of nuclei to the nuclear shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Molecular bond stabilization in the strong-field dissociation of O 2 +

We theoretically examine the rotational and vibrational dynamics of O 2 + molecular ions exposed to intense, short laser pulses for conditions realized in contemporary pump-probe experiments. We solve the time-dependent Schrödinger equation within the Born-Oppenheimer approximation for an initial distribution of randomly aligned molecular ions. For fixed peak intensities, our numerical results show that total, angle-integrated O 2 + → O( 3 P) + O + ( 4 S 0 ) dissociation yields do not monotonically increase with increasing infrared-probe pulse duration. We find this pulse-duration-dependent stabilization to be consistent with the transient trapping of nuclear probability density in a light-induced (bond-hardening) potential-energy surface and robust against rotational excitation. We analyze this stabilization effect and its underlying bond-hardening mechanism (i) in the time domain, by following the evolution of partial nuclear probability densities associated with the dipole-coupled O 2 + (a 4 Π u ) and O 2 + (f 4 Π g ) cationic states, and (ii) in the frequency domain, by examining rovibrational quantum-beat spectra for the evolution of the partial nuclear probability densities associated with these states. In this work, our analysis reveals the characteristic timescale for the bond-hardening mechanism in O 2 + and explains the onset of bond stabilization for sufficiently long pulse durations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Strong-field physics in three-dimensional topological insulators

We investigate theoretically the strong-field regime of light-matter interactions in the topological-insulator class of quantum materials. In particular, we focus on the process of nonperturbative high-order harmonic generation from the paradigmatic three-dimensional topological insulator bismuth selenide (Bi 2 Se 3 ) subjected to intense midinfrared laser fields. We analyze the contributions from the spin-orbit-coupled bulk states and the topological surface bands separately and reveal a major difference in how their harmonic yields depend on the ellipticity of the laser field. Furthermore, bulk harmonics show a monotonic decrease in their yield as the ellipticity increases, in a manner reminiscent of high harmonic generation in gaseous media. However, the surface contribution exhibits a highly nontrivial dependence, culminating with a maximum for circularly polarized fields. We attribute the observed anomalous behavior to (i) the enhanced amplitude and the circular pattern of the interband dipole and the Berry connections in the vicinity of the Dirac point and (ii) the influence of the higher-order, hexagonal warping terms in the Hamiltonian, which are responsible for the hexagonal deformation of the energy surface at higher momenta. The latter are associated directly with spin-orbit-coupling parameters. Our results thus establish the sensitivity of strong-field-driven high harmonic emission to the topology of the band structure as well as to the manifestations of spin-orbit interaction.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Photoionization branching ratios of spin-orbit doublets far above thresholds: Interchannel and relativistic effects in the noble gases

Photoionization branching ratios of all of the spin-orbit doublets in the noble gas atoms Ne, Ar, Kr, and Xe have been investigated theoretically at energies well above the thresholds. The results confirm in all cases that the general behavior of the branching ratios is to decrease monotonically below their statistical value with increasing photon energy. Here, this effect gets more pronounced with increasing Z. In addition, the branching ratios can be strongly affected by interchannel coupling with inner shell photoionization channels. Although these effects are strongest in the neighborhood of the inner shell thresholds, they can persist over a broad range of energies,

74 ATOMIC AND MOLECULAR PHYSICS↗