Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum harmonic oscillator”

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.

90 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations

Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Jarzynski Equality for Conditional Stochastic Work

It has been established that the inclusive work for classical, Hamiltonian dynamics is equivalent to the two-time energy measurement paradigm in isolated quantum systems. However, a plethora of other notions of quantum work has emerged, and thus the natural question arises whether any other quantum notion can provide motivation for purely classical considerations. In the present analysis, we propose the conditional stochastic work for classical, Hamiltonian dynamics, which is inspired by the one-time measurement approach. This novel notion is built upon the change of expectation value of the energy conditioned on the initial energy surface. As main results, we obtain a generalized Jarzynski equality and a sharper maximum work theorem, which account for how non-adiabatic the process is. Our findings are illustrated with the parametric harmonic oscillator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Entanglement Effect and Angular Momentum Conservation in a Nonseparable Tunneling Treatment

The important, and often dominant, role of tunneling in low temperature kinetics has resulted in numerous theoretical explorations into the methodology for predicting it. Nevertheless, there are still key aspects of the derivations that are lacking, particularly for nonseparable systems in the low temperature regime, and further explorations of the physical factors affecting the tunneling rate are warranted. In this work we obtain a closed-form rate expression for the tunneling rate constant that is a direct analog of the rigid-rotor-harmonic-oscillator expression. This expression introduces a novel "entanglement factor" that modulates the reaction rate. Furthermore, we are able to extend this expression, which is valid for nonseparable systems at low temperatures, to properly account for the conservation of angular momentum. In contrast, previous calculations have considered only vibrational transverse modes and so effectively employ a decoupled rotational partition function for the orientational modes. We also suggest a simple theoretical model to describe the tunneling effects in the vicinity of the crossover temperature (the temperature where tunneling becomes the dominating mechanism). This model allows one to naturally classify, interpret, and predict experimental data. Among other things, it quantitatively explains in simple terms the so-called "quantum bobsled" effect, also known as the negative centrifugal effect, which is related to curvature of the reaction path. Taken together, the expressions obtained here allow one to predict the thermal and E-resolved rate constants over broad ranges of temperatures and energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Anomalous two-photon Compton scattering

Abstract X-ray free-electron lasers can generate radiation pulses with extreme peak intensities at short wavelengths. This enables the investigation of laser–matter interactions in a regime of high fields, yet at a non-relativistic ponderomotive potential, where ordinary rules of light–matter interaction may no longer apply and nonlinear processes are starting to become observable. Despite small cross-sections, first nonlinear effects in the hard x-ray regime have recently been observed in solid targets, including x-ray-optical sum-frequency generation (XSFG), x-ray second harmonic generation (XSHG) and two-photon Compton scattering (2PCS). Nonlinear interactions of bound electrons in the x-ray range are fundamentally different from those dominating at optical frequencies. Whereas in the optical regime nonlinearities are predominantly caused by anharmonicities of the atomic potential in the chemical bonds, x-ray nonlinearities far above atomic resonances are expected to be due to nonlinear oscillations of quasi-free electrons, including inner-shell atomic electrons. While the quasi-free-electron model agrees reasonably well with the experimental data for XSFG and XSHG, 2PCS measurements have led to unexpected results: the energy of the nonlinearly scattered photons from non-relativistic electrons shows a substantial unexpected red shift in addition to the Compton shift that is well beyond that predicted by a nonlinear quantum electrodynamics model for free electrons. A potential explanation for the spectral broadening is based on a previously unexplored scattering process that involves the whole atom rather than just quasi-free electrons. A first simulation that includes the atomic binding potential was successful in describing a broadening of the spectrum of the nonlinearly scattered photons to longer wavelengths for soft x-rays. However, the same model does not show any broadening at hard x-ray wavelengths, which is in agreement with other simulation approaches. To this point no calculation has been able to reproduce the experimentally observed broadening. Here we present further experimental data of 2PCS for an extended parameter range using additional diagnostics. In particular, we present measurements of the electron momentum distribution during the interaction that strongly suggest that the spectral broadening is not caused by an increased plasma temperature. We extend our measurement of the magnitude of the red shift in beryllium to > 1.9 k e V in addition to the Compton shift expected for free electrons and expand the measurement of the angular distribution to include forward scattering angles. We also present first measurements of 2PCS from diamond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fermi surface and Berry phase analysis for Dirac nodal line semimetals: Cautionary tale from SrGa 2 and BaGa 2

A Berry phase of odd multiples of 𝜋 inferred from quantum oscillations (QOs) has often been treated as evidence for nontrivial reciprocal space topology. However, disentangling the Berry phase values from the Zeeman effect and the orbital magnetic moment is often challenging. In centrosymmetric compounds, the case is simpler as the orbital magnetic moment contribution is negligible. Although the Zeeman effect can be significant, it is usually overlooked in most studies of QOs in centrosymmetric compounds. Here, we present a detailed study on the nonmagnetic centrosymmetric SrGa 2 and BaGa 2 , which are predicted to be Dirac nodal line semimetals based on density functional theory (DFT) calculations. Evidence of the nontrivial topology is found in magnetotransport measurements. The Fermi surface topology and band structure are carefully studied through a combination of angle-dependent QOs, angle-resolved photoemission spectroscopy (ARPES), and DFT calculations, where the nodal line is observed in the vicinity of the Fermi level. Strong de Haas–van Alphen fundamental oscillations associated with higher harmonics are observed in both compounds, which are well fitted by the Lifshitz-Kosevich (LK) formula. However, even with the inclusion of higher harmonics in the fitting, we found that the Berry phases cannot be unambiguously determined when the Zeeman effect is included. We revisit the LK formula and analyze the phenomena and outcomes that were associated with the Zeeman effect in previous studies. Our experimental results confirm that SrGa 2 and BaGa 2 are Dirac nodal line semimetals. Additionally, we highlight the often overlooked role of spin-damping terms in Berry phase analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photon-recoil and laser-focusing limits to Rydberg gate fidelity

Limits to Rydberg gate fidelity that arise from the entanglement of internal states of neutral atoms with the motional degrees of freedom due to the momentum kick from photon absorption and re-emission is quantified. This occurs when the atom is in a superposition of internal states but only one of these states is manipulated by visible or UV photons. The Schrödinger equation that describes this situation is presented and two cases are explored. In the first case, the entanglement arises because the spatial wave function shifts due to the separation in time between excitation and stimulated emission. For neutral atoms in a harmonic trap, the decoherence can be expressed within a sudden approximation when the duration of the laser pulses are shorter than the harmonic oscillator period. In this limit, the decoherence is given by simple analytic formulas that account for the momentum of the photon, the temperature of the atoms, the harmonic oscillator frequency, and atomic mass. In the second case, there is a reduction in gate fidelity because the photons causing absorption and stimulated emission are in focused beam modes. This leads to a dependence of the optically induced changes in the internal states on the center of mass atomic position. In the limit where the time between pulses is short, the decoherence can be expressed as a simple analytic formula involving the laser waist, temperature of the atoms, the trap frequency, and the atomic mass. These limits on gate fidelity are studied for the standard π-2π-π Rydberg gate and a protocol based on a single adiabatic pulse with a Gaussian envelope.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Ultrafast Nanoimaging of Spin-Mediated Shear Waves in an Acoustic Cavity

Strong spin–lattice coupling in van der Waals (vdW) magnets shows potential for innovative magneto-mechanical applications. Here, nanoscale and picosecond imaging by ultrafast electron microscopy reveal heterogeneous spin-mediated coherent acoustic phonon dynamics in a thin-film cavity of the vdW antiferromagnet FePS3. The harmonics of the interlayer shear acoustic modes are observed, in which the even and odd harmonics exhibit distinct nanoscopic dynamics. Corroborated by acoustic wave simulation, the role of defects in forming even harmonics is elucidated. Above the Néel temperature (T N ), the interlayer shear acoustic harmonics are suppressed, while the in-plane traveling wave is predominantly excited. The dominant acoustic dynamics shifts from the out-of-plane shear to the in-plane traveling wave across T N , demonstrating that magnetic properties can influence phonon scattering pathways. In conclusion, the spatiotemporally resolved structural characterization provides valuable nanoscopic insights for interlayer-shear-mode-based acoustic cavities, opening up possibilities for magneto-mechanical applications of vdW magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Efficiency-optimized relativistic plasma harmonics for extreme fields

Bright harmonic radiation from relativistically oscillating laser plasmas offers a direct route for generating extreme electromagnetic fields. Theory predicts that under optimized conditions, the plasma medium can support strong spatiotemporal compression of laser energy in a coherent harmonic focus (CHF), delivering intensity boosts many orders of magnitude greater than the incident driving laser pulse. Although diffraction-limited performance (spatial compression) and attosecond phase locking (temporal compression) have been demonstrated experimentally, efficient coupling of relativistically intense laser pulse energy into the emitted harmonic cone has not been realized so far. Here we demonstrate that this highly nonlinear interaction can be tailored to deliver the maximum conversion efficiencies predicted from simulations. By fine-tuning the temporal profile of the driving laser on sub-picosecond (<10 −12 s) timescales, energies >9 mJ between the 12th and 47th harmonics are observed. These results are in agreement with the theoretically expected efficiency dependence on harmonic order, verifying that optimal conditions have been achieved in the generation process. This is the important final element required to achieve the expected intensity boosts from a CHF in experiments. Although obtaining spatiotemporal compression and optimal efficiency simultaneously remains challenging, the path to realizing extreme optical field strengths approaching the critical field of quantum electrodynamics (the Schwinger limit at >10 16 V cm −1 or >10 29 W cm −2 ) is now open, permitting all-optical studies of the quantum vacuum and new frontiers for intense attosecond science.

High-harmonic generation↗
DOE OSTI2021

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↗
DOE OSTI2022

Natural orbitals for the ab initio no-core configuration interaction approach

Ab initio no-core configuration interaction (NCCI) calculations for the nuclear many-body problem have traditionally relied upon an antisymmetrized product (Slater determinant) basis built from harmonic oscillator orbitals. The accuracy of such calculations is limited by the finite dimensions which are computationally feasible for the truncated many-body space. We therefore seek to improve the accuracy obtained for a given basis size by optimizing the choice of single-particle orbitals. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which maximizes the occupation of low-lying orbitals, thus accelerating convergence in a configuration-interaction basis, while also possibly providing physical insight into the single-particle structure of the many-body wave function. We describe the implementation of natural orbitals in the NCCI framework and examine the nature of the natural orbitals thus obtained, the properties of the resulting many-body wave functions, and the convergence of observables. After taking 3 He as an illustrative testbed, we explore aspects of NCCI calculations with natural orbitals for the ground state of the p-shell neutron halo nucleus 6 He .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

State-selective probing of CO 2 autoionizing inner valence Rydberg states with attosecond extreme ultraviolet four-wave-mixing spectroscopy

Nonlinear spectroscopies can disentangle spectra that are congested due to inhomogeneous broadening. Here, in conjunction with theoretical calculations, attosecond extreme ultraviolet (XUV) four-wave-mixing (FWM) spectroscopy is utilized to probe the dynamics of autoionizing inner valence excited Rydberg states of the polyatomic molecule, CO 2 . This tabletop nonlinear technique employs a short attosecond XUV pulse train and two noncollinear, few-cycle near-infrared pulses to generate background-free XUV wave-mixing signals. FWM emission is observed from the n=5-7 states of the Henning sharp ndσ g Rydberg series that converges to the ionic $\widetilde{B}$ 2 Σ$^{+}_{u}$ state. However, these transient emission signals decay with lifetimes of 33 ± 6, 53 ± 2, and 94 ± 2 fs, respectively, which calculations show are consistent with the lifetimes of the short-lived n=6-8 members of the nsσ g character Henning diffuse Rydberg series. The oscillator strengths of transitions between states involved in all possible resonant FWM processes are calculated, verifying that the nonlinear spectra are dominated by pathways described by an initial excitation to the diffuse nsσ g Rydberg series and emission from the sharp ndσ g Rydberg series. The results substantiate not only that attosecond XUV FWM spectroscopy produces rigorous and meaningful measurements of ultrafast dynamics in polyatomic systems, but also that nonlinear spectroscopic techniques are versatile tools to selectively probe dynamics that are otherwise difficult to access.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Multidimensional coherent spectroscopy of light-driven states and their collective modes in multiband superconductors

We present a comprehensive theory of light-controlled multiband superconductivity, and apply it to predict distinctive signatures of light-driven superconducting (SC) states in terahertz multidimensional coherent spectroscopy (THz-MDCS) experiments. We first derive gauge-invariant Maxwell-Bloch equations for multiband BCS superconductors with spatial fluctuations. We consider driving electromagnetic fields determined self-consistently by Maxwell's equations. By calculating the THz-MDCS spectra measured experimentally in the clean SC limit, we identify unique signatures of finite-momentum Cooper-pairing states that live longer than the laser pulse. They are controlled by a pair of THz laser pulses with well-defined relative phase (pulse pair). The pseudospin oscillators that describe the properties of these SC states are parametrically driven by both finite-momentum Cooper pairing and by time oscillations of the order-parameter relative phase. We show that such strong parametric driving leads to drastic changes in the THz-MDCS spectral shape from the predictions of third-order nonlinear susceptibility calculations. These spectral changes strongly depend on the interband-to-intraband interaction ratio and on the collective modes of the light-driven state. For negligible interband interaction, the spectra show a transition with increasing field, from traditional pump-probe, four-wave-mixing, and third-harmonic generation peaks determined by the laser frequency to sidebands determined by the excitations of the driven system. These sidebands emerge from difference-frequency Raman processes in the nonequilibrium SC state. For interband couplings weaker than the intraband pairing, we show that the Leggett phase collective mode leads to harmonic sidebands around the traditional pump-probe peaks. Additional Higgs collective mode peaks result from light-induced inversion-symmetry breaking in a thin-film geometry. For strong interband coupling, we find a transition from a nonequilibrium finite Cooper-pair momentum state characterized by hybrid-Higgs amplitude mode peaks in THz-MDCS spectra to a driven state identified experimentally by the emergence of Floquet-type sidebands at bi-Higgs frequencies. Those dominant bi-Higgs-frequency satellites are manifestations of a new order parameter relative phase collective mode that characterizes the nonequilibrium SC state. The predicted interaction- and field-dependent transitions in the spectral profile allow us to propose THz-MDCS experiments for quantum tomography of light-driven superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Tidally Excited Modes and δ Scuti Pulsations in the Eclipsing Triple Star IM Persei

IM Per is a triple star system whose eclipsing pair masses and radii are accurate to within 1%. The Transiting Exoplanet Survey Satellite (TESS) light curve of the program target exhibits partial eclipses and multiple oscillations with millimagnitude-level amplitudes. It is found that the oscillations affect eclipse timing measurements. Binary modeling of the high-quality TESS data indicates that the eclipsing components of the triple system are twin dwarfs with parameters of M {sub 2}/M {sub 1} = 0.995, R {sub 2}/R {sub 1} = 0.901, and Δ(T{sub eff,1}–T{sub eff,2}) = 12 K in an eccentric (e = 0.049), detached configuration. The third light of l {sub 3} = 0.054 may mostly come from a G-type tertiary companion. Our predicted parallax of 1.52 ± 0.09 mas is concurrent with the Gaia measurement of 1.52 ± 0.05 mas. Multifrequency analysis of the outside-eclipse residuals reveals 22 significant pulsation signals: 4 in the gravity-mode region (0.03−2.22 day{sup −1}) and 18 in the pressure-mode region (9.19−25.12 day{sup −1}). Of the low frequencies, f {sub 11} and f {sub 14} are orbital harmonics that can be identified as tidally excited modes. The pulsation periods and constants for the high frequencies, and the position in the Cepheid instability strip demonstrate that the pulsating component of IM Per is a δ Sct variable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Comprehensive anisotropic linear optical properties of the Weyl semimetals TaAs and NbAs

We report that TaAs and NbAs are two of the earliest identified Weyl semimetals that possess many intriguing optical properties, such as chirality-dependent optical excitations and giant second harmonic generation (SHG). Linear and nonlinear optics have been employed as tools to probe the Weyl physics in these crystals. Here, we extend these studies to address two important aspects: determining the complete anisotropic dielectric response and exploring if and how they can reveal essential Weyl physics. We determine the complete anisotropic dielectric functions of TaAs and NbAs by combining spectroscopic ellipsometry and density functional theory (DFT). Parameterized Lorentz oscillators are reported from 1.2 - 6 eV (experiment) and 0 - 6 eV (DFT), and good agreement is shown between them. Both linear and nonlinear optical properties have been previously reported to reveal Weyl physics. We suggest that strong optical resonances from trivial bands are the likely origin of the large optical SHG previously reported at these energies. Furthermore, by comparing the contribution of a small k -space centered around the Weyl cones with the total linear dielectric function, we find that these contributions are highly anisotropic and are < 25% of the total dielectric function below 0.5 eV; abov1 eV, these contributions are negligible. Thus, the study of Weyl physics using optical techniques requires very low energies, and even there, a careful assessment is required in distinguishing the much smaller contributions of the Weyl bands from the dominant contributions of the trivial bands and Drude response to the total dielectric function.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗